A normal distribution for tensor-valued random variables: Applications to diffusion tensor MRI

被引:87
作者
Basser, PJ
Pajevic, S
机构
[1] NICHD, STBB, LIMB, NIH, Bethesda, MD 20892 USA
[2] NIH, MSCL, CIT, Bethesda, MD 20892 USA
关键词
covariance; distribution; experimental design; fourth-order; Gaussian; normal; precision; probability; random variable; second-order; strain-energy; tensor;
D O I
10.1109/TMI.2003.815059
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
Diffusion tensor magnetic resonance imaging (DT-MRI) provides a statistical estimate of a symmetric, second-order diffusion tensor of water, D, in each voxel within an imaging volume. We propose a new normal distribution, p(D) proportional to exp(-1/2 D : A : D), which describes the variability of D in an ideal DT-MRI experiment. The scalar invariant, D : A : D, is the contraction of a positive definite symmetric, fourth-order precision tensor, A, and D. A correspondence is established between D: A: D and the elastic strain energy density function in continuum mechanics-specifically between D and the second-order infinitesimal strain tensor, and between A and the fourth-order tensor of elastic coefficients. We show that A can be further classified according to different classical elastic symmetries (i.e., isotropy, transverse isotropy, orthotropy, planar symmetry, and anisotropy). When A is an isotropic fourth-order tensor, we derive an explicit analytic expression for p(D), and for the distribution of the three eigenvalues of D, p(gamma(1), gamma(2), gamma(3)), which are confirmed by Monte Carlo simulations. We show how A can be estimated from either real or synthetic DT-MRI data for any given experimental design. Here we propose a new criterion for an optimal experimental design: that A be an isotropic fourth-order tensor. This condition ensures that the statistical properties of D (and quantities derived from it) are rotationally invariant. We also investigate the degree of isotropy of several DT-MRI experimental designs. Finally, we show that the univariate and multivariate distributions are special cases of the more general tensor-variate normal distribution, and suggest how to generalize p(D) to treat normal random tensor variables that are of third- (or higher) order. We expect that this new distribution, p(D), should be useful in feature extraction; in developing a hypothesis testing framework for segmenting and classifying noisy, discrete tensor data; and in designing experiments to measure tensor quantities.
引用
收藏
页码:785 / 794
页数:10
相关论文
共 30 条
  • [1] Spatial transformations of diffusion tensor magnetic resonance images
    Alexander, DC
    Pierpaoli, C
    Basser, PJ
    Gee, JC
    [J]. IEEE TRANSACTIONS ON MEDICAL IMAGING, 2001, 20 (11) : 1131 - 1139
  • [2] Anderson T., 1984, INTRO MULTIVARIATE S
  • [3] Arfken G. B., 2000, MATH METHODS PHYSICI
  • [4] MR DIFFUSION TENSOR SPECTROSCOPY AND IMAGING
    BASSER, PJ
    MATTIELLO, J
    LEBIHAN, D
    [J]. BIOPHYSICAL JOURNAL, 1994, 66 (01) : 259 - 267
  • [5] ESTIMATION OF THE EFFECTIVE SELF-DIFFUSION TENSOR FROM THE NMR SPIN-ECHO
    BASSER, PJ
    MATTIELLO, J
    LEBIHAN, D
    [J]. JOURNAL OF MAGNETIC RESONANCE SERIES B, 1994, 103 (03): : 247 - 254
  • [6] BATCHELOR P, 2002, WORKSH DIFF MRI BIOP
  • [7] BATCHELOR PG, 2002, 10 ISMRM C HON HI
  • [8] Fukunaga K., 1972, ELECT SCI
  • [9] Green A. E., 1992, Theoretical Elasticity
  • [10] MEASUREMENT OF SIGNAL INTENSITIES IN THE PRESENCE OF NOISE IN MR IMAGES
    HENKELMAN, RM
    [J]. MEDICAL PHYSICS, 1985, 12 (02) : 232 - 233