[A simulation method for studying scintillation camera collimators]

Study of the point spread function (PSF) for123I 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.

Relevance of accurate Monte Carlo modeling in nuclear medical imaging

Monte Carlo techniques have become popular in different areas of medical physics with advantage of powerful computing systems. In particular, they have been extensively applied to simulate processes involving random behavior and to quantify physical parameters that are difficult or even impossible to calculate by experimental measurements. Recent nuclear medical imaging innovations such as single-photon emission computed tomography (SPECT), positron emission tomography (PET), and multiple emission tomography (MET) are ideal for Monte Carlo modeling techniques because of the stochastic nature of radiation emission, transport and detection processes. Factors which have contributed to the wider use include improved models of radiation transport processes, the practicality of application with the development of acceleration schemes and the improved speed of computers. In this paper we present a derivation and methodological basis for this approach and critically review their areas of application in nuclear imaging. An overview of existing simulation programs is provided and illustrated with examples of some useful features of such sophisticated tools in connection with common computing facilities and more powerful multiple-processor parallel processing systems. Current and future trends in the field are also discussed.

[Simulation study of the influence of collimator defects on the uniformity of scintigraphic images]

The mechanical tolerances in building collimators for scintillation cameras are studied. A simulation method has been used to quantify the effects of defects in hole inclination and hole diameter on the uniformity of planar and tomographic images. The calculation takes into account the geometry of the hexagonal hole collimator, the camera intrinsic resolution, the object size, the pixel size, the effect of low-pass filtering, as well as the type, size and position of the defect. For instance, a 0.03 mm diameter defect on several holes located in the central region of a very high resolution collimator can result in a 12% uniformity artefact in tomographic imaging of an 18 cm diameter cylinder, using a 3.45 mm resolution camera, 4.5 mm pixel size, and Hamming filtering with a Nyquist frequency cut-off. A 0.17 degree inclination defect of a few holes can result in the same uniformity artefact. These results show that the building of a collimator has to be very precise.

[Characterization of the response of hexagonal parallel-hole collimators of scintillation cameras]

This paper presents a formulation of the frequency response of hexagonal parallel-hole collimator scintillation cameras. To describe this response, we propose an equation determined semi-empirically from a great number of simulations. The utility of the equation is that it enables the simple calculation of the response from collimator characteristics by taking into account the collimator's hexagonal structure. Because the equation does not assume translation invariance, the results can be directly compared with experimental measurements obtained with a point source. It is particularly interesting for collimators with large holes, like the medium-resolution ones used for high-energy radiation. Quality control and physical performance measurements are thus facilitated for this kind of collimator. Also, we present a new parameter that gives a quantitative assessment of the importance of partition penetration. This parameter is measured directly from the collimator frequency response. It has been studied by simulation, taking into account gamma photon attenuation in collimator partitions. The experimental measurements that have been made are in accord with the proposed equations.