DOCSLIB.ORG

Centroidal Voronoi Tessellations: Applications and Algorithms Author(S): Qiang Du, Vance Faber and Max Gunzburger Source: SIAM Review, Vol

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Centroidal Voronoi Tessellations: Applications and Algorithms Author(S): Qiang Du, Vance Faber and Max Gunzburger Source: SIAM Review, Vol Centroidal Voronoi Tessellations: Applications and Algorithms Author(s): Qiang Du, Vance Faber and Max Gunzburger Source: SIAM Review, Vol. 41, No. 4 (Dec., 1999), pp. 637-676 Published by: Society for Industrial and Applied Mathematics Stable URL: http://www.jstor.org/stable/2653198 Accessed: 09-08-2016 14:52 UTC REFERENCES Linked references are available on JSTOR for this article: http://www.jstor.org/stable/2653198?seq=1&cid=pdf-reference#references_tab_contents You may need to log in to JSTOR to access the linked references. Your use of the JSTOR archive indicates your acceptance of the Terms & Conditions of Use, available at http://about.jstor.org/terms JSTOR is a not-for-profit service that helps scholars, researchers, and students discover, use, and build upon a wide range of content in a trusted digital archive. We use information technology and tools to increase productivity and facilitate new forms of scholarship. For more information about JSTOR, please contact [email protected]. Society for Industrial and Applied Mathematics is collaborating with JSTOR to digitize, preserve and extend access to SIAM Review This content downloaded from 144.174.6.137 on Tue, 09 Aug 2016 14:52:00 UTC All use subject to http://about.jstor.org/terms SIAM REVIEW (c) 1999 Society for Industrial and Applied Mathematics Vol. 41, No. 4, pp. 637-676 Centroidal Voronoi Tessellations: Applications and Algorithms* Qiang Dut Vance Fabert Max Gunzburger? Abstract. A centroidal Voronoi tessellation is a Voronoi tessellation whose generating points are the centroids (centers of mass) of the corresponding Voronoi regions. We give some applica- tions of such tessellations to problems in image compression, quadrature, finite difference methods, distribution of resources, cellular biology, statistics, and the territorial behavior of animals. We discuss methods for computing these tessellations, provide some analyses concerning both the tessellations and the methods for their determination, and, finally, present the results of some numerical experiments. Key words. Voronoi tessellations, centroids, vector quantization, data compression, clustering AMS subject classifications. 5202, 52B55, 62H30, 6502, 65D30, 65U05, 65Y25, 68U05, 68U10 Pll. S0036144599352836 1. Introduction. Given an open set Q C RN, the set {1}k is called a tessel- lation of Q if vi n Vj = 0 for i V& j and U1=,Vi = Q. Let denote the Euclidean norm on RN. Given a set of poinits {Zi}i=1 belonging to Q, the Voronoi region Vi corresponding to the point zi is defined by (1.1)~ V V= {x EQ I Ix - z1 < x - zl for j=l,...,)k, j V-i I The points {z%}f are called generators. The set {VJ}k= is a Voronoi tessellation or Voronoi diagram of Q, and each Vi is referred to as the Voronoi region corresponding to Zi. (Depending on the application, there exist many different names for Voronoi re- gions, including Dirichlet regions, area of influence polygons, Meijering cells, Thiessen polygons, and S-mosaics.) The Voronoi regions are polyhedra. These tessellations, and their dual tessellations (in R2, the Delaunay triangulations), are very useful in a variety of applications. For a comprehensive treatment, see [47]. *Received by the editors December 15, 1998; accepted for publication (in revised form) March 3, 1999; published electronically October 20, 1999. http://www.siam.org/journals/sirev/41-4/35283.html tDepartment of Mathematics, Iowa State University, Ames, IA 50011-2064, and Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong ([email protected]). The research of this author was supported in part by National Science Foundation grant DMS- 9796208 and by grants from HKRGC and HKUST. tCommunications and Computing Division, Los Alamos National Laboratory, Los Alamos, NM 87545. Current address: Lizardtech Inc., 1520 Bellevue Avenue, Seattle, WA 98122 (vxf@lizardtech. com). The research of this author was supported in part by the Department of Energy. ?Department of Mathematics, Iowa State University, Ames, IA 50011-2064 (gunzburg@iastate. edu). The research of this author was supported in part by National Science Foundation grant DMS-9806358. 637 This content downloaded from 144.174.6.137 on Tue, 09 Aug 2016 14:52:00 UTC All use subject to http://about.jstor.org/terms 638 QIANG DU, VANCE FABER, AND MAX GUNZBURGER I 0 0~~~~ Fig. 1.1 On the left, the Voronoi regions corresponding to 10 randomly selected points in a square; the density function is a constant. The dots are the Voronoi generators and the circles are the centroids of the corresponding Voronoi regions. Note that the generators and the centroids do not coincide. On the right, a 10-point centroidal Voronoi tessellation. The dots are simultaneously the generators for the Voronoi tessellation and the centroids of the Voronoi regions. Given a region V C RN and a density function p, defined in V, the mass centroid z* of V is defined by /YP(y) dy (1.2) z* J /8p(y) dy Given k points zi, i = 1, ... , k, we can define their associated Voronoi regions Vi, i = 1, ... , k. On the other hand, given the regions Vi, i = 1, . , k, we can define their mass centroids zi, i = 1, ..., k. Here, we are interested in the situation where (1.3) zi = Z* i-=1, ...,)k, i.e., the points zi that serve as generators for the Voronoi regions Vi are themselves the mass centroids of those regions. We call such a tessellation a centroidal Voronoi tessellation. This situation is quite special since, in general, arbitrarily chosen points in RJN are not the centroids of their associated Voronoi regions. See Figure 1.1 for an illustration in two dimensions. One may ask, How does one find centroidal Voronoi tessellations, and are they of any use? In this paper, we review some answers to these questions. Concerning the first question, and to be more precise, consider the following problem: Given a region Q C RjtN a positive integer k, and a density function p, defined for y in Q, find k points zi E Q and k regions Vi that tessellate Q such that simultaneously for each i, This content downloaded from 144.174.6.137 on Tue, 09 Aug 2016 14:52:00 UTC All use subject to http://about.jstor.org/terms CENTROIDAL VORONOI TESSELLATIONS 639 Fig. 1.2 Two centroidal Voronoi tessellations of a square. The points zl and Z2 are the centroids of the rectangles on the left or of the triangles on the right. Vi is the Voronoi region for zi and zi is the mass centroid of Vi. The solution of this problem is in general not unique. For example, consider the case of N = 2, Q C JR2 a square, and p 1. Two solutions are depicted in Figure 1.2; others may be obtained through rotation. Another example is provided by the three regular tessellations of JR2 into squares, triangles, and hexagons. The plan of the remainder of the paper is as follows. First, in the remainder of section 1, we give some remarks about various generalizations of centroidal Voronoi diagrams. In section 2, we consider applications. These are drawn from the worlds of data compression, numerical analysis, biology, statistics, and operations research. In section 3, we discuss properties of these tessellations and their relation to the critical points of an energy (or cost or error) functional. In section 4, we present probabilistic approaches for computing these tessellations. In section 5, we present some deterministic methods. In section 6, we discuss some theoretical issues related to one of these methods. Some numerical experiments are reported in section 7, and, finally, concluding remarks are given in section 8. 1. 1. Centroidal Voronoi Tessellations in Other Metrics. We first consider the case where, instead of a region Q, we are given a discrete set of points W y{yi}=L in RN. A set {V}k 1 is a tessellation of W if vi n Vj 0 for i f j and Uk=1V% = W. Given a set of points {Z%}>kl belonging to RN, Voronoi sets are now defined by (1.4) V?={xeW I Ix-zil< x-zjl for j=1,...,k, j# i, where equality holds only for i < j}. Other tie-breaking rules for points equidistant to two or more generators can also be used. Given a density function p defined in W, the mass centroid z* of a set V c W is now defined by (1.5) E P(Y)jy_Z*j2= inf E p(Y)jy_-z2, yEV yEV where the sums extend over the points belonging to V, and V* can be taken to be V or it can be a larger set like RN. In the statistical and vector quantization literature (see, e.g., [15, 19, 34]) discrete centroidal Voronoi tessellations are often related to optimal k-means clusters and Voronoi regions and centroids are referred to as clusters and cluster centers, respectively. More discussion of this topic is provided in section 2.3. This content downloaded from 144.174.6.137 on Tue, 09 Aug 2016 14:52:00 UTC All use subject to http://about.jstor.org/terms 640 QIANG DU, VANCE FABER, AND MAX GUNZBURGER The notions of Voronoi regions and centroids, and therefore of centroidal Voronoi regions, may be generalized to more abstract spaces and to metrics other than the Euclidean f2-norm. For example, (1.6) d(x, y) = llx - yllefp will be of interest to us. For another example, let d(x, y) denote a distance function induced by a norm that is equivalent to the 12 norm on RRN. Let -< be an order of RN. Then, Voronoi tessellations of Q, with respect to this metric, are defined by (1.7) Vi = x E Q I d(x, zi) < d(x, zj), j = 1,.
Load more

Recommended publications
  • Guannan Zhang's CV
    Guannan Zhang, Ph.D. Senior Research Staff Joint-Faculty Professor Computer Science and Mathematics Division Department of Mathematics and Statistics Oak Ridge National Laboratory Auburn University Oak Ridge, TN 37831-6211 Auburn, AL 36849 Tel: (865) 241-4503, Fax: (865) 241-0381 Email: [email protected] Email: [email protected] Education 2009 – 2012 Ph.D. in Computational Science, Florida State University, Tallahassee, Florida Thesis advisor: Prof. Max D. Gunzburger 2009 – 2011 M.S. in Computational Science, Florida State University, Tallahassee, Florida Thesis advisor: Prof. Max D. Gunzburger 2007 – 2009 M.S. in Mathematics, Shandong University, Jinan, Shandong, China Thesis advisors: Prof. Weidong Zhao and Prof. Shige Peng 2003 – 2007 B.S. in Mathematics, Shandong University, Jinan, Shandong, China Professional experience 2020 – Senior Research Staff, Computer Science and Mathematics Division, Oak Ridge National Labo- ratory, Oak Ridge, Tennessee 2015 – 2020 Research Staff, Computer Science and Mathematics Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 2012 – 2015 Research Associate, Computer Science and Mathematics Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 2012 – 2014 Distinguished Staff Fellow (Householder Fellow), Computer Science and Mathematics Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 2014 – Joint-Faculty Professor, Department of Mathematics and Statistics, Auburn University, Auburn, Alabama Research interests High-dimensional function approximation Uncertainty quantification and related applications Scientific machine learning and artificial intelligence Stochastic differential equations and stochastic optimization Numerical solution of parameterized differential equations 1/9 Grants & constracts Current funding support 2019 – 2021 DOE - Advanced Scientific Computing Research, A Stochastic Optimal Control Framework For Quantifying And Reducing Uncertainties In Deep Learning: $800,000, Principal Investigator.
    [Show full text]
  • PUBLICATIONS of MAX GUNZBURGER Books
    PUBLICATIONS OF MAX GUNZBURGER Books 1. Finite Element Methods for Viscous Incompressible Flows: A Guide to Theory, Practice and Algorithms; Academic, Boston, 1989. 2. Perspectives in Flow Control and Optimization; SIAM, Philadelphia, 2003. 3. Least Squares Finite Element Methods; Springer, Berlin, 2009; with P. Bochev. 4. An Algorithmic Introduction to Numerical Methods for Uncertainty Quantification for PDEs with Random Inputs; to be published by Springer, Berlin, 2016; with C. Webster. Edited books 5. Incompressible Computational Fluid Dynamics: Trends and Advances; Cambridge, Cambridge, 1993; with R. Nico- laides. 6. Optimal Control and Design; Birkha¨user, Boston, 1995; with J. Borggaard, J. Burkardt, and J. Peterson. 7. Flow Control; Springer, New York, 1995. Journal articles 1. Diffraction of shock waves by a moving thin wing; J. Fluid Mech 42, 1970, 585-608; with L. Ting. 2. Motion of decaying vortex rings with non-similar vorticity distributions; J. Engrg. Math. 6, 1972, 53-61. 3. Antisymmetric problems in shock diffraction theory; Z. Angew. Math. Phys. 24, 1973, 83-104; with V. Yang. 4. Long time behavior of a decaying vortex; Z. Angew. Math. Mech. 53, 1973, 751-760. 5. Mathematical formulation for the propagation of sound through a turbulent jet; J. Engrg. Math. 10, 1976, 219-230; with C. Liu, L. Maestrello, and L. Ting. 6. Frequency modulation at a moving material interface and a conservation law for wave number; J. Sound Vib. 48, 1976, 169-178; with G. Kleinstein. 7. On the mathematical conditions for the existence of periodic fluctuations in nonuniform media; J. Sound Vib. 48, 1976, 345-357; with G.
    [Show full text]
  • Curriculum Vitae 1
    Xiaoqiang Wang’s Curriculum Vitae 1 Curriculum Vitae by Xiaoqiang Wang School of Computational Science & Voice: (850)644-5792 Department of Mathematics Fax: (850)644-0098 Florida State University Email: [email protected] Dirac Science Library http://www.scs.fsu.edu/˜xwang Tallahassee, FL 32306-4120 RESEARCH INTERESTS Numerical analysis and applied partial differential equations Mathematical biology Image processing, scientific visualization and data mining High-performance scientific computing EDUCATION Ph.D in Applied Mathematics, Pennsylvania State University, University Park, 2005 Dissertation: Phase field models and simulations of vesicle bio-membranes Advisor: Qiang Du M.S. in Institute of Mathematics, Chinese Academy of Sciences, P.R.China, 1998 B.S. in Mathematics, Wuhan University, Wuhan, Hubei, P.R.China, 1995 ACADEMIC EXPERIENCE Assistant Professor (2006 – present), Department of Mathematics and School of Computational Science, Florida State University, Tallahassee, FL ¾ Instructor of MAC2311 Calculus with Analytic Geometry I (4), Fall 2006 IMA Industrial Postdoctoral Associate (2005 – 2), Institute for Mathematics and its Applications (IMA), University of Minnesota ¾ Mathematical Biology: cell membrane related research ¾ Human Brain Imaging Project carried out by the International Neuroimaging Consortium (INC) Research Assistant (2003 – 2005), Pennsylvania State University, University Park ¾ Numerical simulation of vesicle membranes via Phase Field approaches ¾ Vector fields visualization and segmentation ¾ Tessellation
    [Show full text]
  • Notices of the American Mathematical Society ABCD Springer.Com
    ISSN 0002-9920 Notices of the American Mathematical Society ABCD springer.com Highlights in Springer’s eBook Collection of the American Mathematical Society September 2009 Volume 56, Number 8 Bôcher, Osgood, ND and the Ascendance NEW NEW 2 EDITION EDITION of American Mathematics forms bridges between From the reviews of the first edition The theory of elliptic curves is Mathematics at knowledge, tradition, and contemporary 7 Chorin and Hald provide excellent distinguished by the diversity of the Harvard life. The continuous development and explanations with considerable insight methods used in its study. This book growth of its many branches permeates and deep mathematical understanding. treats the arithmetic theory of elliptic page 916 all aspects of applied science and 7 SIAM Review curves in its modern formulation, technology, and so has a vital impact on through the use of basic algebraic our society. The book will focus on these 2nd ed. 2009. X, 162 p. 7 illus. number theory and algebraic geometry. aspects and will benefit from the (Surveys and Tutorials in the Applied contribution of world-famous scientists. Mathematical Sciences) Softcover 2nd ed. 2009. XVIII, 514 p. 14 illus. Invariant Theory ISBN 978-1-4419-1001-1 (Graduate Texts in Mathematics, 2009. XI, 263 p. (Modeling, Simulation & 7 Approx. $39.95 Volume 106) Hardcover of Tensor Product Applications, Volume 3) Hardcover ISBN 978-0-387-09493-9 7 $59.95 Decompositions of ISBN 978-88-470-1121-2 7 $59.95 U(N) and Generalized Casimir Operators For access check with your librarian page 931 Stochastic Partial Linear Optimization Mathematica in Action Riverside Meeting Differential Equations The Simplex Workbook The Power of Visualization A Modeling, White Noise Approach G.
    [Show full text]
  • Positions Education Awards and Honors
    MAX D. GUNZBURGER Francis Eppes Eminent Professor and Chair Department of Scientific Computing Florida State University Tallahassee FL 32306-4120 850-644-7060 850-644-0098 (fax) [email protected] http://www.sc.fsu.edu/egunzburg Positions Full-time positions Florida State University - Frances Eppes Eminent Professor of Scientific Computing and Founding Chair, Department of Scientific Computing, 2008-present; Director, School of Computational Science, 2005-2008; Frances Eppes Eminent Professor of Mathematics, 2002-2008 Iowa State University - Distinguished Professor of Mathematics, 2001-2003; Professor and Chair of Mathematics, 1995-2001 Virginia Polytechnic Institute and State University - Professor of Mathematics, 1987-1997 Los Alamos National Laboratories - Visiting Scientist, 1988-1989 Carnegie Mellon University - Professor of Mathematics, 1981-1989 University of Tennessee - Associate Professor and Professor of Mathematics, 1976-1982 Institute for Computer Applications in Science and Engineering - Staff Scientist, 1973-1976 Naval Ordnance Laboratory - National Research Council Postdoctorate Associate, 1971-1973 New York University - Research Scientist and Assistant Professor of Mathematics, 1969-1971 Adjunct, courtesy, and part-time positions Florida State University - Professor of Mathematics, 2008-present; Professor of Mechanical Engineering, 2006-present Peking University - Guest Professor of Computational Mathematics, 2000-2002 Virginia Polytechnic Institute and State University - Professor of Computer Sciences, 1992-1997; Professor of
    [Show full text]
  • Final Program
    Final Program The SIAM Conference on Control and Its Applications is sponsored by the SIAM Activity Group on Control and Systems Theory (SIAG/CST) The SIAM Activity Group on Control and Systems Theory fosters collaboration and interaction among mathematicians, engineers, and other scientists in those areas of research related to the theory of systems and their control. It seeks to promote the development of theory and methods related to modeling, control, estimation, and approximation of complex biological, physical, and engineering systems. The SIAG organizes a biennial conference, sponsors minisymposia at SIAM meetings and periodic conferences, and maintains a member directory and an electronic discussion group. Every two years, the activity group also awards the SIAG/Control and Systems Theory Prize to a young researcher for outstanding research contributions to mathematical control or systems theory and the SIAG/CST Best SICON Paper Prize to the author(s) of the two most outstanding papers published in the SIAM Journal on Control and Optimization (SICON). Society for Industrial and Applied Mathematics 3600 Market Street, 6th Floor Philadelphia, PA 19104-2688 USA Telephone: +1-215-382-9800 Fax: +1-215- 386-7999 SIAM 2013 Mobile App Conference E-mail: [email protected] Scan the QR code to the right with any QR reader and download the TripBuilder ® EventMobile app Conference Web: www.siam.org/meetings/ to your iPHONE, iPAD, iTOUCH, ANDROID or Membership and Customer Service: (800) 447-7426 (US & BLACKBERRY mobile device. You can also visit Canada) or +1-215-382-9800 (worldwide) www.tripbuilder.net/mobileweb/apps/siam2013 2 2013 SIAM Annual Meeting • SIAM Conference on Control & Its Applications General Information Table of Contents William M.
    [Show full text]
  • Notices of the American Mathematical Society ABCD Springer.Com
    ISSN 0002-9920 Notices of the American Mathematical Society ABCD springer.com Highlights in Springer’s eBook Collection of the American Mathematical Society June/July 2009 Volume 56, Number 6 NEW NEW Using representation theory and This self-contained text off ers a host of Editor-in-Chief: Steven G. Krantz invariant theory to analyze the new mathematical tools and strategies The Journal of Geometric Analysis is a symmetries arising from group actions, which develop a connection between high-quality journal devoted to and with emphasis on the geometry analysis and other mathematical publishing important new results at the and basic theory of Lie groups and Lie disciplines, such as physics and interface of analysis, geometry and algebras, this book reworks an earlier engineering. A broad view of math- partial diff erential equations. Founded highly acclaimed work by the author. ematics is presented throughout; the 17 years ago by its current Editor-in- This comprehensive introduction to Lie text is excellent for the classroom or Chief the journal has maintained theory, representation theory, invariant self-study. standards of innovation and excellence. theory and algebraic groups is more 2009. XX, 452 p. 10 illus. Softcover accessible to students and includes a ISSN 1050-6926 (print ) ISBN 978-0-387-77378-0 7 $79.95 broader ranger of applications. ISSN 1559-002X (online) 2009. XX, 716 p. 10 illus. (Graduate Texts in Mathematics, Volume 255) Hardcover ISBN 978-0-387-79851-6 7 $69.95 For access check with your librarian Fractional Diff erentiation Risk and Asset Allocation Groups and Symmetries Inequalities A.
    [Show full text]
  • Mathematical Association of America Southeastern Section 85Th Annual Meeting
    Mathematical Association of America Southeastern Section 85th Annual Meeting Society for Industrial and Applied Mathematics Southeastern Atlantic Section 30th Annual Meeting 31 March – 1 April 2006 Auburn University Auburn, Alabama 1 MAA-SE Officers Chair, Charlotte Knotts-Zides, Wofford College Vice Chair, Virginia Parks, Georgia Perimeter College-Clarkston Program Chair, Tom Barr, Rhodes College Secretary-Treasurer, Martha Abell, Georgia Southern University Section Governor, Ben Klein, Davidson College Section Lecturer (2005-2006), David Anderson, University of Tennessee Newsletter Editor, Virgina Watson, Kennesaw State University State Directors: Bruce Atkinson, Samford University (Alabama), Curtis Herink, Mercer University (Georgia), Rob Harger (High Point University (North Carolina), John Harris, Furman University (South Carolina), Jeff Knisley, East Tennessee State University (Tennessee) SIAM-SEAS Officers President, Debra Knisley, East Tennessee State University Secretary-Treasurer and President-elect, Amnon J. Meir, Auburn University 2 3 Joint Conference Mathematical Association of America, Southeastern Section http://www.maa.org/southeastern Eighty-fifthAnnual Meeting Society for Industrial and Applied Mathematics, Southeastern Atlantic Section http://www.mathcs.citadel.edu/~chenm/SEASarch.html Thirtieth Annual Meeting Auburn University Auburn, Alabama 31 March – 1 April 2006 Sessions will be held in Parker Hall, 371 W. Roosevelt Drive and the Sciences Center Auditorium (adjacent to Parker Hall). The Friday evening reception and
    [Show full text]
  • Xiaoming He Curriculum Vitae
    Xiaoming He Curriculum Vitae Department of Mathematics and Statistics Phone: 573-341-6208 Missouri University of Science and Technology Email: [email protected] Rolla, MO, 65409 Web: http://web.mst.edu/ hex/ Education Ph.D. in Mathematics; Virginia Tech (08/05|05/09). Advisor: Prof. Tao Lin. M. S. in Mathematics; Virginia Tech (08/05|05/07). Advisor: Prof. Tao Lin. M. S. in Computational Mathematics; Sichuan University (08/02|06/05). Advisor: Prof. Tao L¨u. B. S. in Computational Mathematics; Sichuan University (08/98|06/02). Thesis Advisor: Prof. Xiaoping Xie. Positions Professor; Department of Mathematics and Statistics, Missouri University of Science and Tech- nology (09/21{present). Faculty Fellow; Office of Vice Chancellor of Research, Missouri University of Science and Technology (07/21{present). Associate Professor; Department of Mathematics and Statistics, Missouri University of Science and Technology (09/16{08/21). Assistant Professor; Department of Mathematics and Statistics, Missouri University of Science and Technology (09/10{08/16). Postdoctoral research associate (supervised by Prof. Max Gunzburger); Department of Scien- tific Computing, Florida State University (05/09|08/10). Research Interests • Interface problems, multi-phase problems, computational fluid dynamics, data assimilation, stochastic PDEs, boundary integral equations, feedback control • Finite element methods, domain decomposition methods, lattice Boltzmann methods, extrap- olations • Computational plasma physics (Particle-In-Cell method) 1 Awards and Fellowships • Faculty Research Award, Missouri University of Science and Technology, 2019. • Affordable and Open Educational Resources Initiative Award, University of Missouri, 2019. • Humboldt Research Fellowship for Experienced Researchers, Alexander von Humboldt Foun- dation, 2018. • Faculty Research Award, Missouri University of Science and Technology, 2014.
    [Show full text]