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

Grant T. Gullberg, Ronald H. Huesman, Dilip N Ghosh Roy, Jinyi Qi, Bryan W. Reutter

Research output: Chapter in Book/Report/Conference proceedingConference contribution

3 Citations (Scopus)

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 languageEnglish (US)
Title of host publicationIEEE Nuclear Science Symposium Conference Record
EditorsS.D. Metzler
Pages3019-3023
Number of pages5
Volume5
StatePublished - 2003
Externally publishedYes
Event2003 IEEE Nuclear Science Symposium Conference Record - Nuclear Science Symposium, Medical Imaging Conference - Portland, OR, United States
Duration: Oct 19 2003Oct 25 2003

Other

Other2003 IEEE Nuclear Science Symposium Conference Record - Nuclear Science Symposium, Medical Imaging Conference
CountryUnited States
CityPortland, OR
Period10/19/0310/25/03

Fingerprint

Covariance matrix
Kinetic parameters
Blood
Tissue
Rate constants
Volume fraction
Data acquisition
Numerical methods
Physics
Statistics
Derivatives
Kinetics

ASJC Scopus subject areas

  • Computer Vision and Pattern Recognition
  • Industrial and Manufacturing Engineering

Cite this

Gullberg, G. T., Huesman, R. H., Roy, D. N. G., Qi, J., & Reutter, B. W. (2003). Estimation of the parameter covariance matrix for a one-compartment cardiac perfusion model estimated from a dynamic sequence reconstructed using MAP iterative reconstruction algorithms. In S. D. Metzler (Ed.), 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.

IEEE Nuclear Science Symposium Conference Record. ed. / S.D. Metzler. Vol. 5 2003. p. 3019-3023 M13-1.

Research output: Chapter in Book/Report/Conference proceedingConference contribution

Gullberg, GT, Huesman, RH, Roy, DNG, Qi, J & Reutter, BW 2003, Estimation of the parameter covariance matrix for a one-compartment cardiac perfusion model estimated from a dynamic sequence reconstructed using MAP iterative reconstruction algorithms. in SD Metzler (ed.), 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.
Gullberg GT, Huesman RH, Roy DNG, Qi J, Reutter BW. Estimation of the parameter covariance matrix for a one-compartment cardiac perfusion model estimated from a dynamic sequence reconstructed using MAP iterative reconstruction algorithms. In Metzler SD, editor, IEEE Nuclear Science Symposium Conference Record. Vol. 5. 2003. p. 3019-3023. M13-1
Gullberg, Grant T. ; Huesman, Ronald H. ; Roy, Dilip N Ghosh ; Qi, Jinyi ; Reutter, Bryan W. / Estimation of the parameter covariance matrix for a one-compartment cardiac perfusion model estimated from a dynamic sequence reconstructed using MAP iterative reconstruction algorithms. IEEE Nuclear Science Symposium Conference Record. editor / S.D. Metzler. Vol. 5 2003. pp. 3019-3023
@inproceedings{92a6635d3d1a4278b0bd1b4de5ec47f2,
title = "Estimation of the parameter covariance matrix for a one-compartment cardiac perfusion model estimated from a dynamic sequence reconstructed using MAP iterative reconstruction algorithms",
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.",
author = "Gullberg, {Grant T.} and Huesman, {Ronald H.} and Roy, {Dilip N Ghosh} and Jinyi Qi and Reutter, {Bryan W.}",
year = "2003",
language = "English (US)",
volume = "5",
pages = "3019--3023",
editor = "S.D. Metzler",
booktitle = "IEEE Nuclear Science Symposium Conference Record",

}

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

AN - SCOPUS:11944253937

VL - 5

SP - 3019

EP - 3023

BT - IEEE Nuclear Science Symposium Conference Record

A2 - Metzler, S.D.

ER -