• Journal of the Korean Physical Society
• /
• v.73 no.9
• /
• pp.1393-1398
• /
• 2018

In conventional single-photon emission computed tomography (SPECT), only the photoelectric events in the detectors are used for image reconstruction. However, if the $^{131}I$ isotope, which emits high-energy radiations (364, 637, and 723 keV), is used in nuclear medicine, both photoelectric and Compton scattering events can be used for image reconstruction. The purpose of our work is to perform simulations for Compton SPECT by using the Geant4 application for tomographic emission (GATE). The performance of Compton SPECT is evaluated and compared with that of conventional SPECT. The Compton SPECT unit has an area of $12cm{\times}12cm$ with four gantry heads. Each head is composed of a 2-cm tungsten collimator and a $40{\times}40$ array of CdZnTe (CZT) crystals with a $3{\times}3mm^2$ area and a 6-mm thickness. Compton SPECT can use not only the photoelectric effect but also the Compton scattering effect for image reconstruction. The correct sequential order of the interactions used for image reconstruction is determined using the angular resolution measurement (ARM) method and the energies deposited in each detector. In all the results of simulations using spherical volume sources of various diameters, the reconstructed images of Compton SPECT show higher signal-to-noise ratios (SNRs) without degradation of the image resolution when compared to those of conventional SPECT because the effective count for image reconstruction is higher. For a Derenzo-like phantom, the reconstructed images for different modalities are compared by visual inspection and by using their projected histograms in the X-direction of the reconstructed images.

• Biomedical Engineering Letters
• /
• v.8 no.4
• /
• pp.383-392
• /
• 2018

For prompt gamma ray imaging for biomedical applications and environmental radiation monitoring, we propose herein a multiple-scattering Compton camera (MSCC). MSCC consists of three or more semiconductor layers with good energy resolution, and has potential for simultaneous detection and differentiation of multiple radio-isotopes based on the measured energies, as well as three-dimensional (3D) imaging of the radio-isotope distribution. In this study, we developed an analytic simulator and a 3D image generator for a MSCC, including the physical models of the radiation source emission and detection processes that can be utilized for geometry and performance prediction prior to the construction of a real system. The analytic simulator for a MSCC records coincidence detections of successive interactions in multiple detector layers. In the successive interaction processes, the emission direction of the incident gamma ray, the scattering angle, and the changed traveling path after the Compton scattering interaction in each detector, were determined by a conical surface uniform random number generator (RNG), and by a Klein-Nishina RNG. The 3D image generator has two functions: the recovery of the initial source energy spectrum and the 3D spatial distribution of the source. We evaluated the analytic simulator and image generator with two different energetic point radiation sources (Cs-137 and Co-60) and with an MSCC comprising three detector layers. The recovered initial energies of the incident radiations were well differentiated from the generated MSCC events. Correspondingly, we could obtain a multi-tracer image that combined the two differentiated images. The developed analytic simulator in this study emulated the randomness of the detection process of a multiple-scattering Compton camera, including the inherent degradation factors of the detectors, such as the limited spatial and energy resolutions. The Doppler-broadening effect owing to the momentum distribution of electrons in Compton scattering was not considered in the detection process because most interested isotopes for biomedical and environmental applications have high energies that are less sensitive to Doppler broadening. The analytic simulator and image generator for MSCC can be utilized to determine the optimal geometrical parameters, such as the distances between detectors and detector size, thus affecting the imaging performance of the Compton camera prior to the development of a real system.

• Journal of Biomedical Engineering Research
• /
• v.32 no.1
• /
• pp.15-24
• /
• 2011

While Compton imaging is recognized as a valuable 3-D technique in nuclear medicine, reconstructing an image from Compton scattered data has been of a difficult problem due to its computational complexity. The most complex and time-consuming computation in Compton camera reconstruction is to perform the conical projection and backprojection operations. To alleviate the computational burden imposed by these operations, we investigate a rebinning method which can convert conical projections into parallel projections. The use of parallel projections allows to directly apply the existing deterministic reconstruction methods, which have been useful for conventional emission tomography, to Compton camera reconstruction. To convert conical projections into parallel projections, a cone surface is sampled with a number of lines. Each line is projected onto an imaginary plane that is mostly perpendicular to the line. The projection data rebinned in each imaginary plane can then be treated as the standard parallel projection data. To validate the rebinning method, we tested with the representative deterministic algorithms, such as the filtered backprojection method and the algebraic reconstruction technique. Our experimental results indicate that the rebinning method can be useful when the direct application of existing deterministic methods is needed for Compton camera reconstruction.

• Nuclear Engineering and Technology
• /
• v.52 no.10
• /
• pp.2361-2369
• /
• 2020

Compton cameras using large scintillators have been developed for high imaging sensitivity. These scintillator-based Compton cameras, however, mainly due to relatively low energy resolution, suffer from undesired background-radiation signals, especially when radioactive materials' activity is very low or their location is far from the Compton camera. To alleviate this problem for a large-size Compton camera, in the present study, a hybrid-type shielding system was designed that combines an active shield with a veto detector and a passive shield that surrounds the active shield. Then, the performance of the hybrid shielding system was predicted, by Monte Carlo radiation transport simulation using Geant4, in terms of minimum detectable activity (MDA), signal-to-noise ratio (SNR), and image resolution. Our simulation results show that, for the most cases, the hybrid shielding system significantly improves the performance of the large-size Compton camera. For the cases investigated in the present study, the use of the shielding system decreased the MDA by about 1.4, 1.6, and 1.3 times, increased the SNR by 1.2-1.9, 1.1-1.7, and 1.3-2.1 times, and improved the image resolution (i.e., reduced the FWHM) by 7-8, 1-6, and 3-5% for ¹³⁷Cs, ⁶⁰Co, and ¹³¹I point source located at 1-5 m from the imaging system, respectively.

• Progress in Medical Physics
• /
• v.20 no.2
• /
• pp.51-61
• /
• 2009

A Compton camera, which is based on the geometrical interpretation of Compton scattering, is a very promising gamma-ray imaging device considering its several advantages over the conventional gamma-ray imaging devices: high imaging sensitivity, 3-D imaging capability from a fixed position, multi-tracing functionality, and almost no limitation in photon energy. In the present study, a Monte Carlo-based, user-friendly Compton imaging simulator was developed in the form of a graphical user interface (GUI) based on Geant4 and $MATLAB^{TM}$. The simulator was tested against the experimental result of the double-scattering Compton camera, which is under development at Hanyang University in Korea. The imaging resolution of the simulated Compton image well agreed with that of the measured image. The imaging sensitivity of the measured data was 2~3 times higher than that of the simulated data, which is due to the fact that the measured data contains the random coincidence events. The performance of a stacking-structure type Compton camera was evaluated by using the simulator. The result shows that the Compton camera shows its highest performance when it uses 4 layers of scatterer detectors.

• Proceedings of the Korean Society of Medical Physics Conference
• /
• 2006.08a
• /
• pp.125-125
• /
• 2006

• Journal of Biomedical Engineering Research
• /
• v.28 no.4
• /
• pp.520-529
• /
• 2007

An efficient method for implementing image reconstruction algorithms for Compton cameras is presented. Since Compton scattering formula establishes a cone surface from which the incident photon must have originated, it is crucial to implement a computationally efficient cone-surface integration method for image reconstruction. In this paper we assume that a cone is made up of a series of ellipses (or circles) stacked up one on top of the other. In order to reduce computational burden for tracing ellipses formed by the intersection of a cone and an image plane, we propose a new method using a series of imaginary planes perpendicular to the cone axis so that each plane contains a circle, not an ellipse. In this case the cone surface integral can be performed by simply accumulating the circles along the cone axis. To reduce the computational cost of tracing circles, only one of the circles in the cone is traced and the rest are determined by using simple trigonometric ratios. For our experiments, we used the three different schemes for tracing ellipses; (i) using the samples generated by the ellipse equation, (ii) using the fixed number of samples along a circle on the imaginary plane, and (iii) using the fixed sampling interval along a circle on the imaginary plane. We then compared performance of the above three methods by applying them to the two reconstruction algorithms - the simple back-projection method and the expectation-maximization algorithm. The experimental results demonstrate that our proposed methods (ii) and (iii) using imaginary planes significantly improve reconstruction accuracy as well as computational efficiency.

• Journal of Biomedical Engineering Research
• /
• v.30 no.1
• /
• pp.56-65
• /
• 2009

Compton imaging is often recognized as a potentially more valuable 3-D technique in nuclear medicine than conventional emission tomography. Due to inherent computational limitations, however, it has been of a difficult problem to reconstruct images with good accuracy. In this work we show that the row-action maximum likelihood algorithm (RAMLA), which have proven useful for conventional tomographic reconstruction, can also be applied to the problem of 3-D reconstruction of cone-beam projections from Compton scattered data. The major advantage of RAMLA is that it converges to a true maximum likelihood solution at an order of magnitude faster than the standard expectation maximiation (EM) algorithm. For our simulations, we first model a Compton camera system consisting of the three pairs of scatterer and absorber detectors placed at x-, y- and z-axes, and generate conical projection data using a software phantom. We then compare the quantitative performance of RAMLA and EM reconstructions in terms of the percentage error. The net conclusion based on our experimental results is that the RAMLA applied to Compton camera reconstruction significantly outperforms the EM algorithm in convergence rate; while computational costs of one iteration of RAMLA and EM are about the same, one iteration of RAMLA performs as well as 128 iterations of EM.

• The Journal of Engineering Research
• /
• v.8 no.1
• /
• pp.51-57
• /
• 2006

We propose efficient methods for implementing 3-D reconstruction algorithms for Compton camera. Since reconstructing Compton scattered data involves the surface integral over the cone associated with the measurement bin, it is crucial to develop a computationally efficient surface integration method. In this work we assume that a cone is made up of a series of ellipses (or circles) stacked up one o top of the other. In order to reduce computational burden for tracing ellipses formed by the intersection of a cone and an image plane, we construct a series of imaginary planes perpendicular to the cone axis so that each plane contains a circle, not an ellipse. In this case the surface integral can be performed by adding uniform samples along each circle. The experimental results demonstrate that our method using imaginary planes significantly improves computational efficiency while keeping reconstruction accuracy.

• Journal of Radiation Protection and Research
• /
• v.34 no.3
• /
• pp.107-114
• /
• 2009

A Compton camera, which is based on Compton kinematics, is a very promising gamma-ray imaging device in that it could overcome the limitations of the conventional gamma-ray imaging devices. In the present study, the image quality of a rotating Compton camera was evaluated by using 4-D Monte Carlo simulation technique and the applicability to nuclear industrial applications was examined. It was found that Compton images were significantly improved when the Compton camera rotates around a gamma-ray source. It was also found that the 3-D imaging capability of a Compton camera could enable us to accurately determine the 3-D location of radioactive contamination in a concrete wall for decommissioning purpose of nuclear facilities. The 4-D Monte Carlo simulation technique, which was applied to the Compton camera fields for the first time, could be also used to model the time-dependent geometry for various applications.