Sequential and simultaneous dual-isotope brain SPECT: comparison with PET for estimation and discrimination tasks in early Parkinson disease

Parkinson disease (PD) is the second most frequently occurring cerebral degenerative disease, after Alzheimer disease. Treatments are available, but their efficacy is diminished unless they are administered in the early stages. Therefore, early identification of PD is crucial. In addition to providing perfectly registered studies, simultaneous 99mTc/123I imaging makes possible the assessment of pre- and postsynaptic neurotransmission functions under identical physiological conditions, while doubling the number of counts for the same total imaging time. These advantages are limited, however, by cross talk between the two radionuclides due to the close emission energies of 99mTc (140 keV) and 123I (159 keV). PET, on the other hand, provides good temporal and spatial resolution and sensitivity but usually requires the use of a single radionuclide. In the present work, the authors compared brain PET with sequential and simultaneous dual-isotope SPECT for the task of estimating striatal activity concentration and striatal size for a normal brain and two stages of early PD. Realistic Monte Carlo simulations of a time-of-flight PET scanner and gamma cameras were performed while modeling all interactions in the brain, collimator (gamma camera) and crystal (detector block in PET), as well as population biological variability of pre- and postsynaptic uptake. For SPECT imaging, we considered two values of system energy resolution and scanners with two and three camera heads. The authors used the Cramer-Rao bound, as a surrogate for the best theoretical performance, to optimize the SPECT acquisition energy windows and objectively compare PET and SPECT. The authors determined the discrimination performance between 500 simulated subjects in every disease stage as measured by the area under the ROC curve (AUC). The discrimination accuracy between a normal subject and a subject in the prodromal disease stage was AUC = 0.924 with PET, compared to 0.863 and 0.831 with simultaneous and sequential SPECT, respectively. The significant improvement in performance obtained with simultaneous dual-isotope SPECT compared to sequential imaging (p = 0.019) was due primarily to the increased number of counts detected and resulted in comparable performance when performing simultaneous SPECT on a two-head camera with 9.2% energy resolution to that obtained with sequential SPECT on a three-head camera with 6.2% energy resolution.     

Performance of coincidence imaging with long-lived positron emitters as an alternative to dedicated PET and SPECT

An important application of quantitative imaging in nuclear medicine is the estimation of absorbed doses in radionuclide therapy. Depending on the radionuclide used for therapy, quantitative imaging of the kinetics of the therapeutic radiopharmaceutical could be done using planar imaging, SPECT or PET. Since many nuclear medicine departments have a gamma camera system that is also suitable for coincidence imaging, the performance of these systems with respect to quantitative imaging of PET isotopes that could be of use in radionuclide dosimetry is of interest. We investigated the performance of a gamma camera with coincidence imaging capabilities with 99mTc, 111In, 18F and 76Br and a dedicated PET system with 18F and 76Br, using a single standard set of phantom measurements. Here, 76Br was taken as a typical example of prompt gamma-emitting PET isotopes that are applicable in radionuclide therapy dosimetry such as 86Y and 124I. Image quality measurements show comparable image contrasts for 76Br coincidence imaging and 111In SPECT. Although the spatial resolution of coincidence imaging is better than single photon imaging, the contrast obtained with 76Br is not better than that with 99mTc or 111In because of the prompt gamma involved. Additional improvements are necessary to allow for quantitative coincidence imaging of long-lived, prompt gamma producing positron emitters.     

Monte Carlo simulations of compact gamma cameras based on avalanche photodiodes

Avalanche photodiodes (APDs), and in particular position-sensitive avalanche photodiodes (PSAPDs), are an attractive alternative to photomultiplier tubes (PMTs) for reading out scintillators for PET and SPECT. These solid-state devices offer high gain and quantum efficiency, and can potentially lead to more compact and robust imaging systems with improved spatial and energy resolution. In order to evaluate this performance improvement, we have conducted Monte Carlo simulations of gamma cameras based on avalanche photodiodes. Specifically, we investigated the relative merit of discrete and PSAPDs in a simple continuous crystal gamma camera. The simulated camera was composed of either a 4 x 4 array of four channels 8 x 8 mm2 PSAPDs or an 8 x 8 array of 4 x 4 mm2 discrete APDs. These configurations, requiring 64 channels readout each, were used to read the scintillation light from a 6 mm thick continuous CsI:Tl crystal covering the entire 3.6 x 3.6 cm2 photodiode array. The simulations, conducted with GEANT4, accounted for the optical properties of the materials, the noise characteristics of the photodiodes and the nonlinear charge division in PSAPDs. The performance of the simulated camera was evaluated in terms of spatial resolution, energy resolution and spatial uniformity at 99mTc (140 keV) and 125I ( approximately 30 keV) energies. Intrinsic spatial resolutions of 1.0 and 0.9 mm were obtained for the APD- and PSAPD-based cameras respectively for 99mTc, and corresponding values of 1.2 and 1.3 mm FWHM for 125I. The simulations yielded maximal energy resolutions of 7% and 23% for 99mTc and 125I, respectively. PSAPDs also provided better spatial uniformity than APDs in the simple system studied. These results suggest that APDs constitute an attractive technology especially suitable to build compact, small field of view gamma cameras dedicated, for example, to small animal or organ imaging.     

Image quality assessment of LaBr3-based whole-body 3D PET scanners: a Monte Carlo evaluation

The main thrust for this work is the investigation and design of a whole-body PET scanner based on new lanthanum bromide scintillators. We use Monte Carlo simulations to generate data for a 3D PET scanner based on LaBr3 detectors, and to assess the count-rate capability and the reconstructed image quality of phantoms with hot and cold spheres using contrast and noise parameters. Previously we have shown that LaBr3 has very high light output, excellent energy resolution and fast timing properties which can lead to the design of a time-of-flight (TOF) whole-body PET camera. The data presented here illustrate the performance of LaBr3 without the additional benefit of TOF information, although our intention is to develop a scanner with TOF measurement capability. The only drawbacks of LaBr3 are the lower stopping power and photo-fraction which affect both sensitivity and spatial resolution. However, in 3D PET imaging where energy resolution is very important for reducing scattered coincidences in the reconstructed image, the image quality attained in a non-TOF LaBr3 scanner can potentially equal or surpass that achieved with other high sensitivity scanners. Our results show that there is a gain in NEC arising from the reduced scatter and random fractions in a LaBr3 scanner. The reconstructed image resolution is slightly worse than a high-Z scintillator, but at increased count-rates, reduced pulse pileup leads to an image resolution similar to that of LSO. Image quality simulations predict reduced contrast for small hot spheres compared to an LSO scanner, but improved noise characteristics at similar clinical activity levels.     

Monte Carlo simulations of a scintillation camera using GATE: validation and application modelling

Geant4 application for tomographic emission (GATE) is a recently developed simulation platform based on Geant4, specifically designed for PET and SPECT studies. In this paper we present validation results of GATE based on the comparison of simulations against experimental data, acquired with a standard SPECT camera. The most important components of the scintillation camera were modelled. The photoelectric effect. Compton and Rayleigh scatter are included in the gamma transport process. Special attention was paid to the processes involved in the collimator: scatter, penetration and lead fluorescence. A LEHR and a MEGP collimator were modelled as closely as possible to their shape and dimensions. In the validation study, we compared the simulated and measured energy spectra of different isotopes: 99mTc, 22Na, 57Co and 67Ga. The sensitivity was evaluated by using sources at varying distances from the detector surface. Scatter component analysis was performed in different energy windows at different distances from the detector and for different attenuation geometries. Spatial resolution was evaluated using a 99mTc source at various distances. Overall results showed very good agreement between the acquisitions and the simulations. The clinical usefulness of GATE depends on its ability to use voxelized datasets. Therefore, a clinical extension was written so that digital patient data can be read in by the simulator as a source distribution or as an attenuating geometry. Following this validation we modelled two additional camera designs: the Beacon transmission device for attenuation correction and the Solstice scanner prototype with a rotating collimator. For the first setup a scatter analysis was performed and for the latter design. the simulated sensitivity results were compared against theoretical predictions. Both case studies demonstrated the flexibility and accuracy of GATE and exemplified its potential benefits in protocol optimization and in system design.     

Changes in the energy response of a dedicated gamma camera after exposure to a high-flux irradiation.

This work reports the effects of the gain variation of the photomultiplier tubes (PMTs) observed on a cardiac dedicated gamma camera after accidental high-flux irradiation. One detector of this dual-headed 90 degrees-fixed gamma camera was accidentally left uncollimated during a quality assurance procedure on the other detector with a 57Co flood source (259 MBq) and received a non-uniform high flux of 1.9-0.6 Mcps over 25000 mm2 areas for about 30 min. To evaluate the severity and the duration of the perturbation effect on the energy response of the detector, the photopeak position was monitored for about 1 month with a 99mTc point source. The 140 keV photopeak shifted to 158 keV soon after irradiation, reached the correct position after 9 days and moved to a stable value of 132 keV after 15 days. Afterwards, a new energy calibration reset the photopeak position at 140 keV and the correct energy response of the gamma camera. This experience suggests that particular care should be taken to avoid exposures to high radiation fluxes that induce persistent gain shifts on the PMTs of this system.     

Optimizing Compton camera geometries

Compton cameras promise to improve the characteristics of nuclear medicine imaging, wherein mechanical collimation is replaced with electronic collimation. This leads to huge gains in sensitivity and, consequently, a reduction in the radiation dosage that needs to be administered to the patient. Design modifications that improve the sensitivity invariably compromise resolution. The scope of the current project was to determine an optimal design and configuration of a Compton camera that strikes a balance between these two properties. Transport of the photon flux from the source to the detectors was simulated with the camera geometry serving as the parameter to be optimized. Two variations of the Boltzmann photon transport equation, with and without photon polarization, were employed to model the flux. Doppler broadening of the energy spectra was also included. The simulation was done in a Monte Carlo framework using GEANT4. Two clinically relevant energies, 140 keV and 511 keV, corresponding to 99mTc and 18F were simulated. The gain in the sensitivity for the Compton camera over the conventional camera was 100 fold. Neither Doppler broadening nor polarization had any significant effect on the sensitivity of the camera. However, the spatial resolution of the camera was affected by these processes. Doppler broadening had a deleterious effect on the spatial resolution, but polarization improved the resolution when accounted for in the reconstruction algorithm.     

The PETRRA positron camera: design, characterization and results of a physical evaluation

The PETRRA positron camera is a large-area (600 mm x 400 mm sensitive area) prototype system that has been developed through a collaboration between the Rutherford Appleton Laboratory and the Institute of Cancer Research/Royal Marsden Hospital. The camera uses novel technology involving the coupling of 10 mm thick barium fluoride scintillating crystals to multi-wire proportional chambers filled with a photosensitive gas. The performance of the camera is reported here and shows that the present system has a 3D spatial resolution of approximately 7.5 mm full-width-half-maximum (FWHM), a timing resolution of approximately 3.5 ns (FWHM), a total coincidence count-rate performance of at least 80-90 kcps and a randoms-corrected sensitivity of approximately 8-10 kcps kBq(-1) ml. For an average concentration of 3 kBq ml(-1) as expected in a patient it is shown that, for the present prototype, approximately 20% of the data would be true events. The count-rate performance is presently limited by the obsolete off-camera read-out electronics and computer system and the sensitivity by the use of thin (10 mm thick) crystals. The prototype camera has limited scatter rejection and no intrinsic shielding and is, therefore, susceptible to high levels of scatter and out-of-field activity when imaging patients. All these factors are being addressed to improve the performance of the camera. The large axial field-of-view of 400 mm makes the camera ideally suited to whole-body PET imaging. We present examples of preliminary clinical images taken with the prototype camera. Overall, the results show the potential for this alternative technology justifying further development.     

Improvement of the fluorine-18 fluorodeoxyglucose images in simultaneous F-18 FDG/Tc-99m collimated spect imaging.

Collimated F-18 FDG SPECT imaging has been shown to be an acceptable alternative to F-18 FDG PET imaging for evaluating injured but viable myocardium. Ultra-energy (UHE) imaging is usually performed in simultaneous F-18 FDG/Tc-99m MIBI studies. The main limitations of this technique are degradation of the Tc-99m MIBI images due to F-18 downscatter to the Tc-99m window, and loss of resolution in Tc-99m images caused by using a UHE rather than a low-energy collimator. The quality of F-18 images has not been addressed. In our clinical and phantom studies we have found that F-18 images are inferior to simultaneously acquired Tc-99m MIBI images. This paper compares two correction methods for F-18 FDG images in a realistic cardiac phantom study. One approach is subtractive scatter correction, which employs a third 410 keV energy window image to estimate scatter. The other approach is based on restoration. The phantom acquisition was performed with 7.2 MBq of F-18 and 22.2 MBq of Tc-99m injected into the left ventricular (LV) wall. Three inserts, 3 cm, 2 cm, and 1 cm in diameter, were placed in the LV wall to simulate infarcts. Circumferential profiles were drawn from three successive short-axis slices and compared with true phantom data. The differences were calculated as root-mean-square error (rmse). Scatter correction improved rmse only 4.5 +/- 0.3%, while restoration improved rmse 16.1 +/- 0.4%, when compared with raw data. The same differences, measured as rmse, were 9.5 +/- 0.5, 6.8 +/- 0.4, and 5.1 +/- 0.5 for raw, scatter corrected, and restored F-18 data, respectively, when compared with Tc-99m window 140 keV data. The amount of noise, measured as root-mean-square % (rms%) was 5.3 +/- 0.5% for the Tc-99m image, 4.9 +/- 0.7% for the F-18 restored image, 6.2 +/- 0.6% for the raw F-18 image, and 6.5 +/- 0.9% for the scatter corrected F-18 image. The contrast measured for 2 cm and 3 cm inserts was 0.17 +/- 0.07 and 0.26 +/- 0.06 for F-18 raw data, 0.19 +/- 0.08 and 0.29 +/- 0.06 for the scatter corrected F-18 image, and 0.28 +/- 0.06 and 0.43 +/- 0.07 for the restored F-18 image. The contrast was 0.20 +/- 0.07 and 0.46 +/- 0.05 for the Tc-99m 140 keV window image. The restoration approach provided F-18 images of better contrast and detectibility than uncorrected or scatter corrected F-18 images. Restored F-18 images match better with the simultaneously acquired Tc-99m images.     

Comparison of correction techniques for simultaneous 201Tl/99mTc myocardial perfusion SPECT imaging: a dog study

We compared two correction methods for simultaneous 201Tl/99mTc dual-isotope single-photon emission computed tomography (SPECT). Both approaches use the information from the third energy window placed between the photopeak windows of the 201Tl and 99mTc. The first approach, described by Moore et al, corrects only for the contribution of the 99mTc to the 201Tl primary 70 keV window. We developed the three-window transformation dual isotope correction method, which is a simultaneous cross-talk correction. The two correction methods were compared in a simultaneous 201Tl/99mTc sestamibi cardiac dog study. Three separate acquisitions were performed in this dog study: two single-isotope and one dual-isotope acquisition. The 201Tl single-isotope images were used as references. The total number of counts, and the contrast between the left ventricular cavity (LVC) and the myocardium, were used in 70 keV short axis slices as parameters for evaluating the results of the dual-isotope correction methods. Three consecutive short-axis slices were used to calculate averaged contrast and the averaged total number of counts. The total number of the counts was 667000+/-500 and 414500+/-400 counts for the dual isotope (201Tl+/-99mTc) and single-isotope (201Tl-only) 70 keV images, respectively. The corrected dual-isotope images had 514700+/-700 and 368000+/-600 counts for Moore's correction and our approach, respectively. Moore's method improved contrast in the dual isotope 70 keV image to 0.14+/-0.03 from 0.11+/-0.02, which was the value in the 70 keV non-corrected dual-isotope image. Our method improved the same contrast to 0.22+/-0.03. The contrast in the 201Tl single-isotope 70 keV image was 0.28+/-0.02. Both methods improved the 70 keV dual-isotope images. However, our approach provided slightly better images than Moore's correction when compared with 201Tl-only 70 keV images.     

An electronically collimated gamma camera for single photon emission computed tomography. Part I: Theoretical considerations and design criteria

The detection and imaging characteristics of a new type of gamma camera for single photon emission computed tomography have been investigated. Unlike conventional gamma cameras which use mechanical collimation, the new gamma camera utilizes electronic collimation which is obtained from a sequential interaction of gamma radiation with a dual position-and-energy sensitive detection system. Coincident counting between the two detectors provides localization of activity upon a multitude of conical surfaces throughout the object, wherefrom the three-dimensional activity distribution can be reconstructed. Not only does electronic collimation provide simultaneous multiple views of the object, but a large gain in sensitivity is also indicated over a conventionally collimated gamma camera under conditions of similar spatial resolution. Detector optimization studies have been performed to design a prototype system comprising a 33 X 33 array of high-purity germanium detectors coupled to an uncollimated conventional scintillation camera. The cumulative signal-to-noise ratio in projection images obtained with this system is expected to be about a factor of 4 higher (sensitivity about a factor of 15 higher) than that obtained in a corresponding projection image with a conventional gamma camera for imaging a uniformly distributed Tc-99m source in a 20-cm-diam X 20-cm-tall cylinder. A similar gain is expected in the tomographic images.     

Performance comparison of two dual-head coincidence cameras of the first and latest generation

This study compares the performance and image quality of two gamma camera based PET systems of the first and latest generation. We investigated two dual head coincidence gamma cameras (PRISM 2000XP and AXIS, manufactured in 1997 and 2001 by Picker/Philips) predominantly in accordance with the NEMA NU2-1994 and NU2-2001 protocols. All performance parameters except for spatial resolution and image quality were determined after measuring a standard cylinder over several half-life periods. Scatter and random fractions were evaluated with the sinogram technique. In order to determine spatial resolution and image quality we used phantoms as described in the NEMA NU2-2001 protocol. The efficiency of the new system was found to be increased. True count rate at activity levels used in clinical conditions is improved and scatter fraction is decreased substantially. Accordingly, improved spatial resolution and image quality were observed with the new system. Altogether, the AXIS represents a further approach to the performance of dedicated positron emission tomographs.     

Dual-window scatter correction and energy window setting in cerebral blood flow SPECT: a Monte Carlo study

The image quality in SPECT studies of the regional cerebral blood flow (rCBF) performed with 99mTc-HMPAO is degraded by scattered photons. The finite energy resolution of the gamma camera makes the detection of scattered photons unavoidable, and this is observed in the image as an impaired contrast between grey and white matter structures. In this work, a Monte Carlo simulated SPECT study of a realistic voxel-based brain phantom was used to evaluate the resulting contrast-to-noise ratio for a number of energy window settings, with and without the dual-window scatter correction. Values of the scaling factor k, used to obtain the fraction of scattered photons in the photopeak window, were estimated for each energy window. The use of a narrower, asymmetric, energy discrimination window improved the contrast, with a subsequent increase in statistical noise due to the lower number of counts. The photopeak-window setting giving the best contrast-to-noise ratio was found to be the same whether or not scatter correction was applied. Its value was 17% centred at 142 keV. At the optimum photopeak-window setting, the contrast was improved by using scatter correction, but the contrast-to-noise ratio was made worse.     

Design and evaluation of an LSO PET detector for breast cancer imaging

Functional imaging with positron emission tomography (PET) may be a promising technique in conjunction with x-ray mammography for breast cancer patient management. Conventional whole body PET scanners provide metabolic images of breast cancer patients with several shortcomings related to the general-purpose nature of these systems. In whole body scanners, the detectors are typically 20-30 cm away from the breast or axilla, reducing sensitivity, and these scanners have relatively large detector elements (> 4 mm), limiting spatial resolution. Dedicated PET systems for breast imaging aim to overcome these limitations and improve the overall diagnostic quality of the images by bringing the detectors closer to the area to be imaged, thereby improving sensitivity, and by using smaller detector elements to improve the spatial resolution. We have designed and developed a modular PET detector that is composed of a 9x9 array of 3x3x20 mm3 lutetium oxyorthosilicate (LSO) scintillator crystals coupled to an optical fiber taper, which in turn is coupled to a Hamamatsu R5900-C8 position-sensitive photomultiplier tube. These detectors can be tiled together without gaps to construct large area detector arrays to form a dedicated PET breast cancer imaging system. Two complete detector modules have been built and tested. All detector elements are clearly visualized upon flood irradiation of the module. The intrinsic spatial resolution (full-width at half-maximum) was measured to be 2.26 mm (range 1.8-2.6 mm). The average energy resolution was 19.5% (range 17%-24%) at 511 keV. The coincidence time resolution was measured to be 2.4 ns. The detector efficiency for 511 keV gamma rays was 53% using a 350 keV energy threshold. These promising results support the feasibility of developing a high resolution, high sensitivity dedicated PET scanner for breast cancer applications.     

Experimental comparison of high-density scintillators for EMCCD-based gamma ray imaging

Detection of x-rays and gamma rays with high spatial resolution can be achieved with scintillators that are optically coupled to electron-multiplying charge-coupled devices (EMCCDs). These can be operated at typical frame rates of 50 Hz with low noise. In such a set-up, scintillation light within each frame is integrated after which the frame is analyzed for the presence of scintillation events. This method allows for the use of scintillator materials with relatively long decay times of a few milliseconds, not previously considered for use in photon-counting gamma cameras, opening up an unexplored range of dense scintillators. In this paper, we test CdWO? and transparent polycrystalline ceramics of Lu?O?:Eu and (Gd,Lu)?O?:Eu as alternatives to currently used CsI:Tl in order to improve the performance of EMCCD-based gamma cameras. The tested scintillators were selected for their significantly larger cross-sections at 140 keV ((99m)Tc) compared to CsI:Tl combined with moderate to good light yield. A performance comparison based on gamma camera spatial and energy resolution was done with all tested scintillators having equal (66%) interaction probability at 140 keV. CdWO?, Lu?O?:Eu and (Gd,Lu)?O?:Eu all result in a significantly improved spatial resolution over CsI:Tl, albeit at the cost of reduced energy resolution. Lu?O?:Eu transparent ceramic gives the best spatial resolution: 65 μm full-width-at-half-maximum (FWHM) compared to 147 μm FWHM for CsI:Tl. In conclusion, these 'slow' dense scintillators open up new possibilities for improving the spatial resolution of EMCCD-based scintillation cameras.     

Generalized five-dimensional dynamic and spectral factor analysis

We have generalized the spectral factor analysis and the factor analysis of dynamic sequences (FADS) in SPECT imaging to a five-dimensional general factor analysis model (5D-GFA), where the five dimensions are the three spatial dimensions, photon energy, and time. The generalized model yields a significant advantage in terms of the ratio of the number of equations to that of unknowns in the factor analysis problem in dynamic SPECT studies. We solved the 5D model using a least-squares approach. In addition to the traditional non-negativity constraints, we constrained the solution using a priori knowledge of both time and energy, assuming that primary factors (spectra) are Gaussian-shaped with full-width at half-maximum equal to gamma camera energy resolution. 5D-GFA was validated in a simultaneous pre-/post-synaptic dual isotope dynamic phantom study where 99mTc and 123I activities were used to model early Parkinson disease studies. 5D-GFA was also applied to simultaneous perfusion/dopamine transporter (DAT) dynamic SPECT in rhesus monkeys. In the striatal phantom, 5D-GFA yielded significantly more accurate and precise estimates of both primary 99mTc (bias=6.4 % +/- 4.3 %) and 1231 (-1.7% +/- 6.9%) time activity curves (TAC) compared to conventional FADS (biases = 15.5% +/- 10.6% in 99mTc and 8.3% +/- 12.7% in 123I, p < 0.05). Our technique was also validated in two primate dynamic dual isotope perfusion/DAT transporter studies. Biases of 99mTc-HMPAO and 123I-DAT activity estimates with respect to estimates obtained in the presence of only one radionuclide (sequential imaging) were significantly lower with 5D-GFA (9.4% +/- 4.3% for 99mTc-HMPAO and 8.7% +/-4.1% for 123I-DAT) compared to biases greater than 15% for volumes of interest (VOI) over the reconstructed volumes (p < 0.05). 5D-GFA is a novel and promising approach in dynamic SPECT imaging that can also be used in other modalities. It allows accurate and precise dynamic analysis while compensating for Compton scatter and cross-talk.     

Characterization of a high-purity germanium detector for small-animal SPECT

We present an initial evaluation of a mechanically cooled, high-purity germanium double-sided strip detector as a potential gamma camera for small-animal SPECT. It is 90 mm in diameter and 10 mm thick with two sets of 16 orthogonal strips that have a 4.5 mm width with a 5 mm pitch. We found an energy resolution of 0.96% at 140 keV, an intrinsic efficiency of 43.3% at 122 keV and a FWHM spatial resolution of approximately 1.5 mm. We demonstrated depth-of-interaction estimation capability through comparison of pinhole acquisitions with a point source on and off axes. Finally, a flood-corrected flood image exhibited a strip-level uniformity of less than 1%. This high-purity germanium offers many desirable properties for small-animal SPECT.     

Experimentally measured scatter fractions and energy spectra as a test of Monte Carlo simulations

A method for the validation of Monte Carlo photon transport calculations is presented, with particular emphasis on the scatter component of such calculations. The method is based on a quantitative comparison of calculated and experimental scatter fractions. In addition, the method includes a qualitative comparison of point spread functions and energy spectra. An application of the method is demonstrated by comparing the results of an existing Monte Carlo code with experimental results obtained with a gamma camera viewing a point source of 99Tcm (140 keV gamma rays) centred within a water-filled cylinder. The results of the comparisons show good agreement between experiment and calculation. These results allow the code to be used with increased confidence in a variety of situations, and they define more precisely the region of applicability of the code. In addition, the determination of scatter fractions and energy spectra is useful for other applications. For example, scatter fractions can be a useful parameter for evaluating possible techniques for scatter compensation.     

Validation of the GATE Monte Carlo simulation platform for modelling a CsI(Tl) scintillation camera dedicated to small-animal imaging

Monte Carlo simulations are increasingly used in scintigraphic imaging to model imaging systems and to develop and assess tomographic reconstruction algorithms and correction methods for improved image quantitation. GATE (GEANT4 application for tomographic emission) is a new Monte Carlo simulation platform based on GEANT4 dedicated to nuclear imaging applications. This paper describes the GATE simulation of a prototype of scintillation camera dedicated to small-animal imaging and consisting of a CsI(Tl) crystal array coupled to a position-sensitive photomultiplier tube. The relevance of GATE to model the camera prototype was assessed by comparing simulated 99mTc point spread functions, energy spectra, sensitivities, scatter fractions and image of a capillary phantom with the corresponding experimental measurements. Results showed an excellent agreement between simulated and experimental data: experimental spatial resolutions were predicted with an error less than 100 microns. The difference between experimental and simulated system sensitivities for different source-to-collimator distances was within 2%. Simulated and experimental scatter fractions in a [98-182 keV] energy window differed by less than 2% for sources located in water. Simulated and experimental energy spectra agreed very well between 40 and 180 keV. These results demonstrate the ability and flexibility of GATE for simulating original detector designs. The main weakness of GATE concerns the long computation time it requires: this issue is currently under investigation by the GEANT4 and the GATE collaborations.     

Study of the point spread function (PSF) for 123I SPECT imaging using Monte Carlo simulation

The iterative reconstruction algorithms employed in brain single-photon emission computed tomography (SPECT) allow some quantitative parameters of the image to be improved. These algorithms require accurate modelling of the so-called point spread function (PSF). Nowadays, most in vivo neurotransmitter SPECT studies employ pharmaceuticals radiolabelled with 123I. In addition to an intense line at 159 keV, the decay scheme of this radioisotope includes some higher energy gammas which may have a non-negligible contribution to the PSF. The aim of this work is to study this contribution for two low-energy high-resolution collimator configurations, namely, the parallel and the fan beam. The transport of radiation through the material system is simulated with the Monte Carlo code PENELOPE. We have developed a main program that deals with the intricacies associated with tracking photon trajectories through the geometry of the collimator and detection systems. The simulated PSFs are partly validated with a set of experimental measurements that use the 511 keV annihilation photons emitted by a 18F source. Sensitivity and spatial resolution have been studied, showing that a significant fraction of the detection events in the energy window centred at 159 keV (up to approximately 49% for the parallel collimator) are originated by higher energy gamma rays, which contribute to the spatial profile of the PSF mostly outside the 'geometrical' region dominated by the low-energy photons. Therefore, these high-energy counts are to be considered as noise, a fact that should be taken into account when modelling PSFs for reconstruction algorithms. We also show that the fan beam collimator gives higher signal-to-noise ratios than the parallel collimator for all the source positions analysed.