The performance of a hybrid analytical-Monte Carlo system response matrix in pinhole SPECT reconstruction.

It is well-known that in pinhole SPECT (single-photon-emission computed tomography), iterative reconstruction methods including accurate estimations of the system response matrix can lead to submillimeter spatial resolution. There are two different methods for obtaining the system response matrix: those that model the system analytically using an approach including an experimental characterization of the detector response, and those that make use of Monte Carlo simulations. Methods based on analytical approaches are faster and handle the statistical noise better than those based on Monte Carlo simulations, but they require tedious experimental measurements of the detector response. One suggested approach for avoiding an experimental characterization, circumventing the problem of statistical noise introduced by Monte Carlo simulations, is to perform an analytical computation of the system response matrix combined with a Monte Carlo characterization of the detector response. Our findings showed that this approach can achieve high spatial resolution similar to that obtained when the system response matrix computation includes an experimental characterization. Furthermore, we have shown that using simulated detector responses has the advantage of yielding a precise estimate of the shift between the point of entry of the photon beam into the detector and the point of interaction inside the detector. Considering this, it was possible to slightly improve the spatial resolution in the edge of the field of view.

A parallel computational model for GATE simulations.

GATE/Geant4 Monte Carlo simulations are computationally demanding applications, requiring thousands of processor hours to produce realistic results. The classical strategy of distributing the simulation of individual events does not apply efficiently for Positron Emission Tomography (PET) experiments, because it requires a centralized coincidence processing and large communication overheads. We propose a parallel computational model for GATE that handles event generation and coincidence processing in a simple and efficient way by decentralizing event generation and processing but maintaining a centralized event and time coordinator. The model is implemented with the inclusion of a new set of factory classes that can run the same executable in sequential or parallel mode. A Mann-Whitney test shows that the output produced by this parallel model in terms of number of tallies is equivalent (but not equal) to its sequential counterpart. Computational performance evaluation shows that the software is scalable and well balanced.

Performance Characteristics of BGO Detectors for a Low Cost Preclinical PET Scanner.

PETbox is a low-cost benchtop PET scanner dedicated to high throughput preclinical imaging that is currently under development at our institute. This paper presents the design and characterization of the detectors that are used in the PETbox system. In this work, bismuth germanate scintillator was used for the detector, taking advantage of its high stopping power, high photoelectric event fraction, lack of intrinsic background radiation and low cost. The detector block was segmented into a pixelated array consisting of 20 × 44 elements, with a crystal pitch of 2.2 mm and a crystal cross section of 2 mm × 2 mm. The effective area of the array was 44 mm × 96.8 mm. The array was coupled to two Hamamatsu H8500 position sensitive photomultiplier tubes, forming a flat-panel type detector head with a sensitive area large enough to cover the whole body of a typical laboratory mouse. Two such detector heads were constructed and their performance was characterized. For one detector head, the energy resolution ranged from 16.1% to 38.5% full width at half maximum (FWHM), with a mean of 20.1%; for the other detector head, the energy resolution ranged from 15.5% to 42.7% FWHM, with a mean of 19.6%. The intrinsic spatial resolution was measured to range from 1.55 mm to 2.39 mm FWHM along the detector short axis and from 1.48 mm to 2.33 mm FWHM along the detector long axis, with an average of 1.78 mm. Coincidence timing resolution for the detector pair was measured to be 4.1 ns FWHM. These measurement results show that the detectors are suitable for our specific application.

GATE: a simulation toolkit for PET and SPECT.

Monte Carlo simulation is an essential tool in emission tomography that can assist in the design of new medical imaging devices, the optimization of acquisition protocols and the development or assessment of image reconstruction algorithms and correction techniques. GATE, the Geant4 Application for Tomographic Emission, encapsulates the Geant4 libraries to achieve a modular, versatile, scripted simulation toolkit adapted to the field of nuclear medicine. In particular, GATE allows the description of time-dependent phenomena such as source or detector movement, and source decay kinetics. This feature makes it possible to simulate time curves under realistic acquisition conditions and to test dynamic reconstruction algorithms. This paper gives a detailed description of the design and development of GATE by the OpenGATE collaboration, whose continuing objective is to improve, document and validate GATE by simulating commercially available imaging systems for PET and SPECT. Large effort is also invested in the ability and the flexibility to model novel detection systems or systems still under design. A public release of GATE licensed under the GNU Lesser General Public License can be downloaded at http:/www-lphe.epfl.ch/GATE/. Two benchmarks developed for PET and SPECT to test the installation of GATE and to serve as a tutorial for the users are presented. Extensive validation of the GATE simulation platform has been started, comparing simulations and measurements on commercially available acquisition systems. References to those results are listed. The future prospects towards the gridification of GATE and its extension to other domains such as dosimetry are also discussed.