Computational Analysis and Improvement of SIRT

Iterative X-ray computed tomography (CT) algorithms have the potential for producing high-quality images but are computationally very demanding, especially when applied to high-resolution problems. Focusing on simultaneous iterative reconstruction technique (SIRT), we provide an eigenvalue based scheme for automatically determining a near-optimal value of the relaxation parameter. This accelerates the convergence rate of SIRT to the point where only half the number of iterations normally required is needed. We also modify the way SIRT uses preconditioning to solve a weighted least squares problem. The resulting algorithm, which we call PSIRT, is associated with a smaller memory footprint and calls for less data to be communicated in a distributed-memory implementation. Experimental residual norm and timing results are provided based on cone-beam micro-CT mouse data, including for an ordered subsets study.

[1]  Jens Gregor,et al.  Distributed iterative image reconstruction for micro-CT with ordered-subsets and focus of attention problem reduction , 2004 .

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

[3]  T Nielsen,et al.  Cardiac cone-beam CT volume reconstruction using ART. , 2005, Medical physics.

[4]  H. V. D. Vorst,et al.  SIRT- and CG-type methods for the iterative solution of sparse linear least-squares problems , 1990 .

[5]  Patrick Dupont,et al.  Reduction of metal streak artifacts in X-ray computed tomography using a transmission maximum a posteriori algorithm , 1999 .

[6]  Avinash C. Kak,et al.  Principles of computerized tomographic imaging , 2001, Classics in applied mathematics.

[7]  Jens Gregor,et al.  Iterative reconstruction of cone-beam CT data on a cluster , 2007, Electronic Imaging.

[8]  George Bosilca,et al.  Open MPI: Goals, Concept, and Design of a Next Generation MPI Implementation , 2004, PVM/MPI.

[9]  F J Beekman,et al.  Parallel statistical image reconstruction for cone-beam x-ray CT on a shared memory computation platform. , 2005 .

[10]  J. R. BILBAO-CASTRO,et al.  Performance of Parallel 3 D Iterative Reconstruction Algorithms , 2004 .

[11]  Charles R. Johnson,et al.  Matrix analysis , 1985, Statistical Inference for Engineers and Data Scientists.

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

[13]  Lubomir M. Hadjiiski,et al.  A comparative study of limited-angle cone-beam reconstruction methods for breast tomosynthesis. , 2006, Medical physics.

[14]  Jeffrey A. Fessler,et al.  Statistical image reconstruction for polyenergetic X-ray computed tomography , 2002, IEEE Transactions on Medical Imaging.

[15]  P. Gilbert.  Iterative methods for the three-dimensional reconstruction of an object from projections. , 1972, Journal of theoretical biology.

[16]  R. J. Gaudette,et al.  A comparison study of linear reconstruction techniques for diffuse optical tomographic imaging of absorption coefficient. , 2000, Physics in medicine and biology.

[17]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[18]  M. Knaup,et al.  Statistical Cone-Beam CT Image Reconstruction using the Cell Broadband Engine , 2006, 2006 IEEE Nuclear Science Symposium Conference Record.

[19]  Y. Saad,et al.  Iterative solution of linear systems in the 20th century , 2000 .

[20]  Jean-Baptiste Thibault,et al.  A three-dimensional statistical approach to improved image quality for multislice helical CT. , 2007, Medical physics.

[21]  Tinsu Pan,et al.  Acceleration of Landweber-type algorithms by suppression of projection on the maximum singular vector , 1991 .

[22]  J Gregor,et al.  Three-dimensional focus of attention for iterative cone-beam micro-CT reconstruction. , 2006, Physics in medicine and biology.

[23]  L. Feldkamp,et al.  Practical cone-beam algorithm , 1984 .

[24]  M. Defrise,et al.  Iterative reconstruction for helical CT: a simulation study. , 1998, Physics in medicine and biology.

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