Abstract
Monte Carlo simulation has been used to produce projections from a voxel-based brain phantom, simulating a 99mTc-HMPAO single photon emission computed tomography (SPECT) brain investigation. For comparison, projections free from the effects of attenuation and scattering were also simulated, giving ideal transaxial images after reconstruction. Three methods of attenuation correction were studied: (a) a pre-processing method, (b) a post-processing uniform method and (c) a post-processing non-uniform method using a density map. The accuracy of these methods was estimated by comparison of the reconstructed images with the ideal images using the normalized mean square error, NMSE, and quantitative values of the regional cerebral blood flow, rCBF. A minimum NMSE was achieved for the effective linear attenuation coefficient µeff = 0.07 (0.09) cm-1 for the uniformpre method, the effective mass attenuation coefficient µeff/ρ = 0.08 (0.10) cm2 g-1 for the uniformpost method and µeff/ρ = 0.12 (0.13) cm2 g-1 for the non-uniformpost method. Values in parentheses represent the case of dual-window scatter correction. The non-uniformpost method performed better, as measured by the NMSE, both with and without scatter correction. Furthermore, the non-uniformpost method gave, on average, more accurate rCBF values. Although the difference in rCBF accuracy was small between the various methods, the same method should be used for patient studies as for the reference material.