Modified Catmull-Clark Methods for Modelling, Reparameterization and Grid Generation
|
|
- Joy Moody
- 5 years ago
- Views:
Transcription
1 Modified Catmull-Clark Methods for Modelling, Reparameterization and Grid Generation Karl-Heinz Brakhage RWTH Aachen, 55 Aachen, Deutschland, Abstract In this paper we demonstrate that subdivision schemes can be used for modelling and grid generation in aeronautical projects such as that in the collaborative research center SFB at the RWTH Aachen: F low Modulation and Fluid-Structure Interaction at Airplane Wings. To make the schemes usable for this purposes their rules have to be modified in such a way that the limit surfaces have the desired properties. Furthermore fast algorithms for evaluation and approximation of existing surfaces have to be developed. Introduction We start with a short survey on the history of subdivision. At least fifty years ago the first papers on subdivision were published. In 97 de Rham s first articles on corner cutting appeared []. He suggested to do corner cutting with ratio : : compare Figure ). From the start polygon with vertices P, P,..., P n all the vertices P i are split the following way the indices of the old vertices P i have to be computed modulo n) Qn = P P + P n )/ Q = P + P )/ P P i.. Q = P + P )/ Q = P + P )/ Qi = P i + P i )/ Q i = P i + P i+ )/ ) Since the type of the refinement rule for this task is a vertex split we call it a vertex split scheme. The limit curve of this process is only GC- this is the same as C-, the curve is continuous). Later de Rham generalized the scheme by introducing a weight ω and uses the ratios ω : ω : ω. Furthermore he showed that the limit curves of this corner cutting are GC- this is C- continuously differentiable) after reparameterization) if and only if ω = /. The relevance of such algorithms to curve generation became clear with the paper of the graphics artist G. Chaikin [] in 97. At the 97 CAGD conference at the University of Utah, he invented an algorithm for high speed curve generation which coincides with the de Rham method for ω = /. Already at that conference R. Riesenfeld and M. Sabin opined that Chaikins algorithm is an iterative way to generate uniform quadratic B-spline curves. An analysis of the method had the result that the limit curve is a uniform B-spline curve whose control points are the vertices of the starting polygon. Further investigations showed that uniform B-spline curves of any degree would have such a subdivision construction. This is shown for the cubic P Q Q R R R R n- Q n- R R n- R R 5 P n- Q n- Q 9 Q P Q Q 7 Q 5 P P Q Q Q P P n- P Q n- R n- R n- Q n- Figure : Cubic edge splitting scheme two steps). In every step all edges are split into two parts. The new points on the edges are constructed by ratio : of the adjacent vertices. The old vertices are recomputed with the ratio : :. Q R R Q P R R R R 5 Q Q Figure : Corner cutting of de Rham 97: The outer polygon P, P,..., P n ) is subdivided two times with the edge ratio ::. The result is the inner curve polygon R, R,..., R n ). We use the sequence P, Q, R,... instead of P, P, P,... to reduce the number of indices. Afterwards the outer parts corners) are cut away. The resulting polygon is the input for the next same step. P Q Q 5 P case in Figure. In this method the edges are split into two pieces and hence it is called edge split scheme. We have new vertices on the edges and recomputed old ones. The formulas for this scheme are Q = P n + P + P )/ Q = P + P )/. Q i = P i + P i + P i+ )/ Q i+ = P i + P i+ )/ ) The first subdivision schemes for surfaces were introduced in 97. Doo and Sabin [] established a quadratic
2 subdivision scheme which is a direct extension of the vertex split method of de Rham / Chaikin. Furthermore they analyzed the behavior of such schemes. This was done by writing the subdivision process in form of matrix multiplications and studying the eigenvalues and eigenvectors of the subdivision matrices. face point edge point WN)= NN-) Figure : Subdivision rules for the Catmull-Clark Method valence N = 5). The new) face points are the barycenters of the vertices of each face. This is denoted by averaging with the weights for every face vertex. Next the edge points are computed using the new face points according to the given weights. In the last vertex step the old) vertices are recomputed. Again the weights are given in the figure. Notice that in every case we have to divide by the sum of the weights. Catmull and Clark [] published their descriptions of quadratic and cubic subdivision surfaces in 97 too. Both papers appeared in the same journal on consequent pages.) The latter one regards to the curves generated according to ). Contrary to tensor-product splines their scheme can be applied to meshes that are not regular rectangular grids see Figure ). The number of edges that meet at a vertex is called the valence of that point. We denote it with N. After the vertex steps face vertex, edge vertex and vertex vertex) illustrated in Figure the new edges are build by splitting the old ones edge edge) and connecting the new face points with the new edge points face edge). The new faces are the faces inside the old ones face face). The right plot of Figure can be thought of as an example for a pentagon. The centroid is the new face vertex, the vertices with weight are the new edge vertices and the the vertices with weight are the recomputed vertices. In a regular rectangular grid the valence of all interior vertices is N =. Subdivision schemes that have the same coefficients/rules in every level are called stationary. In our modelling and grid generation concepts see [5]) we are interested in ending up with smooth untrimmed B-spline patches. Thus for our purposes the Catmull- Clark scheme is the method of choice. We summarize its crucial properties: The surfaces can be of arbitrary genus since the subdivision rules can be carried out on a mesh of arbitrary topological type. After one subdivision step all faces are quadrilaterals. Except at extraordinary vertices vertices of valence N ) the limiting surface converges to uniform bi-cubic B-Spline patches. Hence the surface is C except at extraordinary vertices. The number of extraordinary vertices is fixed after the first subdivision step. After two subdivision steps all faces have at most one extraordinary vertex. Near an extraordinary vertex the surface can be shown to have a well defined tangent plane at the limit point, but the curvature there is generally not well defined. The subdivision rules can be modified in such a way that they generate infinitely sharp creases as well as semi-sharp creases, i.e. creases whose sharpness can vary from zero meaning smooth) to infinite. See de Rose et. all [] 99) for more details. Because at least after one subdivision step there are only quadrilaterals we show the rules for this case in Figure. The correspondence to cubic B-splines can easily be seen. Notice that for the regular case A) =. AN) Figure : Subdivision rules for the Catmull-Clark-Method for quadrilaterals only. The weight for the recomputed extraordinary vertex of valence N is AN) = N 7 N. Again we have to divide by the sum of the weights in every case. Analysis of subdivision schemes The ideas of Doo and Sabin [] using eigenanalysis of the subdivision matrix were advanced by Ball and Storry from 9 to 9 in [7], [] and [9]) and later by Stam [] in 99. Stam additionaly gave an algorithm for evaluating the Catmull-Clark-Scheme and its derivatives) at arbitrary points. He used a choice of ordering for the control vertices that the main part of the subdivision matrix has a cyclical structure. Hence the discrete Fourier transform can be used to compute its eigenstructure. For our practical examples we can restrict the valence N of the extraordinary vertices to N, 5, }. For these cases we use a more intuitive numbering and pre-compute the eigenstructure with numerical methods. Our aim is Table : Dimension of subdivision matrices at extra ordinary points valence rows columns regular N N N + + 7
3 to solve interpolation and approximation problems for this methods. Thus we do not only need values at certain points but also need the coefficients of the involved mesh points. For this reason we additionally build up tables for the interpolation and approximation schemes for valence N,, 5, }. More details can be found in Appendix A CD-version of this paper only). For surfaces and even for valences N the subdivision matrices are of regarding analysis high dimension see Table ). Therefore we explain the ideas for surfaces and give the analysis for curves. The main difference between curves and surfaces is the existence of extraordinary points. The polygon for the curve subdivision can be evaluated everywhere at each step by well known B- spline algorithms whereas the surface mesh can not. The control vertex structure near an extraordinary point is not a simple rectangular grid, thus all faces that contain extraordinary vertices cannot be evaluated as uniform B-splines. Since after one subdivision step all faces are quadrilaterals and after two steps all faces have at most one extraordinary vertex, we assume for our mesh that each face is a quadrilateral and contains at most one extraordinary vertex. Figure 5 shows that the region Figure 5: Behavior of the Catmull-Clark-Method near an extraordinary point of valence N =. The non regular region dark) is scaled down with every step. where the surface can not be evaluated with standard methods is reduced in every subdivision step. Since we can evaluate the surface at extraordinary vertices the remaining problem is to demonstrate only how to evaluate a patch corresponding to a face with just one extraordinary vertex, such as the dark region shown in Figure 5. Analogue to curves compare Figure ) we introduce parameter values and define a surface patch xu, v) over the unit square [, ] [, ] such that the point x, ) corresponds to the extraordinary vertex. We can evaluate the surface at those vertices x, )) as a linear combination of the circumfluent vertices the same ones we need for the re-computation step vertex vertex during subdivision, but with different weights). Additionally we can evaluate xu, ) for u [, ] and x, v) for v [, ] as regular B-spline part. The remaining problem is the evaluation xu, v) in the rest of the unit square. This problem is solved in the following way. First we do just enough subdivision steps such that u, v) corresponds to a regular part at that stage. We then pick the right vertices and do the evaluation as regular B-spline. Q P P Q C Q Q C Figure : Cubic edge splitting scheme one step and limit curve xt), t [, ]). The points P, P, P, P determine the curve between C = P + P + P )/ and C = P + P + P )/. We say that C = x) corresponds to P and C = x) to P After one subdivision step C = x.5) can be computed as C = Q + Q + Q )/ = P + P + P + P )/. We will give the analysis for curves now see Figure ). Imagining P as an extraordinary vertex we can not compute xt), t [, ] by standard B-Spline methods cause we have to cross P and for this thought experiment we have an analogon to the surface case. The subdivision matrix A corresponding to P is given by A = and Q Q Q C P Q P = A P P P. ) The eigenvalues of A are, / and /. For every of our subdivision schemes there is a unique largest eigenvalue, all eigenvalues are real and the matrices are diagonalizable equal algebraic and geometric multiplicities for all eigenvalues). One matrix that has as columns the eigenvectors of A is V = V and thus = ) V A V = diag,, } = Λ. 5) On the other hand this yields A = V Λ V. Furthermore we have A n = V Λ V ) n = V Λ n V and lim n Λ n = diag,, } =: Λ. Now we can conclude A := lim n An = V Λ V =. ) If we apply the subdivision for the curve case on P, P, P this yields C = x) = P + P + P 7)
4 and we have proven the consistency to uniform B-spline curves. Due to the fact that λ = > λ = / > λ = / we can already conclude from the given inverse in ) that x) = P + P + P x ) = α P + P ) x ) = α P + P P ) ) If we want to do the same computation for C = x.5) we have to enlarge the subdivision matrix A to A with ) A A = = A A. 9) The new rows directly follow from the usual B-spline knot-insertion rules for curves and surfaces. We make use of the block structure to compute the eigenvalues and eigenvectors of A. First it induces ) ) V Λ V = and Λ = ) V V Λ as structure for the equation A V = ΛV. Inserting ) and expanding the last equation yields ) ) A V V Λ = A V + A V A V V Λ V Λ and we can compute the remaining eigenvalues and the eigenvector block V with A. For surfaces A is a 7 7 matrix and independent of the valence N. Hence it can be computed in advance and used for any valence. The remaining part is the extension V of the eigenvectors in V. V is a 7 N + matrix and can be computed by solving the N + linear systems A V + A V = V Λ ) which in more detail is A Λ j) I 7 ) V j) = A V. In the last equation we have used Λ j) for the j th diagonal element of Λ, I 7 for 7 7 identity matrix and V j) for the j th column of V. In our thought experiment for curves see 9)) we have A = /) = Λ and V = ). Solving equation ) for V, which is ) + V = V, we get V = ) and finally V =. ) For the computation of the inverse of V we can make use of the block structure again. We omit the details and only state the result V =. ) The same computations as above lead to V A V = diag,,, } = Λ. ) and A n = V Λ V ) n = V Λ n V. 5) Repeating the above process we can compute x.5) etc. But if we want to compute values in.5, ) we have to compute Q and thus need another subdivision step given by an usual B-spline knot-insertion rule. Notice that we do not need further vertices as input. Our matrix for such a step is A Â = A A = A A. ) Let us summarize the results from above. The first steps of the subdivision processes in matrix form are Q = A P and ˆQ = Â P. 7) We have collected the vertices in P old level) and in Q and ˆQ new level). If we build the product of the mask, /, /, /) and A we get the coefficient vector c = /, /, /, /) for P, the vector of vertices compare Figure ). The same result can be obtained by evaluating the basis functions of 9) for s = /. Since A is a square matrix we can iterate with it. Thus if we want to compute values inside t [/ i+, / i ] we have to do i subdivision steps with A and than one with Â. This results in a general subdivision matrix S = Â Ai = Â V Λi V. ) After that we pick the points numbered from to see Figure ) to compute the desired value by normal B- spline evaluation. The same way we can compute coefficients for the the vertices of the starting polyhedron everywhere cause we know the basis functions for cubic uniform B-splines. On [, ] they are b s) = s), b s) = + s) + s) s, b s) = + s + s s), b s) = s. 9)
5 We only have to transform the value t [/ i+, / i ] to s [, ] by the formula st) = t / i+ ) i+. In the event of surfaces we make use of the tensor product structure. More details can be found in Appendix A CDversion of this paper only). We now describe how to build the matrices for the approximation of given surfaces with our Catmull-Clark meshes. We can compute the values L i of the limit surface and the corresponding coefficients c j of the involved vertices for arbitrary wanted points. L i can be written as L i = c T i Pi where we have collected the c j in the vector c i and the associated vertices in P i. For a given surface s we can project L i onto it: L i L s i. Know we have the desired equation of form c T i P i = L s i. ) Sampling enough of these collocation points we end up in a sparse linear system for interpolation or an overdetermined sparse linear) system for approximation. Normally we use the latter one and solve it with CGLS, a Conjugate Gradient method for linear Least Squares also called CGNR in []). To apply CGLS for solving the normal equations we only need effective methods to multiply the system matrix M and its transpose M T with vectors. It should be noted that our vectors are vectors of d-vectors. In ) we can see how the sparse multiplication with M is done for one row. From Table we know that for valence N we have at most entries in each row of M. This number is independent of the number of control vertices #V). For the multiplication with M T we have to construct coefficient vectors for each row of M to get an analogue form to ). Again the number of non-zero entries is independent of #V but it depends on the ratio number of sample points)/#v. Figure : Simplified wing tip: Start polyhedron for the Catmull-Clark-Method and polyhedron after three subdivision steps =N+7 = N+ =N 5 II 9 I IV III 9 =N+ 7 =N+ 7=N+ 5 = N+ Figure 7: Area of effect for the Catmull-Clark method valence N =. For the computation of the limit point corresponding to vertex we need the N + = 7 vertices numbered from to N =. For area I we need the additional 7 this number is independent of N) vertices from N + = 7 to N + 7 =. For the computation inside the areas II, III and IV we need another 9 vertices from N + = to N + =. For valence N > we have a vertex 7 which does not exist in the case N = where we have to identify it with vertex. Examples In this section we give some examples from our aeronautical project. All examples show the vertex meshes and Figure 9: Wing-fuselage constellation: Start polyhedron for the Catmull-Clark-Method and constellation after three subdivision steps not the limit surfaces. Furthermore we do not project the subdivision surface onto our given constellation because we want to demonstrate the nice smoothing properties of the subdivision surfaces on a larger scale. For real applications an approximation step is done at the state of the lower parts of the figures and for flow calculations an exact evaluation of the limit surface is done. Figure shows the reparameterization of a wing tip. Some details of Figure 9 can be seen, especially the high grid quality at the extraordinary points valence ) of the wing tip. Furthermore it can be seen that already a very rough approximation leads to a high quality grid. Figure 9 shows a part of a wing-fuselage configuration. We use a non stationary subdivision process with tagged edges, in other words: some edges are treated specially for a few steps. For the first steps for instance we treat
6 Aachen, University of Technology, Aachen, Germany. Figure : Approximation of an engine: Start polyhedron for the Catmull-Clark-Method and polyhedron after three subdivision steps the transition curve from the wing to the fuselage with the masks of ) as a curve. Only in the last step we use the surface masks. This leads to semi sharp part in that region. Tagged edges and their influence on the limit surface where analyzed in [], [] and []. We use a modified version of the latter on. Figure shows the approximation of an engine. The scalings of the different areas have been slightly modified to see some details, namely the sharp edges and the overall high grid quality. Conclusion We have presented the basic framework for modelling, interpolation, approximation and grid generation with Catmull-Clark methods in technical projects. The result is a watertight geometry of untrimmed surface patches. The isolines of them yield high quality block structured grids that are well suited for adaptive flow solvers. In this contents adaptation is simply function evaluation. Our approach is based on the precalculation of coefficient masks for standard or modified Catmull-Clark methods. These coefficients can directly be used to efficiently solve the sparse linear systems. Acknowledgments This work has been performed with funding by the Deutsche Forschungsgemeinschaft in the Collaborative Research Center SFB Flow Modulation and Fluid- Structure Interaction at Airplane Wings of the RWTH References [] G. de Rham. Un peu de mathématiques á propos d une courbe plane. Elemente der Mathematik, ):9, 97. [] G. Chaikin. An algorithm for high speed curve generation. Computer Graphics and Image Processing, : 9, 97. [] M.A. Sabin and D. Doo. Behaviour of recursive division surfaces near extraordinary points. CAD, ):5, 97. [] J. Clark and E. Catmull. Recursively generated B-spline surfaces on arbitrary topological meshes. CAD, ):5 55, 97. [5] K.-H. Brakhage and Ph. Lamby. Application of B-spline techniques to the modeling of airplane wings and numerical grid generation. submitted to CAGD, 7, preprint at: [] T. DeRose, M. Kass, and T. Truong. Subdivision surfaces in character animation. In M. Cohen, editor, Proceedings of the 5th annual conference on Computer graphics and interactive techniques, ACM SIGGRAPH, pages 5 9. Addison Wesley, 99. [7] A.A. Ball and D.J.T. Storry. Recursively generated B-spline surfaces. In Proceedings of CAD, pages 9, Brighton, Sussex, U.K., April [] D.J.T. Storry. B-spline surfaces over an irregular topology by recursive subdivision. Ph.d. thesis, Loughborough University, 9. [9] A.A. Ball and D.J.T. Storry. Conditions for tangent plane continuity over recursively generated B-spline surfaces. ACM Transactions on Graphics, 7):, 9. [] J. Stam. Exact evaluation of catmull-clark subdivision surfaces at arbitrary parameter values. In Proceedings of SIGGRAPH, pages 95, revised version on: [] Y. Saad. Iterative Methods for Sparse Linear Systems. SIAM, second edition,. [] H. Hoppe, T. DeRose, T. Duchamp, M. Halstead, H. Jin, J. McDonald, J. Schweitzer, and W. Stuetzle. Piecewise smooth surface reconsruction. In Computer Graphics Proceedings, Annual Conference Series, ACM SIGGRAPH, page 95. Addison Wesley, 99. [] H. Biermann, A. Levin, and D. Zorin. Piecewisesmooth subdivision surfaces with normal control. In Kurt Akeley, editor, Siggraph, Computer Graphics Proceedings, ACM SIGGRAPH, pages. Addison Wesley Longman,.
Grid Generation and Grid Conversion by Subdivision Schemes
Grid Generation and Grid Conversion by Subdivision Schemes Karl Heinz Brakhage Institute for Geometry and Applied Mathematics RWTH Aachen University D-55 Aachen brakhage@igpm.rwth-aachen.de Abstract In
More informationVolume Mesh Generation for Numerical Flow Simulations Using Catmull Clark and Surface Approximation Methods
Volume Mesh Generation for Numerical Flow Simulations Using Catmull Clark and Surface Approximation Methods Michael Rom 1,2 and Karl Heinz Brakhage 1 Bericht Nr. 334 Dezember 2011 Key words: volume mesh
More informationUNIVERSITY OF CALGARY. Subdivision Surfaces. Advanced Geometric Modeling Faramarz Samavati
Subdivision Surfaces Surfaces Having arbitrary Topologies Tensor Product Surfaces Non Tensor Surfaces We can t find u-curves and v-curves in general surfaces General Subdivision Coarse mesh Subdivision
More informationNon-Uniform Recursive Doo-Sabin Surfaces (NURDSes)
Non-Uniform Recursive Doo-Sabin Surfaces Zhangjin Huang 1 Guoping Wang 2 1 University of Science and Technology of China 2 Peking University, China SIAM Conference on Geometric and Physical Modeling Doo-Sabin
More informationUsing Semi-Regular 4 8 Meshes for Subdivision Surfaces
Using Semi-Regular 8 Meshes for Subdivision Surfaces Luiz Velho IMPA Instituto de Matemática Pura e Aplicada Abstract. Semi-regular 8 meshes are refinable triangulated quadrangulations. They provide a
More informationSubdivision Surfaces
Subdivision Surfaces 1 Geometric Modeling Sometimes need more than polygon meshes Smooth surfaces Traditional geometric modeling used NURBS Non uniform rational B-Spline Demo 2 Problems with NURBS A single
More informationu 0+u 2 new boundary vertex
Combined Subdivision Schemes for the design of surfaces satisfying boundary conditions Adi Levin School of Mathematical Sciences, Tel-Aviv University, Tel-Aviv 69978, Israel. Email:fadilev@math.tau.ac.ilg
More informationSubdivision overview
Subdivision overview CS4620 Lecture 16 2018 Steve Marschner 1 Introduction: corner cutting Piecewise linear curve too jagged for you? Lop off the corners! results in a curve with twice as many corners
More informationRecursive Subdivision Surfaces for Geometric Modeling
Recursive Subdivision Surfaces for Geometric Modeling Weiyin Ma City University of Hong Kong, Dept. of Manufacturing Engineering & Engineering Management Ahmad Nasri American University of Beirut, Dept.
More informationPolar Embedded Catmull-Clark Subdivision Surface
Polar Embedded Catmull-Clark Subdivision Surface Anonymous submission Abstract In this paper, a new subdivision scheme with Polar embedded Catmull-Clark mesh structure is presented. In this new subdivision
More informationCurves and Surfaces 2
Curves and Surfaces 2 Computer Graphics Lecture 17 Taku Komura Today More about Bezier and Bsplines de Casteljau s algorithm BSpline : General form de Boor s algorithm Knot insertion NURBS Subdivision
More informationSubdivision curves and surfaces. Brian Curless CSE 557 Fall 2015
Subdivision curves and surfaces Brian Curless CSE 557 Fall 2015 1 Reading Recommended: Stollnitz, DeRose, and Salesin. Wavelets for Computer Graphics: Theory and Applications, 1996, section 6.1-6.3, 10.2,
More informationSubdivision Curves and Surfaces
Subdivision Surfaces or How to Generate a Smooth Mesh?? Subdivision Curves and Surfaces Subdivision given polyline(2d)/mesh(3d) recursively modify & add vertices to achieve smooth curve/surface Each iteration
More informationNon-Uniform Recursive Doo-Sabin Surfaces
Non-Uniform Recursive Doo-Sabin Surfaces Zhangjin Huang a,b,c,, Guoping Wang d,e a School of Computer Science and Technology, University of Science and Technology of China, PR China b Key Laboratory of
More informationSubdivision Curves and Surfaces: An Introduction
Subdivision Curves and Surfaces: An Introduction Corner Cutting De Casteljau s and de Boor s algorithms all use corner-cutting procedures. Corner cutting can be local or non-local. A cut is local if it
More informationTaxonomy of interpolation constraints on recursive subdivision curves
1 Introduction Taxonomy of interpolation constraints on recursive subdivision curves Ahmad H. Nasri 1, Malcolm A. Sabin 2 1 Dept. of Math. & Computer Science, American University of Beirut, Beirut, PO
More informationSmooth Patching of Refined Triangulations
Smooth Patching of Refined Triangulations Jörg Peters July, 200 Abstract This paper presents a simple algorithm for associating a smooth, low degree polynomial surface with triangulations whose extraordinary
More informationInterpolatory 3-Subdivision
EUROGRAPHICS 2000 / M. Gross and F.R.A. Hopgood (Guest Editors) Volume 19 (2000), Number 3 Interpolatory 3-Subdivision U. Labsik G. Greiner Computer Graphics Group University of Erlangen-Nuremberg Am Weichselgarten
More informationExact Evaluation Of Catmull-Clark Subdivision Surfaces At Arbitrary Parameter Values
Exact Evaluation Of Catmull-Clark Subdivision Surfaces At Arbitrary Parameter Values Jos Stam Alias wavefront Inc Abstract In this paper we disprove the belief widespread within the computer graphics community
More informationCurve Corner Cutting
Subdivision ision Techniqueses Spring 2010 1 Curve Corner Cutting Take two points on different edges of a polygon and join them with a line segment. Then, use this line segment to replace all vertices
More informationLocal Modification of Subdivision Surfaces Based on Curved Mesh
Local Modification of Subdivision Surfaces Based on Curved Mesh Yoshimasa Tokuyama Tokyo Polytechnic University tokuyama@image.t-kougei.ac.jp Kouichi Konno Iwate University konno@cis.iwate-u.ac.jp Junji
More informationTaxonomy of interpolation. constraints on recursive subdivision. Ahmad H. Nasri 1, Malcolm A. Sabin 2. 1 Introduction
1 Introduction Taxonomy of interpolation constraints on recursive subdivision surfaces Ahmad H. Nasri 1, Malcolm A. Sabin 2 1 Department of Mathematics and Computer Science, American University of Beirut,
More informationSubdivision Surfaces. Homework 1: Questions on Homework? Last Time? Today. Tensor Product. What s an illegal edge collapse?
Homework 1: Questions/Comments? Subdivision Surfaces Questions on Homework? Last Time? What s an illegal edge collapse? Curves & Surfaces Continuity Definitions 2 3 C0, G1, C1, C 1 a b 4 Interpolation
More informationApproximating Catmull-Clark Subdivision Surfaces with Bicubic Patches
Approximating Catmull-Clark Subdivision Surfaces with Bicubic Patches Charles Loop Microsoft Research Scott Schaefer Texas A&M University April 24, 2007 Technical Report MSR-TR-2007-44 Microsoft Research
More informationSubdivision Surfaces. Homework 1: Questions/Comments?
Subdivision Surfaces Homework 1: Questions/Comments? 1 Questions on Homework? What s an illegal edge collapse? 1 2 3 a b 4 7 To be legal, the ring of vertex neighbors must be unique (have no duplicates)!
More informationNormals of subdivision surfaces and their control polyhedra
Computer Aided Geometric Design 24 (27 112 116 www.elsevier.com/locate/cagd Normals of subdivision surfaces and their control polyhedra I. Ginkel a,j.peters b,,g.umlauf a a University of Kaiserslautern,
More informationHoneycomb Subdivision
Honeycomb Subdivision Ergun Akleman and Vinod Srinivasan Visualization Sciences Program, Texas A&M University Abstract In this paper, we introduce a new subdivision scheme which we call honeycomb subdivision.
More informationSubdivision Surfaces. Homework 1: Last Time? Today. Bilinear Patch. Tensor Product. Spline Surfaces / Patches
Homework 1: Questions/Comments? Subdivision Surfaces Last Time? Curves & Surfaces Continuity Definitions Spline Surfaces / Patches Tensor Product Bilinear Patches Bezier Patches Trimming Curves C0, G1,
More informationSubdivision. Outline. Key Questions. Subdivision Surfaces. Advanced Computer Graphics (Spring 2013) Video: Geri s Game (outside link)
Advanced Computer Graphics (Spring 03) CS 83, Lecture 7: Subdivision Ravi Ramamoorthi http://inst.eecs.berkeley.edu/~cs83/sp3 Slides courtesy of Szymon Rusinkiewicz, James O Brien with material from Denis
More informationSubdivision Surfaces
Subdivision Surfaces 1 Geometric Modeling Sometimes need more than polygon meshes Smooth surfaces Traditional geometric modeling used NURBS Non uniform rational B-Spline Demo 2 Problems with NURBS A single
More informationA subdivision scheme for hexahedral meshes
A subdivision scheme for hexahedral meshes Chandrajit Bajaj Department of Computer Sciences, University of Texas Scott Schaefer Department of Computer Science, Rice University Joe Warren Department of
More informationJoe Warren, Scott Schaefer Rice University
Joe Warren, Scott Schaefer Rice University Polygons are a ubiquitous modeling primitive in computer graphics. Their popularity is such that special purpose graphics hardware designed to render polygons
More informationNormals of subdivision surfaces and their control polyhedra
Normals of subdivision surfaces and their control polyhedra I. Ginkel, a, J. Peters b, and G. Umlauf a, a University of Kaiserslautern, Germany b University of Florida, Gainesville, FL, USA Abstract For
More informationSubdivision Surfaces
Subdivision Surfaces CS 4620 Lecture 31 Cornell CS4620 Fall 2015 1 Administration A5 due on Friday Dreamworks visiting Thu/Fri Rest of class Surfaces, Animation, Rendering w/ prior instructor Steve Marschner
More informationExample: Loop Scheme. Example: Loop Scheme. What makes a good scheme? recursive application leads to a smooth surface.
Example: Loop Scheme What makes a good scheme? recursive application leads to a smooth surface 200, Denis Zorin Example: Loop Scheme Refinement rule 200, Denis Zorin Example: Loop Scheme Two geometric
More informationSubdivision based Interpolation with Shape Control
Subdivision based Interpolation with Shape Control Fengtao Fan University of Kentucky Deparment of Computer Science Lexington, KY 40506, USA ffan2@uky.edu Fuhua (Frank) Cheng University of Kentucky Deparment
More informationEvaluation of Loop Subdivision Surfaces
Evaluation of Loop Subdivision Surfaces Jos Stam Alias wavefront, Inc. 8 Third Ave, 8th Floor, Seattle, WA 980, U.S.A. jstam@aw.sgi.com Abstract This paper describes a technique to evaluate Loop subdivision
More informationG 2 Interpolation for Polar Surfaces
1 G 2 Interpolation for Polar Surfaces Jianzhong Wang 1, Fuhua Cheng 2,3 1 University of Kentucky, jwangf@uky.edu 2 University of Kentucky, cheng@cs.uky.edu 3 National Tsinhua University ABSTRACT In this
More informationCS354 Computer Graphics Surface Representation IV. Qixing Huang March 7th 2018
CS354 Computer Graphics Surface Representation IV Qixing Huang March 7th 2018 Today s Topic Subdivision surfaces Implicit surface representation Subdivision Surfaces Building complex models We can extend
More informationSubdivision surfaces. University of Texas at Austin CS384G - Computer Graphics Fall 2010 Don Fussell
Subdivision surfaces University of Texas at Austin CS384G - Computer Graphics Fall 2010 Don Fussell Reading Recommended: Stollnitz, DeRose, and Salesin. Wavelets for Computer Graphics: Theory and Applications,
More informationCS354 Computer Graphics Surface Representation III. Qixing Huang March 5th 2018
CS354 Computer Graphics Surface Representation III Qixing Huang March 5th 2018 Today s Topic Bspline curve operations (Brief) Knot Insertion/Deletion Subdivision (Focus) Subdivision curves Subdivision
More informationTernary Butterfly Subdivision
Ternary Butterfly Subdivision Ruotian Ling a,b Xiaonan Luo b Zhongxian Chen b,c a Department of Computer Science, The University of Hong Kong b Computer Application Institute, Sun Yat-sen University c
More informationGenerating Airplane Wings for Numerical Simulation and Manufacturing
Generating Airplane Wings for Numerical Simulation and Manufacturing Karl-Heinz Brakhage Philipp Lamby Institute of Geometry and Applied Mathematics RWTH Aachen, University of Technology D-52056 Aachen,
More informationRemoving Polar Rendering Artifacts in Subdivision Surfaces
This is an electronic version of an article published in Journal of Graphics, GPU, and Game Tools, Volume 14, Issue 2 pp. 61-76, DOI: 10.1080/2151237X.2009.10129278. The Journal of Graphics, GPU, and Game
More informationTechnical Report. Removing polar rendering artifacts in subdivision surfaces. Ursula H. Augsdörfer, Neil A. Dodgson, Malcolm A. Sabin.
Technical Report UCAM-CL-TR-689 ISSN 1476-2986 Number 689 Computer Laboratory Removing polar rendering artifacts in subdivision surfaces Ursula H. Augsdörfer, Neil A. Dodgson, Malcolm A. Sabin June 2007
More informationG 2 Bezier Crust on Quad Subdivision Surfaces
Pacific Graphics (2013) B. Levy, X. Tong, and K. Yin (Editors) Short Papers G 2 Bezier Crust on Quad Subdivision Surfaces paper 1348 Figure 1: Two examples of Bezier crust applied on Catmull-Clark subdivision
More informationUntil now we have worked with flat entities such as lines and flat polygons. Fit well with graphics hardware Mathematically simple
Curves and surfaces Escaping Flatland Until now we have worked with flat entities such as lines and flat polygons Fit well with graphics hardware Mathematically simple But the world is not composed of
More informationPiecewise Smooth Subdivision Surfaces with Normal Control
Piecewise Smooth Subdivision Surfaces with Normal Control Henning Biermann New York University Adi Levin Tel Aviv University Denis Zorin New York University Abstract In this paper we introduce improved
More informationTriangle Mesh Subdivision with Bounded Curvature and the Convex Hull Property
Triangle Mesh Subdivision with Bounded Curvature and the Convex Hull Property Charles Loop cloop@microsoft.com February 1, 2001 Technical Report MSR-TR-2001-24 The masks for Loop s triangle subdivision
More informationINF3320 Computer Graphics and Discrete Geometry
INF3320 Computer Graphics and Discrete Geometry More smooth Curves and Surfaces Christopher Dyken, Michael Floater and Martin Reimers 10.11.2010 Page 1 More smooth Curves and Surfaces Akenine-Möller, Haines
More information09 - Designing Surfaces. CSCI-GA Computer Graphics - Fall 16 - Daniele Panozzo
9 - Designing Surfaces Triangular surfaces A surface can be discretized by a collection of points and triangles Each triangle is a subset of a plane Every point on the surface can be expressed as an affine
More informationSubdivision curves. University of Texas at Austin CS384G - Computer Graphics
Subdivision curves University of Texas at Austin CS384G - Computer Graphics Fall 2010 Don Fussell Reading Recommended: Stollnitz, DeRose, and Salesin. Wavelets for Computer Graphics: Theory and Applications,
More informationQUADRATIC UNIFORM B-SPLINE CURVE REFINEMENT
On-Line Geometric Modeling Notes QUADRATIC UNIFORM B-SPLINE CURVE REFINEMENT Kenneth I. Joy Visualization and Graphics Research Group Department of Computer Science University of California, Davis Overview
More informationSubdivision on Arbitrary Meshes: Algorithms and Theory
Subdivision on Arbitrary Meshes: Algorithms and Theory Denis Zorin New York University 719 Broadway, 12th floor, New York, USA E-mail: dzorin@mrl.nyu.edu Subdivision surfaces have become a standard geometric
More informationA subdivision scheme for hexahedral meshes
A subdivision scheme for hexahedral meshes Chandrajit Bajaj Department of Computer Sciences, University of Texas Scott Schaefer Department of Computer Science, Rice University Joe Warren Department of
More informationComputergrafik. Matthias Zwicker Universität Bern Herbst 2016
Computergrafik Matthias Zwicker Universität Bern Herbst 2016 Today Curves NURBS Surfaces Parametric surfaces Bilinear patch Bicubic Bézier patch Advanced surface modeling 2 Piecewise Bézier curves Each
More informationAdvanced Modeling 2. Katja Bühler, Andrej Varchola, Eduard Gröller. March 24, x(t) z(t)
Advanced Modeling 2 Katja Bühler, Andrej Varchola, Eduard Gröller March 24, 2014 1 Parametric Representations A parametric curve in E 3 is given by x(t) c : c(t) = y(t) ; t I = [a, b] R z(t) where x(t),
More informationAn interpolating 4-point C 2 ternary stationary subdivision scheme
Computer Aided Geometric Design 9 (2002) 8 www.elsevier.com/locate/comaid An interpolating 4-point C 2 ternary stationary subdivision scheme M.F Hassan a,, I.P. Ivrissimitzis a, N.A. Dodgson a,m.a.sabin
More informationFINAL REPORT. Tessellation, Fairing, Shape Design, and Trimming Techniques for Subdivision Surface based Modeling
FINAL REPORT Tessellation, Fairing, Shape Design, and Trimming Techniques for Subdivision Surface based Modeling (DMS-0422126) PI: Fuhua (Frank) Cheng Department of Computer Science College of Engineering
More informationSmooth Surfaces from 4-sided Facets
Smooth Surfaces from -sided Facets T. L. Ni, Y. Yeo, A. Myles, V. Goel and J. Peters Abstract We present a fast algorithm for converting quad meshes on the GPU to smooth surfaces. Meshes with 1,000 input
More informationAdvanced Graphics. Subdivision Surfaces. Alex Benton, University of Cambridge Supported in part by Google UK, Ltd
Advanced Graphics Subdivision Surfaces Alex Benton, University of Cambridge A.Benton@damtp.cam.ac.uk Supported in part by Google UK, Ltd NURBS patches aren t the greatest NURBS patches are nxm, forming
More informationComputergrafik. Matthias Zwicker. Herbst 2010
Computergrafik Matthias Zwicker Universität Bern Herbst 2010 Today Curves NURBS Surfaces Parametric surfaces Bilinear patch Bicubic Bézier patch Advanced surface modeling Piecewise Bézier curves Each segment
More information3D Modeling techniques
3D Modeling techniques 0. Reconstruction From real data (not covered) 1. Procedural modeling Automatic modeling of a self-similar objects or scenes 2. Interactive modeling Provide tools to computer artists
More informationSurfaces for CAGD. FSP Tutorial. FSP-Seminar, Graz, November
Surfaces for CAGD FSP Tutorial FSP-Seminar, Graz, November 2005 1 Tensor Product Surfaces Given: two curve schemes (Bézier curves or B splines): I: x(u) = m i=0 F i(u)b i, u [a, b], II: x(v) = n j=0 G
More informationTrimming for Subdivision Surfaces
Trimming for Subdivision Surfaces Nathan Litke a,1 Adi Levin b,2 Peter Schröder a,3 a Caltech, Pasadena, CA 91125, USA b Tel Aviv University, Tel Aviv 69978, Israel Abstract Trimming is an important primitive
More informationOn Smooth Bicubic Surfaces from Quad Meshes
On Smooth Bicubic Surfaces from Quad Meshes Jianhua Fan and Jörg Peters Dept CISE, University of Florida Abstract. Determining the least m such that one m m bi-cubic macropatch per quadrilateral offers
More informationSubdivision Depth Computation for Extra-Ordinary Catmull-Clark Subdivision Surface Patches
Subdivision Depth Computation for Extra-Ordinary Catmull-Clark Subdivision Surface Patches Fuhua Frank Cheng,GangChen, and Jun-Hai Yong University of Kentucky, Lexington, KY, USA Tsinghua University, Beijing,
More informationApproximating Catmull-Clark Subdivision Surfaces with Bicubic Patches
Approximating Catmull-Clark Subdivision Surfaces with Bicubic Patches CHARLES LOOP Microsoft Research and SCOTT SCHAEFER Texas A&M University We present a simple and computationally efficient algorithm
More informationA Continuous 3-D Medial Shape Model with Branching
A Continuous 3-D Medial Shape Model with Branching Timothy B. Terriberry Guido Gerig Outline Introduction The Generic 3-D Medial Axis Review of Subdivision Surfaces Boundary Reconstruction Edge Curves
More informationNear-Optimum Adaptive Tessellation of General Catmull-Clark Subdivision Surfaces
Near-Optimum Adaptive Tessellation of General Catmull-Clark Subdivision Surfaces Shuhua Lai and Fuhua (Frank) Cheng (University of Kentucky) Graphics & Geometric Modeling Lab, Department of Computer Science,
More informationDesign by Subdivision
Bridges 2010: Mathematics, Music, Art, Architecture, Culture Design by Subdivision Michael Hansmeyer Department for CAAD - Institute for Technology in Architecture Swiss Federal Institute of Technology
More informationComplexity Reduction of Catmull-Clark/Loop Subdivision Surfaces
EUROGRAPHICS 2001 / Jonathan C. Roberts Short Presentations Complexity Reduction of Catmull-Clark/Loop Subdivision Surfaces Eskil Steenberg The Interactive Institute, P.O. Box 24081, SE 104 50 Stockholm,
More informationSubdivision of Curves and Surfaces: An Overview
Subdivision of Curves and Surfaces: An Overview Ben Herbst, Karin M Hunter, Emile Rossouw Applied Mathematics, Department of Mathematical Sciences, University of Stellenbosch, Private Bag X1, Matieland,
More informationTo appear in Computer-Aided Design Revised June 18, J-splines
To appear in Computer-Aided Design Revised June 18, 2008 J-splines Jarek Rossignac School of Interactive Computing, College of Computing, Georgia Institute of Technology, Atlanta, GA http://www.gvu.gatech.edu/~jarek
More informationVolume Enclosed by Example Subdivision Surfaces
Volume Enclosed by Example Subdivision Surfaces by Jan Hakenberg - May 5th, this document is available at vixra.org and hakenberg.de Abstract Simple meshes such as the cube, tetrahedron, and tripod frequently
More informationComputer Graphics Curves and Surfaces. Matthias Teschner
Computer Graphics Curves and Surfaces Matthias Teschner Outline Introduction Polynomial curves Bézier curves Matrix notation Curve subdivision Differential curve properties Piecewise polynomial curves
More informationSmooth Subdivision of Tetrahedral Meshes
Eurographics Symposium on Geometry Processing (2004) R. Scopigno, D. Zorin, (Editors) Smooth Subdivision of Tetrahedral Meshes S. Schaefer J. Hakenberg J. Warren Rice University Abstract We describe a
More informationParametrization of General Catmull-Clark Subdivision Surfaces and its Applications
Parametrization of General Catmull-Clark Subdivision Surfaces and its pplications Shuhua Lai and Fuhua (Frank) Cheng Graphics & Geometric Modeling Lab, Department of Computer Science University of Kentucky,
More informationSubdivision Scheme Tuning Around Extraordinary Vertices
Subdivision Scheme Tuning Around Extraordinary Vertices Loïc Barthe Leif Kobbelt Computer Graphics Group, RWTH Aachen Ahornstrasse 55, 52074 Aachen, Germany Abstract In this paper we extend the standard
More informationSURFACE FAIRING FOR SHIP HULL DESIGN
SURFACE FAIRING FOR SHIP HULL DESIGN Xoán A. Leiceaga Eva Soto GED, Universidad de Vigo, Vigo, España leiceaga@uvigo.es Oscar E. Ruiz Carlos A. Vanegas Laboratorio CAD/CAM/CAE, Universidad EAFIT, Medellín,
More information3D Modeling Parametric Curves & Surfaces. Shandong University Spring 2013
3D Modeling Parametric Curves & Surfaces Shandong University Spring 2013 3D Object Representations Raw data Point cloud Range image Polygon soup Surfaces Mesh Subdivision Parametric Implicit Solids Voxels
More informationPostprocessing of Compressed 3D Graphic Data
Journal of Visual Communication and Image Representation 11, 80 92 (2000) doi:10.1006/jvci.1999.0430, available online at http://www.idealibrary.com on Postprocessing of Compressed 3D Graphic Data Ka Man
More information12.3 Subdivision Surfaces. What is subdivision based representation? Subdivision Surfaces
2.3 Subdivision Surfaces What is subdivision based representation? Subdivision Surfaces Multi-resolution (Scalability) One piece representation (arbitrary topology) What is so special? Numerical stability
More informationPhysically-Based Modeling and Animation. University of Missouri at Columbia
Overview of Geometric Modeling Overview 3D Shape Primitives: Points Vertices. Curves Lines, polylines, curves. Surfaces Triangle meshes, splines, subdivision surfaces, implicit surfaces, particles. Solids
More informationSurface Quality Assessment of Subdivision Surfaces on Programmable Graphics Hardware
Sur Quality Assessment of Subdivision Surs on Programmable Graphics Hardware Yusuke Yasui Takashi Kanai Keio University SFC Faculty of Environmental Information 53 Endo, Fujisawa, Kanagawa, 5-850, JAPAN.
More information3D Modeling Parametric Curves & Surfaces
3D Modeling Parametric Curves & Surfaces Shandong University Spring 2012 3D Object Representations Raw data Point cloud Range image Polygon soup Solids Voxels BSP tree CSG Sweep Surfaces Mesh Subdivision
More information1. Introduction. 2. Parametrization of General CCSSs. 3. One-Piece through Interpolation. 4. One-Piece through Boolean Operations
Subdivision Surface based One-Piece Representation Shuhua Lai Department of Computer Science, University of Kentucky Outline. Introduction. Parametrization of General CCSSs 3. One-Piece through Interpolation
More informationAn Efficient Data Structure for Representing Trilateral/Quadrilateral Subdivision Surfaces
BULGARIAN ACADEMY OF SCIENCES CYBERNETICS AND INFORMATION TECHNOLOGIES Volume 3, No 3 Sofia 203 Print ISSN: 3-9702; Online ISSN: 34-408 DOI: 0.2478/cait-203-0023 An Efficient Data Structure for Representing
More informationGeneralizing the C 4 Four-directional Box Spline to Surfaces of Arbitrary Topology Luiz Velho Abstract. In this paper we introduce a new scheme that g
Generalizing the C 4 Four-directional Box Spline to Surfaces of Arbitrary Topology Luiz Velho Abstract. In this paper we introduce a new scheme that generalizes the four-directional box spline of class
More informationDYNAMIC GEOMETRY AND WEBSITE SETUP BY AUTOMATIC OBJECT RECOGNITION OF FREE-HAND DRAWINGS AND SCANS
16 TH INTERNATIONAL CONFERENCE ON GEOMETRY AND GRAPHICS 2014 ISGG 4 8 AUGUST, 2014, INNSBRUCK, AUSTRIA DYNAMIC GEOMETRY AND WEBSITE SETUP BY AUTOMATIC OBJECT RECOGNITION OF FREE-HAND DRAWINGS AND SCANS
More informationGeometric modeling 1
Geometric Modeling 1 Look around the room. To make a 3D model of a room requires modeling every single object you can see. Leaving out smaller objects (clutter) makes the room seem sterile and unrealistic
More informationAn introduction to interpolation and splines
An introduction to interpolation and splines Kenneth H. Carpenter, EECE KSU November 22, 1999 revised November 20, 2001, April 24, 2002, April 14, 2004 1 Introduction Suppose one wishes to draw a curve
More informationFast Rendering of Subdivision Surfaces
Fast Rendering of Subdivision Surfaces Kari Pulli (Univ. of Washington, Seattle, WA) Mark Segal (SGI) Abstract Subdivision surfaces provide a curved surface representation that is useful in a number of
More informationSmooth Surface Reconstruction using Doo-Sabin Subdivision Surfaces
Smooth Surface Reconstruction using Doo-Sabin Subdivision Surfaces Fuhua (Frank) Cheng, Fengtao Fan, Conglin Huang, Jiaxi Wang Department of Computer Science, University of Kentucky, Lexington, KY 40506,
More informationNonmanifold Subdivision
Nonmanifold Subdivision Lexing Ying, Denis Zorin NYU Abstract Commonly used subdivision schemes require manifold control meshes and produce surfaces which are manifolds, that is, each point on the surface
More informationMA 323 Geometric Modelling Course Notes: Day 36 Subdivision Surfaces
MA 323 Geometric Modelling Course Notes: Day 36 Subdivision Surfaces David L. Finn Today, we continue our discussion of subdivision surfaces, by first looking in more detail at the midpoint method and
More informationSpace Filling Curves and Hierarchical Basis. Klaus Speer
Space Filling Curves and Hierarchical Basis Klaus Speer Abstract Real world phenomena can be best described using differential equations. After linearisation we have to deal with huge linear systems of
More informationShape optimization of smooth surfaces with arbitrary topology
International conference on Innovative Methods in Product Design June 15 th 17 th, 2011, Venice, Italy Shape optimization of smooth surfaces with arbitrary topology Przemysław Kiciak (a) (a) Institut Matematyki
More informationFrom curves to surfaces. Parametric surfaces and solid modeling. Extrusions. Surfaces of revolution. So far have discussed spline curves in 2D
From curves to surfaces Parametric surfaces and solid modeling CS 465 Lecture 12 2007 Doug James & Steve Marschner 1 So far have discussed spline curves in 2D it turns out that this already provides of
More information852 QIN Kaihuai, CHANG Zhengyi et al Vol7 are generated by the recursive subdivision Each subdivision step will bring new vertices, new edges, as well
Vol7 No6 J Comput Sci & Technol Nov 22 Physics-ased Loop Surface Modeling QIN Kaihuai ( ΠΞ), CHANG Zhengyi ( ffω), WANG Huawei (ΦΛΨ) and LI Denggao (± ) Department of Computer Science and Technology, Tsinghua
More information