A model-constrained Monte Carlo method for blind arterial input function estimation in dynamic contrast-enhanced MRI: I. simulations

Matthias Schabel, Jacob U. Fluckiger, Edward V R Dibella

Research output: Contribution to journalArticle

29 Citations (Scopus)

Abstract

Widespread adoption of quantitative pharmacokinetic modeling methods in conjunction with dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) has led to increased recognition of the importance of obtaining accurate patient-specific arterial input function (AIF) measurements. Ideally, DCE-MRI studies use an AIF directly measured in an artery local to the tissue of interest, along with measured tissue concentration curves, to quantitatively determine pharmacokinetic parameters. However, the numerous technical and practical difficulties associated with AIF measurement have made the use of population-averaged AIF data a popular, if sub-optimal, alternative to AIF measurement. In this work, we present and characterize a new algorithm for determining the AIF solely from the measured tissue concentration curves. This Monte Carlo blind estimation (MCBE) algorithm estimates the AIF from the subsets of D concentration-time curves drawn from a larger pool of M candidate curves via nonlinear optimization, doing so for multiple (Q) subsets and statistically averaging these repeated estimates. The MCBE algorithm can be viewed as a generalization of previously published methods that employ clustering of concentration-time curves and only estimate the AIF once. Extensive computer simulations were performed over physiologically and experimentally realistic ranges of imaging and tissue parameters, and the impact of choosing different values of D and Q was investigated. We found the algorithm to be robust, computationally efficient and capable of accurately estimating the AIF even for relatively high noise levels, long sampling intervals and low diversity of tissue curves. With the incorporation of bootstrapping initialization, we further demonstrated the ability to blindly estimate AIFs that deviate substantially in shape from the population-averaged initial guess. Pharmacokinetic parameter estimates for Ktrans, kep, vp and ve all showed relative biases and uncertainties of less than 10% for measurements having a temporal sampling rate of 4 s and a concentration measurement noise level of σ = 0.04 mM. A companion paper discusses the application of the MCBE algorithm to DCE-MRI data acquired in eight patients with malignant brain tumors.

Original languageEnglish (US)
Pages (from-to)4783-4806
Number of pages24
JournalPhysics in Medicine and Biology
Volume55
Issue number16
DOIs
StatePublished - Aug 21 2010
Externally publishedYes

Fingerprint

Monte Carlo Method
Pharmacokinetics
Magnetic Resonance Imaging
Noise
Brain Neoplasms
Computer Simulation
Population
Uncertainty
Cluster Analysis
Arteries

ASJC Scopus subject areas

  • Radiology Nuclear Medicine and imaging
  • Radiological and Ultrasound Technology

Cite this

A model-constrained Monte Carlo method for blind arterial input function estimation in dynamic contrast-enhanced MRI : I. simulations. / Schabel, Matthias; Fluckiger, Jacob U.; Dibella, Edward V R.

In: Physics in Medicine and Biology, Vol. 55, No. 16, 21.08.2010, p. 4783-4806.

Research output: Contribution to journalArticle

@article{da1a655944a54dbb8f1f65b3e394b68e,
title = "A model-constrained Monte Carlo method for blind arterial input function estimation in dynamic contrast-enhanced MRI: I. simulations",
abstract = "Widespread adoption of quantitative pharmacokinetic modeling methods in conjunction with dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) has led to increased recognition of the importance of obtaining accurate patient-specific arterial input function (AIF) measurements. Ideally, DCE-MRI studies use an AIF directly measured in an artery local to the tissue of interest, along with measured tissue concentration curves, to quantitatively determine pharmacokinetic parameters. However, the numerous technical and practical difficulties associated with AIF measurement have made the use of population-averaged AIF data a popular, if sub-optimal, alternative to AIF measurement. In this work, we present and characterize a new algorithm for determining the AIF solely from the measured tissue concentration curves. This Monte Carlo blind estimation (MCBE) algorithm estimates the AIF from the subsets of D concentration-time curves drawn from a larger pool of M candidate curves via nonlinear optimization, doing so for multiple (Q) subsets and statistically averaging these repeated estimates. The MCBE algorithm can be viewed as a generalization of previously published methods that employ clustering of concentration-time curves and only estimate the AIF once. Extensive computer simulations were performed over physiologically and experimentally realistic ranges of imaging and tissue parameters, and the impact of choosing different values of D and Q was investigated. We found the algorithm to be robust, computationally efficient and capable of accurately estimating the AIF even for relatively high noise levels, long sampling intervals and low diversity of tissue curves. With the incorporation of bootstrapping initialization, we further demonstrated the ability to blindly estimate AIFs that deviate substantially in shape from the population-averaged initial guess. Pharmacokinetic parameter estimates for Ktrans, kep, vp and ve all showed relative biases and uncertainties of less than 10{\%} for measurements having a temporal sampling rate of 4 s and a concentration measurement noise level of σ = 0.04 mM. A companion paper discusses the application of the MCBE algorithm to DCE-MRI data acquired in eight patients with malignant brain tumors.",
author = "Matthias Schabel and Fluckiger, {Jacob U.} and Dibella, {Edward V R}",
year = "2010",
month = "8",
day = "21",
doi = "10.1088/0031-9155/55/16/011",
language = "English (US)",
volume = "55",
pages = "4783--4806",
journal = "Physics in Medicine and Biology",
issn = "0031-9155",
publisher = "IOP Publishing Ltd.",
number = "16",

}

TY - JOUR

T1 - A model-constrained Monte Carlo method for blind arterial input function estimation in dynamic contrast-enhanced MRI

T2 - I. simulations

AU - Schabel, Matthias

AU - Fluckiger, Jacob U.

AU - Dibella, Edward V R

PY - 2010/8/21

Y1 - 2010/8/21

N2 - Widespread adoption of quantitative pharmacokinetic modeling methods in conjunction with dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) has led to increased recognition of the importance of obtaining accurate patient-specific arterial input function (AIF) measurements. Ideally, DCE-MRI studies use an AIF directly measured in an artery local to the tissue of interest, along with measured tissue concentration curves, to quantitatively determine pharmacokinetic parameters. However, the numerous technical and practical difficulties associated with AIF measurement have made the use of population-averaged AIF data a popular, if sub-optimal, alternative to AIF measurement. In this work, we present and characterize a new algorithm for determining the AIF solely from the measured tissue concentration curves. This Monte Carlo blind estimation (MCBE) algorithm estimates the AIF from the subsets of D concentration-time curves drawn from a larger pool of M candidate curves via nonlinear optimization, doing so for multiple (Q) subsets and statistically averaging these repeated estimates. The MCBE algorithm can be viewed as a generalization of previously published methods that employ clustering of concentration-time curves and only estimate the AIF once. Extensive computer simulations were performed over physiologically and experimentally realistic ranges of imaging and tissue parameters, and the impact of choosing different values of D and Q was investigated. We found the algorithm to be robust, computationally efficient and capable of accurately estimating the AIF even for relatively high noise levels, long sampling intervals and low diversity of tissue curves. With the incorporation of bootstrapping initialization, we further demonstrated the ability to blindly estimate AIFs that deviate substantially in shape from the population-averaged initial guess. Pharmacokinetic parameter estimates for Ktrans, kep, vp and ve all showed relative biases and uncertainties of less than 10% for measurements having a temporal sampling rate of 4 s and a concentration measurement noise level of σ = 0.04 mM. A companion paper discusses the application of the MCBE algorithm to DCE-MRI data acquired in eight patients with malignant brain tumors.

AB - Widespread adoption of quantitative pharmacokinetic modeling methods in conjunction with dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) has led to increased recognition of the importance of obtaining accurate patient-specific arterial input function (AIF) measurements. Ideally, DCE-MRI studies use an AIF directly measured in an artery local to the tissue of interest, along with measured tissue concentration curves, to quantitatively determine pharmacokinetic parameters. However, the numerous technical and practical difficulties associated with AIF measurement have made the use of population-averaged AIF data a popular, if sub-optimal, alternative to AIF measurement. In this work, we present and characterize a new algorithm for determining the AIF solely from the measured tissue concentration curves. This Monte Carlo blind estimation (MCBE) algorithm estimates the AIF from the subsets of D concentration-time curves drawn from a larger pool of M candidate curves via nonlinear optimization, doing so for multiple (Q) subsets and statistically averaging these repeated estimates. The MCBE algorithm can be viewed as a generalization of previously published methods that employ clustering of concentration-time curves and only estimate the AIF once. Extensive computer simulations were performed over physiologically and experimentally realistic ranges of imaging and tissue parameters, and the impact of choosing different values of D and Q was investigated. We found the algorithm to be robust, computationally efficient and capable of accurately estimating the AIF even for relatively high noise levels, long sampling intervals and low diversity of tissue curves. With the incorporation of bootstrapping initialization, we further demonstrated the ability to blindly estimate AIFs that deviate substantially in shape from the population-averaged initial guess. Pharmacokinetic parameter estimates for Ktrans, kep, vp and ve all showed relative biases and uncertainties of less than 10% for measurements having a temporal sampling rate of 4 s and a concentration measurement noise level of σ = 0.04 mM. A companion paper discusses the application of the MCBE algorithm to DCE-MRI data acquired in eight patients with malignant brain tumors.

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

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

U2 - 10.1088/0031-9155/55/16/011

DO - 10.1088/0031-9155/55/16/011

M3 - Article

C2 - 20679691

AN - SCOPUS:78149354861

VL - 55

SP - 4783

EP - 4806

JO - Physics in Medicine and Biology

JF - Physics in Medicine and Biology

SN - 0031-9155

IS - 16

ER -