Physically realistic Monte Carlo simulation of source, collimator and tomographic data acquisition for emission computed tomography

A sophisticated simulation package has been developed permitting full tomographic acquisition of nuclear medicine data from physically realistic, non-uniform and asymmetric source and scattering objects. The simulation package is based on MCNP (Monte Carlo for neutron-photon transport), a Monte Carlo code developed at the Los Alamos Scientific Laboratory. The MCNP code has been extensively modified with features that allow direct manipulation of the geometric and physical parameters of the simulation experiment. Collimators of various types have been modelled realistically; collimator properties such as hole dimensions, number of holes and packing geometry can be defined. Other parameters such as energy windows and the number of tomographic views can be set, as can spatial and energy resolution. Any 3D photon source and scattering object can be specified with any degree of detail, allowing for extremely realistic patient models. A post-processing algorithm allows the display and analysis of the simulated images as they are being generated. Using three-dimensional, realistic patient or phantom geometries, this simulation package represents a valuable tool for examining many of the physical aspects of the SPECT imaging situation.

[1]  Briesmeister MCNP: a general Monte Carlo code for neutron and photon transport. Version 3A. Revision 2 , 1986 .

[2]  Glenn F. Knoll,et al.  Results of Scattering in Radioisotope Imaging , 1973 .

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

[4]  A. Todd-Pokropek Image Processing in Nuclear Medicine , 1980, IEEE Transactions on Nuclear Science.

[5]  Grant T. Gullberg,et al.  Users Manual: Donner Algorithms for Reconstruction Tomography , 1977 .

[6]  J. S. Hendricks,et al.  Monte Carlo photon benchmark problems , 1991 .

[7]  C E Floyd,et al.  Energy and spatial distribution of multiple order Compton scatter in SPECT: a Monte Carlo investigation. , 1984, Physics in medicine and biology.

[8]  Ronald J. Jaszczak,et al.  Revised scatter fraction results for SPECT , 1987 .

[9]  M. S. Rosenthal,et al.  TECHNICAL NOTE: Scattering in uniform media (for correction of nuclear medicine images) , 1990 .

[10]  H O ANGER,et al.  GAMMA-RAY-DETECTION EFFICIENCY AND IMAGE RESOLUTION IN SODIUM IODIDE. UCRL-11184. , 1964, UCRL [reports]. U.S. Atomic Energy Commission.

[11]  P N Goodwin,et al.  SPECT instrumentation: performance, lesion detection, and recent innovations. , 1987, Seminars in nuclear medicine.

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

[13]  Ronald J. Jaszczak,et al.  Analysis of SPECT including Scatter and Attenuation Using Sophisticated Monte Carlo Modeling Methods , 1982, IEEE Transactions on Nuclear Science.