Riemannian Median and Its Estimation

In this paper, we define the geometric median of a probability measure on a Riemannian manifold, give its characterization and a natural condition to ensure its uniqueness. In order to calculate the median in practical cases, we also propose a subgradient algorithm and prove its convergence as well as estimating the error of approximation and the rate of convergence. The convergence property of this subgradient algorithm, which is a generalization of the classical Weiszfeld algorithm in Euclidean spaces to the context of Riemannian manifolds, also answers a recent question in P. T. Fletcher et al. [13]

[1]  I. Holopainen.  Riemannian Geometry , 1927, Nature.

[2]  Harold W. Kuhn,et al.  A note on Fermat's problem , 1973, Math. Program..

[3]  H. Karcher.  Riemannian center of mass and mollifier smoothing , 1977 .

[4]  Lawrence M. Ostresh.  On the Convergence of a Class of Iterative Methods for Solving the Weber Location Problem , 1978, Oper. Res..

[5]  L. Cooper,et al.  Optimal location on a sphere , 1980 .

[6]  Igor Nikolaev,et al.  Generalized Riemannian spaces , 1986 .

[7]  W. Kendall.  Probability, Convexity, and Harmonic Maps with Small Image I: Uniqueness and Fine Existence , 1990 .

[8]  Séminaire de probabilités,et al.  Séminaire de probabilités XXV , 1991 .

[9]  M. Shirosaki.  Another proof of the defect relation for moving targets , 1991 .

[10]  M. Émery,et al.  Sur le barycentre d'une probabilité dans une variété , 1991 .

[11]  W. Kendall.  Convexity and the Hemisphere , 1991 .

[12]  Claude Lemaréchal,et al.  Convergence of some algorithms for convex minimization , 1993, Math. Program..

[13]  Marc Arnaudon.  Espérances conditionnelles et C-martingales dans les variétés , 1994 .

[14]  J. Eschenburg.  Comparison Theorems in Riemannian Geometry , 1994 .

[15]  J. Picard.  Barycentres et martingales sur une variété , 1994 .

[16]  Séminaire de probabilités,et al.  Séminaire de probabilités XXVIII , 1994 .

[17]  C. Udriste,et al.  Convex Functions and Optimization Methods on Riemannian Manifolds , 1994 .

[18]  J. Jost.  Riemannian geometry and geometric analysis , 1995 .

[19]  Michel Émery,et al.  Séminaire de probabilités XXIX , 1995 .

[20]  C. Small.  Multidimensional medians arising from geodesics on graphs , 1997 .

[21]  Espérances d'une variable aléatoire à valeurs dans un espace métrique , 1998 .

[22]  O. P. Ferreira,et al.  Subgradient Algorithm on Riemannian Manifolds , 1998 .

[23]  M. Bridson,et al.  Metric Spaces of Non-Positive Curvature , 1999 .

[24]  D. Bertsekas,et al.  Incremental subgradient methods for nondifferentiable optimization , 1999, Proceedings of the 38th IEEE Conference on Decision and Control (Cat. No.99CH36304).

[25]  D. Bertsekas,et al.  Convergen e Rate of In remental Subgradient Algorithms , 2000 .

[26]  Dimitri P. Bertsekas,et al.  Incremental Subgradient Methods for Nondifferentiable Optimization , 2001, SIAM J. Optim..

[27]  S. Uryasev,et al.  Stochastic optimization : Algorithms and Applications , 2001 .

[28]  P. Thomas Fletcher,et al.  Statistics of shape via principal geodesic analysis on Lie groups , 2003, 2003 IEEE Computer Society Conference on Computer Vision and Pattern Recognition, 2003. Proceedings..

[29]  H. Le,et al.  Estimation of Riemannian Barycentres , 2004 .

[30]  P. Thomas Fletcher,et al.  Principal Geodesic Analysis on Symmetric Spaces: Statistics of Diffusion Tensors , 2004, ECCV Workshops CVAMIA and MMBIA.

[31]  Marc Arnaudon,et al.  Barycenters of measures transported by stochastic flows , 2005 .

[32]  Xavier Pennec,et al.  Intrinsic Statistics on Riemannian Manifolds: Basic Tools for Geometric Measurements , 2006, Journal of Mathematical Imaging and Vision.

[33]  F. Barbaresco.  Innovative tools for radar signal processing Based on Cartan’s geometry of SPD matrices & Information Geometry , 2008, 2008 IEEE Radar Conference.

[34]  Suresh Venkatasubramanian,et al.  The geometric median on Riemannian manifolds with application to robust atlas estimation , 2009, NeuroImage.