Abstract:We propose several methods that address the problem of fitting a C1\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$C^1$$\end{document} curve γ\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\gamma $$\end{document} to time-labeled data points on a manifold. The methods have a parameter, λ\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\lambda $$\end{document}, to adjust the relative importance of the two goals that the curve should meet: being “straight enough” while fitting the data “closely enough.” The methods are designed for ease of use: they only require to compute Riemannian exponentials and logarithms, they represent the curve by means of a number of tangent vectors that grows linearly with the number of data points, and, once the representation is computed, evaluating γ(t)\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\gamma (t)$$\end{document} at any t requires a small number of exponentials and logarithms (independent of the number of data points). Among the proposed methods, the blended cubic spline technique combines the additional properties of interpolating the data when λ→∞\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\lambda \rightarrow \infty $$\end{document} and reducing to the well-known cubic smoothing spline when the manifold is Euclidean. The methods are illustrated on synthetic and real data.

Computing regularized splines in the Riemannian manifold of probability measures

A New Proof of Some Polynomial Inequalities Related to Pseudo-Splines

Cubic Hermite interpolators on the space of probability measures

A note on the regularity and the existence of Riemannian splines

Variational splines on Riemannian manifolds with applications to integral geometry

Approximation of Splines in Wasserstein Spaces

Multivariate Hermite Interpolation on Riemannian Manifolds

Smoothing splines on Riemannian manifolds, with applications to 3D shape space

Bézier curves and C2 interpolation in Riemannian Symmetric Spaces

Fitting Smoothing Splines to Time-Indexed, Noisy Points on Nonlinear Manifolds

A gradient flow method for smooth splines versus least-squares fitting on Riemannian manifolds

A Cubic Regularized Newton's Method over Riemannian Manifolds

Fast Bivariate Penalized Splines: the Sandwich Smoother

Hessian Based Smoothing Splines for Manifold Learning

Generalized Moving Least-Squares for Solving Vector-valued PDEs on Unknown Manifolds

Splines on Manifolds: a Survey

Inexact Adaptive Cubic Regularization Algorithms on Riemannian Manifolds and Application

Nonconforming Spline Collocation Methods in Irregular Domains II: Error Analysis

Data Fitting on Manifolds with Composite Bézier-Like Curves and Blended Cubic Splines

Sobolev Metrics on Spaces of Discrete Regular Curves