PET Image Reconstruction: A Robust State Space Approach

Statistical iterative reconstruction algorithms have shown improved image quality over conventional nonstatistical methods in PET by using accurate system response models and measurement noise models. Strictly speaking, however, PET measurements, pre-corrected for accidental coincidences, are neither Poisson nor Gaussian distributed and thus do not meet basic assumptions of these algorithms. In addition, the difficulty in determining the proper system response model also greatly affects the quality of the reconstructed images. In this paper, we explore the usage of state space principles for the estimation of activity map in tomographic PET imaging. The proposed strategy formulates the organ activity distribution through tracer kinetics models, and the photon-counting measurements through observation equations, thus makes it possible to unify the dynamic reconstruction problem and static reconstruction problem into a general framework. Further, it coherently treats the uncertainties of the statistical model of the imaging system and the noisy nature of measurement data. Since H(infinity) filter seeks minimummaximum-error estimates without any assumptions on the system and data noise statistics, it is particular suited for PET image reconstruction where the statistical properties of measurement data and the system model are very complicated. The performance of the proposed framework is evaluated using Shepp-Logan simulated phantom data and real phantom data with favorable results.

[1]  Zhengrong Liang,et al.  A theoretically based pre-reconstructing filter for spatio-temporal noise reduction in gated cardiac SPECT , 2000, 2000 IEEE Nuclear Science Symposium. Conference Record (Cat. No.00CH37149).

[2]  Patrick Dupont,et al.  Maximum-likelihood expectation-maximization reconstruction of sinograms with arbitrary noise distribution using NEC-transformations , 2001, IEEE Transactions on Medical Imaging.

[3]  M. Rafecas,et al.  Effect of noise in the probability matrix used for statistical reconstruction of PET data , 2002, 2002 IEEE Nuclear Science Symposium Conference Record.

[4]  E. Hoffman,et al.  Quantitation in Positron Emission Computed Tomography: 4. Effect of Accidental Coincidences , 1981, Journal of computer assisted tomography.

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

[6]  Robert M. Lewitt,et al.  Overview of methods for image reconstruction from projections in emission computed tomography , 2003, Proc. IEEE.

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

[8]  Xuemin Shen,et al.  A dynamic system approach to speech enhancement using the H∞ filtering algorithm , 1999, IEEE Trans. Speech Audio Process..

[9]  E. Veklerov,et al.  MLE reconstruction of a brain phantom using a Monte Carlo transition matrix and a statistical stopping rule , 1988 .

[10]  Steve R. Gunn,et al.  Tracer kinetic modeling via basis pursuit , 2002 .

[11]  Jeffrey A. Fessler,et al.  New Statistical Models for Randoms-Precorrected PET Scans , 1997, IPMI.

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

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

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

[15]  R. Lecomte,et al.  Detector response models for statistical iterative image reconstruction in high resolution PET , 1998, 1998 IEEE Nuclear Science Symposium Conference Record. 1998 IEEE Nuclear Science Symposium and Medical Imaging Conference (Cat. No.98CH36255).

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