title: Convolution-Based Forced Detection Monte Carlo Simulation Incorporating Septal Penetration Modeling.
abstract: In SPECT imaging, photon transport effects such as scatter, attenuation and septal penetration can negatively affect the quality of the reconstructed image and the accuracy of quantitation estimation. As such, it is useful to model these effects as carefully as possible during the image reconstruction process. Many of these effects can be included in Monte Carlo (MC) based image reconstruction using convolution-based forced detection (CFD). With CFD Monte Carlo (CFD-MC), often only the geometric response of the collimator is modeled, thereby making the assumption that the collimator materials are thick enough to completely absorb photons. However, in order to retain high collimator sensitivity and high spatial resolution, it is required that the septa be as thin as possible, thus resulting in a significant amount of septal penetration for high energy radionuclides. A method for modeling the effects of both collimator septal penetration and geometric response using ray tracing (RT) techniques has been performed and included into a CFD-MC program. Two look-up tables are pre-calculated based on the specific collimator parameters and radionuclides, and subsequently incorporated into the SIMIND MC program. One table consists of the cumulative septal thickness between any point on the collimator and the center location of the collimator. The other table presents the resultant collimator response for a point source at different distances from the collimator and for various energies. A series of RT simulations have been compared to experimental data for different radionuclides and collimators. Results of the RT technique matches experimental data of collimator response very well, producing correlation coefficients higher than 0.995. Reasonable values of the parameters in the lookup table and computation speed are discussed in order to achieve high accuracy while using minimal storage space for the look-up tables. In order to achieve noise-free projection images from MC, it is seen that the inclusion of the RT implementation for septal penetration increases the speed of the simulation by a factor of about 7,500 compared to the conventional SIMIND MC program.
authors.affiliation: Department of Electrical and Computer Engineering, McMaster University, Hamilton, ON L8S 4L8, Canada.
authors.forename: Shaoying
authors.initials: S
authors.lastname: Liu
authors.forename: Michael A
authors.initials: MA
authors.lastname: King
authors.forename: Aaron B
authors.initials: AB
authors.lastname: Brill
authors.forename: Michael G
authors.initials: MG
authors.lastname: Stabin
authors.forename: Troy H
authors.initials: TH
authors.lastname: Farncombe
references.citation: Phys Med Biol. 2005 Apr 21;50(8):1791-804
references.pmid: 15815096
references.citation: Phys Med Biol. 2007 Jun 7;52(11):3027-43
references.pmid: 17505087
references.citation: IEEE Trans Nucl Sci. 2008 Jan 1;55(1):560-567
references.pmid: 20811587
references.citation: Phys Med Biol. 1973 Nov;18(6):855-62
references.pmid: 4778448