Gaussian Process Landmarking on Manifolds

by   Tingran Gao, et al.

As a means of improving analysis of biological shapes, we propose a greedy algorithm for sampling a Riemannian manifold based on the uncertainty of a Gaussian process. This is known to produce a near optimal experimental design with the manifold as the domain, and appears to outperform the use of user-placed landmarks in representing geometry of biological objects. We provide an asymptotic analysis for the decay of the maximum conditional variance, which is frequently employed as a greedy criterion for similar variance- or uncertainty-based sequential experimental design strategies, to our knowledge this is the first result of this type for experimental design. The key observation is to link the greedy algorithm with reduced basis methods in the context of model reduction for partial differential equations. We apply the proposed landmarking algorithm to geometric morphometrics, a branch of evolutionary biology focusing on the analysis and comparisons of anatomical shapes, and compare the automatically sampled landmarks with the "ground truth" landmarks manually placed by evolutionary anthropologists, the results suggest that Gaussian process landmarks perform equally well or better, in terms of both spatial coverage and downstream statistical analysis. We expect this approach will find additional applications in other fields of research.


page 11

page 12

page 14

page 15

page 23

page 26

page 28


Gaussian Process Landmarking for Three-Dimensional Geometric Morphometrics

We demonstrate applications of the Gaussian process-based landmarking al...

Automatic Detection and Uncertainty Quantification of Landmarks on Elastic Curves

A population quantity of interest in statistical shape analysis is the l...

On the Inference of Applying Gaussian Process Modeling to a Deterministic Function

The Gaussian process modeling is a standard tool for building emulators ...

A supermartingale approach to Gaussian process based sequential design of experiments

Gaussian process (GP) models have become a well-established frameworkfor...

Optimal Design of Experiments on Riemannian Manifolds

Traditional optimal design of experiment theory is developed on Euclidea...

Gaussian Process Manifold Interpolation for Probabilistic Atrial Activation Maps and Uncertain Conduction Velocity

In patients with atrial fibrillation, local activation time (LAT) maps a...

Algorithms to automatically quantify the geometric similarity of anatomical surfaces

We describe new approaches for distances between pairs of 2-dimensional ...

Please sign up or login with your details

Forgot password? Click here to reset