High accuracy geometric Hermite interpolation

We describe a parametric cubic spline interpolation scheme for planar curves which is based on an idea of Sabin for the construction of C¹ bicubic parametric spline surfaces. The method is a natural generalization of standard] Hermite interpolation. In addition to position and tangent, the curvature is prescribed at each knot. This ensures that the resulting interpolating piecewise cubic curve is twice continuously differentiable with respect to arclength and can be constructed locally. Moreover, under appropriate assumptions, the interpolant preserves convexity and is 6-th order accurate.