Statistical approaches in quantitative positron emission tomography

Positron emission tomography is a medical imaging modality for producing 3D images of the spatial distribution of biochemical tracers within the human body. The images are reconstructed from data formed through detection of radiation resulting from the emission of positrons from radioisotopes tagged onto the tracer of interest. These measurements are approximate line integrals from which the image can be reconstructed using analytical inversion formulae. However these direct methods do not allow accurate modeling either of the detector system or of the inherent statistical fluctuations in the data. Here we review recent progress in developing statistical approaches to image estimation that can overcome these limitations. We describe the various components of the physical model and review different formulations of the inverse problem. The wide range of numerical procedures for solving these problems are then reviewed. Finally, we describe recent work aimed at quantifying the quality of the resulting images, both in terms of classical measures of estimator bias and variance, and also using measures that are of more direct clinical relevance.

[1]  H H Barrett,et al.  Objective assessment of image quality: effects of quantum noise and object variability. , 1990, Journal of the Optical Society of America. A, Optics and image science.

[2]  Anand Rangarajan,et al.  Bayesian reconstruction of functional images using anatomical information as priors , 1993, IEEE Trans. Medical Imaging.

[3]  James F. Young,et al.  MicroPET: a high resolution PET scanner for imaging small animals , 1996, IEEE Nuclear Science Symposium Conference Record.

[4]  R. Leahy,et al.  Fully 3 D Bayesian Image Reconstruction for the ECAT EXACT HR + 1 , 1998 .

[5]  J. M. Ollinger,et al.  Positron Emission Tomography , 2018, Handbook of Small Animal Imaging.

[6]  Jie Yao,et al.  Predicting human performance by a channelized Hotelling observer model , 1992, Optics & Photonics.

[7]  H. Malcolm Hudson,et al.  Accelerated image reconstruction using ordered subsets of projection data , 1994, IEEE Trans. Medical Imaging.

[8]  Michael I. Miller,et al.  The Use of Sieves to Stabilize Images Produced with the EM Algorithm for Emission Tomography , 1985, IEEE Transactions on Nuclear Science.

[9]  Sridhar Lakshmanan,et al.  Simultaneous Parameter Estimation and Segmentation of Gibbs Random Fields Using Simulated Annealing , 1989, IEEE Trans. Pattern Anal. Mach. Intell..

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

[11]  Richard M. Leahy,et al.  A theoretical study of the contrast recovery and variance of MAP reconstructions with applications to the selection of smoothing parameters , 1998, 1998 IEEE Nuclear Science Symposium Conference Record. 1998 IEEE Nuclear Science Symposium and Medical Imaging Conference (Cat. No.98CH36255).

[12]  Craig K. Abbey,et al.  Observer signal-to-noise ratios for the ML-EM algorithm , 1996, Medical Imaging.

[13]  Paul Kinahan,et al.  Analytic 3D image reconstruction using all detected events , 1989 .

[14]  E U Mumcuoğlu,et al.  Bayesian reconstruction of PET images: methodology and performance analysis. , 1996, Physics in medicine and biology.

[15]  Scott T. Grafton,et al.  Results of a clinical receiver operating characteristic study comparing filtered backprojection and maximum likelihood estimator images in FDG PET studies. , 1993, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[16]  Jeffrey A. Fessler,et al.  Statistical image reconstruction methods for randoms-precorrected PET scans , 1998, Medical Image Anal..

[17]  J. Mazziotta,et al.  Brain Mapping: The Methods , 2002 .

[18]  R. Leahy,et al.  Accurate geometric and physical response modelling for statistical image reconstruction in high resolution PET , 1996, 1996 IEEE Nuclear Science Symposium. Conference Record.

[19]  M E Phelps,et al.  Quantitation in Positron Emission Computed Tomography: 6. Effect of Nonuniform Resolution , 1982, Journal of computer assisted tomography.

[20]  Peter Herscovitch,et al.  An approximation formula for the variance of PET region-of-interest values , 1993, IEEE Trans. Medical Imaging.

[21]  Ken D. Sauer,et al.  ML parameter estimation for Markov random fields with applications to Bayesian tomography , 1998, IEEE Trans. Image Process..

[22]  D. Snyder,et al.  Corrections for accidental coincidences and attenuation in maximum-likelihood image reconstruction for positron-emission tomography. , 1991, IEEE transactions on medical imaging.

[23]  Y. Censor Finite series-expansion reconstruction methods , 1983, Proceedings of the IEEE.

[24]  Alvaro R. De Pierro,et al.  A row-action alternative to the EM algorithm for maximizing likelihood in emission tomography , 1996, IEEE Trans. Medical Imaging.

[25]  G. L. Brownell,et al.  Estimation of the Local Statistical Noise in Emission Computed Tomography , 1982, IEEE Transactions on Medical Imaging.

[26]  I. Weir Fully Bayesian Reconstructions from Single-Photon Emission Computed Tomography Data , 1997 .

[27]  Stuart Geman,et al.  Statistical methods for tomographic image reconstruction , 1987 .

[28]  S C Strother,et al.  The convergence of object dependent resolution in maximum likelihood based tomographic image reconstruction. , 1993, Physics in medicine and biology.

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

[30]  J M Hoffman,et al.  Semiquantitative and visual analysis of FDG-PET images in pulmonary abnormalities. , 1994, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[31]  Jeffrey A. Fessler Penalized weighted least-squares image reconstruction for positron emission tomography , 1994, IEEE Trans. Medical Imaging.

[32]  Ken D. Sauer,et al.  A local update strategy for iterative reconstruction from projections , 1993, IEEE Trans. Signal Process..

[33]  Robert M. Lewitt,et al.  Practical considerations for 3-D image reconstruction using spherically symmetric volume elements , 1996, IEEE Trans. Medical Imaging.

[34]  G. McLachlan,et al.  The EM algorithm and extensions , 1996 .

[35]  V J Lowe,et al.  Detection of primary and recurrent lung cancer by means of F-18 fluorodeoxyglucose positron emission tomography (FDG PET). , 1995, The Journal of thoracic and cardiovascular surgery.

[36]  Kimberlee J. Kearfott A Statistical Model for Positron Emission Tomography: Comment , 1985 .

[37]  Jeffrey A. Fessler Mean and variance of implicitly defined biased estimators (such as penalized maximum likelihood): applications to tomography , 1996, IEEE Trans. Image Process..

[38]  Jorge Llacer,et al.  Optimization of Maximum Likelihood Estimator Images for PET: II. Statistical Analysis of Human Brain FDG Studies , 1991 .

[39]  Kimberlee J. Kearfott Comment: Practical Considerations , 1985 .

[40]  V. Johnson A Model for Segmentation and Analysis of Noisy Images , 1994 .

[41]  Ken D. Sauer,et al.  A generalized Gaussian image model for edge-preserving MAP estimation , 1993, IEEE Trans. Image Process..

[42]  Simon R. Cherry,et al.  Fully 3D Bayesian image reconstruction for the ECAT EXACT HR , 1997 .

[43]  Kenneth F. Koral,et al.  Object-dependent performance comparison of two iterative reconstruction algorithms , 1988 .

[44]  Richard E. Carson,et al.  Precision and accuracy of regional radioactivity quantitation using the maximum likelihood EM reconstruction algorithm , 1994, IEEE Trans. Medical Imaging.

[45]  Michel Defrise,et al.  Exact and approximate rebinning algorithms for 3-D PET data , 1997, IEEE Transactions on Medical Imaging.

[46]  Samuel Matej,et al.  A comparison of transform and iterative reconstruction techniques for a volume-imaging PET scanner with a large axial acceptance angle , 1994, Proceedings of 1994 IEEE Nuclear Science Symposium - NSS'94.

[47]  Simon R. Cherry,et al.  Fast gradient-based methods for Bayesian reconstruction of transmission and emission PET images , 1994, IEEE Trans. Medical Imaging.

[48]  Richard M. Leahy,et al.  Incorporation of Anatomical MR Data for Improved Dunctional Imaging with PET , 1991, IPMI.

[49]  Richard M. Leahy,et al.  MAP image reconstruction using intensity and line processes for emission tomography data , 1991, Electronic Imaging.

[50]  Valen E. Johnson,et al.  A note on stopping rules in EM-ML reconstructions of ECT images , 1994, IEEE Trans. Medical Imaging.

[51]  K. Lange,et al.  A Theoretical Study of Some Maximum Likelihood Algorithms for Emission and Transmission Tomography , 1987, IEEE Transactions on Medical Imaging.

[52]  G. Wahba Smoothing noisy data with spline functions , 1975 .

[53]  Á. R. De Pierro,et al.  A modified expectation maximization algorithm for penalized likelihood estimation in emission tomography. , 1995, IEEE transactions on medical imaging.

[54]  Anand Rangarajan,et al.  Bayesian Reconstruction of Functional Images Using Registered Anatomical Images as Priors , 1991, IPMI.

[55]  Charles L. Byrne,et al.  Convergent block-iterative algorithms for image reconstruction from inconsistent data , 1997, IEEE Trans. Image Process..

[56]  J. Varah Pitfalls in the Numerical Solution of Linear Ill-Posed Problems , 1981 .

[57]  Anand Rangarajan,et al.  Bayesian image reconstruction in SPECT using higher order mechanical models as priors , 1995, IEEE Trans. Medical Imaging.

[58]  J. Fessler,et al.  Objective functions for tomographic reconstruction from randoms-precorrected PET scans , 1996, 1996 IEEE Nuclear Science Symposium. Conference Record.

[59]  Sylvia Richardson,et al.  Markov Chain Monte Carlo in Practice , 1997 .

[60]  J.A. Fessler,et al.  Fully 3D PET image reconstruction using a Fourier preconditioned conjugate-gradient algorithm , 1996, 1996 IEEE Nuclear Science Symposium. Conference Record.

[61]  C T Chen,et al.  Bayesian restoration of PET images using Gibbs priors. , 1991, Progress in clinical and biological research.

[62]  Jun Zhang,et al.  Maximum-likelihood parameter estimation for unsupervised stochastic model-based image segmentation , 1994, IEEE Trans. Image Process..

[63]  Linda Kaufman,et al.  Implementing and Accelerating the EM Algorithm for Positron Emission Tomography , 1987, IEEE Transactions on Medical Imaging.

[64]  Jeffrey A. Fessler,et al.  Ieee Transactions on Image Processing: to Appear Hybrid Poisson/polynomial Objective Functions for Tomographic Image Reconstruction from Transmission Scans , 2022 .

[65]  Alfred O. Hero,et al.  Exploring estimator bias-variance tradeoffs using the uniform CR bound , 1996, IEEE Trans. Signal Process..

[66]  J L Sunshine,et al.  Evaluation of breast masses and axillary lymph nodes with [F-18] 2-deoxy-2-fluoro-D-glucose PET. , 1993, Radiology.

[67]  Albert Macovski,et al.  A Maximum Likelihood Approach to Emission Image Reconstruction from Projections , 1976, IEEE Transactions on Nuclear Science.

[68]  T. Hebert,et al.  A generalized EM algorithm for 3-D Bayesian reconstruction from Poisson data using Gibbs priors. , 1989, IEEE transactions on medical imaging.

[69]  Jeffrey A. Fessler,et al.  Grouped-coordinate ascent algorithms for penalized-likelihood transmission image reconstruction , 1997, IEEE Transactions on Medical Imaging.

[70]  R.M. Leahy,et al.  Fully 3D Bayesian image reconstruction for the ECAT EXACT HR+ , 1997, 1997 IEEE Nuclear Science Symposium Conference Record.

[71]  E. Veklerov,et al.  Stopping Rule for the MLE Algorithm Based on Statistical Hypothesis Testing , 1987, IEEE Transactions on Medical Imaging.

[72]  J. Besag,et al.  Bayesian Computation and Stochastic Systems , 1995 .

[73]  Ronald J. Jaszczak,et al.  Bayesian reconstruction and use of anatomical a priori information for emission tomography , 1996, IEEE Trans. Medical Imaging.

[74]  G. Gindi,et al.  Noise analysis of MAP-EM algorithms for emission tomography. , 1997, Physics in medicine and biology.

[75]  S. Holte,et al.  Iterative image reconstruction for positron emission tomography: a study of convergence and quantitation problems , 1990 .

[76]  Alfred O. Hero,et al.  Ieee Transactions on Image Processing: to Appear Penalized Maximum-likelihood Image Reconstruction Using Space-alternating Generalized Em Algorithms , 2022 .

[77]  G. Krishna,et al.  Vector-extrapolated fast maximum likelihood estimation algorithms for emission tomography , 1992, IEEE Trans. Medical Imaging.

[78]  Per Christian Hansen,et al.  Analysis of Discrete Ill-Posed Problems by Means of the L-Curve , 1992, SIAM Rev..

[79]  Robert M. Lewitt,et al.  Accelerated Iterative Reconstruction for Positron Emission Tomography Based on the EM Algorithm for Maximum Likelihood Estimation , 1986, IEEE Transactions on Medical Imaging.

[80]  Yuchen Yan,et al.  A system for the 3D reconstruction of retracted-septa PET data using the EM algorithm , 1994, Proceedings of 1994 IEEE Nuclear Science Symposium - NSS'94.

[81]  C. Geyer,et al.  Constrained Monte Carlo Maximum Likelihood for Dependent Data , 1992 .

[82]  P. Green Bayesian reconstructions from emission tomography data using a modified EM algorithm. , 1990, IEEE transactions on medical imaging.

[83]  Donald W. Wilson,et al.  Noise properties of the EM algorithm. I. Theory , 1994 .

[84]  J. Besag On the Statistical Analysis of Dirty Pictures , 1986 .

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

[86]  Peter Craven,et al.  Smoothing noisy data with spline functions , 1978 .

[87]  R. Kessler,et al.  Analysis of emission tomographic scan data: limitations imposed by resolution and background. , 1984, Journal of computer assisted tomography.

[88]  H. Mayberg,et al.  Correction of PET Data for Partial Volume Effects in Human Cerebral Cortex by MR Imaging , 1990, Journal of computer assisted tomography.

[89]  Alvaro R. De Pierro,et al.  A modified expectation maximization algorithm for penalized likelihood estimation in emission tomography , 1995, IEEE Trans. Medical Imaging.

[90]  K. Lange Convergence of EM image reconstruction algorithms with Gibbs smoothing. , 1990, IEEE transactions on medical imaging.

[91]  George M. Nixon,et al.  Practical Considerations of , 1945 .

[92]  Grace Wahba,et al.  Spline Models for Observational Data , 1990 .

[93]  Jeffrey A. Fessler,et al.  Spatial resolution properties of penalized-likelihood image reconstruction: space-invariant tomographs , 1996, IEEE Trans. Image Process..

[94]  B. Tsui,et al.  Noise properties of the EM algorithm: II. Monte Carlo simulations. , 1994, Physics in medicine and biology.

[95]  Richard M. Leahy,et al.  A theoretical study of the contrast recovery and variance of MAP reconstructions from PET data , 1999, IEEE Transactions on Medical Imaging.

[96]  D. Newport,et al.  A Single Scatter Simulation Technique for Scatter Correction in 3D PET , 1996 .

[97]  Zhenyu Zhou,et al.  Approximate maximum likelihood hyperparameter estimation for Gibbs priors , 1995, Proceedings., International Conference on Image Processing.

[98]  Walter R. Gilks,et al.  MCMC in image analysis , 1995 .

[99]  R. Huesman A new fast algorithm for the evaluation of regions of interest and statistical uncertainty in computed tomography. , 1984, Physics in medicine and biology.

[100]  Linda Kaufman,et al.  Maximum likelihood, least squares, and penalized least squares for PET , 1993, IEEE Trans. Medical Imaging.

[101]  K. Coakley A cross-validation procedure for stopping the EM algorithm and deconvolution of neutron depth profiling spectra , 1991 .

[102]  R. Leahy,et al.  High-resolution 3D Bayesian image reconstruction using the microPET small-animal scanner. , 1998, Physics in medicine and biology.

[103]  C E Metz,et al.  An evaluation of maximum likelihood-expectation maximization reconstruction for SPECT by ROC analysis. , 1992, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[104]  Ken D. Sauer,et al.  A unified approach to statistical tomography using coordinate descent optimization , 1996, IEEE Trans. Image Process..

[105]  B. F. Logan,et al.  The Fourier reconstruction of a head section , 1974 .

[106]  J. R. Baker,et al.  Spatially variant tomographic imaging: estimation, identification, and optimization , 1992 .

[107]  R. Leahy,et al.  Comparing lesion detection performance for PET image reconstruction algorithms: A case study , 1996, 1996 IEEE Nuclear Science Symposium. Conference Record.

[108]  Zang-Hee Cho,et al.  Parallelization of the EM algorithm for 3-D PET image reconstruction , 1991 .

[109]  Richard M. Leahy,et al.  Statistic-based MAP image-reconstruction from Poisson data using Gibbs priors , 1992, IEEE Trans. Signal Process..

[110]  J. M. Ollinger,et al.  Model based scatter correction in three dimensions (positron emission tomography) , 1992, IEEE Conference on Nuclear Science Symposium and Medical Imaging.

[111]  Michael A. King,et al.  Comparison of the channelized Hotelling and human observers for lesion detection in hepatic SPECT imaging , 1997, Medical Imaging.

[112]  S. Derenzo,et al.  Application of mathematical removal of positron range blurring in positron emission tomography , 1990 .

[113]  Donald Geman,et al.  Stochastic relaxation, Gibbs distributions, and the Bayesian restoration of images , 1984 .

[114]  Simon R. Cherry,et al.  18 – Imaging Brain Function with Positron Emission Tomography , 2002 .

[115]  H. Herzog,et al.  High resolution and better quantification by tube of response modelling in 3D PET reconstruction , 1996, 1996 IEEE Nuclear Science Symposium. Conference Record.

[116]  Gabor T. Herman,et al.  Bayesian Image Reconstruction Using a High-Order Interacting MRF Model , 1995, ICIAP.

[117]  Wesley E. Snyder,et al.  Mean field annealing: a formalism for constructing GNC-like algorithms , 1992, IEEE Trans. Neural Networks.

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

[119]  A. Chatziioannou,et al.  Comparing lesion detection performance for PET image reconstruction algorithms: A case study , 1996, 1996 IEEE Nuclear Science Symposium. Conference Record.

[120]  F. O’Sullivan,et al.  Variability Assessment in Positron Emission Tomography and Related Generalized Deconvolution Models , 1998 .

[121]  L G Strauss,et al.  The applications of PET in clinical oncology. , 1991, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[122]  Ronald J. Jaszczak,et al.  Fully Bayesian estimation of Gibbs hyperparameters for emission computed tomography data , 1997, IEEE Transactions on Medical Imaging.

[123]  J. D. Wilson,et al.  A smoothed EM approach to indirect estimation problems, with particular reference to stereology and emission tomography , 1990 .

[124]  D. M. Titterington,et al.  A Study of Methods of Choosing the Smoothing Parameter in Image Restoration by Regularization , 1991, IEEE Trans. Pattern Anal. Mach. Intell..

[125]  S. Strother,et al.  Practical tradeoffs between noise, quantitation, and number of iterations for maximum likelihood-based reconstructions. , 1991, IEEE transactions on medical imaging.

[126]  Richard M. Leahy,et al.  Fast computation of the covariance of MAP reconstructions of PET images , 1999, Medical Imaging.

[127]  M. P. Beddoes,et al.  Effects of Detector Wobble Motion on Image Noise in Positron Emission Tomography , 1985, IEEE Transactions on Medical Imaging.