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.

A revised geometric representative model of the lower part of the colon, including the rectum, the urinary bladder and prostate, is proposed for use in the calculation of absorbed dose from injected radiopharmaceuticals. The lower segment of the sigmoid colon as described in the 1987 Oak Ridge National Laboratory mathematical phantoms does not accurately represent the combined urinary bladder/rectal/prostate geometry. In the revised model in this study, the lower part of the abdomen includes an explicitly defined rectum. The shape of sigmoid colon is more anatomically structured, and the diameters of the descending colon are modified to better approximate their true anatomic dimensions. To avoid organ overlap and for more accurate representation of the urinary bladder and the prostate gland (in the male), these organs are shifted from their originally defined positions. The insertion of the rectum and the shifting of the urinary bladder will not overlap with or displace the female phantom's ovaries or the uterus. In the adult male phantom, the prostatic urethra and seminal duct are also included explicitly in the model. The relevant structures are defined for the newborn and 1-, 5-, 10- and 15-y-old (adult female) and adult male phantoms. METHODS: Values of the specific absorbed fractions and radionuclide S values were calculated with the SIMDOS dosimetry package. Results for 99mTc and other radionuclides are compared with previously reported values. RESULTS: The new model was used to calculate S values that may be crucial to calculations of the effective dose equivalent. For 131I, the S (prostate<--urinary bladder contents) and S (lower large intestine [LLI] wall<--urinary bladder contents) are 6.7 x 10(-6) and 3.41 x 10(-6) mGy/MBq x s, respectively. Corresponding values given by the MIRDOSE3 computer program are 6.23 x 10(-6) and 1.53 x 10(-6) mGy/MBq x s, respectively. The value of S (rectum wall<--urinary bladder contents) is 4.84 x 10(-5) mGy/MBq x s. For 99mTc, we report S (testes<--prostate) and S (LLI wall<--prostate) of 9.41 x 10(-7) and 1.53 x 10(-7) mGy/MBq x s versus 1.33 x 10(-6) and 7.57 x 10(-6) mGy/MBq x s given by MIRDOSE3, respectively. The value of S (rectum wall<--prostate) for 99mTc is given as 4.05 x 10(-6) mGy/MBq x s in the present model. CONCLUSION: The new revised rectal model describes an anatomically realistic lower abdomen region, thus giving improved estimates of absorbed dose. Due to shifting the prostate gland, a 30%-45% reduction in the testes dose and the insertion of the rectum leads to 48%-55% increase in the LLI wall dose when the prostate is the source organ.