TY - JOUR
T1 - Combination of MCNP and SimSET for Monte Carlo simulation of SPECT with medium- and high-energy photons
AU - Du, Y.
AU - Frey, E. C.
AU - Wang, W. T.
AU - Tocharoenchai, C.
AU - Baird, W. H.
AU - Tsui, B. M.W.
N1 - Funding Information:
Manuscript received November 25, 2001; revised February 24, 2002. This work was supported in part by the National Institutes of Health under Grants R01-HL61616 and R01-CA39463.
PY - 2002/6
Y1 - 2002/6
N2 - Monte Carlo (MC) simulation has become an important tool for characterizing and evaluating imaging systems and reconstruction algorithms. In this work, we combined two MC codes-Monte Carlo N-particle transport code system (MCNP) and simulation system for emission tomography (SimSET)-using the advantages of each to produce a simulation tool that is efficient enough to generate single photon emission computed tomography (SPECT) data. The new SimSET-MCNP method allows the use of voxelized three-dimensional (3-D) phantoms and models photon propagation inside collimator-detector systems. This combination provides a tool for evaluating compensation methods applied to imaging of agents where medium- and high-energy photons are important. To validate the new tool, we compared simulated projections of a sphere containing I-123 water solution embedded in a cylindrical water-filled phantom with experimentally measured projections and simulations using SimSET and MCNP alone. Using these data, we compared profiles through the projection data, energy spectra, and relative number of photons in four projection views. The agreement with experiment was good, with disagreements of the order of a few percent. In addition, for the MC simulations, we classified detected photons based on whether they scattered in the phantom, whether they passed through the collimator holes, penetrated the septa, or scattered in the collimator, and whether they resulted from 159-keV or high-energy photons. For all these classes of photons, there was excellent agreement between SimSET-MCNP and MCNP. Finally, we evaluated the new combination in terms of simulation time and found it significantly more efficient than MCNP alone. We conclude that the new simulation tool works and allows the generation of SPECT data using voxelized phantoms for cases when medium- and high-energy photons are important.
AB - Monte Carlo (MC) simulation has become an important tool for characterizing and evaluating imaging systems and reconstruction algorithms. In this work, we combined two MC codes-Monte Carlo N-particle transport code system (MCNP) and simulation system for emission tomography (SimSET)-using the advantages of each to produce a simulation tool that is efficient enough to generate single photon emission computed tomography (SPECT) data. The new SimSET-MCNP method allows the use of voxelized three-dimensional (3-D) phantoms and models photon propagation inside collimator-detector systems. This combination provides a tool for evaluating compensation methods applied to imaging of agents where medium- and high-energy photons are important. To validate the new tool, we compared simulated projections of a sphere containing I-123 water solution embedded in a cylindrical water-filled phantom with experimentally measured projections and simulations using SimSET and MCNP alone. Using these data, we compared profiles through the projection data, energy spectra, and relative number of photons in four projection views. The agreement with experiment was good, with disagreements of the order of a few percent. In addition, for the MC simulations, we classified detected photons based on whether they scattered in the phantom, whether they passed through the collimator holes, penetrated the septa, or scattered in the collimator, and whether they resulted from 159-keV or high-energy photons. For all these classes of photons, there was excellent agreement between SimSET-MCNP and MCNP. Finally, we evaluated the new combination in terms of simulation time and found it significantly more efficient than MCNP alone. We conclude that the new simulation tool works and allows the generation of SPECT data using voxelized phantoms for cases when medium- and high-energy photons are important.
KW - Collimator scatter
KW - Monte Carlo (MC) simulation
KW - Septal penetration
KW - Single photon emission computed tomography (SPECT)
UR - http://www.scopus.com/inward/record.url?scp=0036624601&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=0036624601&partnerID=8YFLogxK
U2 - 10.1109/TNS.2002.1039547
DO - 10.1109/TNS.2002.1039547
M3 - Article
AN - SCOPUS:0036624601
SN - 0018-9499
VL - 49 I
SP - 668
EP - 674
JO - IEEE Transactions on Nuclear Science
JF - IEEE Transactions on Nuclear Science
IS - 3
ER -