Construction and Analysis of Binary Subdivision Schemes for Curves and Surfaces Originated from Chaikin Points

We present a new variant of Lane-Riesenfeld algorithm for curves and surfaces both. Our refining operator is the modification of Chaikin/Doo-Sabin subdivision operator, while each smoothing operator is the weighted average of the four/sixteen adjacent points. Our refining operator depends on two parameters (shape and smoothing parameters). So we get new families of univariate and bivariate approximating subdivision schemes with two parameters. The bivariate schemes are the nontensor product schemes for quadrilateral meshes. Moreover, we also present analysis of our families of schemes. Furthermore, our schemes give cubic polynomial reproduction for a specific value of the shape parameter. The nonuniform setting of our univariate and bivariate schemes gives better performance than that of the uniform schemes.


Introduction and Related Work
Subdivision is a process of generating curves/surfaces by iteratively refining a set of control points according to some specific refinement rules.The set of these specific rules are called the subdivision scheme.For stationary, linear, and uniform subdivision schemes, these refinement rules are same at each refinement level.So the set of new control points at next refinement level can be generated by computing the affine combination of the set of control points at previous refinement level.A general compact form of linear, uniform, and stationary binary -variate ( = 1,2) subdivision scheme   which maps a polygon  −1 = { −1  ,  ∈ Z  } to a refined polygon   = {   ,  ∈ Z  } is defined as The symbol of above subdivision scheme is given by the Laurent polynomial where  = {  ,  ∈ Z  } is called the mask of subdivision scheme.The detailed information about refinement rules and Laurent polynomial can be found in [1][2][3].Surface modeling via subdivision is very important topic in computer graphics and computer aided geometric design and it has been studied by several authors; see surveys [4][5][6][7][8][9][10] and references therein.
The construction of univariate subdivision schemes by different variants of Lane-Riesenfeld algorithm has been studied by few authors in [8,[11][12][13] and construction of bivariate subdivision schemes by a variant of Lane-Riesenfeld algorithm has been studied by Romani [8].In fact, the Lane-Riesenfeld algorithm is the combination of two operators.One operator refines the initial points and other operator smooths these refining points for a fixed number of times.Therefore it is also named as Refining-Smooth algorithm (RSalgorithm).Moreover, more smoothing steps provide limit curves/surfaces of wider support as well as of higher smoothness.
Lane and Riesenfeld [7] proposed an algorithm for generating (+1)-degree B-spline curve, which is the simplest form of the RS-algorithm.In Lane-Riesenfeld algorithm, they combined the symbols of linear B-spline scheme and odd 2 International Journal of Analysis (a) Our refining stage stencil of linear B-spline scheme as refining and smoothing operators, respectively.Therefore the symbol of th family member of Lane-Riesenfeld schemes gives symbol of ( + 1)degree B-spline scheme.Their family gives   continuity and linear polynomial reproduction.
In literature, there are few variants which are applied on the Lane-Riesenfeld algorithm.For example, Cashman et al. [12] also proposed a family of univariate subdivision schemes by using the RS-algorithm which is based on Dubuc Deslauriers 4-point interpolatory scheme [14].They combined the symbols of 4-point interpolatory subdivision scheme and odd stencil of 4-point interpolatory subdivision scheme as refining and smoothing operators, respectively.Their univariate family gives cubic polynomial reproduction but by increasing smoothing stages continuity of their family may or may not be increased.
Ashraf et al. [11] proposed a family of univariate subdivision schemes by using similar technique to that used by Cashman et al. [12] on Dubuc Deslauriers 6-point interpolatory scheme [14].They also combined the symbol and odd stencil's symbol of Dubuc Deslauriers 6-point interpolatory scheme [14] as refining and smoothing operators, respectively.Their univariate family gives quintic polynomial reproduction but level of continuity of their family does not increase in general by increasing the smoothing stages.
Mustafa et al. [13] also proposed a family of univariate subdivision schemes by using the symbol of 4-point interpolatory subdivision scheme [15] as the refining operator and the symbol of even stencil of the 4-point approximating scheme [16] as the smoothing operator.Their univariate family gives cubic polynomial reproduction but level of continuity of their family also does not increase in general by increasing the smoothing stages.
Romani [8] proposed families of univariate and bivariate subdivision schemes by using RS-algorithm in which the refining operator is based on a perturbation of Chaikin's corner cutting subdivision scheme [17] and the smoothing operator takes the average of two adjacent vertices as in Lane-Riesenfeld algorithm.Their univariate family gives linear polynomial reproduction and continuity of their th family member is   .Moreover, if  = 1/16, their th family member becomes ( + 3)th member of Hormann and Sabin's [16] family and it gives cubic polynomial reproduction.
There are few variants of Lane-Riesenfeld algorithm but all those algorithms started with two binary schemes, by taking the symbol of first scheme as the refining operator and the odd or even symbol of second scheme as the smoothing operator.Our approach is different in the sense that we do not use symbol of any existing subdivision scheme as refining operator and also we choose appropriate smoothing operator by ourselves.
In this paper, we propose a new RS-algorithm which is also a variant of Lane-Riesenfeld algorithm that compared with the one proposed by Romani; we modify refining and smoothing operators both.There are two main differences between our RS-algorithm and Romani's RS-algorithm, which are as follows: (i) Both refining operators depend on Chaikin points, computed around initial points, but the difference between our refining operator and Romani's refining operator is that Romani draws two vectors at every initial control point and changes position of both Chaikin points, which have been computed around corresponding initial point, by using these two vectors, while we draw one vector at every Chaikin point and change the position of each Chaikin point by using each corresponding vector which we have drawn on it (see Figure 1).
(ii) Difference between both smoothing operators is that Romani's smoothing operator takes average of two adjacent points while our smoothing operator takes weighted average of the four adjacent points.
Our refining and smoothing operators are different from the operators proposed by Romani.The framework proposed in this paper provides families of schemes where all members are approximating while the first member of Romani's family is interpolatory and all the other members are approximating.Moreover, the flexibility of our shape parameter allows us to adjust the value of shape parameter locally according to the geometry of initial polygon.
The article is organized as follows.In Section 2, we present the framework for the construction of univariate and bivariate families of subdivision schemes and their graphical and mathematical comparisons with Romani's families.In Section 3, we present some basic properties of these families of subdivision schemes that include discussion on basic limit functions, analysis of continuity, polynomial generation, and reproduction of the schemes.In Section 4, applications of some members of the family and comparisons with Romani schemes are presented.Conclusions are given in Section 5.

Construction and Comparison
The Refining-Smooth subdivision scheme is the scheme which consists of two operators.One operator refines the initial points while the other operator smooths the refined points for the required number of times.If we denote refining operator by , smoothing operator by , and the subdivision operator by    , then where In fact, refining operator  maps a polygon

𝛼
,  ∈ Z  }, where  = 1, 2; that is, and smoothing operator  maps a polygon The symbols of refining and smoothing operators are with In RS-algorithm, combination of one refining operator and  smoothing operators is called the subdivision operator.Moreover, we denote proposed univariate and bivariate families of subdivision schemes by    and    , respectively.The construction procedure for our families of schemes is given below.by taking the linear combinations of two vectors, as shown in Figure 1(a): where in Figure 1(b) vectors are computed as with  ∈ N.
In our framework, we compute refining point by adding each Chaikin point to the corresponding vector; that is, International Journal of Analysis  where  ∈ R. While in Romani's framework refining points are computed as Now by substituting ( 9) in ( 11), we get the following subdivision scheme which will refine the course polygon  −1,0 into refined polygon  −1,1 : where Now the second step in our frame work is to smooth the refined polygon  −1,1 = { −1,1  } ∈Z by applying -times smoothing operator.That is, smoothing operator maps by using the following rule: where  = 1, 2, . . ., .
By applying one time refining operator and -times smoothing operator, we get the new family of subdivision schemes with symbol where ⌈ ⌉ is a ceiling function.
The first member of our proposed family is a primal 5point relaxed subdivision scheme with two refinement rules.The brief description of these refinement rules is given below.
(i) Edge splitting rule: for every edge in the coarse polygon, a new edge point is calculated by using the affine combination of 4 neighboring control points as shown in Figure 2(a).
(ii) Vertex updating rule: for every vertex in the coarse polygon, a new vertex point is calculated by using the affine combination of 5 neighboring control points as shown in Figure 2(b).
Remark 1.For  = 0, symbol (16) takes the form Now we shall construct family of bivariate nontensor product subdivision schemes for quadrilateral meshes by choosing Doo-Sabin points as origin.and  −1,1 2+1,2+1 that are computed from the initial points  −1 , ,  −1 +1, ,  −1 ,+1 , and  −1 +1,+1 as shown in Figure 3.In this figure, black solid circles show initial points at ( − 1)th subdivision level, blue cross symbols show Doo-Sabin points, and blue solid circles show refining points at ( − 1)th level; red, gray, green, and orange vectors are involved in the insertion of the refining points
Similarly to the curve case, we compute refining points as where  ∈ R. Now by substituting ( 18) in ( 19), we get the subdivision scheme which will refine the course mesh f k−1,0 into refined mesh f k−1,1 . where } ,∈Z by applying -time smoothing operator.That is, we map , by using the following smoothing operator: where  = 1, 2, . . ., .By applying one time refining operator  and -times smoothing operator , we get the family of bivariate primal/dual subdivision schemes for quadrilateral meshes with symbol Remark 2. The bivariate family for  = 0 corresponding to the symbol (22) reduces to the tensor product univariate family.This tensor product family of schemes can also be obtained from the symbol defined by ( 16).If we apply one time refining operator (20) for  = 1 and one time smoothing operator (21), we get the primal approximating subdivision scheme with four refinement rules.The schematic overview of these refinement rules is given in Figure 4.

Remark 3.
For  1 = 1 and  2 =  or  1 =  and  2 = 1, our family of bivariate schemes reduce to the family of univariate schemes with the symbol This family of schemes is different from our proposed univariate family of schemes.

Analysis of Univariate and Bivariate Schemes
This section deals with the analysis of proposed families of schemes.The analysis includes support width and smoothness/continuity via Laurent polynomial method.The capacity of generating and reproducing polynomials has also been discussed mathematically.
3.1.Support of Basic Limit Function.Support of basic limit function of a subdivision scheme is equal to the number of spans of the curve affected when one control point is moved from its initial place.
A convergent subdivision scheme    defines a basic limit function   =  ∞     with   being the initial data such that Support of basic limit function of our univariate family of schemes depends on parameter .
Basic limit functions for the first three family members of univariate family for three different values of shape parameters are presented in Figure 5, whereas basic limit functions for the first family member of bivariate family for three different values of shape parameter are presented in Figure 6.From these figures, we see that the support width for positive value is less than the support width for negative value of the parameter involved in the bivariate family.

Continuity Analysis.
Here we present the smoothness of proposed families of schemes by using the Laurent polynomial method.For this purpose, we have used the methods that are given ([3], Theorems 4.11, 4.13, 4.29, and 4.30).

Theorem 5. (i) The subdivision scheme 𝑆
When the initial control points are taken at the same intervals from some polynomial and if the new data points after applying subdivision scheme are located on the polynomial of same degree, then we say that the subdivision scheme generates polynomials up to that certain degree (degree of polynomial).But if the new data points after applying subdivision scheme are located on the same polynomial, then we say that the subdivision scheme reproduces polynomials up to that certain degree.See [16] for details.We will discuss the polynomial generation and reproduction of our families in the upcoming subsections.

Generation Degree.
Let Ω  and Ω   be the range of  for which the subdivision schemes    and    are convergent, respectively.
Here we compute the polynomial generation of our proposed family of schemes by using the conditions given in [1] (Conditions  1 and   ).Theorem 8.The subdivision scheme    generates polynomials (i) up to degree  for  ∈ Ω  , (ii) up to degree  + 2 for  = 0.
Proof.Proof is same as Theorem 9.
Theorem 12.The subdivision scheme   1 reproduces linear polynomials for  ∈ Ω  1 . Proof.Since then Now by combing these results with Theorem 9, the required result has been proved.

Proof.
Since then Also, for  = 23/48, we have Now by combing these results with Theorem 10, the required result has been proved.
Remark 14.If we take  = 0, the symbols ( 16) and ( 22) reduce to the symbol of "univariate family of subdivision schemes" and "bivariate family of subdivision schemes."These families produce  +2 limiting curves and surfaces, respectively.Also these families generate polynomials up to degree  + 2, while polynomial reproduction of these families is linear.

Applications and Comparisons
In this section, we have shown the capacity of parameter  to control shapes in different ways.The optical presentation of our proposed families of univariate and bivariate (nontensor product for quadrilateral meshes) subdivision schemes also has been shown.In Figure 7, we have shown limit curves that are obtained by our proposed family of univariate subdivision schemes at five values of shape parameter.In Figures 8 and 9, we have also shown limit surfaces.These surfaces are obtained by   1 (i.e., scheme with one time smoothing operator) and   2 (i.e., scheme with two times smoothing operator) on quadrilateral meshes, respectively.From these figures, we observe that shapes of limit curves/surfaces are affected by shape parameter and number of applied smoothing stages both.
Tables 1 and 2 show comparison between our univariate and bivariate families with Romani's [8] univariate and bivariate families, respectively.From these tables it is clear that our schemes give better smoothness than Romani's schemes for  = 0; however in this case the support size is larger than Romani's schemes and the reproduction degree is only one.
The geometrical construction of our refining operator allows us to extend our subdivision schemes to the nonuniform setting according to the geometry of the initial polygon/mesh by defining local parameter   .Sometimes choosing a variable shape parameter   for each edge of the initial polygon/mesh gives better result than choosing a constant shape parameter  for an entire initial polygon/mesh.Figures 10(a) and 10(b) demonstrate the visual effects of the proposed scheme   2 and the scheme proposed by Romani   2, , respectively, for  = 1/8.Similarly Figures 11(a) and 11(b) display the visual performance of the proposed scheme   2 and Romani's scheme   2, , respectively, for  = 1/8, while Figures 10(c) and 11(c) show the performance of nonuniform setting of the schemes   2 and   2 , respectively.Also Figures 10(d and   2,  , respectively.These figures show that the limit curves/surfaces generated by nonuniform schemes are better than the limit curves/surfaces obtained from uniform schemes.

Conclusion
In this paper, we have presented univariate and bivariate families of approximating subdivision schemes depending on two parameters.These families are called the Refining-Smooth families where refining operator is the operator in International Journal of Analysis    , ,  ∈ Z, in surface case, where  = 1, 2, . . ., .We also have analyzed some basic properties of our proposed families.Our families of schemes give cubic polynomial reproduction for specific values of shape parameter .Polynomial generation and continuity of our schemes are  for specific range of parameter , while for  = 0 they are  + 2. However for  = 0, our schemes reproduce only linear

Figure 1 :
Figure 1: The geometrical difference between our refining stage and Romani's refining stage.

2. 1 .
Univariate Family.Here we denote the Chaikin points by  −1,1 2 and  −1,1 2+1 that are computed around the point  −1,0  =  −1  as shown in Figure 1.In this figure, black solid circles show control points at (−1)th subdivision level, blue crosses show the Chaikin points, and blue solid circles show refining points.Red vectors are involved in the insertion of even point  −1,1 2 and green vectors are involved in the insertion of odd point  −1,1 2+1 .In this figure, we have shown the refinement of only two Chaikin points which have been inserted around point  −1  .In our framework, we compute vectors  −1,1 2 and  −1,1 2+1

Figure 2 :
Figure 2: Schematic overview of the refinement rules when we apply one time refining operator for  = 1 and one time smoothing operator.Solid circles show the points at ( − 1)th subdivision level while circles show the points at th subdivision level after applying one refining and one smoothing stage.

Figure 3 :
Figure 3: The geometrical construction of the refining stage for quadrilateral meshes.

Figure 4 :
Figure 4: Schematic overview of the refinement rules for  = 1.Solid circles show points at ( − 1)th subdivision level where solid squares show the new points at level .

Figure 6 :
Figure 6: Basic limit functions of subdivision scheme   1 for three values of shape parameter .

Figure 8 :
Figure 8: (b)-(f) are limit surfaces obtained by applying one time our refining and one time smoothing operators for different values of shape parameter  after 3 iterations.

Figure 9 :Theorem 13 .
Figure 9: (b)-(f) are limit surfaces obtained by applying one time our refining and two times our smoothing operators for different values of shape parameter  after 3 iterations.

Figure 10 :
Figure 10: Solid lines represent initial polygons while dash lines represent limit curves obtained by subdivision schemes proposed ((a) and (c)) by us and ((b) and (d)) by Romani.

Figure 11 :
Figure 11: Solid red lines represent initial mesh, while red and blue bullets show initial points.
Similarly for  = 23/216 and  = 2,    reduces to the subdivision scheme with mask

Table 1 :
Comparison between our univariate family    and Romani's univariate family   , via different properties.Let S, SS, RD, and GD be the integer smoothness, support size, reproduction degree, and generation degree, respectively.

Table 2 :
Comparison between our bivariate family    and Romani's bivariate family   , via different properties.Let S, SS, RD, and GD be the integer smoothness, support size, reproduction degree, and generation degree, respectively.polynomials.Nonuniform setting of the proposed family members gives better performance than that of uniform setting of the schemes.