Probabilistic models and numerical calculation of system matrix and sensitivity in list-mode MLEM 3D reconstruction of Compton camera images

This paper addresses the problem of evaluating the system matrix and the sensitivity for iterative reconstruction in Compton camera imaging. Proposed models and numerical calculation strategies are compared through the influence they have on the three-dimensional reconstructed images. The study attempts to address four questions. First, it proposes an analytic model for the system matrix. Second, it suggests a method for its numerical validation with Monte Carlo simulated data. Third, it compares analytical models of the sensitivity factors with Monte Carlo simulated values. Finally, it shows how the system matrix and the sensitivity calculation strategies influence the quality of the reconstructed images.

[1]  Wei Chang,et al.  Dependence of angular uncertainties on the energy resolution of Compton cameras , 1997, 1997 IEEE Nuclear Science Symposium Conference Record.

[2]  Bruce Smith,et al.  Reconstruction methods and completeness conditions for two Compton data models. , 2005, Journal of the Optical Society of America. A, Optics, image science, and vision.

[3]  Wolfgang Enghardt,et al.  Test of Compton camera components for prompt gamma imaging at the ELBE bremsstrahlung beam , 2014 .

[4]  A. Zoglauer.  First light for the next generation of compton and pair telescopes : development of new techniques for the data analysis of combined compton and pair telescopes and their application to the MEGA prototype , 2005 .

[5]  K. Lange,et al.  EM reconstruction algorithms for emission and transmission tomography. , 1984, Journal of computer assisted tomography.

[6]  Jeffrey A. Fessler,et al.  List-mode maximum likelihood reconstruction of Compton scatter camera images in nuclear medicine , 1998, 1998 IEEE Nuclear Science Symposium Conference Record. 1998 IEEE Nuclear Science Symposium and Medical Imaging Conference (Cat. No.98CH36255).

[7]  N. Clinthorne,et al.  Modeling of System Response in List Mode EM Reconstruction of Compton Scatter Camera Images , 2000 .

[8]  Liji Cao,et al.  List-mode maximum-likelihood reconstruction for the ClearPEM system , 2011, 2011 IEEE Nuclear Science Symposium Conference Record.

[9]  Jeffrey A. Fessler,et al.  Improved modeling of system response in list mode EM reconstruction of Compton scatter camera images , 1999 .

[10]  Voichita Maxim.  Filtered Backprojection Reconstruction and Redundancy in Compton Camera Imaging , 2014, IEEE Transactions on Image Processing.

[11]  J. Sullivan,et al.  Extended radiation source imaging with a prototype Compton imager. , 2009, Applied radiation and isotopes : including data, instrumentation and methods for use in agriculture, industry and medicine.

[12]  G T Gullberg,et al.  Application of spherical harmonics to image reconstruction for the Compton camera. , 1998, Physics in medicine and biology.

[13]  A. Dell'Acqua,et al.  Geant4 - A simulation toolkit , 2003 .

[14]  C. Lacasta,et al.  A Compton imaging algorithm for on-line monitoring in hadron therapy , 2011, Medical Imaging.

[15]  L. Parra,et al.  Reconstruction of cone-beam projections from Compton scattered data , 1999, 1999 IEEE Nuclear Science Symposium. Conference Record. 1999 Nuclear Science Symposium and Medical Imaging Conference (Cat. No.99CH37019).

[16]  Georg Weidenspointner,et al.  Cosima — The cosmic simulator of MEGAlib , 2009, 2009 IEEE Nuclear Science Symposium Conference Record (NSS/MIC).

[17]  L. Shepp,et al.  A Statistical Model for Positron Emission Tomography , 1985 .

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

[19]  Denis Dauvergne,et al.  Development of a Compton camera for medical applications based on silicon strip and scintillation detectors , 2015 .

[20]  Sam Beddar,et al.  Evaluation of a stochastic reconstruction algorithm for use in Compton camera imaging and beam range verification from secondary gamma emission during proton therapy , 2012, Physics in medicine and biology.

[21]  Robert Andritschke,et al.  MEGAlib – The Medium Energy Gamma-ray Astronomy Library , 2006 .

[22]  Alfred O. Hero,et al.  3D image reconstruction for a Compton SPECT camera model , 1999 .

[23]  Zhong He,et al.  Gamma-ray energy-imaging integrated deconvolution , 2005, IEEE Nuclear Science Symposium Conference Record, 2005.

[24]  T. Tomitani,et al.  Image reconstruction from limited angle Compton camera data , 2002, Physics in medicine and biology.

[25]  Zhong He,et al.  Gamma-ray energy-imaging integrated spectral deconvolution , 2007 .

[26]  W. L. Rogers,et al.  Fast algorithm for list mode back-projection of Compton scatter camera data , 1997, 1997 IEEE Nuclear Science Symposium Conference Record.

[27]  T. Tomitani,et al.  An analytical image reconstruction algorithm to compensate for scattering angle broadening in Compton cameras. , 2003, Physics in medicine and biology.

[28]  D. Rubin,et al.  Maximum likelihood from incomplete data via the EM - algorithm plus discussions on the paper , 1977 .

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

[30]  Denis Dauvergne,et al.  Low Statistics Reconstruction of the Compton Camera Point Spread Function in 3D Prompt-$\gamma$ Imaging of Ion Beam Therapy , 2013, IEEE Transactions on Nuclear Science.

[31]  Yves Goussard,et al.  Three-dimensional edge-preserving image enhancement for computed tomography , 2003, IEEE Transactions on Medical Imaging.

[32]  L. Shepp,et al.  Maximum Likelihood Reconstruction for Emission Tomography , 1983, IEEE Transactions on Medical Imaging.

[33]  Rémy Prost,et al.  Analytical inversion of the Compton transform using the full set of available projections , 2009 .

[34]  D. McGinnis,et al.  Beam decelerations with variable momentum compaction in the Fermilab Antiproton Accumulator , 2003 .

[35]  P Solevi,et al.  Noise evaluation of Compton camera imaging for proton therapy , 2015, Physics in medicine and biology.

[36]  Andreas Zoglauer,et al.  Design, implementation, and optimization of MEGAlib's image reconstruction tool Mimrec , 2011 .

[37]  Van-Giang Nguyen,et al.  Three-dimensional edge-preserving regularization for Compton camera reconstruction , 2008, 2008 IEEE Nuclear Science Symposium Conference Record.

[38]  G L Zeng,et al.  Two-dimensional iterative region-of-interest (ROI) reconstruction from truncated projection data. , 2007, Medical physics.

[39]  A. Zoglauer,et al.  A Tracking Compton-Scattering Imaging System for Hadron Therapy Monitoring , 2010, IEEE Transactions on Nuclear Science.

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

[41]  Philip J. Bones,et al.  Towards direct reconstruction from a gamma camera based on Compton scattering , 1994, IEEE Trans. Medical Imaging.

[42]  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.