### Abstract

In dynamic cardiac SPECT estimates of kinetic parameters of a one-compartment perfusion model are usually obtained in a two step process: 1) first a MAP iterative algorithm, which properly models the Poisson statistics and the physics of the data acquisition, reconstructs a sequence of dynamic reconstructions, 2) then kinetic parameters are estimated from time activity curves generated from the dynamic reconstructions. This paper provides a method for calculating the covariance matrix of the kinetic parameters, which are determined using weighted least squares fitting that incorporates the estimated variance and covariance of the dynamic reconstructions. Sequential tomographic projections are reconstructed into a sequence of transaxial reconstructions for each transaxial slice using for each reconstruction in the time sequence the fixed-point solution to the MAP reconstruction. Time-activity curves for a sum of activity in a blood region inside the left ventricle and a sum in a cardiac tissue region, for the variance of the two estimates of the sum, and for the covariance between the two ROI estimates are generated at convergence. A one-compartment model is fit to the tissue activity curves assuming a noisy blood input function to give weighted least squares estimates of blood volume fraction, wash-in and wash-out rate constants specifying the kinetics for the left ventricular myocardium. Numerical methods are used to calculate the second derivative of the chi-square criterion to obtain estimates of the covariance matrix for the weighted least square parameter estimates. Even though the method requires one matrix inverse for each time interval of tomographic acquisition, efficient estimates of the tissue kinetic parameters in a dynamic cardiac SPECT study can be obtained with present day desk-top computers.

Original language | English (US) |
---|---|

Title of host publication | IEEE Nuclear Science Symposium Conference Record |

Editors | S.D. Metzler |

Pages | 3019-3023 |

Number of pages | 5 |

Volume | 5 |

State | Published - 2003 |

Externally published | Yes |

Event | 2003 IEEE Nuclear Science Symposium Conference Record - Nuclear Science Symposium, Medical Imaging Conference - Portland, OR, United States Duration: Oct 19 2003 → Oct 25 2003 |

### Other

Other | 2003 IEEE Nuclear Science Symposium Conference Record - Nuclear Science Symposium, Medical Imaging Conference |
---|---|

Country | United States |

City | Portland, OR |

Period | 10/19/03 → 10/25/03 |

### Fingerprint

### ASJC Scopus subject areas

- Computer Vision and Pattern Recognition
- Industrial and Manufacturing Engineering

### Cite this

*IEEE Nuclear Science Symposium Conference Record*(Vol. 5, pp. 3019-3023). [M13-1]

**Estimation of the parameter covariance matrix for a one-compartment cardiac perfusion model estimated from a dynamic sequence reconstructed using MAP iterative reconstruction algorithms.** / Gullberg, Grant T.; Huesman, Ronald H.; Roy, Dilip N Ghosh; Qi, Jinyi; Reutter, Bryan W.

Research output: Chapter in Book/Report/Conference proceeding › Conference contribution

*IEEE Nuclear Science Symposium Conference Record.*vol. 5, M13-1, pp. 3019-3023, 2003 IEEE Nuclear Science Symposium Conference Record - Nuclear Science Symposium, Medical Imaging Conference, Portland, OR, United States, 10/19/03.

}

TY - GEN

T1 - Estimation of the parameter covariance matrix for a one-compartment cardiac perfusion model estimated from a dynamic sequence reconstructed using MAP iterative reconstruction algorithms

AU - Gullberg, Grant T.

AU - Huesman, Ronald H.

AU - Roy, Dilip N Ghosh

AU - Qi, Jinyi

AU - Reutter, Bryan W.

PY - 2003

Y1 - 2003

N2 - In dynamic cardiac SPECT estimates of kinetic parameters of a one-compartment perfusion model are usually obtained in a two step process: 1) first a MAP iterative algorithm, which properly models the Poisson statistics and the physics of the data acquisition, reconstructs a sequence of dynamic reconstructions, 2) then kinetic parameters are estimated from time activity curves generated from the dynamic reconstructions. This paper provides a method for calculating the covariance matrix of the kinetic parameters, which are determined using weighted least squares fitting that incorporates the estimated variance and covariance of the dynamic reconstructions. Sequential tomographic projections are reconstructed into a sequence of transaxial reconstructions for each transaxial slice using for each reconstruction in the time sequence the fixed-point solution to the MAP reconstruction. Time-activity curves for a sum of activity in a blood region inside the left ventricle and a sum in a cardiac tissue region, for the variance of the two estimates of the sum, and for the covariance between the two ROI estimates are generated at convergence. A one-compartment model is fit to the tissue activity curves assuming a noisy blood input function to give weighted least squares estimates of blood volume fraction, wash-in and wash-out rate constants specifying the kinetics for the left ventricular myocardium. Numerical methods are used to calculate the second derivative of the chi-square criterion to obtain estimates of the covariance matrix for the weighted least square parameter estimates. Even though the method requires one matrix inverse for each time interval of tomographic acquisition, efficient estimates of the tissue kinetic parameters in a dynamic cardiac SPECT study can be obtained with present day desk-top computers.

AB - In dynamic cardiac SPECT estimates of kinetic parameters of a one-compartment perfusion model are usually obtained in a two step process: 1) first a MAP iterative algorithm, which properly models the Poisson statistics and the physics of the data acquisition, reconstructs a sequence of dynamic reconstructions, 2) then kinetic parameters are estimated from time activity curves generated from the dynamic reconstructions. This paper provides a method for calculating the covariance matrix of the kinetic parameters, which are determined using weighted least squares fitting that incorporates the estimated variance and covariance of the dynamic reconstructions. Sequential tomographic projections are reconstructed into a sequence of transaxial reconstructions for each transaxial slice using for each reconstruction in the time sequence the fixed-point solution to the MAP reconstruction. Time-activity curves for a sum of activity in a blood region inside the left ventricle and a sum in a cardiac tissue region, for the variance of the two estimates of the sum, and for the covariance between the two ROI estimates are generated at convergence. A one-compartment model is fit to the tissue activity curves assuming a noisy blood input function to give weighted least squares estimates of blood volume fraction, wash-in and wash-out rate constants specifying the kinetics for the left ventricular myocardium. Numerical methods are used to calculate the second derivative of the chi-square criterion to obtain estimates of the covariance matrix for the weighted least square parameter estimates. Even though the method requires one matrix inverse for each time interval of tomographic acquisition, efficient estimates of the tissue kinetic parameters in a dynamic cardiac SPECT study can be obtained with present day desk-top computers.

UR - http://www.scopus.com/inward/record.url?scp=11944253937&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=11944253937&partnerID=8YFLogxK

M3 - Conference contribution

VL - 5

SP - 3019

EP - 3023

BT - IEEE Nuclear Science Symposium Conference Record

A2 - Metzler, S.D.

ER -