用python2.7.3运行以下代码非常好import numpy as np
class VolumeViewer(object):
def __init__(self, raw_dicom_stack, nframes, ds):
self.raw_dicom_stack = raw_dicom_stack
self.nframes = nframes
self.summed_image = np.zeros((self.raw_dicom_stack[0].shape[0],self.raw_dicom_stack[0].shape[1]))
self.frame = 0
self.frames_delete = 0
self.ds = ds
self.reset_stack = raw_dicom_stack
# Plot the first slice of the image
self.fig, self.ax = plt.subplots()
self.im = self.ax.imshow(np.array(raw_dicom_stack[0]), cmap = cm.gray)
def update(self, value):
self.frame = int(self.nframes * value)
# Update the image data
dat = np.array(self.raw_dicom_stack[self.frame])
self.im.set_data(dat)
self.im.set_clim([dat.min(), dat.max()])
# Redraw the plot
self.fig.canvas.draw()
def add(self,event):
if self.query_add == True:
print "Image series has already summed"
else:
for x in range(self.nframes):
self.summed_image += self.raw_dicom_stack[x]
del self.raw_dicom_stack[:]
self.raw_dicom_stack.append(self.summed_image[:,:])
self.nframes = len(self.raw_dicom_stack)
# Update the image data
dat = np.array(self.raw_dicom_stack[0])
self.im.set_data(dat)
self.slider.reset()
self.im.set_clim([dat.min(), dat.max()])
# Redraw the plot
self.fig.canvas.draw()
self.query_add = True
#rezero summed_image
#summed_image = np.zeros((self.summed_image.shape[0],self.summed_image.shape[1]))
def show(self):
plt.show()
def reset(self, event):
#self.summed_image = np.zeros((self.raw_dicom_stack[0].shape[0],self.raw_dicom_stack[0].shape[1]))
#self.query_register = False
#self.query_add = False
self.__init__(self.reset_stack, self.nframes, self.ds)
def quit(self,event):
quit()
if __name__ == '__main__':
viewer = VolumeViewer(np.zeros((2,10,10)), 4, None)
viewer.reset()