Optimization on Spaces of Curves

This thesis is concerned with computational and theoretical aspects of Riemannian metrics on spaces of regular curves, and their applications. It was recently proved that second order constant coefficient Sobolev metrics on curves are geodesically complete. We extend this result to the case of Sobolev metrics with coefficient functions depending on the length of the curve.

We show how to apply this result to analyse a wide range of metrics on the submanifold of unit and constant speed curves. We present a numerical discretization of second order Sobolev metrics on the space of regular curves in Rd, and methods to solve the initial and boundary value problem for geodesics allowing us to compute the Karcher mean and principal components analysis of data of curves.

We apply the methods to study shape variation in synthetic data in the Kimia shape database, in HeLa cell nuclei and cycles of cardiac deformations. Finally we investigate a new application of Riemannian shape analysis in shape optimization. We setup a simple elliptic model problem, and describe how to apply shape calculus to obtain directional derivatives in the manifold of planar curves.

We present an implementation based on parametrization of immersions by B-splines, which ties in naturally with Isogeometric Analysis to solve the PDE. We give numerical examples of solutions, and compare the Riemannian optimization algorithms with different choices of metrics to a naive unregularized discretize-first approach.