n Introduction Bezier curves B spline curves and
n Introduction Bezier curves, B spline curves and subdivision curves are all based upon the input of a control polygon and the specification of an algorithmic method that contructs a curve from this sequence of points. Fundamental to these methods is the concept of a refinement.
n A refinement process is a scheme which defines a sequence of control polygons What is a Refinement Scheme? n where for any k>0 , each can be written as:
What is a Refinement Scheme? n For each fixed j and k the sequence is frequently called a mask. n It covers the cases where the number of control points in each successive polygon is allowed to increase (Chaikin´s Curves and Doo Sabin´s subdivision surfaces are examples of this), or must decrease (de Casteljau's algorithm for generating Bézier Curves is an example of this. (
A Matrix Method for Refinement n n The equation can be written in matrix form as
A Matrix Method for Refinement (cont. ) and overall Where Sk is the refinement matrix n This matrix as being sparse (i. e. most of the entries being zero) with non zero entries clustered along the diagonal.
SUBDIVISION CURVES n n A new set of curve generation methods is now becoming popular which utilize a control polygon, as in the Bézier or B spline case, but instead of using analytic methods to directly calculate points on the curve, these methods successively refine the control polygons into an sequence of control polygons that, in the limit, converge to a curve. The curves are commonly called subdivision curves as the methods are based upon the binary subdivision of the uniform B spline curves.
CHAIKIN'S ALGORITHMS n In 1974, George Chaikin gave a lecture at the University of Utah in which he specified a novel procedure for generating curves from a limited number of points. This algorithm is interesting as it was one of the first corner cutting or refinement algorithms specified to generate a curve from a set of control points, or control polygon.
CHAIKIN'S ALGORITHMS n corner cutting or refinement algorithms: the algorithm generates a new control polygon by cutting the corners off the original one.
CHAIKIN'S ALGORITHMS n Given a control polygon , we refine this control polygon by generating a new sequence of control points: n where each new pair of points Qi, Ri are to be taken to be at a ratio of and between the endpoints of the line segment.
CHAIKIN'S ALGORITHMS n These 2 n new points can be considered a new control polygon a refinement of the original control polygon.
CHAIKIN'S ALGORITHMS n Example: How Chaikin's Algorithm works? n consider the following control polygon:
CHAIKIN'S ALGORITHMS n Chaikin's algorithm generates the points Qi and Ri and uses these points to refine the curve and obtain the control polygon shown in the figure below
CHAIKIN'S ALGORITHMS n These points are in turn utilized to generate a new refinement.
CHAIKIN'S ALGORITHMS n The following illustration shows the initial control polygon, the third control polygon in the sequence, and the final curve.
Subdivision Curves n Example: A Closed Curves consider the following figure
Subdivision Curves We can still apply Chaikin's method to this figure, obtaining
Subdivision Curves This new control polygon can then be utilize to obtain a second refinement as in the following figure
Subdivision Curves n The initial control polygon and the second refinement are shown in the following figure along with the resulting curve.
Subdivision Methods for Quadratic B Spline Curves
n Given a set of. Methods control points P 0 , P 1, P 2. the. B Spline quadratic Curves Subdivision for Quadratic uniform B spline curve : For and where
n We can perform a binaryfor subdivision of the curve, by Curves Subdivision Methods Quadratic B Spline applying one of two splitting matrices to the control polygon: n n When applied to the control polygon SL gives the first half of the curve, and SR gives the second half. several of the control points for the two subdivided components are the same. Thus, we can combine these matrices
n n apply it to a control polygon Subdivision Methods for Quadratic B Spline Curves i. e. at and along each of the lines of the control polygon. These are the same points as are developed in Chaikin's method.
n The General Refinement Procedure: Given a control polygon, we can Subdivision Methods forset Quadratic generate a refinement of this of points by B Spline constructing. Curves new points along each edge of the original polygon at a distance of and between the endpoints of the edge. The general idea behind subdivision curves is to assemble these points into a new control polygon which can then be used as input to another refinement operation, generating a new set of points and another control polygon and then continue this process until a refinement is reached that accurately represents the curve to a desired resolution.
Subdivision Methods for Cubic B Spline Curves
n Given a set of control points P 0 Cubic , P 1, P 2 and P 3. The Curves Cubic Subdivision Methods for B Spline uniform B Spline curve : For and where
n n We can perform. Methods a binary subdivision applying one of Subdivision for Cubicby. B Spline Curves the two splitting matrices. We can combine these matrices
n Apply it to a control polygon Subdivision Methods for Cubic B Spline Curves
n n The General Refinement Subdivision Methods. Procedure for Cubic B Spline Curves In the case of the quadratic curve we were able to state exactly a single procedure for the points of the refinement. In the case of the cubic curve, it is not so easy. However, if we examine the rows of matrix used in the refinement, we see that they have two distinct forms. This motivates us to classify the points of the refinement as vertex and edge points. This classification makes the description of the refinement process quite straightforward.
VERTEX POINTS AND EDGE POINTS
n n VERTEX POINTS AND EDGE POINTS Classifying the Refinement Points Suppose we are given a control polygon
VERTEX POINTS AND EDGE POINTS n n we use the refinement method: We note that three of the new control points appear to lie at the midpoints of the three respective line segments. These points will be classified as ``edge points''. The other points all lie close to on of the vertices of the original control polygon, and will be called ``vertex points''.
n VERTEX POINTS AND EDGE POINTS we denote the new control polygon generated by the binary subdivision method as n By combining and applying the splitting matrices that generate the binary subdivision of the curve we obtain
n VERTEX POINTS The edge points E 0 , E 1 AND and EEDGE by: 2 are calculated
n The V 0 and POINTS V 1 are calculated VERTEX AND by: EDGE POINTS
n An Example. POINTS of the Refinement Algorithm VERTEX AND EDGE POINTS
n Calculate the edge and AND vertex EDGE points producing a VERTEX POINTS refinement of the original control polygon
n Assemble these points into. AND a new control polygon and VERTEX POINTS EDGE POINTS utilize it as input again to the refinement process producing new edge and vertex points from this set of points.
Summary: VERTEX POINTS AND EDGE POINTS n In the case of cubic uniform B spline curves we can consider the refinement procedure to consist of the generation of edge points and vertex points from a given set of control points each of which can be generated by taking only the midpoints of line segments.
Subdivision Surfaces
Subdivision Surfaces n n Subdivision surfaces are based upon the binary subdivision of the uniform B spline curve/surface. In general, they are defined by a initial polygonal mesh, along with a subdivision (or refinement) operation which, given a polygonal mesh, will generate a new mesh that has a greater number of polygonal elements, and is hopefully ``closer'' to some resulting surface. By repetitively applying the subdivision procedure to the initial mesh, we generate a sequence of meshes that (hopefully) converges to a resulting surface.
Subdivision Methods for Quadratic B Spline Surfaces
Subdivision Methods for Quadratic B Spline Surfaces n n Consider the biquadratic uniform B spline surface P(u, v) defined by the array of control points M is the matrix:
n This patch can be subdivided into four. B Spline subpatches, which Subdivision Methods for Quadratic Surfaces can be generated from 16 unique sub control points. We focus on the subdivision scheme for only one of the four (the subpatch corresponding to ), as the others will follow by symmetry.
we consider the reparameterization of the surface by Subdivision Methods for Quadratic B Spline Surfaces and define this new surface as n
Subdivision Methods for Quadratic B Spline Surfaces n n Through this process, we have written the surface as The matrix S is typically called the “splitting matrix”
Subdivision Methods for Quadratic B Spline Surfaces
n Subdivision Methods for Quadratic B Spline Surfaces where the can be written as
Subdivision Methods for Quadratic B Spline Surfaces n These equations can be looked at in two ways: 1. Each of these points utilizes the four points on a certain face of the rectangular mesh, and calculates a new point by weighing the four points in the ratio of 9 3 3 1. Thus, this algorithm can be specified by using subdivision masks, which specify the ratios of the points on a face to generate the new points. In this case, the subdivision masks are as follows
Subdivision Methods for Quadratic B Spline Surfaces 2. Each of these equations is built from weighing the points on an edge in the ratio of 3 1. and then weighing the resulting points in the ratio 3 1. These are exactly the ratios of Chaikin's curve and so this method can be looked upon as a extension of Chaikin's curve to surfaces.
Doo Sabin Surfaces
n n Subdivision Methods for Quadratic B Spline Surfaces Given the subdivision masks generated for subdivision of biquadratic uniform B spline patches Doo and Sabin observed that the points are simply the average of four particular points taken in a polygon the vertex for which the new point is being defined, the two edge points (the midpoints of the edges that are adjacent to this vertex in the polygon), and the face point (average of the vertices of the polygon).
To Subdivision verify this gives the subdivision masks, B Spline let the 4 Surfaces vertex Methods for Quadratic of a face be named A (bottom right), B, C, D (upper left) The 2 edge points adjacent to A are respectively equal and n n D n The face point is equal to: n Finally the new point of the Doo Sabin subdivision is equal to: C B A
n n n Subdivision Methods for Quadratic B Spline Surfaces The Procedure for Meshes of Arbitrary Topology For each vertex of each face of the object, generate a new point as average of the vertex, the two edge points and the face point of the face.
Subdivision Methods for Quadratic B Spline Surfaces n Then for each face, connect the new points that have been generated for each vertex of the face
n n Subdivision Methods for Quadratic B Spline Surfaces For each vertex, connect the new points that have been generated for the faces that are adjacent to this vertex. And then for each edge, connect the new points that have been generated for the faces that are adjacent to this edge.
Subdivision Methods for Quadratic B Spline Surfaces n Three iterations of a Doo Sabin bicycle seat n After the first subdivision, all vertices have valence four which is a characteristic of the Doo Sabin method.
Subdivision Methods for Cubic B Spline Surfaces
n n Subdivision Methods for Cubic B Spline Surfaces Consider the bicubic uniform B spline surface P(u, v) defined by the array of control points M is the matrix:
n n Subdivision Methods for Cubic B Spline Surfaces The bicubic uniform B spline patch can be subdivided into four subpatches, which can be generated from 25 unique sub control points. We focus on the subdivision scheme for only one of the four (the subpatch corresponding to ), as the others will follow by symmetry. It should be noted that the nine ``interior'' control points are utilized by each of the four subpatch components of the subdivision.
we consider the reparameterization of the surface by Subdivision Methods for Cubic B Spline Surfaces and define this new surface as n
Subdivision Methods for Cubic B Spline Surfaces n And we obtain:
Subdivision Methods for Cubic B Spline Surfaces n And we obtain (cont. ):
Subdivision Methods for Cubic B Spline Surfaces These equations can be looked in terms of subdivision masks as follows n n The application of mask A generates a new point corresponding to each vertex. The masks B generate new points corresponding to each edge, and mask C generates a new point corresponding to each face. These new points can be classified into three categories vertex points, edge points and face points depending on each points relationship to the original control point mesh.
CATMULL CLARK SURFACES
CATMULL CLARK SURFACES n n Utilizing the subdivision for bicubic uniform B spline surfaces, Ed Catmull and Jim Clark, following the methodology of Doo and Sabin noted that the subdivision rules expressed for the cubic B spline surface not only work for arbitrary rectangular meshes, but can also be extended to meshes of an arbitrary topology. The Catmull Clarksubdivision process generalizes subdivision masks by calculating the face points, edge points and vertex points from an arbitrary mesh of control point to generate the subdivided control mesh.
CATMULL CLARK SURFACES n The points P’ 0, 0, P’ 0, 2, P’ 2, 0 and P’ 2, 2, which are shown in the figure below are called face points noted Fi, j, and are calculated as the average of the four points that bound the respective face.
CATMULL CLARK SURFACES We obtain: F 0, 0 = P’ 0, 0, F 0, 1 = P’ 0, 2, F 1, 0 = P’ 2, 0 and F 1, 1 = P’ 2, 2 n Note that the face points are simply the average of the 4 points which define the face. These face points Fi, jcan be calculted by: Fi, j=¼ (Pi, j+Pi+1, j+Pi, j+1+Pi+1, j+1).
CATMULL CLARK SURFACES n The points P’ 0, 1, P’ 0, 3, P’ 1, 0 P’ 1, 2, P’ 2, 1, P’ 2, 3, P’ 3, 0 and P’ 3, 2, which are called edge points noted Ei, j, are given as the average of 4 points the two points that define the original edge and the two new face points of the faces sharing the edge.
CATMULL CLARK SURFACES The example of the 2 edge points E 0, 1 and E 1, 0 n P’ 0, 1 can be rewrited as follow E 0, 1 = P’ 0, 1 =¼( (P 0, 0+P 1, 0+P 0, 1+P 1, 1)/4 + P 0, 1 + P 1, 1 + (P 0, 1+P 1, 1+P 0, 2+P 1, 2)/4)) which is also equal to ¼ (F 0, 0+F 0, 1+P 1, 1). So this kind of edge points Ei, j can be calculated as follow Ei, j =¼ (Fi, j 1+Fi, j+Pi+1, j) n P’ 1, 0 can be rewrited as follow E 1, 0 = P’ 1, 0=¼( (P 0, 0+P 1, 0+P 0, 1+P 1, 1)/4 + P 1, 0 + P 1, 1 + (P 1, 0+P 1, 1+P 2, 0+P 2, 1)/4)) which is also equal to ¼ (F 0, 0+F 1, 0+P 1, 1). So this kind of edge points Ei, j can be calculated as follow Ei, j =¼ (Fi 1, j +Fi, j+Pi, j+1) n
CATMULL CLARK SURFACES n Note that the edge points are the average of 4 points : the 2 face points of the 2 faces adjacent to the edge and the 2 points which define the edge. These edge points Ei, j can be calculted either as Ei, j=¼ (Fi, j 1+Fi, j+Pi+1, j) Or Ei, j=¼ (Fi 1, j+Fi, j+Pi, j+1)
CATMULL CLARK SURFACES n The remaining four points P’ 1, 1, P’ 1, 3, P’ 3, 1 and P’ 3, 3, as shown in the figure below are called vertex points noted Vi, j.
CATMULL CLARK SURFACES n n Let’s take the example of the vertex point V 1, 1 P’ 1, 1 can be rewrited as follow (1) V 1, 1 = P’ 1, 1 = 1/16 ( (F 0, 0+F 1, 0+F 0, 1+F 1, 1) + (P 1, 0+P 0, 1+P 2, 1+P 1, 2) + 8 P 1, 1) n Note that the vertex points are then the average of the 4 face points adjacent to the vertex point, with the 4 endpoints of the 4 edges adjacent to the vertex point and with 8 times the corresponding vertex from the original mesh.
CATMULL CLARK SURFACES n By denoting Q = 1/4 (F 0, 0+F 1, 0+F 0, 1+F 1, 1), R = 1/4 ( (P 0, 0+P 1, 1)/2 + (P 1, 0+P 1, 1)/2 + (P 2, 1+P 1, 1)/2 + (P 1, 2+P 1, 1)/2) and S = P 1, 1, we obtain V 1, 1 = P’ 1, 1 = 1/4 (Q + 1/4 (P 0, 0+P 1, 0+P 2, 1+P 1, 2) + P 1, 1 + 2 P 1, 1) = 1/4 (Q + 2 R + S) (2)
CATMULL CLARK SURFACES Note that the vertex points are also equal to 1/4 (Q + 2 R + S) with Q the average of the 4 face points adjacent to the vertex point, R the average of midpoints of the 4 edges incident on the vertex point and S is the corresponding vertex from the original mesh. n Depending on the reference, formula (1) or formula (2) is used to define vertex points.
CATMULL CLARK SURFACES Subdividing Meshes of Arbitrary Topology n Given a mesh of control points with an arbitrary topology, we can generalize the face point, edge point, vertex point specification from the uniform cubic B splinesurface calculations to obtain
CATMULL CLARK SURFACES n n For each face of the mesh, generate the new face points which are the average of all the original points defining the face (We note that faces may have 3, 4, 5, or more points now defining them). Generate the new edge points which are calculated as the average of the two points which define the original edge with the two new face points of the faces adjacent to the edge. Calculate the new vertex points which is equal to 1/n (Q + 2 R + (n 3)S), where n is the order of the vertex, and where Q is the average of the new face points of all faces adjacent to the original face point, R is the average of the midpoints of all original edges incident on the original vertex point, and S is the original vertex point. This calculation rule of the vertex point generates surfaces that exhibit tangent plane continuity at all extraordinary points. With the rule 1/4 (Q + 2 R + S) it was observed that in some arbitrary meshes, tangent plane continuity was not maintained at extraordinary points.
CATMULL CLARK SURFACES n Four steps in the Calmull Clark refinement process
CATMULL CLARK SURFACES
CATMULL CLARK SURFACES
Summarize Doo Sabinand Catmull Clark. Subdivisions Surfaces
Summarize Doo Sabinand Catmull Clark Subdivisions Surfaces n Doo Sabin(respectively Catmull Clark) surface is locally a bi quadratic(respectively bi cubic) B splinesurface, except at a finite number of control points (extraordinary points). Even at those points the Catmull Clarksurface satisfies tangent plane continuity. Both subdivision schemes are uniform and stationary approximating scheme mostly for quadrilateral net but have rules for subdividing n sidedpolygons. The Doo Sabin (respectively Catmull Clark) refinement rule is based on corner cutting – dual (respectively vertex insertion primal). After the first subdivision step of Doo Sabinscheme, all vertices have valence four. After a single subdivision step of Catmull Clark scheme, all the faces in the control net are quadrilaterals.
- Slides: 81