Abstract:Background: In Diffusion Tensor Imaging (DTI), Riemannian framework based on Information Geometry theory has been proposed for processing tensors on estimation, interpolation, smoothing, regularization, segmentation, statistical test and so on. Recently Riemannian framework has been generalized to Orientation Distribution Function (ODF) and it is applicable to any Probability Density Function (PDF) under orthonormal basis representation. Spherical Polar Fourier Imaging (SPFI) was proposed for ODF and Ensemble Average Propagator (EAP) estimation from arbitrary sampled signals without any assumption.Purpose: Tensors only can represent Gaussian EAP and ODF is the radial integration of EAP, while EAP has full information for diffusion process. To our knowledge, so far there is no work on how to process EAP data. In this paper, we present a Riemannian framework as a mathematical tool for such task.Method: We propose a state-of-the-art Riemannian framework for EAPs by representing the square root of EAP, called wavefunction based on quantum mechanics, with the Fourier dual Spherical Polar Fourier (dSPF) basis. In this framework, the exponential map, logarithmic map and geodesic have closed forms, and weighted Riemannian mean and median uniquely exist. We analyze theoretically the similarities and differences between Riemannian frameworks for EAPs and for ODFs and tensors. The Riemannian metric for EAPs is diffeomorphism invariant, which is the natural extension of the affine-invariant metric for tensors. We propose Log-Euclidean framework to fast process EAPs, and Geodesic Anisotropy (GA) to measure the anisotropy of EAPs. With this framework, many important data processing operations, such as interpolation, smoothing, atlas estimation, Principal Geodesic Analysis (PGA), can be performed on EAP data.Results and Conclusions: The proposed Riemannian framework was validated in synthetic data for interpolation, smoothing, PGA and in real data for GA and atlas estimation. Riemannian median is much robust for atlas estimation.

Power Euclidean metrics for covariance matrices with application to diffusion tensor imaging

Statistics on the Manifold of Multivariate Normal Distributions: Theory and Application to Diffusion Tensor MRI Processing

An anisotropy preserving metric for DTI processing

Statistical Inference of Auto-correlated Eigenvalues with Applications to Diffusion Tensor Imaging

A Spatial Bayesian Semiparametric Mixture Model for Positive Definite Matrices with Applications to Diffusion Tensor Imaging

The Measure of Diffusion Skewness and Kurtosis in Magnetic Resonance Imaging

Geostatistical Modeling of Positive Definite Matrices: An Application to Diffusion Tensor Imaging

Distances and Riemannian Metrics for Multivariate Spectral Densities

On optimal transport of matrix-valued measures

Statistical analysis of tensor fields

Diffeomorphism Invariant Riemannian Framework for Ensemble Average Propagator Computing

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

Bures-Wasserstein Barycentric Coordinates with Application to Diffusion Tensor Image Smoothing

The Shape of Data and Probability Measures

Characterizing Distances of Networks on the Tensor Manifold.

Elastic Metrics on Spaces of Euclidean Curves: Theory and Algorithms

Tensor Canonical Correlation Analysis with Convergence and Statistical Guarantees

Varying coefficient model for modeling diffusion tensors along white matter tracts

Direct methods of B(D)-eigenpairs of tensors and application in diffusion kurtosis imaging

Matrix Kendall's tau in High-dimensions: A Robust Statistic for Matrix Factor Model

Metric Space Magnitude for Evaluating the Diversity of Latent Representations