A Riemannian Framework for Orientation Distribution Function Computing

International audienceCompared with Diffusion Tensor Imaging (DTI), High Angular Resolution Imaging (HARDI) can better explore the complex microstructure of white matter. Orientation Distribution Function (ODF) is used to describe the probability of the fiber direction. Fisher information metric has been constructed for probability density family in Information Geometry theory and it has been successfully applied for tensor computing in DTI. In this paper, we present a state of the art Riemannian framework for ODF computing based on Information Geometry and sparse representation of orthonormal bases. In this Riemannian framework, the exponential map, logarithmic map and geodesic have closed forms. And the weighted Frechet mean exists uniquely on this manifold. We also propose a novel scalar measurement, named Geometric Anisotropy (GA), which is the Riemannian geodesic distance between the ODF and the isotropic ODF. The Renyi entropy H_{1/2} of the ODF can be computed from the GA. Moreover, we present an Affine-Euclidean framework and a Log-Euclidean framework so that we can work in an Euclidean space. As an application, Lagrange interpolation on ODF field is proposed based on weighted Frechet mean. We validate our methods on synthetic and real data experiments. Compared with existing Riemannian frameworks on ODF, our framework is model-free. The estimation of the parameters, i.e. Riemannian coordinates, is robust and linear. Moreover it should be noted that our theoretical results can be used for any probability density function (PDF) under an orthonormal basis representation

k-Nearest Neighbor Based Consistent Entropy Estimation for Hyperspherical Distributions

A consistent entropy estimator for hyperspherical data is proposed based on the k-nearest neighbor (knn) approach. The asymptotic unbiasedness and consistency of the estimator are proved. Moreover, cross entropy and Kullback-Leibler (KL) divergence estimators are also discussed. Simulation studies are conducted to assess the performance of the estimators for models including uniform and von Mises-Fisher distributions. The proposed knn entropy estimator is compared with the moment based counterpart via simulations. The results show that these two methods are comparable

Directional functions for orientation distribution estimation

Computing the Orientation Distribution Function (ODF) from High Angular Resolution Diffusion Imaging (HARDI) signals makes it possible to determine the orientation of fiber bundles of the brain. The HARDI signals are samples measured from a spherical shell and thus require processing on the sphere. Past work on ODF estimation involved using the spherical harmonics or spherical radial basis functions. In this work, we propose three novel directional functions able to represent the measured signals in a very compact manner, i.e., they require very few parameters to completely describe the measured signal. Analytical expressions are derived for computing the corresponding ODF. The directional functions can represent diffusion in a particular direction and mixture models can be used to represent multi-fiber orientations. We show how to estimate the parameters of this mixture model and elaborate on the differences between these functions. We also compare this general framework with estimation of ODF using spherical harmonics on some real and synthetic data. The proposed method could be particularly useful in applications such as tractography and segmentation. Details are also given on different ways in which interpolation can be performed using directional functions. In particular, we discuss a complete Euclidean as well as a “hybrid” framework, comprising of the Riemannian as well as Euclidean spaces, to perform interpolation and compute geodesic distances between 2 ODF’s

Computational Representation of White Matter Fiber Orientations

We present a new methodology based on directional data clustering to represent white matter fiber orientations in magnetic resonance analyses for high angular resolution diffusion imaging. A probabilistic methodology is proposed for estimating intravoxel principal fiber directions, based on clustering directional data arising from orientation distribution function (ODF) profiles. ODF reconstructions are used to estimate intravoxel fiber directions using mixtures of von Mises-Fisher distributions. The method focuses on clustering data on the unit sphere, where complexity arises from representing ODF profiles as directional data. The proposed method is validated on synthetic simulations, as well as on a real data experiment. Based on experiments, we show that by clustering profile data using mixtures of von Mises-Fisher distributions it is possible to estimate multiple fiber configurations in a more robust manner than currently used approaches, without recourse to regularization or sharpening procedures. The method holds promise to support robust tractographic methodologies and to build realistic models of white matter tracts in the human brain

Modeling diffusion directions of Corpus Callosum

Diffusion Tensor Imaging (DTI) has been used to study the characteristics of Multiple Sclerosis (MS) in the brain. The von Mises- Fisher distribution (vmf) is a probability distribution for modeling directional data on the unit hypersphere. In this paper we modeled the diffusion directions of the Corpus Callosum (CC) as a mixture of vmf distributions for both MS subjects and healthy controls. Higher diffusion concentration around the mean directions and smaller sum of angles between the mean directions are observed on the normal-appearing CC of the MS subjects as compared to the healthy controls

Von Mises-Fisher models in the total variability subspace for language recognition

Personal use of this material is permitted. Permission from IEEE must be obtained for all other uses, in any current or future media, including reprinting/republishing this material for advertising or promotional purposes, creating new collective works, for resale or redistribution to servers or lists, or reuse of any copyrighted component of this work in other works. I. Lopez-Moreno, D. Ramos, J. Gonzalez-Dominguez, and J. Gonzalez-Rodriguez, "Von Mises-Fisher models in the total variability subspace for language recognition", IEEE Signal Processing Letters, vol. 18, no. 12, pp. 705-708, October 2011This letter proposes a new modeling approach for the Total Variability subspace within a Language Recognition task. Motivated by previous works in directional statistics, von Mises-Fisher distributions are used for assigning language-conditioned probabilities to language data, assumed to be spherically distributed in this subspace. The two proposed methods use Kernel Density Functions or Finite Mixture Models of such distributions. Experiments conducted on NIST LRE 2009 show that the proposed techniques significantly outperform the baseline cosine distance approach in most of the considered experimental conditions, including different speech conditions, durations and the presence of unseen languages.This work was supported by the Ministerio de Ciencia e Innovación under FPI Grant TEC2009-14719-C02-01 and cátedra UAM-Telefónic