### Abstract

We develop statistical methodology for a popular brain imaging technique HARDI based on the high order tensor model by Özarslan and Mareci [10]. We investigate how uncertainty in the imaging procedure propagates through all levels of the model: signals, tensor fields, vector fields, and fibers. We construct asymptotically normal estimators of the integral curves or fibers which allow us to trace the fibers together with confidence ellipsoids. The procedure is computationally intense as it blends linear algebra concepts from high order tensors with asymptotical statistical analysis. The theoretical results are illustrated on simulated and real datasets. This work generalizes the statistical methodology proposed for low angular resolution diffusion tensor imaging by Carmichael and Sakhanenko [3], to several fibers per voxel. It is also a pioneering statistical work on tractography from HARDI data. It avoids all the typical limitations of the deterministic tractography methods and it delivers the same information as probabilistic tractography methods. Our method is computationally cheap and it provides well-founded mathematical and statistical framework where diverse functionals on fibers, directions and tensors can be studied in a systematic and rigorous way.

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

Pages (from-to) | 377-403 |

Number of pages | 27 |

Journal | Linear Algebra and Its Applications |

Volume | 473 |

DOIs | |

State | Published - May 15 2015 |

### Fingerprint

### Keywords

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

### ASJC Scopus subject areas

- Algebra and Number Theory
- Discrete Mathematics and Combinatorics
- Geometry and Topology
- Numerical Analysis

### Cite this

*Linear Algebra and Its Applications*,

*473*, 377-403. https://doi.org/10.1016/j.laa.2014.12.007

**Estimation of integral curves from high angular resolution diffusion imaging (HARDI) data.** / Carmichael, Owen; Sakhanenko, Lyudmila.

Research output: Contribution to journal › Article

*Linear Algebra and Its Applications*, vol. 473, pp. 377-403. https://doi.org/10.1016/j.laa.2014.12.007

}

TY - JOUR

T1 - Estimation of integral curves from high angular resolution diffusion imaging (HARDI) data

AU - Carmichael, Owen

AU - Sakhanenko, Lyudmila

PY - 2015/5/15

Y1 - 2015/5/15

N2 - We develop statistical methodology for a popular brain imaging technique HARDI based on the high order tensor model by Özarslan and Mareci [10]. We investigate how uncertainty in the imaging procedure propagates through all levels of the model: signals, tensor fields, vector fields, and fibers. We construct asymptotically normal estimators of the integral curves or fibers which allow us to trace the fibers together with confidence ellipsoids. The procedure is computationally intense as it blends linear algebra concepts from high order tensors with asymptotical statistical analysis. The theoretical results are illustrated on simulated and real datasets. This work generalizes the statistical methodology proposed for low angular resolution diffusion tensor imaging by Carmichael and Sakhanenko [3], to several fibers per voxel. It is also a pioneering statistical work on tractography from HARDI data. It avoids all the typical limitations of the deterministic tractography methods and it delivers the same information as probabilistic tractography methods. Our method is computationally cheap and it provides well-founded mathematical and statistical framework where diverse functionals on fibers, directions and tensors can be studied in a systematic and rigorous way.

AB - We develop statistical methodology for a popular brain imaging technique HARDI based on the high order tensor model by Özarslan and Mareci [10]. We investigate how uncertainty in the imaging procedure propagates through all levels of the model: signals, tensor fields, vector fields, and fibers. We construct asymptotically normal estimators of the integral curves or fibers which allow us to trace the fibers together with confidence ellipsoids. The procedure is computationally intense as it blends linear algebra concepts from high order tensors with asymptotical statistical analysis. The theoretical results are illustrated on simulated and real datasets. This work generalizes the statistical methodology proposed for low angular resolution diffusion tensor imaging by Carmichael and Sakhanenko [3], to several fibers per voxel. It is also a pioneering statistical work on tractography from HARDI data. It avoids all the typical limitations of the deterministic tractography methods and it delivers the same information as probabilistic tractography methods. Our method is computationally cheap and it provides well-founded mathematical and statistical framework where diverse functionals on fibers, directions and tensors can be studied in a systematic and rigorous way.

KW - Asymptotic normality

KW - Diffusion tensor imaging

KW - Integral curve

KW - Kernel smoothing

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

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

U2 - 10.1016/j.laa.2014.12.007

DO - 10.1016/j.laa.2014.12.007

M3 - Article

VL - 473

SP - 377

EP - 403

JO - Linear Algebra and Its Applications

JF - Linear Algebra and Its Applications

SN - 0024-3795

ER -