Abstract:We explore the connection between two problems that have arisen independently in the signal processing and related fields: the estimation of the geometric mean of a set of symmetric positive definite (SPD) matrices and their approximate joint diagonalization (AJD). Today there is a considerable interest in estimating the geometric mean of a SPD matrix set in the manifold of SPD matrices endowed with the Fisher information metric. The resulting mean has several important invariance properties and has proven very useful in diverse engineering applications such as biomedical and image data processing. While for two SPD matrices the mean has an algebraic closed form solution, for a set of more than two SPD matrices it can only be estimated by iterative algorithms. However, none of the existing iterative algorithms feature at the same time fast convergence, low computational complexity per iteration and guarantee of convergence. For this reason, recently other definitions of geometric mean based on symmetric divergence measures, such as the Bhattacharyya divergence, have been considered. The resulting means, although possibly useful in practice, do not satisfy all desirable invariance properties. In this paper we consider geometric means of covariance matrices estimated on high-dimensional time-series, assuming that the data is generated according to an instantaneous mixing model, which is very common in signal processing. We show that in these circumstances we can approximate the Fisher information geometric mean by employing an efficient AJD algorithm. Our approximation is in general much closer to the Fisher information geometric mean as compared to its competitors and verifies many invariance properties. Furthermore, convergence is guaranteed, the computational complexity is low and the convergence rate is quadratic. The accuracy of this new geometric mean approximation is demonstrated by means of simulations.

Differential Geometry with Extreme Eigenvalues in the Positive Semidefinite Cone

Differential geometry with extreme eigenvalues in the positive semidefinite cone

Geometric statistics with subspace structure preservation for SPD matrices

Kernel Methods on the Riemannian Manifold of Symmetric Positive Definite Matrices

Totally geodesic submanifolds in the manifold SPD of symmetric positive-definite real matrices

Scaling-rotation distance and interpolation of symmetric positive-definite matrices

Extreme eigenvalues of sample covariance matrices under generalized elliptical models with applications

Geometric Characteristics of the Wasserstein Metric on SPD(n) and Its Applications on Data Processing

A New Riemannian Structure in SPD(n)

Positive definite matrices and the S-divergence

Differentially Private Fréchet Mean on the Manifold of Symmetric Positive Definite (SPD) Matrices with log-Euclidean Metric

Extreme Eigenvalues of Large-Dimensional Spiked Fisher Matrices with Application

Dimensionality Reduction on SPD Manifolds: The Emergence of Geometry-Aware Methods

Positive semidefinite kernels that are axially symmetric on the sphere and stationary in time: spectral and semi-spectral theory, and constructive approaches

On Semi-Definiteness and Minimal H-eigenvalue of a Symmetric Space Tensor Using Nonnegative Polynomial Optimization Techniques.

Approximate joint diagonalization and geometric mean of symmetric positive definite matrices

Online Dictionary Learning On Symmetric Positive Definite Manifolds With Vision Applications

Topological analysis for 3D real, symmetric second-order tensor fields using Deviatoric Eigenvalue Wheel

Learning Log-Determinant Divergences for Positive Definite Matrices

Geodesic Convexity of the Symmetric Eigenvalue Problem and Convergence of Steepest Descent

The Geometry of Differential Harnack Estimates