Integral curves from noisy diffusion MRI data with closed-form uncertainty estimates

Owen Carmichael, Lyudmila Sakhanenko

Research output: Contribution to journalArticle

Abstract

We present a method for estimating the trajectories of axon fibers through diffusion tensor MRI (DTI) data that provides theoretically rigorous estimates of trajectory uncertainty. We develop a three-step estimation procedure based on a kernel estimator for a tensor field based on the raw DTI measurements, followed by a plug-in estimator for the leading eigenvectors of the tensors, and a plug-in estimator for integral curves through the resulting vector field. The integral curve estimator is asymptotically normal; the covariance of the limiting Gaussian process allows us to construct confidence ellipsoids for fixed points along the curve. Complete trajectories of fibers are assembled by stopping integral curve tracing at locations with multiple viable leading eigenvector directions and tracing a new curve along each direction. Unlike probabilistic tractography approaches to this problem, we provide a rigorous, theoretically sound model of measurement uncertainty as it propagates from the raw MRI data, to the tensor field, to the vector field, to the integral curves. In addition, trajectory uncertainty is estimated in closed form while probabilistic tractography relies on sampling the space of tensors, vectors, or curves. We show that our estimator provides more realistic trajectory uncertainty estimates than a more simplified prior approach for closed-form trajectory uncertainty estimation due to Koltchinskii et al. (Ann Stat 35:1576–1607, 2007) and a popular probabilistic tractography method due to Behrens et al. (Magn Reson Med 50:1077–1088, 2003) using theory, simulation, and real DTI scans.

Original languageEnglish (US)
JournalStatistical Inference for Stochastic Processes
DOIs
StateAccepted/In press - Aug 27 2015
Externally publishedYes

Fingerprint

Closed-form
Tensor
Uncertainty
Curve
Trajectory
Estimate
Plug-in Estimator
Tracing
Eigenvector
Vector Field
Fiber
Uncertainty Estimation
Estimator
Measurement Uncertainty
Kernel Estimator
Probabilistic Methods
Probabilistic Approach
Ellipsoid
Gaussian Process
Confidence

Keywords

  • Asymptotic normality
  • Diffusion tensor imaging
  • Integral curve
  • Kernel smoothing

ASJC Scopus subject areas

  • Statistics and Probability

Cite this

Integral curves from noisy diffusion MRI data with closed-form uncertainty estimates. / Carmichael, Owen; Sakhanenko, Lyudmila.

In: Statistical Inference for Stochastic Processes, 27.08.2015.

Research output: Contribution to journalArticle

@article{f178d54ad7694a7787ca473708c453a0,
title = "Integral curves from noisy diffusion MRI data with closed-form uncertainty estimates",
abstract = "We present a method for estimating the trajectories of axon fibers through diffusion tensor MRI (DTI) data that provides theoretically rigorous estimates of trajectory uncertainty. We develop a three-step estimation procedure based on a kernel estimator for a tensor field based on the raw DTI measurements, followed by a plug-in estimator for the leading eigenvectors of the tensors, and a plug-in estimator for integral curves through the resulting vector field. The integral curve estimator is asymptotically normal; the covariance of the limiting Gaussian process allows us to construct confidence ellipsoids for fixed points along the curve. Complete trajectories of fibers are assembled by stopping integral curve tracing at locations with multiple viable leading eigenvector directions and tracing a new curve along each direction. Unlike probabilistic tractography approaches to this problem, we provide a rigorous, theoretically sound model of measurement uncertainty as it propagates from the raw MRI data, to the tensor field, to the vector field, to the integral curves. In addition, trajectory uncertainty is estimated in closed form while probabilistic tractography relies on sampling the space of tensors, vectors, or curves. We show that our estimator provides more realistic trajectory uncertainty estimates than a more simplified prior approach for closed-form trajectory uncertainty estimation due to Koltchinskii et al. (Ann Stat 35:1576–1607, 2007) and a popular probabilistic tractography method due to Behrens et al. (Magn Reson Med 50:1077–1088, 2003) using theory, simulation, and real DTI scans.",
keywords = "Asymptotic normality, Diffusion tensor imaging, Integral curve, Kernel smoothing",
author = "Owen Carmichael and Lyudmila Sakhanenko",
year = "2015",
month = "8",
day = "27",
doi = "10.1007/s11203-015-9126-9",
language = "English (US)",
journal = "Statistical Inference for Stochastic Processes",
issn = "1387-0874",
publisher = "Springer Netherlands",

}

TY - JOUR

T1 - Integral curves from noisy diffusion MRI data with closed-form uncertainty estimates

AU - Carmichael, Owen

AU - Sakhanenko, Lyudmila

PY - 2015/8/27

Y1 - 2015/8/27

N2 - We present a method for estimating the trajectories of axon fibers through diffusion tensor MRI (DTI) data that provides theoretically rigorous estimates of trajectory uncertainty. We develop a three-step estimation procedure based on a kernel estimator for a tensor field based on the raw DTI measurements, followed by a plug-in estimator for the leading eigenvectors of the tensors, and a plug-in estimator for integral curves through the resulting vector field. The integral curve estimator is asymptotically normal; the covariance of the limiting Gaussian process allows us to construct confidence ellipsoids for fixed points along the curve. Complete trajectories of fibers are assembled by stopping integral curve tracing at locations with multiple viable leading eigenvector directions and tracing a new curve along each direction. Unlike probabilistic tractography approaches to this problem, we provide a rigorous, theoretically sound model of measurement uncertainty as it propagates from the raw MRI data, to the tensor field, to the vector field, to the integral curves. In addition, trajectory uncertainty is estimated in closed form while probabilistic tractography relies on sampling the space of tensors, vectors, or curves. We show that our estimator provides more realistic trajectory uncertainty estimates than a more simplified prior approach for closed-form trajectory uncertainty estimation due to Koltchinskii et al. (Ann Stat 35:1576–1607, 2007) and a popular probabilistic tractography method due to Behrens et al. (Magn Reson Med 50:1077–1088, 2003) using theory, simulation, and real DTI scans.

AB - We present a method for estimating the trajectories of axon fibers through diffusion tensor MRI (DTI) data that provides theoretically rigorous estimates of trajectory uncertainty. We develop a three-step estimation procedure based on a kernel estimator for a tensor field based on the raw DTI measurements, followed by a plug-in estimator for the leading eigenvectors of the tensors, and a plug-in estimator for integral curves through the resulting vector field. The integral curve estimator is asymptotically normal; the covariance of the limiting Gaussian process allows us to construct confidence ellipsoids for fixed points along the curve. Complete trajectories of fibers are assembled by stopping integral curve tracing at locations with multiple viable leading eigenvector directions and tracing a new curve along each direction. Unlike probabilistic tractography approaches to this problem, we provide a rigorous, theoretically sound model of measurement uncertainty as it propagates from the raw MRI data, to the tensor field, to the vector field, to the integral curves. In addition, trajectory uncertainty is estimated in closed form while probabilistic tractography relies on sampling the space of tensors, vectors, or curves. We show that our estimator provides more realistic trajectory uncertainty estimates than a more simplified prior approach for closed-form trajectory uncertainty estimation due to Koltchinskii et al. (Ann Stat 35:1576–1607, 2007) and a popular probabilistic tractography method due to Behrens et al. (Magn Reson Med 50:1077–1088, 2003) using theory, simulation, and real DTI scans.

KW - Asymptotic normality

KW - Diffusion tensor imaging

KW - Integral curve

KW - Kernel smoothing

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

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

U2 - 10.1007/s11203-015-9126-9

DO - 10.1007/s11203-015-9126-9

M3 - Article

AN - SCOPUS:84940213702

JO - Statistical Inference for Stochastic Processes

JF - Statistical Inference for Stochastic Processes

SN - 1387-0874

ER -