Evaluation of 3D Monte Carlo-based scatter correction for 99mTc cardiac perfusion SPECT.

201Tl-Chloride (201Tl) is a myocardial perfusion SPECT agent with excellent biochemical properties commonly used for assessing tissue viability. However, cardiac 201Tl SPECT images are severely degraded by photons scattered in the thorax. Accurate correction for this scatter is complicated by the nonuniform density and varied sizes of thoraxes, by the additional attenuation and scatter caused by female patients' breasts, and by the energy spectrum of 201Tl. Monte Carlo simulation is a general and accurate method well suited to modeling this scatter. Methods: Statistical reconstruction that includes Monte Carlo modeling of scatter was compared with statistical reconstruction algorithms not corrected for scatter. In the ADS method, corrections for attenuation, detector response, and scatter (Monte Carlo–based) were implemented simultaneously via the dual-matrix ordered-subset expectation maximization algorithm with a Monte Carlo simulator as part of the forward projector. The ADS method was compared with the A method (ordered-subset expectation maximization with attenuation correction) and with the AD method (a method like the A method but with detector response modeling added). A dual-head SPECT system equipped with two 153Gd scanning line sources was used for simultaneously acquiring transmission and emission data. Four clinically realistic phantom configurations (a large thorax and a small thorax, each with and without breasts) with a cardiac insert containing 2 cold defects were used to evaluate the proposed reconstruction algorithms. We compared the performance of the different algorithms in terms of noise properties, contrast-to-noise ratios, the contrast separability of perfusion defects, uniformity, and robustness to anatomic variations. Results: The ADS method provided images with clearly better visual defect contrast than did the other methods. The contrasts achieved with the ADS method were 10%–24% higher than those achieved with the AD method and 11%–37% higher than those achieved with the A method. For a typical contrast level, the ADS method exhibited noise levels around 27% lower than the AD method and 34% lower than the A method. Compared with the other 2 algorithms, the ADS reconstructions were less sensitive to anatomic variations and had better image uniformity in the homogeneously perfused myocardium. Finally, we found that the improvements that can be achieved with Monte Carlo–based scatter correction are stronger for 201Tl than for 99mTc imaging. Conclusion: Our results indicate that Monte Carlo–based scatter correction is suitable for 201Tl cardiac imaging and that such correction simultaneously improves several image-quality metrics.

[1]  Max A. Viergever,et al.  Improved SPECT quantitation using fully three-dimensional iterative spatially variant scatter response compensation , 1996, IEEE Trans. Medical Imaging.

[2]  Michael A. King,et al.  Comparison of scatter compensation strategies for myocardial perfusion imaging using Tc-99m labeled sestamibi , 2001 .

[3]  G Germano,et al.  Quantitative LVEF and qualitative regional function from gated thallium-201 perfusion SPECT. , 1997, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[4]  Max A. Viergever,et al.  Evaluation of OS-EM vs. ML-EM for 1D, 2D and fully 3D SPECT reconstruction , 1996 .

[5]  S. Bergmann,et al.  Assessment of myocardial viability. , 2005, Seminars in nuclear medicine.

[6]  S. Underwood,et al.  Myocardial perfusion scintigraphy in the UK: insights from the British Nuclear Cardiology Society Survey 2000 , 2005, Heart.

[7]  S. Walrand,et al.  Quantitation in SPECT using an effective model of the scattering. , 1994, Physics in medicine and biology.

[8]  H H Barrett,et al.  Objective assessment of image quality: effects of quantum noise and object variability. , 1990, Journal of the Optical Society of America. A, Optics and image science.

[9]  Bruce H. Hasegawa,et al.  A Comparison of 180° and 360° Acquisition for Attenuation-Compensated Thallium-201 SPECT Images , 1998 .

[10]  Daniel S. Berman,et al.  Cardiac SPECT imaging , 1994, Cardiovascular Drugs and Therapy.

[11]  G Germano,et al.  A new algorithm for the quantitation of myocardial perfusion SPECT. I: technical principles and reproducibility. , 2000, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

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

[13]  Wojciech Zbijewski,et al.  Evaluation of 3D Monte Carlo–Based Scatter Correction for 201Tl Cardiac Perfusion SPECT , 2007, Journal of Nuclear Medicine.

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

[15]  C E Floyd,et al.  Experimentally measured scatter fractions and energy spectra as a test of Monte Carlo simulations. , 1987, Physics in medicine and biology.

[16]  Eric C. Frey,et al.  Application of task-based measures of image quality to optimization and evaluation of three-dimensional reconstruction-based compensation methods in myocardial perfusion SPECT , 2002, IEEE Transactions on Medical Imaging.

[17]  Brian F. Hutton,et al.  Transmission-based scatter correction of 180° myocardial single-photon emission tomographic studies , 2005, European Journal of Nuclear Medicine.

[18]  J. Fessler,et al.  Simultaneous transmission/emission myocardial perfusion tomography. Diagnostic accuracy of attenuation-corrected 99mTc-sestamibi single-photon emission computed tomography. , 1996, Circulation.

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

[20]  Eric C. Frey,et al.  A fast projector-backprojector pair modeling the asymmetric, spatially varying scatter response function for scatter compensation in SPECT imaging , 1993 .

[21]  H.W.A.M. de Jong,et al.  Acceleration of Monte Carlo SPECT simulation using convolution-based forced detection , 1999, 1999 IEEE Nuclear Science Symposium. Conference Record. 1999 Nuclear Science Symposium and Medical Imaging Conference (Cat. No.99CH37019).

[22]  Brad Kemp,et al.  A multicenter evaluation of commercial attenuation compensation techniques in cardiac SPECT using phantom models. , 2002, Journal of nuclear cardiology : official publication of the American Society of Nuclear Cardiology.

[23]  W. Niessen,et al.  Selection of task-dependent diffusion filters for the post-processing of SPECT images. , 1998, Physics in medicine and biology.

[24]  Eric C. Frey,et al.  Application of reconstruction-based scatter compensation to thallium-201 SPECT: implementations for reduced reconstructed image noise , 1998, IEEE Transactions on Medical Imaging.

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

[26]  Michael A. King,et al.  AN INVESTIGATION OF THE FILTERING OF TEW SCATTER ESTIMATES USED TO COMPENSATE FOR SCATTER WITH ORDERED SUBSET RECONSTRUCTIONS , 1996 .

[27]  M. King,et al.  Human-observer receiver-operating-characteristic evaluation of attenuation, scatter, and resolution compensation strategies for (99m)Tc myocardial perfusion imaging. , 2003, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[28]  Michael A. King,et al.  A scatter correction method for Tl-201 images: A Monte Carlo investigation , 1992 .

[29]  M A King,et al.  A Monte Carlo investigation of artifacts caused by liver uptake in single-photon emission computed tomography perfusion imaging with technetium 99m-labeled agents , 1996, Journal of nuclear cardiology : official publication of the American Society of Nuclear Cardiology.

[30]  M A Viergever,et al.  SPECT scatter modelling in non-uniform attenuating objects. , 1997, Physics in medicine and biology.

[31]  B F Hutton,et al.  A scanning line source for simultaneous emission and transmission measurements in SPECT. , 1993, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[32]  F. A. Seiler,et al.  Numerical Recipes in C: The Art of Scientific Computing , 1989 .

[33]  William H. Press,et al.  The Art of Scientific Computing Second Edition , 1998 .

[34]  S R Meikle,et al.  Transmission-based scatter correction of 180 degrees myocardial single-photon emission tomographic studies. , 1996, European journal of nuclear medicine.

[35]  E C Frey,et al.  The importance and implementation of accurate 3D compensation methods for quantitative SPECT. , 1994, Physics in medicine and biology.

[36]  Max A. Viergever,et al.  Fast SPECT simulation including object shape dependent scatter , 1995, IEEE Trans. Medical Imaging.

[37]  G V Heller,et al.  Multicenter clinical trial to evaluate the efficacy of correction for photon attenuation and scatter in SPECT myocardial perfusion imaging. , 1999, Circulation.

[38]  D. Bailey,et al.  A transmission-dependent method for scatter correction in SPECT. , 1994, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[39]  E C Frey,et al.  Simultaneous technetium-99m/thallium-201 SPECT imaging with model-based compensation for cross-contaminating effects. , 1999, Physics in medicine and biology.

[40]  T R Miller,et al.  Clinically important characteristics of maximum-likelihood reconstruction. , 1992, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[41]  J. Terry,et al.  Three-dimensional iterative reconstruction algorithms with attenuation and geometric point response correction , 1990 .

[42]  K. Ogawa,et al.  A practical method for position-dependent Compton-scatter correction in single photon emission CT. , 1991, IEEE transactions on medical imaging.

[43]  Freek J. Beekman,et al.  Efficient fully 3-D iterative SPECT reconstruction with Monte Carlo-based scatter compensation , 2002, IEEE Transactions on Medical Imaging.

[44]  Michael A. King,et al.  Evaluation of scatter compensation strategies and their impact on human detection performance in Tc-99m myocardial perfusion imaging , 2002 .

[45]  E C Frey,et al.  Fast implementations of reconstruction-based scatter compensation in fully 3D SPECT image reconstruction. , 1998, Physics in medicine and biology.

[46]  E C Frey,et al.  Scatter compensation methods in 3D iterative SPECT reconstruction: a simulation study. , 1997, Physics in medicine and biology.

[47]  F. Beekman,et al.  Monte Carlo-based statistical SPECT reconstruction: influence of number of photon tracks , 2005, IEEE Transactions on Nuclear Science.

[48]  Michael A. King,et al.  Quantitative myocardial perfusion SPECT , 1998, Journal of nuclear cardiology : official publication of the American Society of Nuclear Cardiology.

[49]  The effect of truncation reduction in fan beam transmission for attenuation correction of cardiac SPECT , 1995 .

[50]  Freek J. Beekman,et al.  Acceleration of Monte Carlo SPECT simulation using convolution-based forced detection , 1999 .