skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: GPU-based prompt gamma ray imaging from boron neutron capture therapy

Abstract

Purpose: The purpose of this research is to perform the fast reconstruction of a prompt gamma ray image using a graphics processing unit (GPU) computation from boron neutron capture therapy (BNCT) simulations. Methods: To evaluate the accuracy of the reconstructed image, a phantom including four boron uptake regions (BURs) was used in the simulation. After the Monte Carlo simulation of the BNCT, the modified ordered subset expectation maximization reconstruction algorithm using the GPU computation was used to reconstruct the images with fewer projections. The computation times for image reconstruction were compared between the GPU and the central processing unit (CPU). Also, the accuracy of the reconstructed image was evaluated by a receiver operating characteristic (ROC) curve analysis. Results: The image reconstruction time using the GPU was 196 times faster than the conventional reconstruction time using the CPU. For the four BURs, the area under curve values from the ROC curve were 0.6726 (A-region), 0.6890 (B-region), 0.7384 (C-region), and 0.8009 (D-region). Conclusions: The tomographic image using the prompt gamma ray event from the BNCT simulation was acquired using the GPU computation in order to perform a fast reconstruction during treatment. The authors verified the feasibility of the prompt gamma ray imagemore » reconstruction using the GPU computation for BNCT simulations.« less

Authors:
; ;  [1];  [2];  [3]
  1. Department of Biomedical Engineering and Research Institute of Biomedical Engineering, College of Medicine, Catholic University of Korea, Seoul 505 137-701 (Korea, Republic of)
  2. Molecular Imaging Program at Stanford (MIPS), Department of Radiology, Stanford University, 300 Pasteur Drive, Stanford, California 94305 (United States)
  3. Department of Radiation Oncology, Stanford University School of Medicine, 875 Blake Wilbur Drive, Stanford, California 94305-5847 (United States)
Publication Date:
OSTI Identifier:
22413373
Resource Type:
Journal Article
Resource Relation:
Journal Name: Medical Physics; Journal Volume: 42; Journal Issue: 1; Other Information: (c) 2015 American Association of Physicists in Medicine; Country of input: International Atomic Energy Agency (IAEA)
Country of Publication:
United States
Language:
English
Subject:
60 APPLIED LIFE SCIENCES; ACCURACY; BORON; COMPUTERIZED SIMULATION; GAMMA RADIATION; IMAGE PROCESSING; MONTE CARLO METHOD; NEUTRON CAPTURE THERAPY; PHANTOMS; UPTAKE

Citation Formats

Yoon, Do-Kun, Jung, Joo-Young, Suk Suh, Tae, E-mail: suhsanta@catholic.ac.kr, Jo Hong, Key, and Sil Lee, Keum. GPU-based prompt gamma ray imaging from boron neutron capture therapy. United States: N. p., 2015. Web. doi:10.1118/1.4903265.
Yoon, Do-Kun, Jung, Joo-Young, Suk Suh, Tae, E-mail: suhsanta@catholic.ac.kr, Jo Hong, Key, & Sil Lee, Keum. GPU-based prompt gamma ray imaging from boron neutron capture therapy. United States. doi:10.1118/1.4903265.
Yoon, Do-Kun, Jung, Joo-Young, Suk Suh, Tae, E-mail: suhsanta@catholic.ac.kr, Jo Hong, Key, and Sil Lee, Keum. Thu . "GPU-based prompt gamma ray imaging from boron neutron capture therapy". United States. doi:10.1118/1.4903265.
@article{osti_22413373,
title = {GPU-based prompt gamma ray imaging from boron neutron capture therapy},
author = {Yoon, Do-Kun and Jung, Joo-Young and Suk Suh, Tae, E-mail: suhsanta@catholic.ac.kr and Jo Hong, Key and Sil Lee, Keum},
abstractNote = {Purpose: The purpose of this research is to perform the fast reconstruction of a prompt gamma ray image using a graphics processing unit (GPU) computation from boron neutron capture therapy (BNCT) simulations. Methods: To evaluate the accuracy of the reconstructed image, a phantom including four boron uptake regions (BURs) was used in the simulation. After the Monte Carlo simulation of the BNCT, the modified ordered subset expectation maximization reconstruction algorithm using the GPU computation was used to reconstruct the images with fewer projections. The computation times for image reconstruction were compared between the GPU and the central processing unit (CPU). Also, the accuracy of the reconstructed image was evaluated by a receiver operating characteristic (ROC) curve analysis. Results: The image reconstruction time using the GPU was 196 times faster than the conventional reconstruction time using the CPU. For the four BURs, the area under curve values from the ROC curve were 0.6726 (A-region), 0.6890 (B-region), 0.7384 (C-region), and 0.8009 (D-region). Conclusions: The tomographic image using the prompt gamma ray event from the BNCT simulation was acquired using the GPU computation in order to perform a fast reconstruction during treatment. The authors verified the feasibility of the prompt gamma ray image reconstruction using the GPU computation for BNCT simulations.},
doi = {10.1118/1.4903265},
journal = {Medical Physics},
number = 1,
volume = 42,
place = {United States},
year = {Thu Jan 15 00:00:00 EST 2015},
month = {Thu Jan 15 00:00:00 EST 2015}
}
  • Purpose: The purpose of this research is to perform the fast reconstruction of a prompt gamma ray image using a graphics processing unit (GPU) computation from boron neutron capture therapy (BNCT) simulations. Methods: To evaluate the accuracy of the reconstructed image, a phantom including four boron uptake regions (BURs) was used in the simulation. After the Monte Carlo simulation of the BNCT, the modified ordered subset expectation maximization reconstruction algorithm using the GPU computation was used to reconstruct the images with fewer projections. The computation times for image reconstruction were compared between the GPU and the central processing unit (CPU).more » Also, the accuracy of the reconstructed image was evaluated by a receiver operating characteristic (ROC) curve analysis. Results: The image reconstruction time using the GPU was 196 times faster than the conventional reconstruction time using the CPU. For the four BURs, the area under curve values from the ROC curve were 0.6726 (A-region), 0.6890 (B-region), 0.7384 (C-region), and 0.8009 (D-region). Conclusion: The tomographic image using the prompt gamma ray event from the BNCT simulation was acquired using the GPU computation in order to perform a fast reconstruction during treatment. The authors verified the feasibility of the prompt gamma ray reconstruction using the GPU computation for BNCT simulations.« less
  • Purpose of paper is to confirm the feasibility of acquisition of three dimensional single photon emission computed tomography image from boron neutron capture therapy using Monte Carlo simulation. Prompt gamma ray (478 keV) was used to reconstruct image with ordered subsets expectation maximization method. From analysis of receiver operating characteristic curve, area under curve values of three boron regions were 0.738, 0.623, and 0.817. The differences between length of centers of two boron regions and distance of maximum count points were 0.3 cm, 1.6 cm, and 1.4 cm.
  • Purpose: Purpose of paper is to confirm the feasibility of acquisition of three dimensional single photon emission computed tomography (SPECT) image from boron neutron capture therapy (BNCT) using Monte Carlo simulation. Methods: In case of simulation, the pixelated SPECT detector, collimator and phantom were simulated using Monte Carlo n particle extended (MCNPX) simulation tool. A thermal neutron source (<1 eV) was used to react with the boron uptake region (BUR) in the phantom. Each geometry had a spherical pattern, and three different BURs (A, B and C region, density: 2.08 g/cm3) were located in the middle of the brain phantom.more » The data from 128 projections for each sorting process were used to achieve image reconstruction. The ordered subset expectation maximization (OSEM) reconstruction algorithm was used to obtain a tomographic image with eight subsets and five iterations. The receiver operating characteristic (ROC) curve analysis was used to evaluate the geometric accuracy of reconstructed image. Results: The OSEM image was compared with the original phantom pattern image. The area under the curve (AUC) was calculated as the gross area under each ROC curve. The three calculated AUC values were 0.738 (A region), 0.623 (B region), and 0.817 (C region). The differences between length of centers of two boron regions and distance of maximum count points were 0.3 cm, 1.6 cm and 1.4 cm. Conclusion: The possibility of extracting a 3D BNCT SPECT image was confirmed using the Monte Carlo simulation and OSEM algorithm. The prospects for obtaining an actual BNCT SPECT image were estimated from the quality of the simulated image and the simulation conditions. When multiple tumor region should be treated using the BNCT, a reasonable model to determine how many useful images can be obtained from the SPECT could be provided to the BNCT facilities. This research was supported by the Leading Foreign Research Institute Recruitment Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Science, Information and Communication Technologies (ICT) and Future Planning (MSIP)(Grant No.200900420) and the Radiation Technology Research and Development program (Grant No.2013043498), Republic of Korea.« less
  • Purpose: The purpose of this study is to verify the acquisition of the three dimensional single photon emission computed tomography (SPECT) image using prompt gamma ray originated from proton boron fusion therapy (PBFT). Methods: The real-time imaging system during the PBFT was simulated to acquire the tomographic image of the prompt gamma ray from treated tumor by the proton boron reaction, using Monte Carlo simulation (MCNPX). We acquired percentage depth dose (PDD) of the proton beam in the water phantom including the boron uptake region (BUR), the energy spectra of the prompt gamma ray and tomographic image which can showmore » treated tumor regions. The prompt gamma ray image was reconstructed using maximum likelihood estimation maximization (MLEM) reconstruction algorithm with 64 projections. In addition, in order to evaluate the reconstructed image, the image profiles between BURs were extracted from the image. Results: In the PDD results, the Bragg-peak was amplified definitely when the proton’s maximum dose level was located at the BUR. This amplification is based on the generation of alpha particles. In addition, the prompt gamma ray peak of 719 keV was observed from the energy spectrum. Through the previous process, the tomographic image of prompt gamma ray from the BUR was reconstructed. The line profile was extracted from image including BURs, and it shows high signal-to-noise ratio. Conclusion: We confirmed that the real-time prompt gamma ray image during the PBFT was successfully deducted, and results of quantitative image analysis show good agreement with the original pattern of the BUR. This study first verified the imaging capability of the PBFT-SPECT system. In conclusion, the PBFT-SPECT system can realize the treated tumor monitoring during the PBFT.« less
  • For the boron neutron capture therapy (NCT) of deep-seated metastatic melanoma, an epithermal (up to a few keV energy) neutron beam from a reactor horizontal facility could be useful if the inherent contamination from fast neutrons and gamma rays could be minimized. Calculations for ANSTO's 10 MW research reactor HIFAR have shown that, even though a filter material such as AlF3 attenuates the fast neutron dose, the beam quality improvement is counteracted by a relative increase in the gamma dose because of the gammas arising from neutron captures in the filter material, particularly the aluminium. The aluminium gammas, most ofmore » which arise from thermal neutron capture, are hard and cannot be attenuated by lead or bismuth without comparable attenuation of the epithermal neutron flux. Addition of an absorber such as 6Li to the AlF3 filter was investigated as a means of reducing the hard gamma dose, but the improvement in beam quality was small and at considerable cost to dose intensity. Dose characteristics calculations confirmed the superiority of a tangential beam over a radial beam with better results from an unfiltered tangential beam than from an AlF3 filter in a radial beam. This study showed conclusively that assessments of filter assemblies based on the effect of individual components on either the neutron or gamma dose in isolation are inadequate. In assessing any epithermal neutron filter, thermal neutron shield, and gamma shield combination, the total effect of each on the neutron, gamma, and boron-10 dose must be considered.« less