Continuous Time Dynamic PET Imaging Using List Mode Data

We describe a method for computing a continuous time estimate of dynamic changes in tracer density using list mode PET data. The tracer density in each voxel is modeled as an inhomogeneous Poisson process whose rate function can be represented using a cubic B-spline basis. An estimate of these rate functions is obtained by maximizing the likelihood of the arrival times of each detected photon pair over the control vertices of the spline. By resorting the list mode data into a standard sinogram plus a "timogram" that retains the arrival times of each of the events, we are able to perform efficient computation that exploits the symmetry inherent in the ordered sinogram. The maximum likelihood estimator uses quadratic temporal and spatial smoothness penalties and an additional penalty term to enforce non-negativity. Corrections for scatter and randoms are described and the results of studies using simulated and human data are included.

[1]  Carl de Boor,et al.  A Practical Guide to Splines , 1978, Applied Mathematical Sciences.

[2]  Lucas C. Parra,et al.  List-Mode Likelihood EM Algorithm and Noise Estimation Demonstrated on 2D-PET , 1998, IEEE Trans. Medical Imaging.

[3]  Chong Gu,et al.  Smoothing spline density estimation: theory , 1993 .

[4]  D. Snyder.  Utilizing Side Information in Emission Tomography , 1984, IEEE Transactions on Nuclear Science.

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

[6]  Donald L. Snyder,et al.  Random Point Processes in Time and Space , 1991 .

[7]  David A. Huffman,et al.  A method for the construction of minimum-redundancy codes , 1952, Proceedings of the IRE.

[8]  Peter Herscovitch,et al.  Brain blood flow measured with intravenous H/sub 2//sup 15/O. I. Theory and error analysis , 1983 .

[9]  Donald L. Snyder,et al.  Parameter Estimation for Dynamic Studies in Emission-Tomography Systems Having List-Mode Data , 1984, IEEE Transactions on Nuclear Science.

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

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

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

[13]  H H Barrett,et al.  List-mode likelihood. , 1997, Journal of the Optical Society of America. A, Optics, image science, and vision.

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

[15]  G. Wahba.  Interpolating Spline Methods for Density Estimation I. Equi-Spaced Knots , 1975 .

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

[17]  B. Barsky,et al.  An Introduction to Splines for Use in Computer Graphics and Geometric Modeling , 1987 .

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

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

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

[21]  M.E. Casey,et al.  Evaluation of simulation-based scatter correction for 3D PET cardiac imaging , 1995, 1995 IEEE Nuclear Science Symposium and Medical Imaging Conference Record.

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

[23]  D. Newport,et al.  Evaluation of simulation-based scatter correction for 3-D PET cardiac imaging , 1995 .

[24]  Finbarr O'Sullivan,et al.  Imaging radiotracer model parameters in PET: a mixture analysis approach , 1993, IEEE Trans. Medical Imaging.

[25]  M. Raichle,et al.  Brain blood flow measured with intravenous H2(15)O. I. Theory and error analysis. , 1983, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.