Evaluation of a stochastic reconstruction algorithm for use in Compton camera imaging and beam range verification from secondary gamma emission during proton therapy

In this paper, we study the feasibility of using the stochastic origin ensemble (SOE) algorithm for reconstructing images of secondary gammas emitted during proton radiotherapy from data measured with a three-stage Compton camera. The purpose of this study was to evaluate the quality of the images of the gamma rays emitted during proton irradiation produced using the SOE algorithm and to measure how well the images reproduce the distal falloff of the beam. For our evaluation, we performed a Monte Carlo simulation of an ideal three-stage Compton camera positioned above and orthogonal to a proton pencil beam irradiating a tissue phantom. Scattering of beam protons with nuclei in the phantom produces secondary gamma rays, which are detected by the Compton camera and used as input to the SOE algorithm. We studied the SOE reconstructed images as a function of the number of iterations, the voxel probability parameter, and the number of detected gammas used by the SOE algorithm. We quantitatively evaluated the capabilities of the SOE algorithm by calculating and comparing the normalized mean square error (NMSE) of SOE reconstructed images. We also studied the ability of the SOE reconstructed images to predict the distal falloff of the secondary gamma production in the irradiated tissue. Our results show that the images produced with the SOE algorithm converge in ~10,000 iterations, with little improvement to the image NMSE for iterations above this number. We found that the statistical noise of the images is inversely proportional to the ratio of the number of gammas detected to the SOE voxel probability parameter value. In our study, the SOE predicted distal falloff of the reconstructed images agrees with the Monte Carlo calculated distal falloff of the gamma emission profile in the phantom to within ±0.6 mm for the positions of maximum emission (100%) and 90%, 50% and 20% distal falloff of the gamma emission profile. We conclude that the SOE algorithm is an effective method for reconstructing images of a proton pencil beam from the data collected by an ideal Compton camera and that these images accurately model the distal falloff of secondary gamma emission during proton irradiation.

[1]  Proton therapy in combination with PET as monitor: a feasibility study , 1992 .

[2]  James D. Kurfess,et al.  Considerations for the next Compton telescope mission , 2000 .

[3]  J. Archambeau,et al.  Visualization and Transport of Positron Emission from Proton Activation in vivo , 1978, Science.

[4]  J. F. Bajema,et al.  On-line monitoring and PET imaging of proton radiotherapy beams , 1992, IEEE Conference on Nuclear Science Symposium and Medical Imaging.

[5]  D Robertson,et al.  Optimizing a three-stage Compton camera for measuring prompt gamma rays emitted during proton radiotherapy , 2010, Physics in medicine and biology.

[6]  Chan Hyeong Kim,et al.  Prompt gamma measurements for locating the dose falloff region in the proton therapy , 2006 .

[7]  K. Vetter,et al.  Erratum: A gamma-ray tracking algorithm for the GRETA spectrometer , 1999 .

[8]  David W. Scott,et al.  Multivariate Density Estimation: Theory, Practice, and Visualization , 1992, Wiley Series in Probability and Statistics.

[9]  K. Vetter,et al.  GRETA: utilizing new concepts in γ-ray detection , 1999 .

[10]  S Beddar,et al.  Measurement and calculation of characteristic prompt gamma ray spectra emitted during proton irradiation , 2009, Physics in medicine and biology.

[11]  E. Parzen.  On Estimation of a Probability Density Function and Mode , 1962 .

[12]  G. W. Phillips.  Gamma-ray imaging with Compton cameras , 1995 .

[13]  R. Siddon.  Fast calculation of the exact radiological path for a three-dimensional CT array. , 1985, Medical physics.

[14]  Van-Giang Nguyen,et al.  GPU-accelerated 3D Bayesian image reconstruction from Compton scattered data. , 2011, Physics in medicine and biology.

[15]  Robert J. Ott,et al.  Gamma ray imaging with silicon detectors — A Compton camera for radionuclide imaging in medicine , 1988 .

[16]  K Parodi,et al.  Potential application of PET in quality assurance of proton therapy. , 2000, Physics in medicine and biology.

[17]  A. Celler,et al.  Fast image reconstruction for Compton camera using stochastic origin ensemble approach. , 2010, Medical physics.

[18]  A Sitek,et al.  Representation of photon limited data in emission tomography using origin ensembles , 2008, Physics in medicine and biology.

[19]  R Leahy,et al.  Three-dimensional maximum-likelihood reconstruction for an electronically collimated single-photon-emission imaging system. , 1990, Journal of the Optical Society of America. A, Optics and image science.

[20]  Uncertainty analysis of a Compton camera imaging system for radiation therapy dose reconstruction. , 2010, Medical physics.

[21]  William A. Wulf,et al.  Three-Compton telescope: theory, simulations, and performance , 2001 .

[22]  H Paganetti,et al.  Quantitative assessment of the physical potential of proton beam range verification with PET/CT , 2008, Physics in medicine and biology.

[23]  Sam Beddar,et al.  Material efficiency studies for a Compton camera designed to measure characteristic prompt gamma rays emitted during proton beam radiotherapy , 2011, Physics in medicine and biology.

[24]  Denis Dauvergne,et al.  Dose profile monitoring with carbon ions by means of prompt-gamma measurements , 2009 .

[25]  M. Rosenblatt.  Remarks on Some Nonparametric Estimates of a Density Function , 1956 .

[26]  F. Rademakers,et al.  ROOT — An object oriented data analysis framework , 1997 .

[27]  Chan Hyeong Kim,et al.  Two approaches to implementing projector–backprojector pairs for 3D reconstruction from Compton scattered data , 2007 .

[28]  H Paganetti,et al.  Monte Carlo patient study on the comparison of prompt gamma and PET imaging for range verification in proton therapy , 2011, Physics in medicine and biology.

[29]  Byung-Hwi Kang,et al.  Monte Carlo Design Study of a Gamma Detector System to Locate Distal Dose Falloff in Proton Therapy , 2009 .

[30]  W. L. Rogers,et al.  Fast algorithm for list mode back-projection of Compton scatter camera data , 1997 .

[31]  G Ciangaru,et al.  Prompt gamma-ray emission from biological tissues during proton irradiation: a preliminary study , 2009, Physics in medicine and biology.

[32]  Pierre L'Ecuyer,et al.  Maximally equidistributed combined Tausworthe generators , 1996, Math. Comput..

[33]  W. Enghardt,et al.  In-beam PET at high-energy photon beams: a feasibility study , 2006, Physics in medicine and biology.