Approximate 3D iterative reconstruction for SPECT.

Compared with slice-by-slice approaches for SPECT reconstruction, three-dimensional iterative methods provide a more accurate physical model and an improved SPECT image. Clinical application of these methods, however, is limited primarily to their computational demands. This paper investigates the methods for approximate 3D iterative reconstruction that greatly reduce this demand by excluding from the reconstruction the smaller magnitude elements of the system matrix. A new method is described which is designed to control the resulting bias in the SPECT image for a given reduction in computation. The approximate methods were compared to fully 3D iterative reconstruction in terms of SPECT image bias and visual quality. All methods were incorporated into the ML-EM algorithm and applied to data from 3D mathematical and experimental brain phantoms. The SPECT images reconstructed by the approximate methods exhibited a positive bias throughout the image that was in general smaller with the new method (in the rage of 2%-6%). The bias was smallest in locally hot regions and largest in locally cold regions. The high quality brain phantom images demonstrated the capability of the new method in realistic imaging contexts. The time per iteration for an entire 3D brain phantom on a modern workstation using the approximate 3D method was 7.0 s.

[1]  Ronald J. Jaszczak,et al.  Bayesian reconstruction and use of anatomical a priori information for emission tomography , 1996, IEEE Trans. Medical Imaging.

[2]  E. Hoffman,et al.  3-D phantom to simulate cerebral blood flow and metabolic images for PET , 1990 .

[3]  T A Riauka,et al.  Photon propagation and detection in single-photon emission computed tomography--an analytical approach. , 1994, Medical physics.

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

[5]  R. Jaszczak,et al.  Inverse Monte Carlo: A Unified Reconstruction Algorithm for SPECT , 1985, IEEE Transactions on Nuclear Science.

[6]  E. Frey,et al.  A practical method for incorporating scatter in a projector-backprojector for accurate scatter compensation in SPECT , 1993 .

[7]  N. Pelc,et al.  An attenuated projector-backprojector for iterative SPECT reconstruction. , 1985, Physics in medicine and biology.

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

[9]  D. Gilland,et al.  Implementation of simultaneous attenuation and detector response correction in SPECT , 1988 .

[10]  T. Budinger,et al.  Three-dimensional reconstruction in nuclear medicine emission imaging , 1974 .

[11]  Lee-Tzuu Chang,et al.  A Method for Attenuation Correction in Radionuclide Computed Tomography , 1978, IEEE Transactions on Nuclear Science.

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

[13]  Max A. Viergever,et al.  Object shape dependent PSF model for SPECT imaging , 1993 .

[14]  B. C. Penney,et al.  A projector, backprojector pair which accounts for the two-dimensional depth and distance dependent blurring in SPECT , 1990 .

[15]  James E. Bowsher,et al.  Out-of-plane photons in SPECT , 1990 .

[16]  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.

[17]  T G Turkington,et al.  A 3D model of non-uniform attenuation and detector response for efficient iterative reconstruction in SPECT. , 1994, Physics in medicine and biology.