Performance comparison of projector-backprojector pairs for iterative tomographic reconstruction

Iterative reconstruction methods, such as the expectation maximization (EM) algorithm and its extended approaches, have played a prominent role in emission computed tomography due to their remarkable advantages over the conventional filtered backprojection method. However, since iterative reconstructions typically are comprised of repeatedly projecting and backprojecting the data, the computational load required for reconstructing an image highly depends on the performance of the projector-backprojector pair used in the algorithm. In this work we compare quantitative performance of representative methods for implementing projector-backprojector pairs-ray-tracing methods, rotation-based methods, and pixel-driven methods. To reduce the overall cost for the projection-backprojection operations for each method, we investigate how previously computed results can be reused so that the number of redundant calculations can be minimized. Our experimental results demonstrate that, while the rotation based methods can be useful for simplifying the correction of important physical factors, the computational cost to achieve good accuracy is considerably higher than that of the ray-tracing methods.

[1]  R. Siddon Fast calculation of the exact radiological path for a three-dimensional CT array. , 1985, Medical physics.

[2]  G. Gullberg,et al.  A three-dimensional ray-driven attenuation, scatter and geometric response correction technique for SPECT in inhomogeneous media. , 2000, Physics in medicine and biology.

[3]  R L Siddon,et al.  Calculation of the radiological depth. , 1985, Medical physics.

[4]  R. W. Schafer,et al.  A comparison of rotation-based methods for iterative reconstruction algorithms , 1995 .

[5]  Alan W. Paeth,et al.  A fast algorithm for general raster rotation , 1986 .

[6]  D P Harrington,et al.  Quantitative reconstruction for myocardial perfusion SPECT: an efficient approach by depth-dependent deconvolution and matrix rotation. , 1994, Physics in medicine and biology.

[7]  Jeffrey A. Fessler,et al.  Grouped-coordinate ascent algorithms for penalized-likelihood transmission image reconstruction , 1997, IEEE Transactions on Medical Imaging.

[8]  G.T. Gullberg,et al.  Frequency domain implementation of the three-dimensional geometric point response correction in SPECT imaging , 1991, Conference Record of the 1991 IEEE Nuclear Science Symposium and Medical Imaging Conference.

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

[10]  Alvaro R. De Pierro,et al.  A modified expectation maximization algorithm for penalized likelihood estimation in emission tomography , 1995, IEEE Trans. Medical Imaging.

[11]  T. Turkington,et al.  Simultaneous compensation for attenuation, scatter and detector response for SPECT reconstruction in three dimensions. , 1992, Physics in medicine and biology.

[12]  Per-Erik Danielsson,et al.  High-accuracy rotation of images , 1992, CVGIP Graph. Model. Image Process..

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