Construction and Application of Subdivision Surface Scheme Using Lagrange Interpolation Polynomial

This paper offers a general formula for surface subdivision rules for quad meshes by using 2-D Lagrange interpolating polynomial [1]. We also see that the result obtained is equivalent to the tensor product of (2N + 4)-point n-ary interpolating curve scheme for N ≥ 0 and n ≥ 2. The simple interpolatory subdivision scheme for quadrilateral nets with arbitrary topology is presented by L. Kobbelt [2], which can be directly calculated from the proposed formula. Furthermore, some characteristics and applications of the proposed work are also discussed.


Introduction
There are two general classes of subdivision schemes, namely, approximating and interpolating schemes.The limit curve of an approximating scheme usually does not pass through the control points of control polygon.As the level of refinement increases, the polygon usually shrinks towards the final limit curve.The interpolating schemes are more attractive than approximating schemes because of their interpolation property.All vertices in the control polygon are located on the limit curve of the interpolation scheme, which facilitates and simplifies the graphics algorithms and engineering designs.
Lian generalized the classical binary 4-point and 6-point interpolatory subdivision schemes to a-ary setting for any integer a ≥ 3.After that, the a-ary 3-point and 5-point interpolatory subdivision schemes for curve design for arbitrary odd integer a ≥ 3 [3,4] were introduced.After that, Lian [5] investigated both the 2m-point, a-ary for any a ≥ 2 and (2m + 1)-point, a-ary for any odd a ≥ 3 interpolatory subdivision schemes for curve design.Ko [6] presented explicitly a new formula for the mask of (2N + 4)-point binary interpolating and approximating subdivision schemes with two parameters.The proposed work presents a new observation about the curve case given by Najma [7].In this work, we avoid finding the mask of subdivision schemes separately, as a result, its approach is simple and avoids complex computation when deriving subdivision rules.
The rest of the paper is organized as follows.Section 2 gives some preliminaries results and a new relation for (2N + 4)-point n-ary interpolating curve scheme for closed and open polygon to access main result.Section 3 presents the construction for general formula of the surface case using Lagrange interpolating polynomial, and some characteristics are also discussed.In Section 4, we also give some numerical examples for the visual performance of the proposed work.This work also provides some special cases of the classical subdivision schemes.

Preliminary Results
Let  be the set of integers and { } i a a i = ∈  a set of constants.The general form of univariaten-ary subdivision scheme which maps a polygon where the set { } i a a i = ∈  of coefficients is called mask of the subdivision scheme.A necessary condition for the uniform convergence of the subdivision scheme is be the space of all polynomials of degree 2 is fundamental Lagrange polynomial corresponding to the nodes { } (2.4) Using all the above mentioned identities Ko [6] presented the general formula for the mask of ( ) 2 4 N + -po- int binary interpolating symmetric subdivision schemes.After that Najma [7] generalized the result for ( ) 2 4 N + -point n-ary interpolating symmetric subdivision scheme and gave the following formula for the mask of n-ary interpolating schemes.
Here, 0 N ≥ with the symmetry condition is for the mask of ( ) (2.5).Following the procedure of binary case, we have derived the following form of ( ) , where 0,1, 2.

Construction of the Schemes for Open Polygon
When dealing with open initial polygon it is not possible to refine the first and last edges by rules (2.10) for interpolating subdivision schemes.However the extension of this strategy to deal with open polygon requires a well-define neighborhood of end points.Since the first and last edges can be treated analogously, it will be sufficient to derive the rules only for one side of the polygon.To this aim define the auxiliary point 0 0 0 0 as extrapolatory rule in the initial polygon 0 f .Then the nonrefined open polygon can be refined by the rules defined below.The formula described in (2.10) for interpolating scheme is not helpful to refine first and last edges of open polygon.Then to refine the open polygon by ( ) 2 4 N + -point interpolating scheme using auxiliary points is defined as following ( ) ( ) where ( ) , where, the weights satisfies the same condition (2.11).
Example: If an open polygon is refined by using the 6-point ternary interpolating subdivision scheme using (2.10), then two auxiliary points 0 has to be defined in the coarsest polygon 0 f .The first two edges 0 1 can be refined by the rules that can be calculated directly by (2.12).Substituting 3, in (2.12), ( ) ( )

Tensor Product of (2N + 4)-Point Interpolating Subdivision Scheme
Given a set of control points , , , , where k is a non-negative integer indicates the subdivision level.n-ary subdivision surface is tensor product of n-ary subdivision curve defined by where, ,r Given initial values 0 , , , , , , is related, in a natural way, with a diadic mesh points , , The process then defines a scheme whereby , (where α and β are not zero at the same time).Labeling of old and new points is shown in Figure 1, which illustrates subdivision schemes (3.1). 1, ,? , , , , are Kroneker delta symbols defined as, , 1, 0, Here, some important results for the formulation of required form of tensor product scheme can be verified using (2.3).That is for each 1, , µ The mask of a subdivision scheme shows the contribution of a single original vertex to each new, subdivided vertex.To find the mask of a scheme, we need to find all ways to get from the origin to each point in the grid.For the tensor product scheme, this is simply the tensor product of the univariate case.
Lemma 3.1.[8] Given initial control polygon 0 ., i j i j = , , i j Z ∈ ，let the values ., k i j p 1 k ≥ be defined recursively by subdivision process (3.1) together with (3.2) then the scheme derived by tensor product naturally get four-sided support region.
It can be loosely say that the support is the tensor product of the supports of the two regions, just as one can loosely say that Doo-Sabin is the generalization of the tensor product of two Chaikin constructions.
Lemma 3.2.[9] Given initial control polygon 0 , , , p k ≥ be defined recursively by subdivision process (3.1) together with (3.2), then if a scheme is derived from a tensor product, then the level of continuity can be determined between pieces by reference to the underlying basis functions, i.e. all the tensor product schemes have the same continuity as their counterparts.
The general formula which generates the mask { } presented by [7] is ,0 2 ,0 2 where Here, n, m stands for n-ary, m-ary subdivision schemes respectively (i.e n, m = 2(binary), 3(ternary), 4(quaternary)···), , 0 where , where , i j b b are the mask of univariate subdivision schemes, then ( ) , ni s mj s ni s mj s a a a The tensor product of ( ) 2 4 N + -point interpolating subdivision scheme is presented as, ( ) , where, ( ) γ σ ≥ and , 2 n m ≥ and symmetry conditions are, . Example: Consider the tensor product of the 4-point DD interpolating subdivision scheme, while DD scheme can be calculated using the result (2.5) mentioned in Section 2. The Laurent polynomial of the scheme is given as ( ) ( ) , a z a z z a z = then, we can obtain the Laurent polynomial of the 4-point tensor product binary interpolating scheme a S .So that the suggested 4-point tensor product binary interpulating scheme is ( ) Using the result obtained above for the tensor product of interpolatory scheme (3.8), tensor product of 4-point DD scheme can be calculated directly.Since the DD scheme has 1  C continuity, then by lemma (3.2) its tensor product has the same continuity.Substituting , Example: A simple interpolatory subdivision scheme for quadrilateral nets with arbitrary topology is presented by L. Kobbelt [2] which generates 1  C surfaces in the limit.In the first step they present the refinement rules derived by the modification of the well-known Dyn et al. [10] 4-point interpolatory subdivision scheme for curve design.The natural way to define refinement operators for quadrilateral nets is therefore to modify a tensor product scheme such that special rules for the vicinity of non regular vertices are found.
The modified form of Dyn scheme can be evaluated by setting the value of n = 2, N =  = 0, 4 0 a = and They used the simple tensor product as the basis for the modification of refinement rules of irregular quadrilateral nets.Since it is interpolating scheme, so 1 2 ,2 , and the edge points

(
) ( ) ( ) ( ) Finally, the face point  γ σ = , the symmetry conditions are then written as, ,  Example: Using the results for the interpolating curve subdivision schemes (2.10) the 4-point interpolatory scheme [3] is obtained.Further here the tensor product of the scheme is evaluated by using the result (3.8).Put ,        (

Numerical Examples
Here, the performances of some of the schemes which are deduced from the proposed formula shown.
comes from the generalized formula

2 4 N
+ -point ternary interpolating subdivision schemes are presented in the following

Figure 1 .
Figure 1.Solid lines show one face of coarse polygons whereas dotted lines are refined polygons.(a)-(c) can be obtained by subdividing one face into four, nine and sixteen new faces by using (3.1) for n = 2,3,4 respectively.
is substituted in (3.8) and (3.9) our requirement is fulfilled, that is the rules (3.10) are obtained.
the tensor product the above rules can be directly obtained by substituting then the formula (3.8) acquires form, ( ) 17 and 4 0 a = for both n and m.After substituting the weights in (3.17) we get the same rules (3.14) and (3.15).
a = for both n and m.After calculating the mask from (3.4) and (3.5) and substituting all the results in equation (3.18) following 4-point ternary interpolating tensor product scheme is obtained

Fig- ure 2
shows the tensor product of 4-point ternary interpolating scheme (3.19), and Figure3gives the performance of the proposed 4-point binary scheme (3.10).