Iterative Computation of Moment Forms for Subdivision Surfaces
|
|
- Cornelia Gordon
- 5 years ago
- Views:
Transcription
1 Iterative Computation of Moment Forms for Subivision Surfaces Jan P. Hakenberg ETH Zürich Ulrich Reif TU Darmstat Figure : Solis boune by Catmull-Clark an Loop subivision surfaces with centrois inicate by ashe lines, an inertia inicate by the principal axes of the ellipsoi with equivalent inertia. The iterative erivation of multilinear forms etermines the moments with unpreceente accuracy. Abstract The erivation of multilinear forms use to compute the moments of sets boune by subivision surfaces requires solving a number of systems of linear equations. As the support of the subivision mask or the egree of the moment grows, the corresponing linear system becomes intractable to construct, let alone to solve by Gaussian elimination. In the paper, we argue that the power iteration an the geometric series are feasible methos to approximate the multilinear forms. The tensor iterations investigate in this work are shown to converge at favorable rates, achieve arbitrary numerical accuracy, an have a small memory footprint. In particular, our approach makes it possible to compute the volume, centroi, an inertia of spatial omains boune by Catmull-Clark an Loop subivision surfaces. Keywors: Stationary subivision, volume, centroi, inertia. Concepts: Theory of computation Computational geometry; Computing methoologies Shape analysis; Introuction Volume, centroi, an inertia are important shape characteristics in 3D moeling an animation. Accurate formulas for the etermination of the moments of solis boune by stationary subivision schemes have only recently become available. The approach is restricte to the computation of moments of low egree an to subivision schemes with weight masks of small support. The iterative methos presente in this publication exten the limits of what is possible. Using the new framework to compute moment forms, solis boune by subivision surfaces such as Catmull- Clark or Loop can be esigne with exact volume, centroi, an inertia up to machine precision. The moment formula for soli shapes in the context of subivision has the following applications: By translation of vertices in the control mesh, limit surfaces can be eforme subject to preservation of volume an centroi of the interior set. If a subivision surface is the contour of an animate entity with constant mass ensity, the formulas help to make the motion physically more accurate. In the absence of applie torques, the inertia (secon moment) is use to etermine motions that preserve angular momentum.. Relate Work [Peters an Nasri 997] approximate the moment of solis boune by [Doo an Sabin 978] an [Catmull an Clark 978] surfaces by increasing the level of refinement aroun non-regular vertices of a mesh. As facets ajacent to non-regular vertices ecrease in size, their moment contribution can be approximate using a planar surface patch. The approach requires regular submeshes to have a polynomial parametrization. Moments efine by the Butterfly scheme by [Dyn et al. 990] are not covere by their framework. [Schwal 999] achieves accurate results for solis boune by the Doo-Sabin algorithm. The moment contribution of a non-regular facet is expresse as a multilinear form in the coorinates of the vertices of the control mesh. In orer to compute the form coefficients, the volume contributions of all spline rings are taken into account. This is achieve by eigen-ecomposition of the subivision matrix. The contributions of infinite sequences of scale eigen-rings form a geometric series. The exact limit of the series is etermine algebraically. [Warren an Weimer 200] pp calculate the exact area form for the 4-point curve subivision scheme by [Dubuc 986]. The erivation is the first application of the homogeneous refinement equation () in the literature. [Hakenberg et al. 204] erive the exact moment forms for facets ajacent to non-regular vertices as well as sharp creases from the general refinement equation (7). The limitation of the approach is the size of the non-sparse linear systems that are solve to etermine the coefficients of the tensors. As the complexity of the subivision scheme, patch topology, an egree of moment grows, the corresponing linear system becomes intractable to store in memory an solve by Gaussian elimination. For instance, the inertia form for a qua facet ajacent to a vertex of valence 5 in a Catmull- Clark mesh has coefficients..2 Contributions We buil upon the work of [Hakenberg et al. 204] by introucing two iterative methos that approximate the coefficients of the moment forms up to arbitrary numerical precision. The first metho is the power iteration that etermines the form as an eigenvector associate with the largest eigenvalue of a linear system. The secon
2 iteration type is a geometric series. Eigen-analysis guarantees that the iterations converge at favorable rates for stanar subivision schemes so that builing the matrix of the linear system becomes unnecessary. Catmull-Clark surfaces are the e facto stanar for 3D moeling in the entertainment inustry. Our approach allows to calculate volume, centroi, an inertia of solis boune by these surfaces with unpreceente accuracy an minimal computational effort. 2 Preliminaries The moment E p,q,r(ω) of a boune set Ω R 3 for non-negative integers p, q, r is efine as the integral E p,q,r(ω) := x p y q z r xyz. () Ω The egree is efine as := p + q + r. Moments of low egree have interpretations as volume(ω) = E 0,0,0(Ω), centroi(ω) = (E,0,0(Ω), E0,,0(Ω), E0,0,(Ω)), volume(ω) inertia(ω) = (E 2,0,0(Ω), E 0,2,0(Ω), E 0,0,2(Ω), E,,0(Ω), E 0,,(Ω), E,0,(Ω)). The mass ensity is assume to be constant across the omain Ω. The ivergence theorem with vector fiel G : R 3 R 3 as G(x, y, z) = ( p+ xp+ y q z r, 0, 0), an iv G = x p y q z r reformulates () to an integral over the bounary Ω E p,q,r(ω) = x p+ y q z r n xa, (2) p + where n x enotes the x-component of the surface normal. 2. Multilinear Forms Let M be an orientable mesh of control points. Accoring to [Schwal 999] an [Hakenberg et al. 204], the moment of the set Ω boune by the subivision limit surface Ω = S (M) is compute as a sum over all facets of the mesh E p,q,r(ω) = E p,q,r(f). (3) Ω f M The contribution of the facet f to the global moment is E p,q,r(f) = p + M τ(f). x. {{.x. y..y. z. {{ {{.z.y.z. (4) p+-times q-times r-times The control points that etermine the surface patch associate to f are arrange into the coorinates vectors x, y, z. For instance, in a Catmull-Clark control mesh, each qua is a facet. The control points in the one ring of the qua etermine the surface patch. The tensor M τ(f) epens on the topology τ(f) of the facet, an the egree of the moment. [Peters an Nasri 997] establish that the coefficients of the tensor M τ(f) have the following integral expression ) i,...,i +,j,k = b i b i+ ( sb j tb k tb j sb k )st. (5) D For qua-base subivision schemes the omain is D = [0, ] 2, an for triangle base schemes D = {(s, t) R 2 0 s, t s+t. The basis function b i : D R is associate to control point i. For surfaces that consist entirely of B-spline, or Bézier-Bernstein patches, [Gonzalez-Ochoa et al. 998] integrate (5) symbolically in orer to yiel the exact tensor coefficients. For facets in subivision meshes that are ajacent to non-regular vertices, the basis functions typically o not have a close-form expression. 2.2 Recursive Formulation We summarize the erivation given in [Hakenberg et al. 204]. One roun of subivision partitions a facet f M into a number of smaller facets f h := S h (f). For stanar surface schemes h {, 2, 3, 4, see Figures 3 an 4, but this is not a requirement. The contribution E p,q,r(f) to the global moment of Ω is ientical to the total contribution of the subivie parts E p,q,r(f) = h E p,q,r(s h (f)). (6) Substituting x S h.x, y S h.y, an z S h.z in (4) to expan the rhs of (6) yiels the general tensor equation inepenent of control points M τ(f) = h M τ(f h) [S h ]. (7) The tensor transformation M[S] is shorthan for the basis transformation along each imension of tensor M with S, i.e. (M[S]) i,...,i m := The matrix S oes not have to be square. j,...,j m M j,...,j m S j i S jm i m. In orer to solve the coefficients of M τ(f), equation (7) is written as a linear system m = c + A.m. Matrix A encoes the linear relationship between the sought coefficients of M τ(f) that make up vector m. In case of homogeneous refinement, i.e. τ(f) = τ(f h ) for all h, it follows that vector c = 0. Then, the problem is reuce to ientifying m as an eigenvector to eigenvalue of matrix A m = A.m. (8) Once establishe, the eigenvector m has to be scale to the proper length in a calibration step using a configuration of control points x, y, z for which the moment E p,q,r(f) of the surface patch is known analytically. The calibration proceure is etaile in [Hakenberg et al. 204] an [Hakenberg an Reif 206b]. In case of non-homogeneous refinement, i.e. τ(f) τ(f h ) for at least one h, vector c aggregates the contribution of the subivie facets for which the moment forms M τ(f h) are alreay known. The resulting linear system is of the form (I A).m = c. (9) For non-regular vertices, the facet topology τ(f) typically appears once in the sum of the rhs of (7), see Figure 4. Then, the system matrix A = +3 S is the ( + 3)-fol Kronecker prouct of the square subivision matrix S. S maps the coorinates of the control points associate to facet f to those of the subivie non-regular facet f h with τ(f) = τ(f h ). The eigenvalues of A follow from the eigenvalues of S. For common subivision schemes, A has
3 eigenvalue with multiplicity, an all other eigenvalues λ with λ <. That means I A is rank eficient by. The nullspace of I A correspons to a symmetric tensor. But m is skew in the last two imensions as seen from (5), that means m follows uniquely from (9). 2.3 Tensor Symmetry When arranging the coefficients of M τ(f) into a vector m, reunant entries stemming from tensor symmetries shoul be omitte in orer to reuce the complexity of the linear system. Two symmetry relations are intrinsic to the inner prouct (5): Invariance uner permutation of the first + inices ) i,...,i +,j,k = an skew symmetry in the last two imensions ) i,...,i +,j,k = ) sort(i,...,i + ),j,k, ) i,...,i +,k,j. Permutations of control points uner which the image of the surface patch is invariant also contribute to tensor symmetries. For instance, rotation ρ of control point inices of a regular facet results in the same tensor coefficient ) i,...,i +3 = ) ρ(i ),...,ρ(i +3 ). Mirroring σ of control point inices of a facet along an axis of mirror symmetry flips the sign ) i,...,i +3 = ) σ(i ),...,σ(i +3 ). For egree = 0, one may aitionally assume that the trilinear form is alternating, i.e. ( ˆM τ(f) 0 ) i,j,k = 6 π S 3 sign(π) 0 ) π(i,j,k), (0) where S 3 is the permutation group of orer 3! = 6. [Hakenberg an Reif 206b] treat alternating volume forms for refinable basis functions b i of which subivision surfaces are a special case. Figure 2: Rotation (pink arrow) an mirror (purple axis) symmetries of facets f. Examples The surface patch associate to a regular Catmull-Clark qua facet f is etermine by 6 control points. The left iagram in Figure 2 illustrates the symmetries. The moment contribution of the facet is invariant up to sign uner rotation an mirroring of the control points. The number of unique coefficients in the centroi form M 4 is boune by The surface patch associate to a Loop triangular facet f with a non-regular vertex of valence τ(f) = 7 is etermine by 3 control points. The moment contribution of the facet is invariant up to sign uner mirroring of the control points across the axis through the non-regular vertex, see the right iagram in Figure 2. The number of unique coefficients in the inertia form M 7 2 is boune by Iteration of Forms Despite taking coefficient symmetries in M τ(f) into account, the linear systems (8) an (9) may still be too large to buil an etermine m. Moreover, the system matrix is typically non-sparse an has a large conition number. Memory constraints an loss of numerical accuracy prevent the application of Gaussian elimination. In this section, we argue that the power iteration an the geometric series are feasible methos to approximate the moment forms. In the subsequent treatment, the egree an facet topology τ(f) are fixe. We rop these two inices to simplify the notation. 3. Homogeneous Refinement Figure 3: Subivision of a regular Loop facet f into f h for h {, 2, 3, 4. S h are matrices of imension 2 2 that map the coorinates of the control points in the -ring of f to the control points of f h. In case of homogeneous refinement, all moment forms in (7) are ientical. The equation simplifies to M = h M[S h ]. () The system of linear equations is of the form m = A.m. An eigenvector m is approximate by the power iteration, if is the largest eigenvalue of matrix A, see [von Mises an Pollaczek-Geiringer 929]. The smaller the subominant eigenvalue of A, the faster the convergence of the iteration. We are not aware of a metho to euce the eigenvalues of A irectly from the collection of subivision matrices S h. For the subivision schemes that we investigate, the largest eigenvalue of matrix A is. For egree = 0, an assuming that M is an alternating trilinear form (0), the multiplicity of eigenvalue is. For {, 2, the multiplicity is 2. For Doo-Sabin, Loop, an Catmull-Clark subivision an {0,, 2 the eigenvalues of A were experimentally etermine to be of the form 2 2n for n {0,, 2,... For the Butterfly scheme with tension parameter w = /6 an = 0, the set of eigenvalues {, /4, 3/32, /6,... contains values outsie of this pattern. The construction of matrix A is not require to perform the power iteration. In tensor notation, the sequence is equivalent to M 0 := π(ranom) M k := M k / M k ( ) M k+ := π M k [S h ], h (2) which is more memory efficient. Coefficient symmetries as erive in Section 2.3 are enforce by a linear projection π in between iteration steps. After the iteration has converge, scaling of the form M = µ lim k M k with an appropriate factor µ R is performe analogous to the calibration of eigenvector m in (8).
4 Figure 4: Subivision of a Catmull-Clark facet f with a vertex of valence τ(f) := 5 into f h for h {, 2, 3, 4. S, S 2, S 3 are 6 8 matrices, S 4 is of imension Non-Homogeneous Refinement For a facet ajacent to a non-regular vertex, (7) simplifies to M = C + M[S]. (3) The tensor M = M τ(f) is unknown. Square matrix S maps the coorinates of the control points that support f to the control points of f h with τ(f) = τ(f h ). C is the constant tensor aggregate from the known terms. The eigenanalysis of the linear system (9) carrie out in Section 2.2 permits to write m as the geometric series m = (I A).c = lim k i=0 k A i.c. In tensor notation, the geometric series is compute as M 0 := C M k+ := C + M k [S] (4) with M = lim k M k. The subominant eigenvalue of A is the subominant eigenvalue λ of S. The smaller λ, the faster the convergence, see Figure 6. 4 Summary The moment forms M τ(f) for a specific subivision scheme have to be erive only once. Subsequently, the tensors apply universally to any close, orientable control mesh M using (3) an (4). The methos to erive the forms are selecte base on the type of basis functions, number of coefficients of M τ(f), an requirements on accuracy. In the following overview, methos 3 solve the homogeneous refinement equation (), whereas methos 4 5 solve the non-homogeneous refinement equation (3). Symbolic Integration The basis functions that parametrize regular surface patches of Doo-Sabin an Catmull-Clark surfaces are polynomials. [Peters an Nasri 997] compute the moment form for arbitrary egree by explicit integration of (5). 2 Nullspace an Calibration For the Loop, an Butterfly subivision schemes, the basis functions of the regular facet are not reaily available. (8) is rewritten to (I A).m = 0. The nullspace of matrix I A can be compute algebraically. For Loop s algorithm, the size of matrix A is 43, 874, 4032 square for = 0,, 2 respectively, see [Hakenberg et al. 204]. For the volume form of the Butterfly scheme, the size of matrix A is 508 square, see [Hakenberg an Reif 206b]. Subsequent calibration of the eigenvector m yiels the exact, rational tensor coefficients. 3 Power Iteration The regular patch of the Butterfly scheme epens on 29 control points. The centroi form is etermine by 2294 coefficients. We estimate that the system matrix A has subominant eigenvalue of λ = /4. We use the power iteration (2) to approximate the tensor. Machine precision in the coefficients is surpasse after few iterations, see Figure Figure 5: Coefficients sorte by magnitue of the centroi form M k for the regular facet of the Butterfly scheme in the power iteration for k =, 2,..., 7. 4 Gaussian Elimination [Hakenberg et al. 204] solve the linear system (9) by Gaussian elimination to obtain the exact moment forms for moerate non-regular vertex valences. If the subivision weights are rational numbers, the tensor coefficients are also rational. For the Doo-Sabin scheme, the egree is limite to {0,, 2; for Loop, {0, ; an for Catmull-Clark, = 0. The construction of matrix A for moments beyon the liste egrees is prohibitive ue to the large number of form coefficients. 5 Geometric Series We use the tensor iteration (4) to approximate the inertia forms for the Loop scheme, as well as the centroi an inertia forms for Catmull-Clark surfaces an moerate vertex valences τ(f). The collection of tensors is compute within a ay using an implementation in Mathematica. The memory footprint of the iteration is low, since only the tensor M k nees to be store. The numerical accuracy of the coefficients increases with each iteration step. Figure 6 shows the rate of convergence Figure 6: Convergence of Catmull-Clark inertia forms for nonregular valences up to 8. The slope is etermine by the subominant eigenvalue λ of the subivision matrix: For vertex valence 3, λ 0.4; for vertex valence 8, λ 0.6. Machine precision in the form coefficients is surpasse after no more than 50 iteration steps. 5 Future Work We plan to investigate the use of the iterative methos to calculate inner proucts ifferent from (5). Such inner proucts arise for instance when moeling a solution space to elliptic partial ifferential equations. Acknowlegements The authors thank Brian Paen for his suggestions that le to the improvement of the article
5 References CATMULL, E., AND CLARK, J Recursively generate b- spline surfaces on arbitrary topological meshes. Computer-Aie Design 0 (Nov.), DOO, D., AND SABIN, M Behaviour of recursive ivision surfaces near extraorinary points. Computer-Aie Design 0, DUBUC, S Interpolation through an iterative scheme. Journal of Mathematical Analysis an Applications 4, DYN, N., LEVIN, D., AND GREGORY, J A butterfly subivision scheme for surface interpolation with tension control. ACM Transactions on Graphics 9, GONZALEZ-OCHOA, C., MCCAMMON, S., AND PETERS, J A butterfly subivision scheme for surface interpolation with tension control. ACM Transactions on Graphics 7, HAKENBERG, J., AND REIF, U., 206. Iterative computation of moment forms for subivision curves an surfaces. youtube.com/watch?v=oeuojofggg. HAKENBERG, J., AND REIF, U On the volume of sets boune by refinable functions. Applie Mathematics an Computation 272, (Jan.), 2 9. HAKENBERG, J., REIF, U., SCHAEFER, S., AND WARREN, J., 204. On moments of sets boune by subivision surfaces. LOOP, C Smooth subivision surfaces base on triangles. Master s thesis, University of Utah. PETERS, J., AND NASRI, A. H Computing volumes of solis enclose by recursive subivision surfaces. Computer Graphics Forum 6, SCHWALD, B Exakte Volumenberechnung von urch Doo- Sabin-Flächen begrenzten Körpern. Master s thesis, Universität Stuttgart. VON MISES, R., AND POLLACZEK-GEIRINGER, H Praktische verfahren er gleichungsauflösung. Zeitschrift für Angewante Mathematik un Mechanik 9, WARREN, J., AND WEIMER, H Subivision Methos for Geometric Design: A Constructive Approach. Morgan Kaufmann.
UNIVERSITY 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 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 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 informationKinematic Analysis of a Family of 3R Manipulators
Kinematic Analysis of a Family of R Manipulators Maher Baili, Philippe Wenger an Damien Chablat Institut e Recherche en Communications et Cybernétique e Nantes, UMR C.N.R.S. 6597 1, rue e la Noë, BP 92101,
More informationTransient analysis of wave propagation in 3D soil by using the scaled boundary finite element method
Southern Cross University epublications@scu 23r Australasian Conference on the Mechanics of Structures an Materials 214 Transient analysis of wave propagation in 3D soil by using the scale bounary finite
More informationA Classification of 3R Orthogonal Manipulators by the Topology of their Workspace
A Classification of R Orthogonal Manipulators by the Topology of their Workspace Maher aili, Philippe Wenger an Damien Chablat Institut e Recherche en Communications et Cybernétique e Nantes, UMR C.N.R.S.
More informationB-Splines and NURBS Week 5, Lecture 9
CS 430/585 Computer Graphics I B-Splines an NURBS Week 5, Lecture 9 Davi Breen, William Regli an Maxim Peysakhov Geometric an Intelligent Computing Laboratory Department of Computer Science Drexel University
More informationExercises of PIV. incomplete draft, version 0.0. October 2009
Exercises of PIV incomplete raft, version 0.0 October 2009 1 Images Images are signals efine in 2D or 3D omains. They can be vector value (e.g., color images), real (monocromatic images), complex or binary
More information5th International Conference on Advanced Design and Manufacturing Engineering (ICADME 2015)
5th International Conference on Avance Design an Manufacturing Engineering (ICADME 25) Research on motion characteristics an application of multi egree of freeom mechanism base on R-W metho Xiao-guang
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 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 informationResearch Article Inviscid Uniform Shear Flow past a Smooth Concave Body
International Engineering Mathematics Volume 04, Article ID 46593, 7 pages http://x.oi.org/0.55/04/46593 Research Article Invisci Uniform Shear Flow past a Smooth Concave Boy Abullah Mura Department of
More informationBIJECTIONS FOR PLANAR MAPS WITH BOUNDARIES
BIJECTIONS FOR PLANAR MAPS WITH BOUNDARIES OLIVIER BERNARDI AND ÉRIC FUSY Abstract. We present bijections for planar maps with bounaries. In particular, we obtain bijections for triangulations an quarangulations
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 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 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 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 informationClassical Mechanics Examples (Lagrange Multipliers)
Classical Mechanics Examples (Lagrange Multipliers) Dipan Kumar Ghosh Physics Department, Inian Institute of Technology Bombay Powai, Mumbai 400076 September 3, 015 1 Introuction We have seen that the
More information6 Gradient Descent. 6.1 Functions
6 Graient Descent In this topic we will iscuss optimizing over general functions f. Typically the function is efine f : R! R; that is its omain is multi-imensional (in this case -imensional) an output
More informationBends, Jogs, And Wiggles for Railroad Tracks and Vehicle Guide Ways
Ben, Jogs, An Wiggles for Railroa Tracks an Vehicle Guie Ways Louis T. Klauer Jr., PhD, PE. Work Soft 833 Galer Dr. Newtown Square, PA 19073 lklauer@wsof.com Preprint, June 4, 00 Copyright 00 by Louis
More informationApproximation with Active B-spline Curves and Surfaces
Approximation with Active B-spline Curves an Surfaces Helmut Pottmann, Stefan Leopolseer, Michael Hofer Institute of Geometry Vienna University of Technology Wiener Hauptstr. 8 10, Vienna, Austria pottmann,leopolseer,hofer
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 informationGeneralized Edge Coloring for Channel Assignment in Wireless Networks
TR-IIS-05-021 Generalize Ege Coloring for Channel Assignment in Wireless Networks Chun-Chen Hsu, Pangfeng Liu, Da-Wei Wang, Jan-Jan Wu December 2005 Technical Report No. TR-IIS-05-021 http://www.iis.sinica.eu.tw/lib/techreport/tr2005/tr05.html
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 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 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 informationLinear First-Order PDEs
MODULE 2: FIRST-ORDER PARTIAL DIFFERENTIAL EQUATIONS 9 Lecture 2 Linear First-Orer PDEs The most general first-orer linear PDE has the form a(x, y)z x + b(x, y)z y + c(x, y)z = (x, y), (1) where a, b,
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 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 informationCurve Subdivision in SE(2)
Curve Subdivision in SE(2) Jan Hakenberg, ETH Zürich 2018-07-26 Figure: A point in the special Euclidean group SE(2) consists of a position in the plane and a heading. The figure shows two rounds of cubic
More informationAdvanced Geometric Modeling CPSC789
Advanced Geometric Modeling CPSC789 Fall 2004 General information about the course CPSC 789 Advanced Geometric Modeling Fall 2004 Lecture Time and Place ENF 334 TR 9:30 10:45 Instructor : Office: MS 618
More informationLearning convex bodies is hard
Learning convex boies is har Navin Goyal Microsoft Research Inia navingo@microsoftcom Luis Raemacher Georgia Tech lraemac@ccgatecheu Abstract We show that learning a convex boy in R, given ranom samples
More informationPairwise alignment using shortest path algorithms, Gunnar Klau, November 29, 2005, 11:
airwise alignment using shortest path algorithms, Gunnar Klau, November 9,, : 3 3 airwise alignment using shortest path algorithms e will iscuss: it graph Dijkstra s algorithm algorithm (GDU) 3. References
More informationAnimated Surface Pasting
Animate Surface Pasting Clara Tsang an Stephen Mann Computing Science Department University of Waterloo 200 University Ave W. Waterloo, Ontario Canaa N2L 3G1 e-mail: clftsang@cgl.uwaterloo.ca, smann@cgl.uwaterloo.ca
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 informationComputer Graphics Chapter 7 Three-Dimensional Viewing Viewing
Computer Graphics Chapter 7 Three-Dimensional Viewing Outline Overview of Three-Dimensional Viewing Concepts The Three-Dimensional Viewing Pipeline Three-Dimensional Viewing-Coorinate Parameters Transformation
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 information1 Surprises in high dimensions
1 Surprises in high imensions Our intuition about space is base on two an three imensions an can often be misleaing in high imensions. It is instructive to analyze the shape an properties of some basic
More informationA C 2 Four-Point Subdivision Scheme with Fourth Order Accuracy and its Extensions
A C 2 Four-Point Subdivision Scheme with Fourth Order Accuracy and its Extensions Nira Dyn School of Mathematical Sciences Tel Aviv University Michael S. Floater Department of Informatics University of
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 informationMultilevel Linear Dimensionality Reduction using Hypergraphs for Data Analysis
Multilevel Linear Dimensionality Reuction using Hypergraphs for Data Analysis Haw-ren Fang Department of Computer Science an Engineering University of Minnesota; Minneapolis, MN 55455 hrfang@csumneu ABSTRACT
More informationRefinement of scene depth from stereo camera ego-motion parameters
Refinement of scene epth from stereo camera ego-motion parameters Piotr Skulimowski, Pawel Strumillo An algorithm for refinement of isparity (epth) map from stereoscopic sequences is propose. The metho
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 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 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 informationGeneralized Edge Coloring for Channel Assignment in Wireless Networks
Generalize Ege Coloring for Channel Assignment in Wireless Networks Chun-Chen Hsu Institute of Information Science Acaemia Sinica Taipei, Taiwan Da-wei Wang Jan-Jan Wu Institute of Information Science
More informationMOTION IDENTIFICATION OF PLANAR FREE-FORM OBJECTS
MOTION IDENTIFICATION OF PLANAR FREE-FORM OBJECTS Mustafa Unel Faculty of Engineering an Natural Sciences Sabanci University Orhanli-Tuzla 34956, Istanbul, Turkey email: munel@sabanciuniv.eu ABSTRACT It
More informationMathematical Tools in Computer Graphics with C# Implementations Table of Contents
Mathematical Tools in Computer Graphics with C# Implementations by Hardy Alexandre, Willi-Hans Steeb, World Scientific Publishing Company, Incorporated, 2008 Table of Contents List of Figures Notation
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 informationAlgebraic transformations of Gauss hypergeometric functions
Algebraic transformations of Gauss hypergeometric functions Raimunas Viūnas Faculty of Mathematics, Kobe University Abstract This article gives a classification scheme of algebraic transformations of Gauss
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 informationOptimal Oblivious Path Selection on the Mesh
Optimal Oblivious Path Selection on the Mesh Costas Busch Malik Magon-Ismail Jing Xi Department of Computer Science Rensselaer Polytechnic Institute Troy, NY 280, USA {buschc,magon,xij2}@cs.rpi.eu Abstract
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 informationNew Geometric Interpretation and Analytic Solution for Quadrilateral Reconstruction
New Geometric Interpretation an Analytic Solution for uarilateral Reconstruction Joo-Haeng Lee Convergence Technology Research Lab ETRI Daejeon, 305 777, KOREA Abstract A new geometric framework, calle
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 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 informationParameterization. Michael S. Floater. November 10, 2011
Parameterization Michael S. Floater November 10, 2011 Triangular meshes are often used to represent surfaces, at least initially, one reason being that meshes are relatively easy to generate from point
More informationHere are a couple of warnings to my students who may be here to get a copy of what happened on a day that you missed.
Preface Here are my online notes for my Calculus I course that I teach here at Lamar University. Despite the fact that these are my class notes, they shoul be accessible to anyone wanting to learn Calculus
More informationSURVIVABLE IP OVER WDM: GUARANTEEEING MINIMUM NETWORK BANDWIDTH
SURVIVABLE IP OVER WDM: GUARANTEEEING MINIMUM NETWORK BANDWIDTH Galen H Sasaki Dept Elec Engg, U Hawaii 2540 Dole Street Honolul HI 96822 USA Ching-Fong Su Fuitsu Laboratories of America 595 Lawrence Expressway
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 informationDistributed Decomposition Over Hyperspherical Domains
Distribute Decomposition Over Hyperspherical Domains Aron Ahmaia 1, Davi Keyes 1, Davi Melville 2, Alan Rosenbluth 2, Kehan Tian 2 1 Department of Applie Physics an Applie Mathematics, Columbia University,
More informationCONSTRUCTION AND ANALYSIS OF INVERSIONS IN S 2 AND H 2. Arunima Ray. Final Paper, MATH 399. Spring 2008 ABSTRACT
CONSTUCTION AN ANALYSIS OF INVESIONS IN S AN H Arunima ay Final Paper, MATH 399 Spring 008 ASTACT The construction use to otain inversions in two-imensional Eucliean space was moifie an applie to otain
More informationParameterization of triangular meshes
Parameterization of triangular meshes Michael S. Floater November 10, 2009 Triangular meshes are often used to represent surfaces, at least initially, one reason being that meshes are relatively easy to
More informationFast Fractal Image Compression using PSO Based Optimization Techniques
Fast Fractal Compression using PSO Base Optimization Techniques A.Krishnamoorthy Visiting faculty Department Of ECE University College of Engineering panruti rishpci89@gmail.com S.Buvaneswari Visiting
More informationFreeform Curves on Spheres of Arbitrary Dimension
Freeform Curves on Spheres of Arbitrary Dimension Scott Schaefer and Ron Goldman Rice University 6100 Main St. Houston, TX 77005 sschaefe@rice.edu and rng@rice.edu Abstract Recursive evaluation procedures
More informationFitting Subdivision Surfaces to Unorganized Point Data Using SDM
Fitting Subivision Surfaces to Unorganize Point Data Using SDM Kin-Shing D. Cheng Wenping Wang Hong Qin Kwan-Yee K. Wong Huaiping Yang Yang Liu Department of Computer Science he University of Hong Kong
More informationApproximating Subdivision Surfaces with Gregory Patches for Hardware Tessellation
Approximating Subdivision Surfaces with Gregory Patches for Hardware Tessellation Charles Loop Microsoft Research Scott Schaefer Texas A&M University Tianyun Ni NVIDIA Ignacio Castaño NVIDIA Goal Real-Time
More informationModified Catmull-Clark Methods for Modelling, Reparameterization and Grid Generation
Modified Catmull-Clark Methods for Modelling, Reparameterization and Grid Generation Karl-Heinz Brakhage RWTH Aachen, 55 Aachen, Deutschland, Email: brakhage@igpm.rwth-aachen.de Abstract In this paper
More informationLarge Monochromatic Components in Two-colored Grids
Large Monochromatic Components in Two-colore Gris Jiří Matoušek 1 Aleš Přívětivý 2 Department of Applie Mathematics 1,2 an Institute of Theoretical Computer Science (ITI) 1 Charles University Malostranské
More informationThe goal is the definition of points with numbers and primitives with equations or functions. The definition of points with numbers requires a
The goal is the definition of points with numbers and primitives with equations or functions. The definition of points with numbers requires a coordinate system and then the measuring of the point with
More informationClassifying Facial Expression with Radial Basis Function Networks, using Gradient Descent and K-means
Classifying Facial Expression with Raial Basis Function Networks, using Graient Descent an K-means Neil Allrin Department of Computer Science University of California, San Diego La Jolla, CA 9237 nallrin@cs.ucs.eu
More informationParticle Swarm Optimization Based on Smoothing Approach for Solving a Class of Bi-Level Multiobjective Programming Problem
BULGARIAN ACADEMY OF SCIENCES CYBERNETICS AND INFORMATION TECHNOLOGIES Volume 17, No 3 Sofia 017 Print ISSN: 1311-970; Online ISSN: 1314-4081 DOI: 10.1515/cait-017-0030 Particle Swarm Optimization Base
More information(a) Find the equation of the plane that passes through the points P, Q, and R.
Math 040 Miterm Exam 1 Spring 014 S o l u t i o n s 1 For given points P (, 0, 1), Q(, 1, 0), R(3, 1, 0) an S(,, 0) (a) Fin the equation of the plane that passes through the points P, Q, an R P Q = 0,
More informationHoly Halved Heaquarters Riddler
Holy Halve Heaquarters Riler Anonymous Philosopher June 206 Laser Larry threatens to imminently zap Riler Heaquarters (which is of regular pentagonal shape with no courtyar or other funny business) with
More informationAn Interpolatory Subdivision for Volumetric Models over Simplicial Complexes
An Interpolatory Subdivision for Volumetric Models over Simplicial Complexes Yu-Sung Chang Kevin T. McDonnell Hong Qin Department of Computer Science State University of New York at Stony Brook {yusung
More informationElement Quality Metrics for Higher-Order Bernstein Bézier Elements
Element Quality Metrics for Higher-Order Bernstein Bézier Elements Luke Engvall and John A. Evans Abstract In this note, we review the interpolation theory for curvilinear finite elements originally derived
More informationCalculation on diffraction aperture of cube corner retroreflector
November 10, 008 / Vol., No. 11 / CHINESE OPTICS LETTERS 8 Calculation on iffraction aperture of cube corner retroreflector Song Li (Ó Ø, Bei Tang (», an Hui Zhou ( ï School of Electronic Information,
More informationA Parameterized Mask Model for Lithography Simulation
A Parameterize Mask Moel for Lithography Simulation henhai hu Caence Research Labs, Berkeley, CA 94704 zhzhu@caence.com 1. ABSTRACT e formulate the mask moeling as a parametric moel orer reuction problem
More informationDepartment of Computer Science, POSTECH, Pohang , Korea. (x 0 (t); y 0 (t)) 6= (0; 0) and N(t) is well dened on the
Comparing Oset Curve Approximation Methos Gershon er +, In-Kwon, an Myung-Soo Kim + Department of Computer Science, Technion, IIT, Haifa 32000, Israel Department of Computer Science, POSTECH, Pohang 790-784,
More informationLoop Scheduling and Partitions for Hiding Memory Latencies
Loop Scheuling an Partitions for Hiing Memory Latencies Fei Chen Ewin Hsing-Mean Sha Dept. of Computer Science an Engineering University of Notre Dame Notre Dame, IN 46556 Email: fchen,esha @cse.n.eu Tel:
More informationImpact of changing the position of the tool point on the moving platform on the dynamic performance of a 3RRR planar parallel manipulator
IOSR Journal of Mechanical an Civil Engineering (IOSR-JMCE) e-issn: 78-84,p-ISSN: 0-4X, Volume, Issue 4 Ver. I (Jul. - Aug. 05), PP 7-8 www.iosrjournals.org Impact of changing the position of the tool
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 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 informationFINDING OPTICAL DISPERSION OF A PRISM WITH APPLICATION OF MINIMUM DEVIATION ANGLE MEASUREMENT METHOD
Warsaw University of Technology Faculty of Physics Physics Laboratory I P Joanna Konwerska-Hrabowska 6 FINDING OPTICAL DISPERSION OF A PRISM WITH APPLICATION OF MINIMUM DEVIATION ANGLE MEASUREMENT METHOD.
More informationUnknown Radial Distortion Centers in Multiple View Geometry Problems
Unknown Raial Distortion Centers in Multiple View Geometry Problems José Henrique Brito 1,2, Rolan Angst 3, Kevin Köser 3, Christopher Zach 4, Pero Branco 2, Manuel João Ferreira 2, Marc Pollefeys 3 1
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 informationA C 2 Four-Point Subdivision Scheme with Fourth Order Accuracy and its Extensions
A C 2 Four-Point Subdivision Scheme with Fourth Order Accuracy and its Extensions Nira Dyn Michael S. Floater Kai Hormann Abstract. We present a new four-point subdivision scheme that generates C 2 curves.
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 informationCOMPUTER AIDED ENGINEERING DESIGN (BFF2612)
COMPUTER AIDED ENGINEERING DESIGN (BFF2612) BASIC MATHEMATICAL CONCEPTS IN CAED by Dr. Mohd Nizar Mhd Razali Faculty of Manufacturing Engineering mnizar@ump.edu.my COORDINATE SYSTEM y+ y+ z+ z+ x+ RIGHT
More informationImage compression predicated on recurrent iterated function systems
2n International Conference on Mathematics & Statistics 16-19 June, 2008, Athens, Greece Image compression preicate on recurrent iterate function systems Chol-Hui Yun *, Metzler W. a an Barski M. a * Faculty
More informationMatrix-valued 4-point Spline and 3-point Non-spline Interpolatory Curve Subdivision Schemes
Matrix-valued 4-point Spline and -point Non-spline Interpolatory Curve Subdivision Schemes Charles K. Chui, Qingtang Jiang Department of Mathematics and Computer Science University of Missouri St. Louis
More informationUsing Vector and Raster-Based Techniques in Categorical Map Generalization
Thir ICA Workshop on Progress in Automate Map Generalization, Ottawa, 12-14 August 1999 1 Using Vector an Raster-Base Techniques in Categorical Map Generalization Beat Peter an Robert Weibel Department
More informationSung-Eui Yoon ( 윤성의 )
CS480: Computer Graphics Curves and Surfaces Sung-Eui Yoon ( 윤성의 ) Course URL: http://jupiter.kaist.ac.kr/~sungeui/cg Today s Topics Surface representations Smooth curves Subdivision 2 Smooth Curves and
More informationOn-Line Geometric Modeling Notes REFINEMENT
On-Line Geometric Modeling Notes REFINEMENT Kenneth I Joy Visualization and Graphics Research Group Department of Computer Science University of California, Davis Overview Bézier curves, B-spline curves
More informationSkyline Community Search in Multi-valued Networks
Syline Community Search in Multi-value Networs Rong-Hua Li Beijing Institute of Technology Beijing, China lironghuascut@gmail.com Jeffrey Xu Yu Chinese University of Hong Kong Hong Kong, China yu@se.cuh.eu.h
More informationEstimation of large-amplitude motion and disparity fields: Application to intermediate view reconstruction
c 2000 SPIE. Personal use of this material is permitte. However, permission to reprint/republish this material for avertising or promotional purposes or for creating new collective works for resale or
More informationNon-homogeneous Generalization in Privacy Preserving Data Publishing
Non-homogeneous Generalization in Privacy Preserving Data Publishing W. K. Wong, Nios Mamoulis an Davi W. Cheung Department of Computer Science, The University of Hong Kong Pofulam Roa, Hong Kong {wwong2,nios,cheung}@cs.hu.h
More information4.2 Implicit Differentiation
6 Chapter 4 More Derivatives 4. Implicit Differentiation What ou will learn about... Implicitl Define Functions Lenses, Tangents, an Normal Lines Derivatives of Higher Orer Rational Powers of Differentiable
More informationUNIT 3 EXPRESSIONS AND EQUATIONS Lesson 3: Creating Quadratic Equations in Two or More Variables
Guided Practice Example 1 Find the y-intercept and vertex of the function f(x) = 2x 2 + x + 3. Determine whether the vertex is a minimum or maximum point on the graph. 1. Determine the y-intercept. The
More information