State-Space Reconstruction of Pet Parametric Maps

The primary goal of dynamic positron emission tomography (PET) is to quantify the physiological and biological processes through tracer kinetics analysis. However, the process is difficult and complicated because of the compromising imaging data quality, i.e. either longer scans with good counting statistics but poor temporal resolution, or noisy shorter scans with good temporal resolution. In this paper, we explore the usage of state space principles for physiological parameter estimation in dynamic PET imaging. The system equation is constructed from particular tracer kinetic models, with the number and relationship between tissue compartments dictated by the physiological and biochemical properties of the process under study. And the observation equation on measurement data is formed based on the specific types of imaging or image-derived data. Once the Poisson distributed PET data are converted to Gaussian ones through the Anscombe transformation, an extended Kalman filter is adopted to estimate the tracer kinetics parameters from the system and observations. More appropriate estimation strategies which better take care of the PET statistics are also under development. The framework is tested on simulated digital phantom data, and the results are of sufficient accuracy and robustness.

[1]  D. Bailey,et al.  The direct calculation of parametric images from dynamic PET data using maximum-likelihood iterative reconstruction. , 1997 .

[2]  M. Nørgaard,et al.  Advances in Derivative-Free State Estimation for Nonlinear Systems , 1998 .

[3]  Ion Chryssoverghi,et al.  Approximate Gradient Projection Method with Runge-Kutta Schemes for Optimal Control Problems , 2004, Comput. Optim. Appl..

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

[5]  Huafeng Liu,et al.  Uncertainty penalized weighted least squares framework for PET reconstruction under uncertain system models , 2005, IEEE International Conference on Image Processing 2005.

[6]  David H. Anderson Compartmental analysis in Biology and Medicine : John A. Jacquez, second edition, The University of Michigan Press, Ann Arbor, 1985, x+560pp., $39.95 , 1986 .

[7]  David Dagan Feng,et al.  Simultaneous estimation of physiological parameters and the input function - in vivo PET data , 2001, IEEE Transactions on Information Technology in Biomedicine.

[8]  Richard E. Carson,et al.  Comment: The EM Parametric Image Reconstruction Algorithm , 1985 .

[9]  Ken D. Sauer,et al.  Direct reconstruction of kinetic parameter images from dynamic PET data , 2003, The Thrity-Seventh Asilomar Conference on Signals, Systems & Computers, 2003.

[10]  Finbarr O'Sullivan,et al.  Use of ridge regression for improved estimation of kinetic constants from PET data , 1999, IEEE Transactions on Medical Imaging.

[11]  Huafeng Liu,et al.  PET Image Reconstruction: A Robust State Space Approach , 2005, IPMI.

[12]  F. J. Anscombe,et al.  THE TRANSFORMATION OF POISSON, BINOMIAL AND NEGATIVE-BINOMIAL DATA , 1948 .

[13]  Claudio Cobelli,et al.  Tracer Kinetics In Biomedical Research , 2000 .

[14]  Ken D. Sauer,et al.  Direct reconstruction of kinetic parameter images from dynamic PET data , 2005, IEEE Transactions on Medical Imaging.

[15]  Richard E. Carson,et al.  The EM parametric image reconstruction algorithm , 1985 .