reconstruction_3d.py

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

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

array = ds.pixel_array

recon_3d = numpy.zeros((64,64,64))

for i in range(0,63):
    recon_3d[i,:,:] = iradon(numpy.transpose(array[:,i,:]), output_size=64)


for i in range(26,39): #coronary view
    pylab.imshow(recon_3d[:,i,:])
    pylab.show()