Geodesic regression on spheres from a numerical optimization viewpoint

Abstract

In the geodesic regression problem it is given a set of data points at known times and the goal is to find a geodesic that best fits the data. This problem corresponds to the generalization of the classical linear regression problem to curved spaces. Here we are interested in the geodesic regression problem on Euclidean spheres. Contrary to the Euclidean situation, the normal equations turn out to be highly nonlinear. To overcome this difficulty, we look at the geodesic regression problem in the unit n-sphere as an optimization problem in the Euclidean space Double-struck capital Rn+1 and use the MATLAB optimization toolbox to solve it numerically.