Reconstruction of Emission Tomography Data Using Origin Ensembles

A new statistical reconstruction method based on origin ensembles (OE) for emission tomography (ET) is examined. Using a probability density function (pdf) derived from first principles, an ensemble expectation of numbers of detected event origins per voxel is determined. These numbers divided by sensitivities of voxels and acquisition time provide OE estimates of the voxel activities. The OE expectations are shown to be the same as expectations calculated using the complete-data space. The properties of the OE estimate are examined. It is shown that OE estimate approximates maximum likelihood (ML) estimate for conditions usually achieved in practical applications in emission tomography. Three numerical experiments with increasing complexity are used to validate theoretical findings and demonstrate similarities of ML and OE estimates. Recommendations for achieving improved accuracy and speed of OE reconstructions are provided.

[1]  R. Casey,et al.  Advances in Pattern Recognition , 1971 .

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

[3]  Gabor T. Herman,et al.  Fundamentals of Computerized Tomography: Image Reconstruction from Projections , 2009, Advances in Pattern Recognition.

[4]  T. K. Narayan,et al.  Evaluation of task-oriented performance of several fully 3D PET reconstruction algorithms. , 1994, Physics in medicine and biology.

[5]  R. Carson A Maximum Likelihood Method for Region-of-Interest Evaluation in Emission Tomography , 1986, Journal of computer assisted tomography.

[6]  Gabor T. Herman,et al.  Image reconstruction from projections : the fundamentals of computerized tomography , 1980 .

[7]  Anna Celler,et al.  Stochastic image reconstruction method for Compton camera , 2009, 2009 IEEE Nuclear Science Symposium Conference Record (NSS/MIC).

[8]  Jinyi Qi,et al.  Iterative reconstruction techniques in emission computed tomography , 2006, Physics in medicine and biology.

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

[10]  Arkadiusz Sitek,et al.  All-count reconstruction: Effects of energy band width , 2009 .

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

[12]  W. K. Hastings,et al.  Monte Carlo Sampling Methods Using Markov Chains and Their Applications , 1970 .

[13]  M. S. Kaplan,et al.  Preliminary Experience With The Photon History Generator Module Of A Public-domain Simulation System For Emission Tomography , 1993, 1993 IEEE Conference Record Nuclear Science Symposium and Medical Imaging Conference.

[14]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[15]  Madhav P. Desai,et al.  On the convergence of reversible Markov chains , 1993 .

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

[17]  Gerard T. Barkema,et al.  Monte Carlo Methods in Statistical Physics , 1999 .

[18]  William H. Press,et al.  Numerical recipes in C , 2002 .

[19]  Hans Herzog,et al.  PET Reconstruction Software Toolkit - PRESTO a novel, universal C++ library for fast, iterative, fully 3D PET image reconstruction using highly compressed, memory-resident system matrices , 2008, 2008 IEEE Nuclear Science Symposium Conference Record.

[20]  Arkadiusz Sitek Image reconstruction in emission tomography using canonical origin ensembles , 2010, IEEE Nuclear Science Symposuim & Medical Imaging Conference.

[21]  Michel Defrise,et al.  REVIEW: Image reconstruction , 2006 .

[22]  M. Defrise,et al.  Image reconstruction. , 2006, Physics in medicine and biology.

[23]  D. Doria,et al.  An electronically collimated gamma camera for single photon emission computed tomography. Part II: Image reconstruction and preliminary experimental measurements. , 1983, Medical physics.

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

[25]  Leon Kaufman,et al.  A general algorithm for oblique image reconstruction , 1990, IEEE Computer Graphics and Applications.

[26]  M. Tamal,et al.  Noise Properties of Four Strategies for Incorporation of Scatter and Attenuation Information in PET Reconstruction Using the EM-ML Algorithm , 2006, IEEE Transactions on Nuclear Science.

[27]  Ronald H. Huesman,et al.  Tomographic reconstruction using an adaptive tetrahedral mesh defined by a point cloud , 2006, IEEE Transactions on Medical Imaging.

[28]  Manbir Singh,et al.  An electronically collimated gamma camera for single photon emission computed tomography. Part I: Theoretical considerations and design criteria , 1983 .