Parallel iterative cone beam CT image reconstruction on a PC cluster

The iterative reconstruction (IR) algorithms are able to generate CT images with higher quality compared with the filtered backprojection method when the projection data is noisy or the radiation dose is low. IR can also be used when the data is incomplete. The major disadvantage of the IR is its high demand on computation and slow reconstruction. To improve the time performance of the IR we parallelized four representative iterative algorithms: EM, SART and their ordered subset (OS) versions on a Linux PC cluster. A micro-forward-back-projector was implemented for increased parallelization. Parameters are cached at the ray level during forward-projection to further reduce time used for back-projection. The speed-up and efficiency factors of our parallel implementation are reported.

[1]  Brian F. Hutton,et al.  A clinical perspective of accelerated statistical reconstruction , 1997, European Journal of Nuclear Medicine.

[2]  M I Miller,et al.  Maximum likelihood SPECT in clinical computation times using mesh-connected parallel computers. , 1991, IEEE transactions on medical imaging.

[3]  M S Atkins,et al.  Use Of Transputers In A 3-d Positron Emission Tomograph , 1990, 1990 IEEE Nuclear Science Symposium Conference Record.

[4]  Giulia Spaletta,et al.  An image reconstruction algorithm in tomography: A version for the CRAY X-MP vector computer , 1989, Comput. Graph..

[5]  L. Shepp,et al.  Maximum Likelihood Reconstruction for Emission Tomography , 1983, IEEE Transactions on Medical Imaging.

[6]  M. S. Atkins,et al.  Use of transputers in a 3-D positron emission tomograph , 1991 .

[7]  A. Kak,et al.  Simultaneous Algebraic Reconstruction Technique (SART): A Superior Implementation of the Art Algorithm , 1984, Ultrasonic imaging.

[8]  Klaus Mueller,et al.  Rapid 3-D cone-beam reconstruction with the simultaneous algebraic reconstruction technique (SART) using 2-D texture mapping hardware , 2000, IEEE Transactions on Medical Imaging.

[9]  William Gropp,et al.  User''s Guide for mpich , 1996 .

[10]  Xiang Li,et al.  A numerical simulator in VC++ on PC for iterative image reconstruction , 2003 .

[11]  Zang-Hee Cho,et al.  A parallel implementation of 3-D CT image reconstruction on hypercube multiprocessor , 1990 .

[12]  H. Malcolm Hudson,et al.  Accelerated image reconstruction using ordered subsets of projection data , 1994, IEEE Trans. Medical Imaging.

[13]  M. King,et al.  Optimization of regularization of attenuation and scatter corrected /sup 99 m/Tc cardiac SPECT studies for defect detection using hybrid images , 2000, 2000 IEEE Nuclear Science Symposium. Conference Record (Cat. No.00CH37149).

[14]  Joseph A. O'Sullivan,et al.  Iterative deblurring for CT metal artifact reduction , 1996, IEEE Trans. Medical Imaging.

[15]  Alessandro Bevilacqua.  A Dynamic Load Balancing Method On A Heterogeneous Cluster Of Workstations , 1999, Informatica.

[16]  N Popa,et al.  Maximum-likelihood reconstruction with ordered subsets in bone SPECT. , 1999, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[17]  K. Lange,et al.  EM reconstruction algorithms for emission and transmission tomography. , 1984, Journal of computer assisted tomography.

[18]  Anthony Skjellum,et al.  A High-Performance, Portable Implementation of the MPI Message Passing Interface Standard , 1996, Parallel Comput..

[19]  Xiang Li,et al.  P2P-enhanced Distributed Computing in EM Medical Image Reconstruction , 2004, PDPTA.

[20]  S. Manglos,et al.  Transmission maximum-likelihood reconstruction with ordered subsets for cone beam CT. , 1995, Physics in medicine and biology.

[21]  Freek J. Beekman,et al.  Accelerated iterative transmission CT reconstruction using an ordered subsets convex algorithm , 1998, IEEE Transactions on Medical Imaging.

[22]  R. Leahy,et al.  Recent developments in iterative image reconstruction for PET and SPECT [Editorial] , 2000, IEEE Transactions on Medical Imaging.

[23]  K. Mueller,et al.  Fast implementations of algebraic methods for three-dimensional reconstruction from cone-beam data , 1999, IEEE Transactions on Medical Imaging.

[24]  M. Casey,et al.  Design of a super fast three-dimensional projection system for positron emission tomography , 1990 .

[25]  R M Leahy,et al.  Internet2-based 3D PET image reconstruction using a PC cluster. , 2002, Physics in medicine and biology.

[26]  Ming Jiang,et al.  Convergence of the simultaneous algebraic reconstruction technique (SART) , 2003, Conference Record of Thirty-Fifth Asilomar Conference on Signals, Systems and Computers (Cat.No.01CH37256).

[27]  Ariela Sofer,et al.  A data-parallel algorithm for iterative tomographic image reconstruction , 1999, Proceedings. Frontiers '99. Seventh Symposium on the Frontiers of Massively Parallel Computation.

[28]  Michael I. Miller,et al.  3-D maximum a posteriori estimation for single photon emission computed tomography on massively-parallel computers , 1993, IEEE Trans. Medical Imaging.

[29]  Albert Macovski,et al.  A maximum likelihood approach to image reconstruction , 1976 .

[30]  Max A. Viergever,et al.  Dual matrix ordered subsets reconstruction for accelerated 3D scatter compensation in single-photon emission tomography , 1997, European Journal of Nuclear Medicine.

[31]  David Jones.  High performance , 1989, Nature.

[32]  Siegfried Benkner,et al.  Web-Based Parallel ML-EM Reconstruction for SPECT on SMP Clusters , 2001 .

[33]  Alessandro Bevilacqua.  Evaluation of a Fully 3-D Bpf Method for Small Animal PET Images on Mimd Architectures , 1999 .

[34]  Ryan J. Murphy,et al.  Deblurring subject to nonnegativity constraints when known functions are present with application to object-constrained computerized tomography , 2001, IEEE Transactions on Medical Imaging.

[35]  M. Defrise,et al.  HeinzelCluster: accelerated reconstruction for FORE and OSEM3D , 2001, 2001 IEEE Nuclear Science Symposium Conference Record (Cat. No.01CH37310).

[36]  SkjellumAnthony,et al.  A high-performance, portable implementation of the MPI message passing interface standard , 1996 .

[37]  Paul Kinahan,et al.  Postinjection single photon transmission tomography with ordered-subset algorithms for whole-body PET imaging , 2002 .

[38]  A. H. Andersen.  Algebraic reconstruction in CT from limited views. , 1989, IEEE transactions on medical imaging.

[39]  José Miguel Fernández-Güell.  Madrid , 2014, Juan Perón.