Abstract
Monte Carlo algorithms have a growing impact on nuclear medicine reconstruction processes. One of
the main limitations of myocardial perfusion imaging (MPI) is the effective mitigation of the scattering component,
which is particularly challenging in Single Photon Emission Computed Tomography (SPECT). In SPECT,
no timing information can be retrieved to locate the primary source photons. Monte Carlo methods allow an
event-by-event simulation of the scattering kinematics, which can be incorporated into a model of the imaging
system response. This approach was adopted in the late Nineties by several authors, and recently took advantage
of the increased computational power made available by high-performance CPUs and GPUs. These recent developments
enable a fast image reconstruction with improved image quality, compared to deterministic approaches.
Deterministic approaches are based on energy-windowing of the detector response, and on the cumulative estimate
and subtraction of the scattering component. In this paper, we review the main strategies and algorithms to
correct the scattering effect in SPECT and focus on Monte Carlo developments, which nowadays allow the threedimensional
reconstruction of SPECT cardiac images in a few seconds.
Keywords:
Monte carlo algorithms, SPECT, MPI, cardiac imaging reconstruction, scattering correction.
[1]
U. D. of Health and H. Services. Summary Health Statistics Tables for US Adults: National Health Interview Survey. 2018.
[4]
J. B. (Ed). Mcnp-a general monte carlo n-particle transport code, version 4c. 2000.
[6]
Battistoni G, Cerutti F, Fassò A, et al. The fluka code: Description and benchmarking. 2007; 896: 31-49.
[8]
Bielajew A, Rogers D. Presta: The parameter reduced electron-step transport algorithm for electron monte carlo transport. Nuclear Inst Methods Phys Res 1986; 18(1-6)
[13]
Jong HD. Accelerated Monte Carlo simulation for scatter correction in SPECT 2001. Eur J Nucl Med Mol Imaging 2004; 31(8): 1173-81.
[14]
Ramachandran GN, Lakshminarayanan AV. Three-Dimensional Reconstruction from Radiographs and Electron Micrographs. Proc Natl Acad Sci USA 68(9): 2236-40.
[20]
Floyd CE, Jaszczak RJ, Greer KL, Coleman RE. Deconvolution of Compton Scatter in SPECT. J Nucl Med 1985; 26(4): 403-8.
[21]
King SJGMA, Coleman M, Penney BC. Active quantitation in SPECT: A study of prereconstruction Metz filtering and use of the scatter degradation factor. Med Phys 1991; 18(2): 184-9.
[29]
Dempster AP, Laird NM, Rubin DB. Maximum Likelihood from Incomplete Data Via the EM Algorithm. 1977; Vol. 39.
[30]
Floyd CE, Jaszczak RJ. Coleman, and R. Edward. “Maximum Likelihood Reconstruction for SPECT with Monte Carlo Modeling. Asymptotic Behavior 1987; 34(1): 285-7.
[36]
De Wit TC, Xiao J, Beekman FJ. Monte carlobased statistical SPECT reconstruction: Influence of number of photon tracks. IEEE Transactions on Nuclear Science 2005; 52(5): 1365-9.