skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Fourier domain closed-form formulas for estimation of kinetic parameters in reversible multi-compartment models

Journal Article · · Biomedical Engineering Online
 [1];  [1];  [2]
  1. Univ. of Utah, Salt Lake City, UT (United States)
  2. Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States)

Background: Compared with static imaging, dynamic emission computed tomographic imaging with compartment modeling can quantify in vivo physiologic processes, providing useful information about molecular disease processes. Dynamic imaging involves estimation of kinetic rate parameters. For multi-compartment models, kinetic parameter estimation can be computationally demanding and problematic with local minima. Methods: This paper offers a new perspective to the compartment model fitting problem where Fourier linear system theory is applied to derive closed-form formulas for estimating kinetic parameters for the two-compartment model. The proposed Fourier domain estimation method provides a unique solution, and offers very different noise response as compared to traditional non-linear chi-squared minimization techniques. Results: The unique feature of the proposed Fourier domain method is that only low frequency components are used for kinetic parameter estimation, where the DC (i.e., the zero frequency) component in the data is treated as the most important information, and high frequency components that tend to be corrupted by statistical noise are discarded. Computer simulations show that the proposed method is robust without having to specify the initial condition. The resultant solution can be fine tuned using the traditional iterative method. Conclusions: The proposed Fourier-domain estimation method has closed-form formulas. The proposed Fourier-domain curve-fitting method does not require an initial condition, it minimizes a quadratic objective function and a closed-form solution can be obtained. The noise is easier to control, simply by discarding the high frequency components, and emphasizing the DC component.

Research Organization:
Lawrence Berkeley National Laboratory (LBNL), Berkeley, CA (United States)
Sponsoring Organization:
Ben B. and Iris M. Margolis Foundation; National Institutes of Health (NIH); USDOE Office of Science (SC), Biological and Environmental Research (BER)
Grant/Contract Number:
AC02-05CH11231; R01 HL108350; R01 CA135556; R01 HL50663; R01 EB007219
OSTI ID:
1626567
Journal Information:
Biomedical Engineering Online, Vol. 11, Issue 1; ISSN 1475-925X
Publisher:
BioMed CentralCopyright Statement
Country of Publication:
United States
Language:
English

References (23)

Weighted Integration Method for Local Cerebral Blood Flow Measurements with Positron Emission Tomography journal April 1986
Models for computer simulation studies of input functions for tracer kinetic modeling with positron emission tomography journal March 1993
On the Construction of Functional Maps in Positron Emission Tomography journal December 1984
Use of ridge regression for improved estimation of kinetic constants from PET data journal January 1999
Optimization algorithms and weighting factors for analysis of dynamic PET studies journal August 2006
An evaluation of the algorithms for determining local cerebral metabolic rates of glucose using positron emission tomography dynamic data journal January 1995
Kinetic modelling using basis functions derived from two-tissue compartmental models with a plasma input function: General principle and application to [18F]fluorodeoxyglucose positron emission tomography journal May 2010
Estimation of input function and kinetic parameters using simulated annealing: application in a flow model journal June 2002
Rapid algorithms for the construction of cerebral blood flow and oxygen utilization images with oxygen-15 and dynamic positron emission tomography journal February 1999
Linear ridge regression with spatial constraint for generation of parametric images in dynamic positron emission tomography studies journal February 2001
Joint estimation of dynamic PET images and temporal basis functions using fully 4D ML-EM journal October 2006
Linear regression with spatial constraint to generate parametric images of ligand-receptor dynamic PET studies with a simplified reference tissue model journal April 2003
A New Approach of Weighted Integration Technique Based on Accumulated Images Using Dynamic PET and H 15 2 O journal May 1991
Generalized linear least squares method for fast generation of myocardial blood flow parametric images with N-13 ammonia PET journal April 1998
GLLS for optimally sampled continuous dynamic system modeling: theory and algorithm journal April 1999
PET kinetic analysis—compartmental model journal November 2006
Improved Parametric Image Generation Using Spatial-Temporal Analysis of Dynamic PET Studies journal March 2002
Closed-form kinetic parameter estimation solution to the truncated data problem journal November 2010
Evaluation of Basis Function and Linear Least Squares Methods for Generating Parametric Blood Flow Images Using 15O-Water and Positron Emission Tomography journal July 2005
Performance evaluation of kinetic parameter estimation methods in dynamic FDG-PET studies journal January 2011
Dynamic single photon emission computed tomography—basic principles and cardiac applications journal September 2010
Positron Emission Tomography Compartmental Models journal June 2001
Fourier domain closed-form formulas for estimation of kinetic parameters in multi-compartment models conference October 2011

Cited By (2)

Modeling of Microvascular Permeability Changes after Electroporation journal March 2015
Voxel-level mapping of tracer kinetics in PET studies: A statistical approach emphasizing tissue life tables text January 2014