Parallelization of iterative reconstruction algorithms in multiple modalities

In this work we have parallelized the Maximum Likelihood Expectation-Maximization (MLEM) and Ordered Subset Expectation Maximization (OSEM) algorithms for improving efficiency of reconstructions of multiple pinholes SPECT, and cone-bean CT data. We implemented the parallelized versions of the algorithms on a General Purpose Graphic Processing Unit (GPGPU): 448 cores of a NVIDIA Tesla M2070 GPU with 6GB RAM per thread of computing. We compared their run times against those from the corresponding CPU implementations running on 8 cores CPU of an AMD Opteron 6128 with 32 GB RAM. We have further shown how an optimization of thread balancing can accelerate the speed of the GPU implementation.

[1]  P. Cochat,et al.  Et al , 2008, Archives de pediatrie : organe officiel de la Societe francaise de pediatrie.

[2]  Bipin Singh,et al.  Ultrafast multipinhole single photon emission computed tomography iterative reconstruction using CUDA , 2011, 2011 IEEE Nuclear Science Symposium Conference Record.

[3]  Robert M. Lewitt,et al.  Fourier correction for spatially variant collimator blurring in SPECT , 1995, IEEE Trans. Medical Imaging.

[4]  Tianyu Ma,et al.  Development of GPU based image reconstruction method for clinical SPECT , 2012, 2012 IEEE Nuclear Science Symposium and Medical Imaging Conference Record (NSS/MIC).

[5]  G Pratx,et al.  Fast List-Mode Reconstruction for Time-of-Flight PET Using Graphics Hardware , 2011, IEEE Transactions on Nuclear Science.