Shape preserving C2C^2 interpolatory subdivision schemes

Stationary interpolatory subdivision schemes which preserve shape properties such as convexity or monotonicity are constructed. The schemes are rational in the data and generate limit functions that are at least $C^2$. The emphasis is on a class of six-point convexity preserving subdivision schemes that generate $C^2$ limit functions. In addition, a class of six-point monotonicity preserving schemes that also leads to $C^2$ limit functions is introduced. As the algebra is far too complicated for an analytical proof of smoothness, validation has been performed by a simple numerical methodology

An interpolatory subdivision algorithm for surfaces over arbitrary triangulations

In this paper, an interpolatory subdivision algorithm for surfaces over ar-bitrary triangulations is introduced and its convergence properties over nonuni-form triangulations studied. The so called Butterfly Scheme (interpolatory) is a special case of this algorithm. In our analysis of the algorithm over uniform triangulations, a matrix approach is employed and the idea, of "Cross Differ-ence of Directional Divided Difference" analysis is presented. This method is a generalization of the technique used by Dyn, Gregory and Levin etc. to analyse univariate subdivision algorithms. While for nonuniform data, an extraordi-nary point analysis is introduced and the local subdivision matrix analysis is presented. It is proved that the algorithm produces smooth surfaces over ar-bitrary triangular networks provided the shape parameters are kept within an appropriate range

Convexity preserving interpolatory subdivision with conic precision

The paper is concerned with the problem of shape preserving interpolatory subdivision. For arbitrarily spaced, planar input data an efficient non-linear subdivision algorithm is presented that results in $G^1$ limit curves, reproduces conic sections and respects the convexity properties of the initial data. Significant numerical examples illustrate the effectiveness of the proposed method

Ellipse-preserving Hermite interpolation and subdivision

We introduce a family of piecewise-exponential functions that have the Hermite interpolation property. Our design is motivated by the search for an effective scheme for the joint interpolation of points and associated tangents on a curve with the ability to perfectly reproduce ellipses. We prove that the proposed Hermite functions form a Riesz basis and that they reproduce prescribed exponential polynomials. We present a method based on Green's functions to unravel their multi-resolution and approximation-theoretic properties. Finally, we derive the corresponding vector and scalar subdivision schemes, which lend themselves to a fast implementation. The proposed vector scheme is interpolatory and level-dependent, but its asymptotic behaviour is the same as the classical cubic Hermite spline algorithm. The same convergence properties---i.e., fourth order of approximation---are hence ensured

The PCHIP subdivision scheme

In this paper we propose and analyze a nonlinear subdivision scheme based on the monotononicity-preserving third order Hermite-type interpolatory technique implemented in the PCHIP package in Matlab. We prove the convergence and the stability of the PCHIP nonlinear subdivision process by employing a novel technique based on the study of the generalized Jacobian of the first difference scheme.MTM2011-2274

Non-uniform interpolatory subdivision schemes with improved smoothness

Subdivision schemes are used to generate smooth curves or surfaces by iteratively refining an initial control polygon or mesh. We focus on univariate, linear, binary subdivision schemes, where the vertices of the refined polygon are computed as linear combinations of the current neighbouring vertices. In the classical stationary setting, there are just two such subdivision rules, which are used throughout all subdivision steps to construct the new vertices with even and odd indices, respectively. These schemes are well understood and many tools have been developed for deriving their properties, including the smoothness of the limit curves. For non-stationary schemes, the subdivision rules are not fixed and can be different in each subdivision step. Non-uniform schemes are even more general, as they allow the subdivision rules to be different for every new vertex that is generated by the scheme. The properties of non-stationary and non-uniform schemes are usually derived by relating the scheme to a corresponding stationary scheme and then exploiting the fact that the properties of the stationary scheme carry over under certain proximity conditions. In particular, this approach can be used to show that the limit curves of a non-stationary or non-uniform scheme are as smooth as those of a corresponding stationary scheme. In this paper we show that non-uniform subdivision schemes have the potential to generate limit curves that are smoother than those of stationary schemes with the same support size of the subdivision rule. For that, we derive interpolatory 2-point and 4-point schemes that generate C-1 and C-2 limit curves, respectively. These values of smoothness exceed the smoothness of classical interpolating schemes with the same support size by one. (C) 2022 The Author(s). Published by Elsevier B.V

Point-Normal Subdivision Curves and Surfaces

This paper proposes to generalize linear subdivision schemes to nonlinear subdivision schemes for curve and surface modeling by refining vertex positions together with refinement of unit control normals at the vertices. For each round of subdivision, new control normals are obtained by projections of linearly subdivided normals onto unit circle or sphere while new vertex positions are obtained by updating linearly subdivided vertices along the directions of the newly subdivided normals. Particularly, the new position of each linearly subdivided vertex is computed by weighted averages of end points of circular or helical arcs that interpolate the positions and normals at the old vertices at one ends and the newly subdivided normal at the other ends. The main features of the proposed subdivision schemes are three folds: (1) The point-normal (PN) subdivision schemes can reproduce circles, circular cylinders and spheres using control points and control normals; (2) PN subdivision schemes generalized from convergent linear subdivision schemes converge and can have the same smoothness orders as the linear schemes; (3) PN $C^2$ subdivision schemes generalizing linear subdivision schemes that generate $C^2$ subdivision surfaces with flat extraordinary points can generate visually $C^2$ subdivision surfaces with non-flat extraordinary points. Experimental examples have been given to show the effectiveness of the proposed techniques for curve and surface modeling.Comment: 30 pages, 17 figures, 22.5M

Non-linear subdivision using local spherical coordinates

In this paper, we present an original non-linear subdivision scheme suitable for univariate data, plane curves and discrete triangulated surfaces, while keeping the complexity acceptable. The proposed technique is compared to linear subdivision methods having an identical support. Numerical criteria are proposed to verify basic properties, such as convergence of the scheme and the regularity of the limit function