Voronoi diagram explained

In mathematics, a Voronoi diagram is a partition of a plane into regions close to each of a given set of objects. It can be classified also as a tessellation. In the simplest case, these objects are just finitely many points in the plane (called seeds, sites, or generators). For each seed there is a corresponding region, called a Voronoi cell, consisting of all points of the plane closer to that seed than to any other. The Voronoi diagram of a set of points is dual to that set's Delaunay triangulation.

The Voronoi diagram is named after mathematician Georgy Voronoy, and is also called a Voronoi tessellation, a Voronoi decomposition, a Voronoi partition, or a Dirichlet tessellation (after Peter Gustav Lejeune Dirichlet). Voronoi cells are also known as Thiessen polygons, after Alfred H. Thiessen.[1] [2] [3] Voronoi diagrams have practical and theoretical applications in many fields, mainly in science and technology, but also in visual art.[4] [5]

The simplest case

In the simplest case, shown in the first picture, we are given a finite set of points

\{p1,...pn\}

in the Euclidean plane. In this case each site

pk

is one of these given points, and its corresponding Voronoi cell

Rk

consists of every point in the Euclidean plane for which

pk

is the nearest site: the distance to

pk

is less than or equal to the minimum distance to any other site

pj

. For one other site

pj

, the points that are closer to

pk

than to

pj

, or equally distant, form a closed half-space, whose boundary is the perpendicular bisector of line segment

pjpk

. Cell

Rk

is the intersection of all of these

n-1

half-spaces, and hence it is a convex polygon.[6] When two cells in the Voronoi diagram share a boundary, it is a line segment, ray, or line, consisting of all the points in the plane that are equidistant to their two nearest sites. The vertices of the diagram, where three or more of these boundaries meet, are the points that have three or more equally distant nearest sites.

Formal definition

Let X be a metric space with distance function d. Let K be a set of indices and let (P_k)_ be a tuple (indexed collection) of nonempty subsets (the sites) in the space X. The Voronoi cell, or Voronoi region, R_k, associated with the site P_k is the set of all points in X whose distance to P_k is not greater than their distance to the other sites P_j, where j is any index different from k. In other words, if d(x,\, A) = \inf\ denotes the distance between the point x and the subset A, then

R_k = \

The Voronoi diagram is simply the tuple of cells (R_k)_ . In principle, some of the sites can intersect and even coincide (an application is described below for sites representing shops), but usually they are assumed to be disjoint. In addition, infinitely many sites are allowed in the definition (this setting has applications in geometry of numbers and crystallography), but again, in many cases only finitely many sites are considered.

In the particular case where the space is a finite-dimensional Euclidean space, each site is a point, there are finitely many points and all of them are different, then the Voronoi cells are convex polytopes and they can be represented in a combinatorial way using their vertices, sides, two-dimensional faces, etc. Sometimes the induced combinatorial structure is referred to as the Voronoi diagram. In general however, the Voronoi cells may not be convex or even connected.

In the usual Euclidean space, we can rewrite the formal definition in usual terms. Each Voronoi polygon R_k is associated with a generator point P_k.Let X be the set of all points in the Euclidean space. Let P_1 be a point that generates its Voronoi region R_1, P_2 that generates R_2, and P_3 that generates R_3, and so on. Then, as expressed by Tran et al,[7] "all locations in the Voronoi polygon are closer to the generator point of that polygon than any other generator point in the Voronoi diagram in Euclidean plane".

Illustration

As a simple illustration, consider a group of shops in a city. Suppose we want to estimate the number of customers of a given shop. With all else being equal (price, products, quality of service, etc.), it is reasonable to assume that customers choose their preferred shop simply by distance considerations: they will go to the shop located nearest to them. In this case the Voronoi cell

Rk

of a given shop

Pk

can be used for giving a rough estimate on the number of potential customers going to this shop (which is modeled by a point in our city).

For most cities, the distance between points can be measured using the familiarEuclidean distance

\ell2=d\left[\left(a1,a2\right),\left(b1,b2\right)\right]=\sqrt{\left(a1-

2
b
1\right)

+\left(a2-

2}
b
2\right)

or the Manhattan distance:

d\left[\left(a1,a2\right),\left(b1,b2\right)\right]=\left|a1-b1\right|+\left|a2-b2\right|

.

The corresponding Voronoi diagrams look different for different distance metrics.

Properties

History and research

Informal use of Voronoi diagrams can be traced back to Descartes in 1644.[10] Peter Gustav Lejeune Dirichlet used two-dimensional and three-dimensional Voronoi diagrams in his study of quadratic forms in 1850.British physician John Snow used a Voronoi-like diagram in 1854 to illustrate how the majority of people who died in the Broad Street cholera outbreak lived closer to the infected Broad Street pump than to any other water pump.

Voronoi diagrams are named after Georgy Feodosievych Voronoy who defined and studied the general n-dimensional case in 1908.[11] Voronoi diagrams that are used in geophysics and meteorology to analyse spatially distributed data are called Thiessen polygons after American meteorologist Alfred H. Thiessen, who used them to estimate rainfall from scattered measurements in 1911. Other equivalent names for this concept (or particular important cases of it): Voronoi polyhedra, Voronoi polygons, domain(s) of influence, Voronoi decomposition, Voronoi tessellation(s), Dirichlet tessellation(s).

Examples

Voronoi tessellations of regular lattices of points in two or three dimensions give rise to many familiar tessellations.

Certain body-centered tetragonal lattices give a tessellation of space with rhombo-hexagonal dodecahedra.

For the set of points (xy) with x in a discrete set X and y in a discrete set Y, we get rectangular tiles with the points not necessarily at their centers.

Higher-order Voronoi diagrams

Although a normal Voronoi cell is defined as the set of points closest to a single point in S, an nth-order Voronoi cell is defined as the set of points having a particular set of n points in S as its n nearest neighbors. Higher-order Voronoi diagrams also subdivide space.

Higher-order Voronoi diagrams can be generated recursively. To generate the nth-order Voronoi diagram from set S, start with the (n − 1)th-order diagram and replace each cell generated by X =  with a Voronoi diagram generated on the set S − X.

Farthest-point Voronoi diagram

For a set of n points the (n − 1)th-order Voronoi diagram is called a farthest-point Voronoi diagram.

For a given set of points S =  the farthest-point Voronoi diagram divides the plane into cells in which the same point of P is the farthest point. A point of P has a cell in the farthest-point Voronoi diagram if and only if it is a vertex of the convex hull of P. Let H =  be the convex hull of P; then the farthest-point Voronoi diagram is a subdivision of the plane into k cells, one for each point in H, with the property that a point q lies in the cell corresponding to a site hi if and only if d(q, hi) > d(q, pj) for each pj ∈ S with hipj, where d(p, q) is the Euclidean distance between two points p and q.[12] [13]

The boundaries of the cells in the farthest-point Voronoi diagram have the structure of a topological tree, with infinite rays as its leaves. Every finite tree is isomorphic to the tree formed in this way from a farthest-point Voronoi diagram.[14]

Generalizations and variations

As implied by the definition, Voronoi cells can be defined for metrics other than Euclidean, such as the Mahalanobis distance or Manhattan distance. However, in these cases the boundaries of the Voronoi cells may be more complicated than in the Euclidean case, since the equidistant locus for two points may fail to be subspace of codimension 1, even in the two-dimensional case.

A weighted Voronoi diagram is the one in which the function of a pair of points to define a Voronoi cell is a distance function modified by multiplicative or additive weights assigned to generator points. In contrast to the case of Voronoi cells defined using a distance which is a metric, in this case some of the Voronoi cells may be empty. A power diagram is a type of Voronoi diagram defined from a set of circles using the power distance; it can also be thought of as a weighted Voronoi diagram in which a weight defined from the radius of each circle is added to the squared Euclidean distance from the circle's center.[15]

The Voronoi diagram of

n

points in

d

-dimensional space can have O(n^) vertices, requiring the same bound for the amount of memory needed to store an explicit description of it. Therefore, Voronoi diagrams are often not feasible for moderate or high dimensions. A more space-efficient alternative is to use approximate Voronoi diagrams.[16]

Voronoi diagrams are also related to other geometric structures such as the medial axis (which has found applications in image segmentation, optical character recognition, and other computational applications), straight skeleton, and zone diagrams.

Applications

Meteorology/Hydrology

It is used in meteorology and engineering hydrology to find the weights for precipitation data of stations over an area (watershed). The points generating the polygons are the various station that record precipitation data. Perpendicular bisectors are drawn to the line joining any two stations. This results in the formation of polygons around the stations. The area

(Ai)

touching station point is known as influence area of the station. The average precipitation is calculated by the formula
\bar{P}=\sumAiPi
\sumAi

Humanities and social sciences

Natural sciences

Health

Engineering

Mathematics

Informatics

Civics and planning

Bakery

Algorithms

Several efficient algorithms are known for constructing Voronoi diagrams, either directly (as the diagram itself) or indirectly by starting with a Delaunay triangulation and then obtaining its dual.Direct algorithms include Fortune's algorithm, an O(n log(n)) algorithm for generating a Voronoi diagram from a set of points in a plane.Bowyer–Watson algorithm, an O(n log(n)) to O(n2) algorithm for generating a Delaunay triangulation in any number of dimensions, can be used in an indirect algorithm for the Voronoi diagram. The Jump Flooding Algorithm can generate approximate Voronoi diagrams in constant time and is suited for use on commodity graphics hardware.[42] [43]

Lloyd's algorithm and its generalization via the Linde–Buzo–Gray algorithm (aka k-means clustering), use the construction of Voronoi diagrams as a subroutine.These methods alternate between steps in which one constructs the Voronoi diagram for a set of seed points, and steps in which the seed points are moved to new locations that are more central within their cells. These methods can be used in spaces of arbitrary dimension to iteratively converge towards a specialized form of the Voronoi diagram, called a Centroidal Voronoi tessellation, where the sites have been moved to points that are also the geometric centers of their cells.

Voronoi in 3D

Voronoi meshes can also be generated in 3D.

See also

References

External links

Notes and References

  1. Book: Peter A. . Burrough . Rachael . McDonnell . Rachael A. . McDonnell . Christopher D. . Lloyd . Principles of Geographical Information Systems . 8.11 Nearest neighbours: Thiessen (Dirichlet/Voroni) polygons . https://books.google.com/books?id=kvoJCAAAQBAJ&pg=PA160 . 2015 . Oxford University Press . 978-0-19-874284-5 . 160–.
  2. Book: Paul A. . Longley . Michael F. . Goodchild . David J. . Maguire . David W. . Rhind . Geographic Information Systems and Science . 14.4.4.1 Thiessen polygons . https://books.google.com/books?id=-FbVI-2tSuYC&pg=PA333 . 2005 . Wiley . 978-0-470-87001-3 . 333–.
  3. Book: Sen, Zekai . Spatial Modeling Principles in Earth Sciences . 2.8.1 Delaney, Varoni, and Thiessen Polygons . https://books.google.com/books?id=6N0yDQAAQBAJ&pg=PA57 . 2016 . Springer . 978-3-319-41758-5 . 57–.
  4. Franz . Aurenhammer . Franz Aurenhammer . 1991 . Voronoi Diagrams – A Survey of a Fundamental Geometric Data Structure . ACM Computing Surveys . 23 . 3 . 345–405 . 10.1145/116873.116880. 4613674 .
  5. Book: Atsuyuki . Okabe . Barry . Boots . Kokichi . Sugihara . Sung Nok . Chiu . 2000 . Spatial Tessellations – Concepts and Applications of Voronoi Diagrams . 2nd . John Wiley . 978-0-471-98635-5.
  6. Book: Boyd . Stephen . Vandenberghe . Lieven . Convex Optimization . 2004 . Cambridge University Press . Exercise 2.9 . 60.
  7. Book: Q. T. . Tran . D. . Tainar . M. . Safar . 2009 . Transactions on Large-Scale Data- and Knowledge-Centered Systems . 357 . Springer . 9783642037214.
  8. .
  9. .
  10. Senechal . Marjorie . 1993-05-21 . Mathematical Structures: Spatial Tessellations . Concepts and Applications of Voronoi Diagrams. Atsuyuki Okabe, Barry Boots, and Kokichi Sugihara. Wiley, New York, 1992. xii, 532 pp., illus. $89.95. Wiley Series in Probability and Mathematical Statistics. . Science . en . 260 . 5111 . 1170–1173 . 10.1126/science.260.5111.1170 . 17806355 . 0036-8075.
  11. and .
  12. Book: 2008 . Computational Geometry . 978-3-540-77974-2 . . Third . Mark . de Berg . Marc . van Kreveld . Mark . Overmars . Otfried . Schwarzkopf . Mark de Berg . Marc van Kreveld . Mark Overmars . Otfried Schwarzkopf. 7.4 Farthest-Point Voronoi Diagrams. Includes a description of the algorithm.
  13. Sven . Skyum . A simple algorithm for computing the smallest enclosing circle . Information Processing Letters . 37 . 3 . 18 February 1991 . 121–125 . 10.1016/0020-0190(91)90030-L., contains a simple algorithm to compute the farthest-point Voronoi diagram.
  14. Biedl . Therese . Therese Biedl . Grimm . Carsten . Palios . Leonidas . Shewchuk . Jonathan . Jonathan Shewchuk . Verdonschot . Sander . Realizing farthest-point Voronoi diagrams . Proceedings of the 28th Canadian Conference on Computational Geometry (CCCG 2016) . 2016.
  15. Book: Edelsbrunner, Herbert . Herbert Edelsbrunner . 13.6 Power Diagrams . 327–328 . Springer-Verlag . EATCS Monographs on Theoretical Computer Science . Algorithms in Combinatorial Geometry . 10 . 1987 . 2012 . 9783642615689.
  16. Book: Space-efficient approximate Voronoi diagrams . 1581134959 . 10.1145/509907.510011 . Proceedings of the thiry-fourth annual ACM symposium on Theory of computing . Sunil . Sunil Arya . Theocharis . Malamatos . David M. . Mount . David Mount . 2002 . 721–730 . 1727373 .
  17. Book: Laver . Michael . Sergenti . Ernest . Party competition : an agent-based model . 2012 . Princeton University Press . Princeton . 978-0-691-13903-6.
  18. Bock . Martin . Tyagi . Amit Kumar . Kreft . Jan-Ulrich . Alt . Wolfgang . Generalized Voronoi Tessellation as a Model of Two-dimensional Cell Tissue Dynamics . Bulletin of Mathematical Biology . 72 . 7 . 1696–1731 . 10.1007/s11538-009-9498-3 . 20082148 . 2009 . 0901.4469v1 . 2009arXiv0901.4469B . 16074264.
  19. Hui Li . Robert . Sitnik . Atilla M . Baskurt . Spatial Modeling of Bone Microarchitecture . Three-Dimensional Image Processing (3Dip) and Applications II . 8290 . 82900P . 2012 . 2012SPIE.8290E..0PL . 10.1117/12.907371 . 1505014.
  20. Sanchez-Gutierrez . D. . Tozluoglu . M. . Barry . J. D. . Pascual . A. . Mao . Y. . Escudero . L. M. . 2016-01-04 . Fundamental physical cellular constraints drive self-organization of tissues . The EMBO Journal . 35 . 1 . 77–88 . 10.15252/embj.201592374 . 4718000 . 26598531.
  21. Book: Feinstein . Joseph . Bionoi: A Voronoi Diagram-Based Representation of Ligand-Binding Sites in Proteins for Machine Learning Applications . 2021 . Protein-Ligand Interactions and Drug Design . 299–312 . Ballante . Flavio . Methods in Molecular Biology . New York, NY . Springer US . en . 10.1007/978-1-0716-1209-5_17 . 978-1-0716-1209-5 . 2021-04-23 . Shi . Wentao . Ramanujam . J. . Brylinski . Michal . 2266 . 33759134 . 232338911.
  22. E pur si muove: Galilean-invariant cosmological hydrodynamical simulations on a moving mesh . Springel . Volker . 2010 . MNRAS . 401 . 2 . 791–851 . 10.1111/j.1365-2966.2009.15715.x . free . 2010MNRAS.401..791S . 0901.4107 . 119241866.
  23. Kasim . Muhammad Firmansyah . 2017-01-01 . Quantitative shadowgraphy and proton radiography for large intensity modulations . Physical Review E . 95 . 2 . 023306 . 10.1103/PhysRevE.95.023306 . 28297858 . 1607.04179 . 2017PhRvE..95b3306K . 13326345.
  24. Book: Steven Johnson . The Ghost Map: The Story of London's Most Terrifying Epidemic — and How It Changed Science, Cities, and the Modern World . 19 October 2006 . 16 October 2017 . Penguin Publishing Group . 978-1-101-15853-1 . 187.
  25. Mulheran . P. A. . Blackman . J. A. . Capture zones and scaling in homogeneous thin-film growth . Physical Review B . 53 . 15 . 1996 . 10261–7 . 10.1103/PhysRevB.53.10261 . 9982595 . 1996PhRvB..5310261M.
  26. Pimpinelli . Alberto . Tumbek . Levent . Winkler . Adolf . Scaling and Exponent Equalities in Island Nucleation: Novel Results and Application to Organic Films . The Journal of Physical Chemistry Letters . 5 . 6 . 2014 . 995–8 . 10.1021/jz500282t . 24660052 . 3962253.
  27. Fanfoni . M. . Placidi . E. . Arciprete . F. . Orsini . E. . Patella . F. . Balzarotti . A. . Sudden nucleation versus scale invariance of InAs quantum dots on GaAs . Physical Review B . 75 . 24 . 245312 . 2007 . 1098-0121 . 10.1103/PhysRevB.75.245312 . 2007PhRvB..75x5312F . 120017577.
  28. Miyamoto . Satoru . Moutanabbir . Oussama . Haller . Eugene E. . Itoh . Kohei M. . Spatial correlation of self-assembled isotopically pure Ge/Si(001) nanoislands . Physical Review B . 79 . 165415 . 165415 . 2009 . 1098-0121 . 10.1103/PhysRevB.79.165415 . 2009PhRvB..79p5415M . 13719907 .
  29. Löbl . Matthias C. . Zhai . Liang . Jahn . Jan-Philipp . Ritzmann . Julian . Huo . Yongheng . Wieck . Andreas D. . Schmidt . Oliver G. . Ludwig . Arne . Rastelli . Armando . Warburton . Richard J. . Correlations between optical properties and Voronoi-cell area of quantum dots . Physical Review B . 100 . 15 . 2019-10-03 . 2469-9950 . 10.1103/physrevb.100.155402 . 155402 . 1902.10145 . 2019PhRvB.100o5402L . 119443529.
  30. Web site: GOLD COAST CULTURAL PRECINCT. ARM Architecture. 2014-04-28. 2016-07-07. https://web.archive.org/web/20160707155535/http://www.a-r-m.com.au/projects_GoldCoastCP.html. dead.
  31. Lopez . C. . Zhao . C.-L. . Magniol . S . Chiabaut . N . Leclercq . L . 28 February 2019 . Microscopic Simulation of Cruising for Parking of Trucks as a Measure to Manage Freight Loading Zone . Sustainability . 11 (5), 1276.
  32. Singh . K. . Sadeghi . F. . Correns . M. . Blass . T. . December 2019 . A microstructure based approach to model effects of surface roughness on tensile fatigue . International Journal of Fatigue . 129 . 105229 . 10.1016/j.ijfatigue.2019.105229 . 202213370.
  33. Niu . Hanlin . Savvaris . Al . Tsourdos . Antonios . Ji . Ze . Voronoi-visibility roadmap-based path planning algorithm for unmanned surface vehicles . The Journal of Navigation . 2019 . 72 . 4 . 850–874 . 10.1017/S0373463318001005. 67908628 .
  34. Cortes . J. . Martinez . S. . Karatas . T. . Bullo . F. . April 2004 . Coverage control for mobile sensing networks . IEEE Transactions on Robotics and Automation . 20 . 2 . 243–255 . 10.1109/TRA.2004.824698 . 2022860 . 2374-958X.
  35. Teruel . Enrique . Aragues . Rosario . López-Nicolás . Gonzalo . April 2021. A Practical Method to Cover Evenly a Dynamic Region With a Swarm . IEEE Robotics and Automation Letters . 6 . 2 . 1359–1366 . 10.1109/LRA.2021.3057568 . 232071627 . 2377-3766.
  36. Pólya, G. On the zeros of the derivatives of a function and its analytic character. Bulletinof the AMS, Volume 49, Issue 3, 178-191, 1943.
  37. Book: Mitchell, Tom M. . Machine Learning . limited . 1997 . McGraw-Hill . International . 978-0-07-042807-2 . 233.
  38. Web site: Shenwai. Tanushree. 2021-11-18. A Novel Deep Learning Technique That Rebuilds Global Fields Without Using Organized Sensor Data. 2021-12-05 . MarkTechPost . en-US.
  39. Archived at Ghostarchive and the Wayback Machine: Web site: Mark DiMarco: User Interface Algorithms [JSConf2014]]. 11 June 2014 . www.youtube.com.
  40. Web site: Find my School . Victorian Government Department of Education . en . 2023-07-25.
  41. Haridy . Rich . 2017-09-06 . Architect turned cake-maker serves up mouth-watering geometric 3D-printed cakes . New Atlas .
  42. Rong . Guodong . Tan . Tiow Seng . Olano . Marc . Séquin . Carlo H. . Jump flooding in GPU with applications to Voronoi diagram and distance transform . https://www.comp.nus.edu.sg/~tants/jfa/i3d06.pdf . 10.1145/1111411.1111431 . 109–116 . ACM . Proceedings of the 2006 Symposium on Interactive 3D Graphics, SI3D 2006, March 14-17, 2006, Redwood City, California, USA . 2006. 1-59593-295-X .
  43. Web site: Shadertoy.