High quality local interpolation by composite parametric surfaces

The proposed method is generalization of the local univariate spline interpolants introduced in Beccari et al. (2013a) and Antonelli et al. (2014), that have arbitrary continuity and arbitrary order of polynomial reproduction. The generated surfaces retain the same smoothness of the underlying class of univariate splines in the regular regions of the mesh (where, locally, all vertices have valence 4). Mesh regions containing vertices of valence other than 4 are covered with suitably defined surface patches joining the neighboring regular ones with G1- or G2-continuity.