Quadratic and cubic b-splines by generalizing higher-order voronoi diagrams
|
|
- Damian Sharp
- 6 years ago
- Views:
Transcription
1 Quadratic and cubic b-splines by generalizing higher-order voronoi diagrams Yuanxin Liu and Jack Snoeyink Joshua Levine April 18, 2007 Computer Science and Engineering, The Ohio State University 1 / 24
2 Background Motivation Desire a way to a define a function of arbitrary degree by piecewise polynomials. Can give lots of useful guarantees, like Smoothness, Partition of unity, and Compact support. Often used to interpolate other functions. Also used in the context of defining curves/surfaces that fit through a set of points. Computer Science and Engineering, The Ohio State University 2 / 24
3 Univariate B-Splines Basic Definition A (univariate) degree k B-spline is a smooth degree k piecewise polynomial function defined over a set K of k + 2 reals called knots. From an arbitrary set K R of knots, we construct a basis for the degree k B-spline by forming a B-spline for each sequence of k + 2 consecutive knots in K. See example at Computer Science and Engineering, The Ohio State University 3 / 24
4 Multivariate Splines A Challenge Univariate B-splines work well, but we also would like to have splines which are multivariate, or have domain R s for some s. Why is this hard? In R, the knots have a total ordering on them, so we can easily select subsets to build the spline bases. We desire an s-variate analogue that has the following properties: 1. Each spline is define as a piecewise polynomial for a partition of the space determined by K. 2. Splines are smooth. 3. (Polynomial Reproducibility) Can reproduce all polynomials of degree k. 4. (Specialization) For s = 1, the formulation reduces to the univariate case. Of these, 3. and 4. are the most challenging. Computer Science and Engineering, The Ohio State University 4 / 24
5 Multivariate Splines One solution: Simplex Splines A degree k simplex-spline is a smooth, degree k, piecewise-polynomial function defined over a set of k + s + 1 points X R s called knots. Geometric interpretation: The shadow of some simplex [Y ] in R k+s, projecting the vertices onto X. How to choose the right simplices? Computer Science and Engineering, The Ohio State University 5 / 24
6 Neamtu s Solution Higher-Order Voronoi Diagrams The Voronoi diagram of order K of the set P is the subdivision of R s into regions, called Voronoi cells (of degree k) such that each cell contains the same k nearest points from K. (Alternatively) let X K be such that X = k. The Voronoi cell, c(x ) is defined as c(x ) := {x R s d(x, X ) d(x, Y ), Y K, Y = k}. Computer Science and Engineering, The Ohio State University 6 / 24
7 Neamtu s Solution Higher-Order Voronoi Diagrams The Voronoi diagram of order K of the set P is the subdivision of R s into regions, called Voronoi cells (of degree k) such that each cell contains the same k nearest points from K. (Alternatively) let X K be such that X = k. The Voronoi cell, c(x ) is defined as c(x ) := {x R s d(x, X ) d(x, Y ), Y K, Y = k}. Computer Science and Engineering, The Ohio State University 6 / 24
8 Neamtu s Solution Delaunay Configurations In a knot set K R s, a k + s + 1-subset X K is a degree k Delaunay configuration iff X can be partitioned into two sets, an s + 1-set t and a k-set I such that the circumsphere of t has, of all points in K, exactly I inside. These are connected to the k-th order Voronoi diagram. Computer Science and Engineering, The Ohio State University 7 / 24
9 Neamtu s Solution Formulation of the spline A set of k + s + 1 knots X R s defines a degree k, s-variate simplex spline M( X ) as follows. Let π( ) : R k+s R s be the vertical projection map. Let Y be any set of points in R k+s such that π(y ) = X. Then the simplex spline is the shadow cast by the simplex [Y ], assumed to be uniformly dense: M(x X ) := vol k{y y Y, π(y) = x}. vol k+s [Y ] Each subset X is called a configuration of degree k w.r.t. K. A set Γ of configurations defines a simplex spline of degree k by span{m( X ) X Γ}. Computer Science and Engineering, The Ohio State University 8 / 24
10 Neamtu s Solution Extension Neamtu s solution is elegant, but restrictive. All simplices must be Delaunay configurations. Why do we need this? One goal of the paper is generalizing his work to a broader class of triangulations, based on two observations: 1. For polynomial reproduction, we only require few properties of the Delaunay configurations: Degree 0 configurations form a triangulation Configurations of degree i and i + 1 match facets. 2. In two dimensions, the algorithm of Lee (for computing higher order Voronoi diagrams) may be generalized to iteratively compute configurations of successively higher degrees Computer Science and Engineering, The Ohio State University 9 / 24
11 Facet-Matching Some Notation In R s an n-dimensional simplex is the convex hull of a set of n + 1 affinely independent points (vertices). Often thought of as an ordered tuple [t 0,..., t n ] For a simplex t, the simplex defined by t \ {v} is the facet of t opposite v. We can drop vertices: i t := ( 1) i [t 0,..., t i 1, t i+1,..., t n ]. And add them: v f := [v, f 0,..., f n 1 ]. Or swap: v i t := [t 0,..., t i 1, v, t i+1,..., t n ]. Remark: And of course the usual combinatorial properties of summing simplices up and taking their boundary still apply. In their notation t := { i t i = 0,..., n}. Computer Science and Engineering, The Ohio State University 10 / 24
12 Facet-Matching Boundary-Interior Configurations In R s, a degree k boundary-interior configuration is a pair (t, I ) where t is a tuple of s + 1 points (representing an s-simplex) and I, the interior set, is a set of k points disjoint from t. A normalized simplex spline over (t, I ) is N( t, I ) := vol(t)m( t I ). Note the union t I. This makes the notion consistent with the Delaunay configuration used before. Much like facets of a simplex are constructed by dropping vertices, we construct facet configurations by dropping vertices, but maintaining the interior. For the configuration (t, I ), v t, I is the facet configuration opposite of v. For a set Γ of configurations, F (Γ) denotes all the facet configurations of Γ. Computer Science and Engineering, The Ohio State University 11 / 24
13 Facet-Matching Facet Matching Property Definition (3.1) Let Γ and Γ be two sets of configurations of degree k and k + 1, resp. Γ and Γ are said to match facets iff 1. for each facet configuration ( u s, I ) F (Γ), (a) there is either another facet configuration ( v t, J) F (Γ) s.t. v t = u s and I v = J u, (b) or there is a facet configuration ( v t, I ) F (Γ ) s.t. v t = u s and I v = I, 2. and for each facet configuration ( v t, I ) F (Γ ), (a) there is either another facet configuration ( u s, J ) F (Γ ) s.t. v t = u s and I = J, (b) or there is another facet configuration ( u s, I ) F (Γ) s.t. v t = u s and I = I v. Computer Science and Engineering, The Ohio State University 12 / 24
14 Facet-Matching Facet Matching Property (Example) Computer Science and Engineering, The Ohio State University 13 / 24
15 Facet-Matching Polynomial Reproduction Theorem Theorem (3.2) Let (Γ 0,..., Γ k ) be sets of configurations of degree 0 to k in R s. If Γ 0 is a triangulation of the space R s, and, for 0 i k 1, sets Γ i and Γ i+1 match facets, then the normalized simplex splines associated with Γ k reproduce all degree k polynomials. For a degree k polynomial p with polar form P, p = P(I )N( t, I ). (t,i ) Γ k The polar form P of a degree n, s-variate polynomial p is the unique function of n vector variables x 1,..., x n R s that is symmetric, affine in each of the n variables, and such that for every x R s the value p(x) is equal to P on the diagonal x 1 =... = x n = x, or p(x) = P(x 1,..., x n ). Computer Science and Engineering, The Ohio State University 14 / 24
16 Definition Preliminaries Theorem 3.2 tells us that all we need is a sequence of configurations (Γ 0,..., Γ k ) which match facets, but how do we find such a sequence? One method is the link triangulation procedure, which is a generalization of the (dual) construction for higher order Voronoi diagrams. Idea: For a configuration (t, I ) Γ i, we should be able to locally construct a configuration (s, J) for Γ i+1 which can ensure facet matching. Computer Science and Engineering, The Ohio State University 15 / 24
17 Definition Facet Matching (Restated) The interior sets I (Γ) := {I (t, I ) Γ}. The vertices V (Γ) := {I {v} v t, (t, I ) Γ}. For a set of k points, I, the neighborhood nb(i, Γ) := {t (t, I ) Γ}. For a set of k + 1 points, J, The link lk(j, Γ) := { v t v J, (t, J \ {v}) Γ}. Lemma (4.1) Let Γ and Γ be two sets of configurations of degree k and k + 1. Then Γ and Γ match facets iff for any set of k + 1 points J, the boundary of the neighborhood of J in Γ is the sum of the edges in the link of J in Γ: nb(j, Γ ) = lk(j, Γ). Computer Science and Engineering, The Ohio State University 16 / 24
18 Definition Algorithm - Link Triangulation Procedure 1 LinkTriangulate(Γ i ) { 2 Γ i FORALL J V (Γ i ) { 5 P lk(j, Γ i ) 6 IF (P ) { 7 T triangulate(p) 8 Γ i+1 Γ i+1 {(t, J)} t T 9 } 10 } RETURN Γ i+1 13 } //end LinkTriangulate Computer Science and Engineering, The Ohio State University 17 / 24
19 Definition Remarks The link triangulation procedure is well defined if the links of the input configurations are simple polygons. Then there is flexibility in the algorithm as to how to triangulate the link, this gives local control over the construction The authors prove that this can be applied 3 times, starting with a triangulation Γ 0. Hence we can create a sequence of up to degree 3 configurations Computer Science and Engineering, The Ohio State University 18 / 24
20 Centroid Triangulations A Dual to order k Voronoi Diagrams By an observation of Aurenhammer, an order k Voronoi diagram has a dual triangulation: For each of the k sites defining a Voronoi cell, compute its centroid of the lifted points. The projection of the convex hull of the centroids is a triangulation that is dual. This triangulation can also be constructed by mapping each Delaunay configuration corresponding to the Voronoi vertices to a set of centroid triangles. We can apply the same construction to configurations made by the link triangulation procedure the result is a planar triangulation. Computer Science and Engineering, The Ohio State University 19 / 24
21 Centroid Triangulations Centroid Triangles For a set of points V R 2, let V denote the centroid of V. A triangle is an order k centroid triangle if 1. its vertices are centroids of k-sets and 2. for any of its edges {U, V }, U V = k 1. The type of an order k centroid triangle [U, V, W ] is the number U V W k. This value is always either 1 or 2. Computer Science and Engineering, The Ohio State University 20 / 24
22 Centroid Triangulations Connection to configurations For a set of knots K, there is a bijection between the set of all possible order k centroid triangles and the set of all possible degree k 1 and k 2 configurations in K. The following average maps map from a degree i configuration (t, I ) to an order i + 1 and i + 2 centroid triangle: A 1 (t, I ) := [I {t 0 }, I {t 1 }, I {t 2 }] A 2 (t, I ) := [I 0 t, I 1 t, I 2 t] Computer Science and Engineering, The Ohio State University 21 / 24
23 Centroid Triangulations Example One simple description of the link triangulation procedure: It triangulates the holes left after placing the triangles A 2 (Γ) on the plane. Computer Science and Engineering, The Ohio State University 22 / 24
24 Centroid Triangulations Example One simple description of the link triangulation procedure: It triangulates the holes left after placing the triangles A 2 (Γ) on the plane. Computer Science and Engineering, The Ohio State University 22 / 24
25 A Coarser Basis Collecting Splines Given Γ a set of degree k configurations, it makes sense to partition the simplex splines based on neighborhoods to form a coarser basis: {B I ( ) := N( t, I )} I I (Γ). t nb(i,γ) This doesn t break the two properties we were originally seeking. Computer Science and Engineering, The Ohio State University 23 / 24
26 Conclusions Conclusions Give a generalization of Neamtu s Delaunay configurations for multivariate B-splines that works for a larger class of triangulations. Able to prove some nice properties, but only for low orders. Experimented with reproducing the ZP-Element, a classic box spline. Computer Science and Engineering, The Ohio State University 24 / 24
Planar Graphs. 1 Graphs and maps. 1.1 Planarity and duality
Planar Graphs In the first half of this book, we consider mostly planar graphs and their geometric representations, mostly in the plane. We start with a survey of basic results on planar graphs. This chapter
More informationMonotone Paths in Geometric Triangulations
Monotone Paths in Geometric Triangulations Adrian Dumitrescu Ritankar Mandal Csaba D. Tóth November 19, 2017 Abstract (I) We prove that the (maximum) number of monotone paths in a geometric triangulation
More informationVoronoi diagram and Delaunay triangulation
Voronoi diagram and Delaunay triangulation Ioannis Emiris & Vissarion Fisikopoulos Dept. of Informatics & Telecommunications, University of Athens Computational Geometry, spring 2015 Outline 1 Voronoi
More informationChapter 8. Voronoi Diagrams. 8.1 Post Oce Problem
Chapter 8 Voronoi Diagrams 8.1 Post Oce Problem Suppose there are n post oces p 1,... p n in a city. Someone who is located at a position q within the city would like to know which post oce is closest
More informationManifolds. Chapter X. 44. Locally Euclidean Spaces
Chapter X Manifolds 44. Locally Euclidean Spaces 44 1. Definition of Locally Euclidean Space Let n be a non-negative integer. A topological space X is called a locally Euclidean space of dimension n if
More informationSurfaces, meshes, and topology
Surfaces from Point Samples Surfaces, meshes, and topology A surface is a 2-manifold embedded in 3- dimensional Euclidean space Such surfaces are often approximated by triangle meshes 2 1 Triangle mesh
More informationCAT(0)-spaces. Münster, June 22, 2004
CAT(0)-spaces Münster, June 22, 2004 CAT(0)-space is a term invented by Gromov. Also, called Hadamard space. Roughly, a space which is nonpositively curved and simply connected. C = Comparison or Cartan
More informationCS 532: 3D Computer Vision 14 th Set of Notes
1 CS 532: 3D Computer Vision 14 th Set of Notes Instructor: Philippos Mordohai Webpage: www.cs.stevens.edu/~mordohai E-mail: Philippos.Mordohai@stevens.edu Office: Lieb 215 Lecture Outline Triangulating
More informationComputational Geometry
Lecture 12: Lecture 12: Motivation: Terrains by interpolation To build a model of the terrain surface, we can start with a number of sample points where we know the height. Lecture 12: Motivation: Terrains
More informationPreferred directions for resolving the non-uniqueness of Delaunay triangulations
Preferred directions for resolving the non-uniqueness of Delaunay triangulations Christopher Dyken and Michael S. Floater Abstract: This note proposes a simple rule to determine a unique triangulation
More informationVoronoi diagrams Delaunay Triangulations. Pierre Alliez Inria
Voronoi diagrams Delaunay Triangulations Pierre Alliez Inria Voronoi Diagram Voronoi Diagram Voronoi Diagram The collection of the non-empty Voronoi regions and their faces, together with their incidence
More informationVoronoi Diagrams in the Plane. Chapter 5 of O Rourke text Chapter 7 and 9 of course text
Voronoi Diagrams in the Plane Chapter 5 of O Rourke text Chapter 7 and 9 of course text Voronoi Diagrams As important as convex hulls Captures the neighborhood (proximity) information of geometric objects
More information6.854J / J Advanced Algorithms Fall 2008
MIT OpenCourseWare http://ocw.mit.edu 6.854J / 18.415J Advanced Algorithms Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 18.415/6.854 Advanced
More information08 - Designing Approximating Curves
08 - Designing Approximating Curves Acknowledgement: Olga Sorkine-Hornung, Alexander Sorkine-Hornung, Ilya Baran Last time Interpolating curves Monomials Lagrange Hermite Different control types Polynomials
More informationPebble Sets in Convex Polygons
2 1 Pebble Sets in Convex Polygons Kevin Iga, Randall Maddox June 15, 2005 Abstract Lukács and András posed the problem of showing the existence of a set of n 2 points in the interior of a convex n-gon
More informationGeometric Modeling in Graphics
Geometric Modeling in Graphics Part 10: Surface reconstruction Martin Samuelčík www.sccg.sk/~samuelcik samuelcik@sccg.sk Curve, surface reconstruction Finding compact connected orientable 2-manifold surface
More informationLecture 0: Reivew of some basic material
Lecture 0: Reivew of some basic material September 12, 2018 1 Background material on the homotopy category We begin with the topological category TOP, whose objects are topological spaces and whose morphisms
More informationImproved Results on Geometric Hitting Set Problems
Improved Results on Geometric Hitting Set Problems Nabil H. Mustafa nabil@lums.edu.pk Saurabh Ray saurabh@cs.uni-sb.de Abstract We consider the problem of computing minimum geometric hitting sets in which,
More informationApproximating Polygonal Objects by Deformable Smooth Surfaces
Approximating Polygonal Objects by Deformable Smooth Surfaces Ho-lun Cheng and Tony Tan School of Computing, National University of Singapore hcheng,tantony@comp.nus.edu.sg Abstract. We propose a method
More informationA Polynomial Time Algorithm for Multivariate Interpolation in Arbitrary Dimension via the Delaunay Triangulation
A Polynomial Time Algorithm for Multivariate Interpolation in Arbitrary Dimension via the Delaunay Triangulation Tyler Chang, Layne Watson, Thomas Lux, Bo Li, Li Xu, Ali Butt, Kirk Cameron, and Yili Hong
More informationParametric Curves. University of Texas at Austin CS384G - Computer Graphics Fall 2010 Don Fussell
Parametric Curves University of Texas at Austin CS384G - Computer Graphics Fall 2010 Don Fussell Parametric Representations 3 basic representation strategies: Explicit: y = mx + b Implicit: ax + by + c
More informationLecture 3 Randomized Algorithms
Lecture 3 Randomized Algorithms Jean-Daniel Boissonnat Winter School on Computational Geometry and Topology University of Nice Sophia Antipolis January 23-27, 2017 Computational Geometry and Topology Randomized
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 informationShellings, the Euler-Poincaré Formula for Polytopes, Dehn-Sommerville Equations, the Upper Bound Theorem
Chapter 8 Shellings, the Euler-Poincaré Formula for Polytopes, Dehn-Sommerville Equations, the Upper Bound Theorem 8.1 Shellings The notion of shellability is motivated by the desire to give an inductive
More informationParametric Curves. University of Texas at Austin CS384G - Computer Graphics
Parametric Curves University of Texas at Austin CS384G - Computer Graphics Fall 2010 Don Fussell Parametric Representations 3 basic representation strategies: Explicit: y = mx + b Implicit: ax + by + c
More informationScientific Computing WS 2018/2019. Lecture 12. Jürgen Fuhrmann Lecture 12 Slide 1
Scientific Computing WS 2018/2019 Lecture 12 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 12 Slide 1 Recap For more discussion of mesh generation, see J.R. Shewchuk: Lecture Notes on Delaunay
More informationSurfaces Beyond Classification
Chapter XII Surfaces Beyond Classification In most of the textbooks which present topological classification of compact surfaces the classification is the top result. However the topology of 2- manifolds
More informationComputation of Voronoi Diagrams and Delaunay Triangulation via Parametric Linear Programming
Computation of Voronoi Diagrams and Delaunay Triangulation via Parametric Linear Programming Saša V. Raković, Pascal Grieder and Colin Jones September 14, 2004 Abstract This note illustrates how Voronoi
More informationBezier Curves, B-Splines, NURBS
Bezier Curves, B-Splines, NURBS Example Application: Font Design and Display Curved objects are everywhere There is always need for: mathematical fidelity high precision artistic freedom and flexibility
More informationThe Borsuk-Ulam theorem- A Combinatorial Proof
The Borsuk-Ulam theorem- A Combinatorial Proof Shreejit Bandyopadhyay April 14, 2015 1 Introduction The Borsuk-Ulam theorem is perhaps among the results in algebraic topology having the greatest number
More informationDelaunay Triangulations
Delaunay Triangulations (slides mostly by Glenn Eguchi) Motivation: Terrains Set of data points A R 2 Height ƒ(p) defined at each point p in A How can we most naturally approximate height of points not
More informationLecture 25: Bezier Subdivision. And he took unto him all these, and divided them in the midst, and laid each piece one against another: Genesis 15:10
Lecture 25: Bezier Subdivision And he took unto him all these, and divided them in the midst, and laid each piece one against another: Genesis 15:10 1. Divide and Conquer If we are going to build useful
More informationSimplicial Complexes: Second Lecture
Simplicial Complexes: Second Lecture 4 Nov, 2010 1 Overview Today we have two main goals: Prove that every continuous map between triangulable spaces can be approximated by a simplicial map. To do this,
More informationPTAS for geometric hitting set problems via Local Search
PTAS for geometric hitting set problems via Local Search Nabil H. Mustafa nabil@lums.edu.pk Saurabh Ray saurabh@cs.uni-sb.de Abstract We consider the problem of computing minimum geometric hitting sets
More informationLecture 1 Discrete Geometric Structures
Lecture 1 Discrete Geometric Structures Jean-Daniel Boissonnat Winter School on Computational Geometry and Topology University of Nice Sophia Antipolis January 23-27, 2017 Computational Geometry and Topology
More informationOutline of the presentation
Surface Reconstruction Petra Surynková Charles University in Prague Faculty of Mathematics and Physics petra.surynkova@mff.cuni.cz Outline of the presentation My work up to now Surfaces of Building Practice
More informationconsisting of compact sets. A spline subdivision scheme generates from such
Spline Subdivision Schemes for Compact Sets with Metric Averages Nira Dyn and Elza Farkhi Abstract. To dene spline subdivision schemes for general compact sets, we use the representation of spline subdivision
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 information9 Bounds for the Knapsack Problem (March 6)
9 Bounds for the Knapsack Problem (March 6) In this lecture, I ll develop both upper and lower bounds in the linear decision tree model for the following version of the (NP-complete) Knapsack 1 problem:
More informationComputational Geometry
Motivation Motivation Polygons and visibility Visibility in polygons Triangulation Proof of the Art gallery theorem Two points in a simple polygon can see each other if their connecting line segment is
More informationIntegral Geometry and the Polynomial Hirsch Conjecture
Integral Geometry and the Polynomial Hirsch Conjecture Jonathan Kelner, MIT Partially based on joint work with Daniel Spielman Introduction n A lot of recent work on Polynomial Hirsch Conjecture has focused
More informationand the crooked shall be made straight, and the rough ways shall be made smooth; Luke 3:5
ecture 8: Knot Insertion Algorithms for B-Spline Curves and Surfaces and the crooked shall be made straight, and the rough ways shall be made smooth; uke 3:5. Motivation B-spline methods have several advantages
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 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 informationTiling Three-Dimensional Space with Simplices. Shankar Krishnan AT&T Labs - Research
Tiling Three-Dimensional Space with Simplices Shankar Krishnan AT&T Labs - Research What is a Tiling? Partition of an infinite space into pieces having a finite number of distinct shapes usually Euclidean
More informationBasics of Combinatorial Topology
Chapter 7 Basics of Combinatorial Topology 7.1 Simplicial and Polyhedral Complexes In order to study and manipulate complex shapes it is convenient to discretize these shapes and to view them as the union
More informationCOMPUTATIONAL GEOMETRY
Thursday, September 20, 2007 (Ming C. Lin) Review on Computational Geometry & Collision Detection for Convex Polytopes COMPUTATIONAL GEOMETRY (Refer to O'Rourke's and Dutch textbook ) 1. Extreme Points
More informationIntro to Curves Week 1, Lecture 2
CS 536 Computer Graphics Intro to Curves Week 1, Lecture 2 David Breen, William Regli and Maxim Peysakhov Department of Computer Science Drexel University Outline Math review Introduction to 2D curves
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 informationGeometric structures on manifolds
CHAPTER 3 Geometric structures on manifolds In this chapter, we give our first examples of hyperbolic manifolds, combining ideas from the previous two chapters. 3.1. Geometric structures 3.1.1. Introductory
More informationVoronoi Diagrams and Delaunay Triangulation slides by Andy Mirzaian (a subset of the original slides are used here)
Voronoi Diagrams and Delaunay Triangulation slides by Andy Mirzaian (a subset of the original slides are used here) Voronoi Diagram & Delaunay Triangualtion Algorithms Divide-&-Conquer Plane Sweep Lifting
More informationComputational Geometry Lecture Delaunay Triangulation
Computational Geometry Lecture Delaunay Triangulation INSTITUTE FOR THEORETICAL INFORMATICS FACULTY OF INFORMATICS 7.12.2015 1 Modelling a Terrain Sample points p = (x p, y p, z p ) Projection π(p) = (p
More informationA Grid-Based Approximation Algorithm for the Minimum Weight Triangulation Problem
A Grid-Based Approximation Algorithm for the Minimum Weight Triangulation Problem arxiv:1706.03263v1 [cs.cg] 10 Jun 2017 Sharath Raghvendra, Mariëtte C. Wessels Abstract Given a set of n points on a plane,
More informationHowever, this is not always true! For example, this fails if both A and B are closed and unbounded (find an example).
98 CHAPTER 3. PROPERTIES OF CONVEX SETS: A GLIMPSE 3.2 Separation Theorems It seems intuitively rather obvious that if A and B are two nonempty disjoint convex sets in A 2, then there is a line, H, separating
More informationCS 536 Computer Graphics Intro to Curves Week 1, Lecture 2
CS 536 Computer Graphics Intro to Curves Week 1, Lecture 2 David Breen, William Regli and Maxim Peysakhov Department of Computer Science Drexel University 1 Outline Math review Introduction to 2D curves
More informationLecture notes for Topology MMA100
Lecture notes for Topology MMA100 J A S, S-11 1 Simplicial Complexes 1.1 Affine independence A collection of points v 0, v 1,..., v n in some Euclidean space R N are affinely independent if the (affine
More informationThe Game of Hex and Brouwer s Fixed-Point Theorem. Michael Cao
The Game of Hex and Brouwer s Fixed-Point Theorem Michael Cao June, 2017 Contents 1 Abstract 2 2 The Game of Hex 3 2.1 Introduction.............................. 3 2.2 Setup and Rules...........................
More information2D Geometry. Pierre Alliez Inria Sophia Antipolis
2D Geometry Pierre Alliez Inria Sophia Antipolis Outline Sample problems Polygons Graphs Convex hull Voronoi diagram Delaunay triangulation Sample Problems Line Segment Intersection Theorem: Segments (p
More informationbe a polytope. has such a representation iff it contains the origin in its interior. For a generic, sort the inequalities so that
( Shelling (Bruggesser-Mani 1971) and Ranking Let be a polytope. has such a representation iff it contains the origin in its interior. For a generic, sort the inequalities so that. a ranking of vertices
More informationThe Graphs of Triangulations of Polygons
The Graphs of Triangulations of Polygons Matthew O Meara Research Experience for Undergraduates Summer 006 Basic Considerations Let Γ(n) be the graph with vertices being the labeled planar triangulation
More informationLecture 2 - Introduction to Polytopes
Lecture 2 - Introduction to Polytopes Optimization and Approximation - ENS M1 Nicolas Bousquet 1 Reminder of Linear Algebra definitions Let x 1,..., x m be points in R n and λ 1,..., λ m be real numbers.
More informationSimplicial Global Optimization
Simplicial Global Optimization Julius Žilinskas Vilnius University, Lithuania September, 7 http://web.vu.lt/mii/j.zilinskas Global optimization Find f = min x A f (x) and x A, f (x ) = f, where A R n.
More informationVoronoi Diagrams and Delaunay Triangulations. O Rourke, Chapter 5
Voronoi Diagrams and Delaunay Triangulations O Rourke, Chapter 5 Outline Preliminaries Properties and Applications Computing the Delaunay Triangulation Preliminaries Given a function f: R 2 R, the tangent
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 informationGeometric structures on 2-orbifolds
Geometric structures on 2-orbifolds Section 1: Manifolds and differentiable structures S. Choi Department of Mathematical Science KAIST, Daejeon, South Korea 2010 Fall, Lectures at KAIST S. Choi (KAIST)
More informationDirichlet Voronoi Diagrams and Delaunay Triangulations
Chapter 9 Dirichlet Voronoi Diagrams and Delaunay Triangulations 9.1 Dirichlet Voronoi Diagrams In this chapter we present very briefly the concepts of a Voronoi diagram and of a Delaunay triangulation.
More informationINTRODUCTION TO FINITE ELEMENT METHODS
INTRODUCTION TO FINITE ELEMENT METHODS LONG CHEN Finite element methods are based on the variational formulation of partial differential equations which only need to compute the gradient of a function.
More informationCAD & Computational Geometry Course plan
Course plan Introduction Segment-Segment intersections Polygon Triangulation Intro to Voronoï Diagrams & Geometric Search Sweeping algorithm for Voronoï Diagrams 1 Voronoi Diagrams Voronoi Diagrams or
More informationPlanarity. 1 Introduction. 2 Topological Results
Planarity 1 Introduction A notion of drawing a graph in the plane has led to some of the most deep results in graph theory. Vaguely speaking by a drawing or embedding of a graph G in the plane we mean
More informationVoronoi Diagram. Xiao-Ming Fu
Voronoi Diagram Xiao-Ming Fu Outlines Introduction Post Office Problem Voronoi Diagram Duality: Delaunay triangulation Centroidal Voronoi tessellations (CVT) Definition Applications Algorithms Outlines
More informationVoronoi Diagrams, Delaunay Triangulations and Polytopes
Voronoi Diagrams, Delaunay Triangulations and Polytopes Jean-Daniel Boissonnat MPRI, Lecture 2 Computational Geometry Learning Voronoi, Delaunay & Polytopes MPRI, Lecture 2 1 / 43 Voronoi diagrams in nature
More informationarxiv: v1 [cs.cg] 7 Oct 2017
A Proof of the Orbit Conjecture for Flipping Edge-Labelled Triangulations Anna Lubiw 1, Zuzana Masárová 2, and Uli Wagner 2 arxiv:1710.02741v1 [cs.cg] 7 Oct 2017 1 School of Computer Science, University
More informationElementary Combinatorial Topology
Elementary Combinatorial Topology Frédéric Meunier Université Paris Est, CERMICS, Ecole des Ponts Paristech, 6-8 avenue Blaise Pascal, 77455 Marne-la-Vallée Cedex E-mail address: frederic.meunier@cermics.enpc.fr
More informationProperties of Blending Functions
Chapter 5 Properties of Blending Functions We have just studied how the Bernstein polynomials serve very nicely as blending functions. We have noted that a degree n Bézier curve always begins at P 0 and
More informationDelaunay Triangulations. Presented by Glenn Eguchi Computational Geometry October 11, 2001
Delaunay Triangulations Presented by Glenn Eguchi 6.838 Computational Geometry October 11, 2001 Motivation: Terrains Set of data points A R 2 Height ƒ(p) defined at each point p in A How can we most naturally
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 informationLecture 17: Continuous Functions
Lecture 17: Continuous Functions 1 Continuous Functions Let (X, T X ) and (Y, T Y ) be topological spaces. Definition 1.1 (Continuous Function). A function f : X Y is said to be continuous if the inverse
More informationPart 3. Topological Manifolds
Part 3 Topological Manifolds This part is devoted to study of the most important topological spaces, the spaces which provide a scene for most of geometric branches in mathematics such as Differential
More informationThe Cyclic Cycle Complex of a Surface
The Cyclic Cycle Complex of a Surface Allen Hatcher A recent paper [BBM] by Bestvina, Bux, and Margalit contains a construction of a cell complex that gives a combinatorial model for the collection of
More information6. Concluding Remarks
[8] K. J. Supowit, The relative neighborhood graph with an application to minimum spanning trees, Tech. Rept., Department of Computer Science, University of Illinois, Urbana-Champaign, August 1980, also
More informationIntroduction to Voronoi Diagrams and Delaunay Triangulations
Introduction to Voronoi Diagrams and Delaunay Triangulations Solomon Boulos Introduction to Voronoi Diagrams and Delaunay Triangulations p.1 Voronoi Diagrams Voronoi region: V (p i ) = {x R n p i x p j
More informationLecture 2 September 3
EE 381V: Large Scale Optimization Fall 2012 Lecture 2 September 3 Lecturer: Caramanis & Sanghavi Scribe: Hongbo Si, Qiaoyang Ye 2.1 Overview of the last Lecture The focus of the last lecture was to give
More informationA PROOF OF THE LOWER BOUND CONJECTURE FOR CONVEX POLYTOPES
PACIFIC JOURNAL OF MATHEMATICS Vol. 46, No. 2, 1973 A PROOF OF THE LOWER BOUND CONJECTURE FOR CONVEX POLYTOPES DAVID BARNETTE A d polytope is defined to be a cz-dimensional set that is the convex hull
More informationFigure 1: A positive crossing
Notes on Link Universality. Rich Schwartz: In his 1991 paper, Ramsey Theorems for Knots, Links, and Spatial Graphs, Seiya Negami proved a beautiful theorem about linearly embedded complete graphs. These
More informationLifting Transform, Voronoi, Delaunay, Convex Hulls
Lifting Transform, Voronoi, Delaunay, Convex Hulls Subhash Suri Department of Computer Science University of California Santa Barbara, CA 93106 1 Lifting Transform (A combination of Pless notes and my
More informationRigidity of ball-polyhedra via truncated Voronoi and Delaunay complexes
!000111! NNNiiinnnttthhh IIInnnttteeerrrnnnaaatttiiiooonnnaaalll SSSyyymmmpppooosssiiiuuummm ooonnn VVVooorrrooonnnoooiii DDDiiiaaagggrrraaammmsss iiinnn SSSccciiieeennnccceee aaannnddd EEEnnngggiiinnneeeeeerrriiinnnggg
More informationCLASSIFICATION OF ORDERABLE AND DEFORMABLE COMPACT COXETER POLYHEDRA IN HYPERBOLIC SPACE
CLASSIFICATION OF ORDERABLE AND DEFORMABLE COMPACT COXETER POLYHEDRA IN HYPERBOLIC SPACE DHRUBAJIT CHOUDHURY, SUHYOUNG CHOI, AND GYE-SEON LEE Abstract. The aim of this work is to investigate properties
More informationReflection groups 4. Mike Davis. May 19, Sao Paulo
Reflection groups 4 Mike Davis Sao Paulo May 19, 2014 https://people.math.osu.edu/davis.12/slides.html 1 2 Exotic fundamental gps Nonsmoothable aspherical manifolds 3 Let (W, S) be a Coxeter system. S
More informationCS130 : Computer Graphics Curves (cont.) Tamar Shinar Computer Science & Engineering UC Riverside
CS130 : Computer Graphics Curves (cont.) Tamar Shinar Computer Science & Engineering UC Riverside Blending Functions Blending functions are more convenient basis than monomial basis canonical form (monomial
More informationSperner and Tucker s lemma
Sperner and Tucker s lemma Frédéric Meunier February 8th 2010 Ecole des Ponts, France Sperner s lemma Brouwer s theorem Let B be a d-dimensional ball and f be a continuous B B map. Then always Theorem
More informationPS Geometric Modeling Homework Assignment Sheet I (Due 20-Oct-2017)
Homework Assignment Sheet I (Due 20-Oct-2017) Assignment 1 Let n N and A be a finite set of cardinality n = A. By definition, a permutation of A is a bijective function from A to A. Prove that there exist
More informationClassifications in Low Dimensions
Chapter XI Classifications in Low Dimensions In different geometric subjects there are different ideas which dimensions are low and which high. In topology of manifolds low dimension means at most 4. However,
More informationarxiv: v1 [cs.cg] 8 Jan 2018
Voronoi Diagrams for a Moderate-Sized Point-Set in a Simple Polygon Eunjin Oh Hee-Kap Ahn arxiv:1801.02292v1 [cs.cg] 8 Jan 2018 Abstract Given a set of sites in a simple polygon, a geodesic Voronoi diagram
More informationSPERNER S LEMMA MOOR XU
SPERNER S LEMMA MOOR XU Abstract. Is it possible to dissect a square into an odd number of triangles of equal area? This question was first answered by Paul Monsky in 970, and the solution requires elements
More informationRemark. Jacobs University Visualization and Computer Graphics Lab : ESM4A - Numerical Methods 331
Remark Reconsidering the motivating example, we observe that the derivatives are typically not given by the problem specification. However, they can be estimated in a pre-processing step. A good estimate
More informationLinear Programming in Small Dimensions
Linear Programming in Small Dimensions Lekcija 7 sergio.cabello@fmf.uni-lj.si FMF Univerza v Ljubljani Edited from slides by Antoine Vigneron Outline linear programming, motivation and definition one dimensional
More informationTopological Data Analysis - I. Afra Zomorodian Department of Computer Science Dartmouth College
Topological Data Analysis - I Afra Zomorodian Department of Computer Science Dartmouth College September 3, 2007 1 Acquisition Vision: Images (2D) GIS: Terrains (3D) Graphics: Surfaces (3D) Medicine: MRI
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 informationFACES OF CONVEX SETS
FACES OF CONVEX SETS VERA ROSHCHINA Abstract. We remind the basic definitions of faces of convex sets and their basic properties. For more details see the classic references [1, 2] and [4] for polytopes.
More informationFinite Element Methods
Chapter 5 Finite Element Methods 5.1 Finite Element Spaces Remark 5.1 Mesh cells, faces, edges, vertices. A mesh cell is a compact polyhedron in R d, d {2,3}, whose interior is not empty. The boundary
More information