Reconstruction of SPECT

reconstruction.py

import dicom
from skimage.transform import iradon
import pylab
import numpy

ds = dicom.read_file("img.dcm") # read a dicom of pre-reconstructed SPECT

array = ds.pixel_array

y36 = numpy.transpose(array[:,36,:]) #important!!

reconstruction_fbp = iradon(y36, output_size=64)

pylab.imshow(reconstruction_fbp)
pylab.show()