Project Description

The diffusion tensor MRI (DT-MRI) data has been used for studying brain structures, their connectivities, and functionalities consisting of fields of diffusion tensor matrices in the imaged volume. The problems involving denoising, interpolation, tractography, and related analyses for such data require a framework for comparing, averaging, and manipulating these tensor structures. In recent years, the increased strengths of MRI scanners have led to the possibility of measuring diffusion orientations beyond the three canonical directions, and have led to HARDI (High Angular Resolution Diffusion Imaging) technology that measures fluid flow at each voxel in numerous directions [1]. The basic unit of HARDI data at each anatomical location is an orientation diffusion function (ODF) that captures the diffusivity as a function of the direction. Using smooth interpolations (e.g. using spherical harmonics) one can obtain an ODF on the full sphere. From those limited observations, an ODF can be viewed as a non-negative function on a sphere. The space of ODFs is, thus, the set of all such functions on a sphere.

Many researchers tend to focus on the shape and orientations of the ODFs by removing the scale variability . They do so by forming a corresponding space of probability density functions (PDFs) by normalizing the space of ODFs. The space of all PDFs forms a Banach manifold, and it can become a Hilbert manifold by choosing a Riemannian structure. The Fisher-Rao Metric [2] has become increasingly popular in recent years. In HARDI data analysis, this metric was first proposed in Chang et al. [3], and later described in several other papers, including [4, 5]. In this paper, we highlight a major disadvantage of using the Fisher-Rao metric directly. The main problem is that in a direct formulation, the shape and the orientation features are intermingled and, as a result, the geodesic paths between PDFs go through PDFs that are not amenable to biological interpretations.

In our work, we propose a novel Riemannian framework for analyzing PDFs in HARDI data sets for use in comparing, interpolating, averaging, and denoising PDFs. This is accomplished by separating the shape and orientation features of PDFs, and then analyzing them separately under their own Riemannian metrics. We formulate the action of the rotation group on the space of PDFs, and define the shape space as the quotient space of PDFs modulo the rotations. In other words, any two PDFs are compared in: (1) shape by rotationally aligning one PDF to another, using the Fisher-Rao distance on the aligned PDFs, and (2) orientation by comparing their rotation matrices. By separating the shapes and the orientations of the PDFs, This allows for a formal Riemannian framework for analyzing PDFs, while generating geodesic paths that are biologically feasible. This framework leads to definitions and efficient computations for the Karcher mean and covariance that provide tools for improved interpolation and denoising.

Fisher-Rao Metric
Joint Shape-Orientation Space

Application: interpolation between given PDFs. In the above figure, we assume that we are given PDFs at the four corners of a square (in blue-green), and our goal is to interpolate PDFs (in grey) at points along a uniform grid in that square. We use unimodal densities with clear orientations (i.e. high diffusivity in a certain direction). Under the Fisher-Rao metric (a), the tract orientations at the intermediate voxels are ambiguous; in the joint shape- orientation space, the orientation is structured and consistent with the direction of the given (corner) PDFs.


  • [1] Tuch, D., Reese, T., Wiegell, M., Makris, N., Belliveau, J., Wedeen, V.: High angular resolution diffusion imaging reveals intravoxel white matter fiber heterogeneity. Magnetic Resonance in Medicine 48, 577{582 (2002)
  • [2] Rao, C.: Information and the accuracy attainable in the estimation of statistical parameters. Bull. Calcutta Match Soc. 37, 81{91 (1945)
  • [3] Chiang, M.C., et al.: Mapping genetic influences on brain fiber architecture with high angular resolution diffusion imaging HARDI. In: ISBI. pp. 871{874 (2008)
  • [4] Goh, A., Lenglet, C., Thompson, P.M., Vidal, R.: R.: A nonparametric Riemannian framework for processing high angular resolution diffusion images HARDI. In: IEEE CVPR (2009)
  • [5] Cheng, J., Ghosh, A., Jiang, T., Deriche, R.: A Riemannian framework for orientation distribution function computing. In: MICCAI. pp. 911{918 (2009)