Abstract
Special collimators used in imaging systems play an important part in obtaining qualified images to improve diagnosis in medicine. Methods: The primary aim of this study was to compare resolution between fan beam and parallel beam collimators using Monte Carlo simulation in the shape of cubic holes. Also, parameters such as geometric efficiency, geometric resolution, scatter, penetration, and full width at half maximum were studied to compare their special characteristics. Results: The simulation results demonstrated that the geometry efficiency of a fan beam collimator increased as the angle of the slant hole increased, and the geometric resolution decreased as the angle of the slant hole increased, at a distinct distance from a monoenergetic source of γ-rays. In contrast, at a distinct angle, geometric resolution increased as the distance between the source and the collimator surface increased. For both collimators, scatter and penetration decreased as the distance increased. These results were in agreement with ADAC company data. Finally, fan beam collimators were found to have better resolution than parallel beam collimators with a cubic hole shape in a wire mesh design. Conclusion: Estimation of the fan beam by parallel beam parameters as cubic holes can be suitable in collimator design to improve resolution and efficiency.
Collimation of low-energy photons, allowing detection of only those photons propagating in the appropriate direction, plays a key role in obtaining a suitable map in SPECT imaging. Some collimators, such as the parallel beam and fan beam types, are widely used in different applications and for organs centered within the field of view. Therefore, it is important to study collimator geometric properties to improve diagnosis in medicine. Fan beam collimators are a special type of converging collimator with the holes focusing toward a focal line parallel to the axis of rotation of the camera (1).
For optimization of collimator design, photon interactions in γ-camera collimators have been simulated by the Monte Carlo N-particle code (MCNP5) (2–4). In this code, source geometry, collimators, and detectors need to be defined in the input file—a cumbersome process (5,6). Recently, researchers have clinically compared parallel beam and fan beam collimators and simulated the parallel beam collimator by this code (7–22). In this study, fan beam parameters were deliberated by this code.
To improve resolution and efficiency, it is helpful to model collimators of various hole shapes, edge effects, septal materials, and geometric configurations. The aim of this study was to investigate resolution between fan beam and parallel beam collimators using MCNP5 simulation in the shape of cubic holes. Also, distinctive parameters such as geometric efficiency, scatter, penetration, and full width at half maximum (FWHM) have been compared, with assessment of the fan beam and parallel beam collimator responses.
MATERIALS AND METHODS
Computation of Cubic Hole Angles
In collimator design, static and dynamic parameters are considered. Static parameters include material and geometry, such as the shape and size of the holes and the degree of the angle. Dynamic parameters include the position of the source and the location of the collimator. In this study, the geometric characteristic of low-energy parallel beam and fan beam collimators used at ADAC were applied, as indicated in Table 1 (23). In view of research by Gunter et al. showing that parallel beam parameters can be chosen locally for fan beam collimators, we estimated the fan beam geometric dimensions by parallel beam collimator parameters (24).
The MCNP5 code and Matlab software (MathWorks, Inc.) were used to simulate collimators with cubic holes as shown in Figure 1. On the central z-axis, both the fan beam collimator and the parallel beam collimator have the same hole length. On the x-axis, as the offset distance changes, the hole length of the fan beam collimator increases because of the various angles and the hole size.
According to ADAC company data, the focal length of the fan beam collimator is 50 cm and the slant hole cross section is variable (23). A simulated collimator was defined by a cubic region because the unit cell, when replicated, would produce the desired hole pattern in parallel beam mode, as indicated in Figure 1A. In fan beam mode, each hole must be considered, with different angles, hole sizes, and hole lengths and constant septa in a similar volume. The cross section of the holes is of a parallelepiped shape, and the hole angles on the fan beam collimator can be computed given the geometric parameters. As the angle increases, the cross section varies according to a row of holes focused on a point. Also, there is another row of holes centered on another focal point in the same direction at an equal distance from the collimator surface. The focal line is formed by the focal points. Figure 2 represents this objective with dissimilar diameters for 2 holes beside each other and can be defined as:
Monte Carlo Simulations
Each camera system has a constant intrinsic resolution, depending on detector thickness, area, and substance density. This study used a NaI detector with a thickness of 0.95 cm in a 40 × 54 cm area, and a scintillation camera with an intrinsic resolution of 0.39 cm was simulated by MCNP5. Several parameters were evaluated, such as FWHM, geometric resolution, and geometric efficiency for 89,600 slant holes and 89,600 cross sections.
In fan beam collimator simulation, angles were calculated for 89,600 slant holes in a 40 × 54 × 3.28 cm region with a septal thickness of 0.0152 cm and various cross sections, in which the holes focused to a line. Metz et al. (22) derived a closed algebraic expression to determine the geometric component of the parallel beam collimators and found that if penetration is negligible, the spatially geometric component is effectively independent of hole shape for a given collimator thickness and a given ratio of hole area to septal area. Table 1 indicates the geometric properties of the collimators used for the simulation.
The point source, technetium (140 keV), was placed in air and medium (water phantom) to obtain adequate scatter, penetration, and point spread. This source was located at z = 10 and 15 cm in front of the collimator surface, and energy spectra with a 20% photopeak window were obtained (25). In each case, 500 million photons were detected by the NaI detector.
The phantom was considered a cylinder with a length of 20 cm and a radius of 7 cm to study source depth in tissue. The depths were assumed to be 2 cm (z = 10 cm) and 7 cm (z = 15 cm), as shown in Figure 3. Meanwhile, for calculation of scatter and penetration, simulations were performed with and without interactions in collimators of the same sensitivity.
Calculation and Assessment of Collimator Parameters
The ratio of hole area (Aopen) to septal area (Aunit) can be defined in the geometric efficiency of the parallel beam collimator (GPB) as follows:
RESULTS
The effective and total resolutions increased when the distance between the source and the parallel beam collimator increased, as shown in Table 2. The effective and total resolutions reported by ADAC at a 10-cm distance for the parallel beam collimator were 0.63 and 0.74 cm, respectively. These values agree with our results, with errors of 3.2% and 2.7%, respectively. Also, the geometric efficiency of the parallel beam collimator and effective hole length were calculated to be 7.69387E−05 and 3.189 cm, respectively.
The geometric efficiency of a fan beam collimator is related to angle and distance and, for a distinct angle, could be increased by increasing the distance, as indicated in Table 3. Geometric efficiency could also be decreased by increasing the slant hole angle, as shown in Figure 4.
The effective and total resolutions of a fan beam collimator with 89,600 holes, for all angles, were calculated to be 0.59 and 0.68 cm, respectively, at a 10-cm distance. The reported effective and total resolutions from ADAC under the same conditions were 0.60 and 0.71 cm, respectively. At a distance of 15 cm, the respective calculated values were 0.69 and 0.74 cm, as shown in Tables 4 and 5 and Figure 5. The values increased with an increase in the hole angle at a distinct distance, as well as with an increase in distance at a distinct angle.
For both collimators, the FWHM using cubic holes was investigated by radial point spread functions. FWHM increased as the distance increased, as shown in Table 6. With cubic holes, the FWHM of the fan beam collimator is less than that of the parallel beam collimator and results in a better resolution. Also, the radial point spread functions were simulated in a cubic shape at a 10-cm distance in air by MCNP5. The fan beam collimator has better resolution than the parallel beam collimator, as shown in Figure 6A. The FWHM of collimators in the water phantom was increased because of Compton scattering, and the resolution was decreased, as indicated in Figure 6B.
Table 7 presents relative differences between the simulated and the reported resolutions at a 10-cm distance from collimators including 89,600 holes: 6% and 3% for the parallel beam and fan beam collimators, respectively. Table 8 indicates the scatter and penetration values in air and water for a 20% photopeak window at various distances. For both collimators, scatter and penetration decreased as distance increased. At a distinct distance, scatter and penetration were found to be higher for the fan beam collimator than for the parallel beam collimator in both air and water because of the hole shape, geometric pattern, and edge effect.
DISCUSSION
This study investigated the resolution and geometric efficiency of parallel beam and fan beam collimators with cubic holes. The spatial resolution of a scintillation camera stems from a combination of the type of collimator and the intrinsic resolution characteristics of the crystal and electronic detection system. When the angle of the hole is changed, efficiency and resolution are altered, and vice versa. Thus, collimator parameters are designed with regard to the necessary trade-off between efficiency and resolution.
Cot et al. (28) indicated the square holes give rise to sensitivity higher than the hexagonal holes without considering the effective length and efficiency in actual extensive collimators. In this study, the simulations included all components plus photoelectric and electron x-rays with higher accuracy (<0.02).
To study the effects of interactions, the simulations included penetration, Compton and coherent scatter, photoelectric effects, and emission of characteristic lead x-rays. The programs also included interactions in a NaI scintillation crystal, and backscatter from a compartment representing detector components behind the crystal. The results of the simulation demonstrated that with an increase in distance there is a decrease in scatter and penetration. The fan beam had higher scatter and penetration, and the parallel beam had lower values, with a 1.04 ratio in air and 1.76 in phantom at z = 10 cm. Scatter and penetration depend on the hole shape, geometric pattern, and edge effect. When scatter and penetration are not negligible, the point spread function will have a long tail, which will not affect FWHM but can significantly affect the contrast resolution of the collimator. Although system resolution was higher than geometric resolution, the fan beam, in comparison with the parallel beam with cubic holes, had the lowest system and geometric resolutions for a distinct distance.
CONCLUSION
This comparison demonstrated that simulated responses were close to the experimental data provided by ADAC. Compared with the parallel beam collimator, the fan beam collimator with a cubic hole shape resulted in a smaller full width at half maximum and better resolution. Modeling of collimators can be commercially beneficial to the manufacturing industry. In addition, modeling can be clinically relevant in allowing for the design of collimators with better qualities and features.
Acknowledgments
No potential conflict of interest relevant to this article was reported.
Footnotes
Published online Jan. 11, 2012.
REFERENCES
- Received for publication March 3, 2011.
- Accepted for publication October 27, 2011.