-
Learning from landmarks, curves, surfaces, and shapes in Geomstats
Authors:
Luís F. Pereira,
Alice Le Brigant,
Adele Myers,
Emmanuel Hartman,
Amil Khan,
Malik Tuerkoen,
Trey Dold,
Mengyang Gu,
Pablo Suárez-Serrato,
Nina Miolane
Abstract:
We introduce the shape module of the Python package Geomstats to analyze shapes of objects represented as landmarks, curves and surfaces across fields of natural sciences and engineering. The shape module first implements widely used shape spaces, such as the Kendall shape space, as well as elastic spaces of discrete curves and surfaces. The shape module further implements the abstract mathematica…
▽ More
We introduce the shape module of the Python package Geomstats to analyze shapes of objects represented as landmarks, curves and surfaces across fields of natural sciences and engineering. The shape module first implements widely used shape spaces, such as the Kendall shape space, as well as elastic spaces of discrete curves and surfaces. The shape module further implements the abstract mathematical structures of group actions, fiber bundles, quotient spaces and associated Riemannian metrics which allow users to build their own shape spaces. The Riemannian geometry tools enable users to compare, average, interpolate between shapes inside a given shape space. These essential operations can then be leveraged to perform statistics and machine learning on shape data. We present the object-oriented implementation of the shape module along with illustrative examples and show how it can be used to perform statistics and machine learning on shape spaces.
△ Less
Submitted 14 June, 2024;
originally announced June 2024.
-
Optimal matching between curves in a manifold
Authors:
Alice Le Brigant,
Marc Arnaudon,
Frédéric Barbaresco
Abstract:
This paper is concerned with the computation of an optimal matching between two manifold-valued curves. Curves are seen as elements of an infinite-dimensional manifold and compared using a Riemannian metric that is invariant under the action of the reparameterization group. This group induces a quotient structure classically interpreted as the ''shape space''. We introduce a simple algorithm allow…
▽ More
This paper is concerned with the computation of an optimal matching between two manifold-valued curves. Curves are seen as elements of an infinite-dimensional manifold and compared using a Riemannian metric that is invariant under the action of the reparameterization group. This group induces a quotient structure classically interpreted as the ''shape space''. We introduce a simple algorithm allowing to compute geodesics of the quotient shape space using a canonical decomposition of a path in the associated principal bundle. We consider the particular case of elastic metrics and show simulations for open curves in the plane, the hyperbolic plane and the sphere.
△ Less
Submitted 9 January, 2024;
originally announced January 2024.
-
The $L^p$-Fisher-Rao metric and Amari-Cencov $α$-connections
Authors:
Martin Bauer,
Alice Le Brigant,
Yuxiu Lu,
Cy Maor
Abstract:
We introduce a family of Finsler metrics, called the $L^p$-Fisher-Rao metrics $F_p$, for $p\in (1,\infty)$, which generalizes the classical Fisher-Rao metric $F_2$, both on the space of densities Dens$_+(M)$ and probability densities Prob$(M)$. We then study their relations to the Amari-uCencov $α$-connections $\nabla^{(α)}$ from information geometry: on Dens$_+(M)$, the geodesic equations of…
▽ More
We introduce a family of Finsler metrics, called the $L^p$-Fisher-Rao metrics $F_p$, for $p\in (1,\infty)$, which generalizes the classical Fisher-Rao metric $F_2$, both on the space of densities Dens$_+(M)$ and probability densities Prob$(M)$. We then study their relations to the Amari-uCencov $α$-connections $\nabla^{(α)}$ from information geometry: on Dens$_+(M)$, the geodesic equations of $F_p$ and $\nabla^{(α)}$ coincide, for $p = 2/(1-α)$. Both are pullbacks of canonical constructions on $L^p(M)$, in which geodesics are simply straight lines. In particular, this gives a new variational interpretation of $α$-geodesics as being energy minimizing curves. On Prob$(M)$, the $F_p$ and $\nabla^{(α)}$ geodesics can still be thought as pullbacks of natural operations on the unit sphere in $L^p(M)$, but in this case they no longer coincide unless $p=2$. Using this transformation, we solve the geodesic equation of the $α$-connection by showing that the geodesic are pullbacks of projections of straight lines onto the unit sphere, and they always cease to exists after finite time when they leave the positive part of the sphere. This unveils the geometric structure of solutions to the generalized Proudman-Johnson equations, and generalizes them to higher dimensions. In addition, we calculate the associate tensors of $F_p$, and study their relation to $\nabla^{(α)}$.
△ Less
Submitted 31 July, 2023; v1 submitted 26 June, 2023;
originally announced June 2023.
-
Conjugate points along Kolmogorov flows on the torus
Authors:
Alice Le Brigant,
Stephen C. Preston
Abstract:
The geodesics in the group of volume-preserving diffeomorphisms (volumorphisms) of a manifold $M$, for a Riemannian metric defined by the kinetic energy, can be used to model the movement of ideal fluids in that manifold. The existence of conjugate points along such geodesics reveal that these cease to be infinitesimally length-minimizing between their endpoints. In this work, we focus on the case…
▽ More
The geodesics in the group of volume-preserving diffeomorphisms (volumorphisms) of a manifold $M$, for a Riemannian metric defined by the kinetic energy, can be used to model the movement of ideal fluids in that manifold. The existence of conjugate points along such geodesics reveal that these cease to be infinitesimally length-minimizing between their endpoints. In this work, we focus on the case of the torus $M=\T^2$ and on geodesics corresponding to steady solutions of the Euler equation generated by stream functions $ψ=-\cos(mx)\cos(ny)$ for integers $m$ and $n$, called Kolmogorov flows. We show the existence of conjugate points along these geodesics for all pairs of strictly positive integers $(m,n)$, thereby completing the characterization of all pairs $(m,n)$ such that the associated Kolmogorov flow generates a geodesic with conjugate points.
△ Less
Submitted 5 December, 2023; v1 submitted 12 April, 2023;
originally announced April 2023.
-
Parametric information geometry with the package Geomstats
Authors:
Alice Le Brigant,
Jules Deschamps,
Antoine Collas,
Nina Miolane
Abstract:
We introduce the information geometry module of the Python package Geomstats. The module first implements Fisher-Rao Riemannian manifolds of widely used parametric families of probability distributions, such as normal, gamma, beta, Dirichlet distributions, and more. The module further gives the Fisher-Rao Riemannian geometry of any parametric family of distributions of interest, given a parameteri…
▽ More
We introduce the information geometry module of the Python package Geomstats. The module first implements Fisher-Rao Riemannian manifolds of widely used parametric families of probability distributions, such as normal, gamma, beta, Dirichlet distributions, and more. The module further gives the Fisher-Rao Riemannian geometry of any parametric family of distributions of interest, given a parameterized probability density function as input. The implemented Riemannian geometry tools allow users to compare, average, interpolate between distributions inside a given family. Importantly, such capabilities open the door to statistics and machine learning on probability distributions. We present the object-oriented implementation of the module along with illustrative examples and show how it can be used to perform learning on manifolds of parametric probability distributions.
△ Less
Submitted 21 November, 2022;
originally announced November 2022.
-
Classifying histograms of medical data using information geometry of beta distributions
Authors:
Alice Le Brigant,
Nicolas Guigui,
Sana Rebbah,
Stéphane Puechmorel
Abstract:
In this paper, we use tools of information geometry to compare, average and classify histograms. Beta distributions are fitted to the histograms and the corresponding Fisher information geometry is used for comparison. We show that this geometry is negatively curved, which guarantees uniqueness of the notion of mean, and makes it suitable to classify histograms through the popular K-means algorith…
▽ More
In this paper, we use tools of information geometry to compare, average and classify histograms. Beta distributions are fitted to the histograms and the corresponding Fisher information geometry is used for comparison. We show that this geometry is negatively curved, which guarantees uniqueness of the notion of mean, and makes it suitable to classify histograms through the popular K-means algorithm. We illustrate the use of these geometric tools in supervised and unsupervised classification procedures of two medical data-sets, cardiac shape deformations for the detection of pulmonary hypertension and brain cortical thickness for the diagnosis of Alzheimer's disease.
△ Less
Submitted 25 August, 2021; v1 submitted 3 June, 2020;
originally announced June 2020.
-
Fisher-Rao geometry of Dirichlet distributions
Authors:
Alice Le Brigant,
Stephen Preston,
Stéphane Puechmorel
Abstract:
In this paper, we study the geometry induced by the Fisher-Rao metric on the parameter space of Dirichlet distributions. We show that this space is geodesically complete and has everywhere negative sectional curvature. An important consequence of this negative curvature for applications is that the Fr{é}chet mean of a set of Dirichlet distributions is uniquely defined in this geometry.
In this paper, we study the geometry induced by the Fisher-Rao metric on the parameter space of Dirichlet distributions. We show that this space is geodesically complete and has everywhere negative sectional curvature. An important consequence of this negative curvature for applications is that the Fr{é}chet mean of a set of Dirichlet distributions is uniquely defined in this geometry.
△ Less
Submitted 19 November, 2020; v1 submitted 12 May, 2020;
originally announced May 2020.
-
Geomstats: A Python Package for Riemannian Geometry in Machine Learning
Authors:
Nina Miolane,
Alice Le Brigant,
Johan Mathe,
Benjamin Hou,
Nicolas Guigui,
Yann Thanwerdas,
Stefan Heyder,
Olivier Peltre,
Niklas Koep,
Hadi Zaatiti,
Hatem Hajri,
Yann Cabanes,
Thomas Gerald,
Paul Chauchat,
Christian Shewmake,
Bernhard Kainz,
Claire Donnat,
Susan Holmes,
Xavier Pennec
Abstract:
We introduce Geomstats, an open-source Python toolbox for computations and statistics on nonlinear manifolds, such as hyperbolic spaces, spaces of symmetric positive definite matrices, Lie groups of transformations, and many more. We provide object-oriented and extensively unit-tested implementations. Among others, manifolds come equipped with families of Riemannian metrics, with associated expone…
▽ More
We introduce Geomstats, an open-source Python toolbox for computations and statistics on nonlinear manifolds, such as hyperbolic spaces, spaces of symmetric positive definite matrices, Lie groups of transformations, and many more. We provide object-oriented and extensively unit-tested implementations. Among others, manifolds come equipped with families of Riemannian metrics, with associated exponential and logarithmic maps, geodesics and parallel transport. Statistics and learning algorithms provide methods for estimation, clustering and dimension reduction on manifolds. All associated operations are vectorized for batch computation and provide support for different execution backends, namely NumPy, PyTorch and TensorFlow, enabling GPU acceleration. This paper presents the package, compares it with related libraries and provides relevant code examples. We show that Geomstats provides reliable building blocks to foster research in differential geometry and statistics, and to democratize the use of Riemannian geometry in machine learning applications. The source code is freely available under the MIT license at \url{geomstats.ai}.
△ Less
Submitted 7 April, 2020;
originally announced April 2020.
-
Intrinsic Riemannian metrics on spaces of curves: theory and computation
Authors:
Martin Bauer,
Nicolas Charon,
Eric Klassen,
Alice Le Brigant
Abstract:
This chapter reviews some past and recent developments in shape comparison and analysis of curves based on the computation of intrinsic Riemannian metrics on the space of curves modulo shape-preserving transformations. We summarize the general construction and theoretical properties of quotient elastic metrics for Euclidean as well as non-Euclidean curves before considering the special case of the…
▽ More
This chapter reviews some past and recent developments in shape comparison and analysis of curves based on the computation of intrinsic Riemannian metrics on the space of curves modulo shape-preserving transformations. We summarize the general construction and theoretical properties of quotient elastic metrics for Euclidean as well as non-Euclidean curves before considering the special case of the square root velocity metric for which the expression of the resulting distance simplifies through a particular transformation. We then examine different numerical approaches that have been proposed to estimate such distances in practice and in particular to quotient out curve reparametrization in the resulting minimization problems.
△ Less
Submitted 21 October, 2020; v1 submitted 11 March, 2020;
originally announced March 2020.
-
The Fisher-Rao geometry of beta distributions applied to the study of canonical moments
Authors:
Alice Le Brigant,
Stéphane Puechmorel
Abstract:
This paper studies the Fisher-Rao geometry on the parameter space of beta distributions. We derive the geodesic equations and the sectional curvature, and prove that it is negative. This leads to uniqueness for the Riemannian centroid in that space. We use this Riemannian structure to study canonical moments, an intrinsic representation of the moments of a probability distribution. Drawing on the…
▽ More
This paper studies the Fisher-Rao geometry on the parameter space of beta distributions. We derive the geodesic equations and the sectional curvature, and prove that it is negative. This leads to uniqueness for the Riemannian centroid in that space. We use this Riemannian structure to study canonical moments, an intrinsic representation of the moments of a probability distribution. Drawing on the fact that a uniform distribution in the regular moment space corresponds to a product of beta distributions in the canonical moment space, we propose a mapping from the space of canonical moments to the product beta manifold, allowing us to use the Fisher-Rao geometry of beta distributions to compare and analyze canonical moments.
△ Less
Submitted 17 April, 2019;
originally announced April 2019.
-
Math in the Black Forest: Workshop on New Directions in Shape Analysis
Authors:
Martin Bauer,
Nicolas Charon,
Philipp Harms,
Boris Khesin,
Alice Le Brigant,
Elodie Maignant,
Stephen Marsland,
Peter Michor,
Xavier Pennec,
Stephen Preston,
Stefan Sommer,
François-Xavier Vialard
Abstract:
These are the proceedings of the workshop "Math in the Black Forest", which brought together researchers in shape analysis to discuss promising new directions. Shape analysis is an inter-disciplinary area of research with theoretical foundations in infinite-dimensional Riemannian geometry, geometric statistics, and geometric stochastics, and with applications in medical imaging, evolutionary devel…
▽ More
These are the proceedings of the workshop "Math in the Black Forest", which brought together researchers in shape analysis to discuss promising new directions. Shape analysis is an inter-disciplinary area of research with theoretical foundations in infinite-dimensional Riemannian geometry, geometric statistics, and geometric stochastics, and with applications in medical imaging, evolutionary development, and fluid dynamics. The workshop is the 6th instance of a series of workshops on the same topic.
△ Less
Submitted 4 November, 2018;
originally announced November 2018.
-
Optimal Riemannian quantization with an application to air traffic analysis
Authors:
Alice Le Brigant,
Stéphane Puechmorel
Abstract:
The goal of optimal quantization is to find the best approximation of a probability distribution by a discrete measure with finite support. When dealing with empirical distributions, this boils down to finding the best summary of the data by a smaller number of points, and automatically yields a K-means-type clustering. In this paper, we introduce Competitive Learning Riemannian Quantization (CLRQ…
▽ More
The goal of optimal quantization is to find the best approximation of a probability distribution by a discrete measure with finite support. When dealing with empirical distributions, this boils down to finding the best summary of the data by a smaller number of points, and automatically yields a K-means-type clustering. In this paper, we introduce Competitive Learning Riemannian Quantization (CLRQ), an online algorithm that computes the optimal summary when the data does not belong to a vector space, but rather a Riemannian manifold. We prove its convergence and show simulated examples on the sphere and the hyperbolic plane. We also provide an application to real data by using CLRQ to create summaries of images of covariance matrices estimated from air traffic images. These summaries are representative of the air traffic complexity and yield clusterings of the airspaces into zones that are homogeneous with respect to that criterion. They can then be compared using discrete optimal transport and be further used as inputs of a machine learning algorithm or as indexes in a traffic database.
△ Less
Submitted 2 July, 2018; v1 submitted 20 June, 2018;
originally announced June 2018.
-
A discrete framework to find the optimal matching between manifold-valued curves
Authors:
Alice Le Brigant
Abstract:
The aim of this paper is to find an optimal matching between manifold-valued curves, and thereby adequately compare their shapes, seen as equivalent classes with respect to the action of reparameterization. Using a canonical decomposition of a path in a principal bundle, we introduce a simple algorithm that finds an optimal matching between two curves by computing the geodesic of the infinite-dime…
▽ More
The aim of this paper is to find an optimal matching between manifold-valued curves, and thereby adequately compare their shapes, seen as equivalent classes with respect to the action of reparameterization. Using a canonical decomposition of a path in a principal bundle, we introduce a simple algorithm that finds an optimal matching between two curves by computing the geodesic of the infinite-dimensional manifold of curves that is at all time horizontal to the fibers of the shape bundle. We focus on the elastic metric studied in the so-called square root velocity framework. The quotient structure of the shape bundle is examined, and in particular horizontality with respect to the fibers. These results are more generally given for any elastic metric. We then introduce a comprehensive discrete framework which correctly approximates the smooth setting when the base manifold has constant sectional curvature. It is itself a Riemannian structure on the product manifold of "discrete curves" given by a finite number of points, and we show its convergence to the continuous model as the size of the discretization goes to infinity. Illustrations of optimal matching between discrete curves are given in the hyperbolic plane, the plane and the sphere, for synthetic and real data, and comparison with dynamic programming is established.
△ Less
Submitted 19 January, 2018; v1 submitted 15 March, 2017;
originally announced March 2017.
-
Computing distances and geodesics between manifold-valued curves in the SRV framework
Authors:
Alice Le Brigant
Abstract:
This paper focuses on the study of open curves in a Riemannian manifold M, and proposes a reparametrization invariant metric on the space of such paths. We use the square root velocity function (SRVF) introduced by Srivastava et al. to define a Riemannian metric on the space of immersions M'=Imm([0,1],M) by pullback of a natural metric on the tangent bundle TM'. This induces a first-order Sobolev…
▽ More
This paper focuses on the study of open curves in a Riemannian manifold M, and proposes a reparametrization invariant metric on the space of such paths. We use the square root velocity function (SRVF) introduced by Srivastava et al. to define a Riemannian metric on the space of immersions M'=Imm([0,1],M) by pullback of a natural metric on the tangent bundle TM'. This induces a first-order Sobolev metric on M' and leads to a distance which takes into account the distance between the origins in M and the L2-distance between the SRV representations of the curves. The geodesic equations for this metric are given and exploited to define an exponential map on M'. The optimal deformation of one curve into another can then be constructed using geodesic shooting, which requires to characterize the Jacobi fields of M'. The particular case of curves lying in the hyperbolic half-plane is considered as an example, in the setting of radar signal processing.
△ Less
Submitted 15 March, 2017; v1 submitted 11 January, 2016;
originally announced January 2016.
-
Reparameterization invariant metric on the space of curves
Authors:
Alice Le Brigant,
Marc Arnaudon,
Frédéric Barbaresco
Abstract:
This paper focuses on the study of open curves in a manifold M, and proposes a reparameterization invariant metric on the space of such paths. We use the square root velocity function (SRVF) introduced by Srivastava et al. in [11] to define a reparameterization invariant metric on the space of immersions M' = Imm([0,1], M) by pullback of a metric on the tangent bundle TM' derived from the Sasaki m…
▽ More
This paper focuses on the study of open curves in a manifold M, and proposes a reparameterization invariant metric on the space of such paths. We use the square root velocity function (SRVF) introduced by Srivastava et al. in [11] to define a reparameterization invariant metric on the space of immersions M' = Imm([0,1], M) by pullback of a metric on the tangent bundle TM' derived from the Sasaki metric. We observe that such a natural choice of Riemannian metric on TM' induces a first-order Sobolev metric on M' with an extra term involving the origins, and leads to a distance which takes into account the distance between the origins and the distance between the SRV representations of the curves. The geodesic equations for this metric are given, as well as an idea of how to compute the exponential map for observed trajectories in applications. This provides a generalized theoretical SRV framework for curves lying in a general manifold M .
△ Less
Submitted 27 October, 2015; v1 submitted 23 July, 2015;
originally announced July 2015.
