Geodesic spline interpolation is a simple and efficient approach for landmark matching by nonambiguous mappings (diffeomorphisms), combining classic spline interpolation and flows of diffeomorphisms. Here, we extend the method to incorporate the estimation of a affine transformation, yielding a consistent and numerically stable algorithm. A theoretical justification is provided by studying the existence of the global minimum of the energy.