Kinetic parameter estimation from SPECT cone-beam projection measurements.

Kinetic parameters are commonly estimated from dynamically acquired nuclear medicine data by first reconstructing a dynamic sequence of images and subsequently fitting the parameters to time-activity curves generated from regions of interest overlaid upon the image sequence. Biased estimates can result from images reconstructed using inconsistent projections of a time-varying distribution of radiopharmaceutical acquired by a rotating SPECT system. If the SPECT data are acquired using cone-beam collimators wherein the gantry rotates so that the focal point of the collimators always remains in a plane, additional biases can arise from images reconstructed using insufficient, as well as truncated, projection samples. To overcome these problems we have investigated the estimation of kinetic parameters directly from SPECT cone-beam projection data by modelling the data acquisition process. To accomplish this it was necessary to parametrize the spatial and temporal distribution of the radiopharmaceutical within the SPECT field of view. In a simulated chest image volume, kinetic parameters were estimated for simple one-compartment models for four myocardial regions of interest. Myocardial uptake and washout parameters estimated by conventional analysis of noiseless simulated cone-beam data had biases ranging between 3-26% and 0-28%, respectively. Parameters estimated directly from the noiseless projection data were unbiased as expected, since the model used for fitting was faithful to the simulation. Statistical uncertainties of parameter estimates for 10,000,000 events ranged between 0.2-9% for the uptake parameters and between 0.3-6% for the washout parameters.

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

[2]  L. Feldkamp,et al.  Practical cone-beam algorithm , 1984 .

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

[4]  Douglas M. Bates,et al.  Nonlinear Regression Analysis and Its Applications , 1988 .

[5]  A Geissbuhler,et al.  Implementation of three-dimensional image reconstruction for multi-ring positron tomographs. , 1990, Physics in medicine and biology.

[6]  Andreas Robert Formiconi,et al.  Least squares algorithm for region-of-interest evaluation in emission tomography , 1993, IEEE Trans. Medical Imaging.

[7]  G. Gullberg,et al.  Kinetic modeling of teboroxime using dynamic SPECT imaging of a canine model. , 1994, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[8]  Ronald H. Huesman,et al.  Using linear time-invariant system theory to estimate kinetic parameters directly from projection measurements , 1994, Proceedings of 1994 IEEE Nuclear Science Symposium - NSS'94.

[9]  Alfred O. Hero,et al.  Model-based estimation for dynamic cardiac studies using ECT , 1994, IEEE Trans. Medical Imaging.

[10]  Alfred O. Hero,et al.  Model-based estimation with boundary side information or boundary regularization [cardiac emission CT] , 1994, IEEE Trans. Medical Imaging.

[11]  J. Borwein,et al.  Direct reconstruction of functional parameters for dynamic SPECT , 1995 .

[12]  G. Zeng,et al.  Using linear time-invariant system theory to estimate kinetic parameters directly from projection measurements , 1995 .

[13]  G. Gullberg,et al.  Experimental verification of technetium 99m-labeled teboroxime kinetic parameters the the myocardium with dynamic single-photon emission computed tomography: Reproducibility, correlation to flow, and susceptibility to extravascular contamination , 1996, Journal of nuclear cardiology : official publication of the American Society of Nuclear Cardiology.

[14]  B. W. Reutter,et al.  Kinetic parameter estimation from SPECT projection measurements , 1997 .