Dobrinat
Dobrinat
Dobrinat
VORONOI DIAGRAMS
ADAM DOBRIN
1. Introduction
where they can be seen and how they are applied. Sections 3 and
First, it should be noted that for any positive integer n, there are
and all of the points in one region are closer to the corresponding site
than to any other site. Where there is not one closest point, there is a
has a boat, with each boat capable of going the same speed. Let every
point in R that can be reached from the boat from island x before any
A REVIEW OF PROPERTIES AND VARIATIONS OF VORONOI DIAGRAMS 3
both within and outside the math world. Voronoi diagrams can be used
ogy) and neolithic clans and tribes (anthropology and archaelogy), and
First, we shall denote the location of a point pi as (xi1 , xi2 ), and the
V (pi ) = {~x ||~x − ~xi || ≤ ||~x − ~xj || ∀j 3 i 6= j}
Sn
the Voronoi diagram of P . Different notation for V is i=1 Vi .
P is the set of generators. Every point on the plane that is not a vertex
e(pi , pj ) has the property that ||~x − ~xi || = ||~x − ~xj ||. An edge can be
denoted as ei , where i is an index for the edges and does not have to
our edges from 1 to n, n being the total number of edges, going top
of Vi .
(or more) nearest generator points on the plane. Vertices are denoted
qi , and are the endpoints of edges. The number of edges that meet at
region is not connected, and in both diagrams, many of the regions are
not convex. Note that the non-convex regions are boundary regions.
diagram on R2 generated by P .
we can see how a Voronoi region is created using the method of finding
T
j∈Z+ ≤n H(pi , pj ). Using this method for all of the points in a Voronoi
8 ADAM DOBRIN
(i) A Voronoi edge e(pi , pj )(6= ∅) is a line segment iff the line connecting
points within its boundary. For each vertex qi ∈ V(P ), there exists a
points. Note that the largest empty circle represented in Figure 4 has
note are the correlations between basic graph theory and Voronoi di-
We find that
(1) nv − ne + n = 1.
ne ≤ 3n − 6 and nv ≤ 2n − 5
1
nv ≥ (n − nc ) + 1 and ne ≤ 3nv − nc − 3
2
10 ADAM DOBRIN
Another interesting theorem is that for any Voronoi diagram, the av-
erage number of Voronoi edges per Voronoi polygon does not exceed 6.
2(3n−6)
More accurately, this number is less than or equal to n
.
[7].
line edges into line segments by capping them with a vertex. All new
vertices added are added to nv . Simultaneously, all but one of the un-
4. Delaunay Tessellations
non-collinearity. This means that for our generator set P , the points in
P are not all on the same line. Given a Voronoi diagram V(P ), join all
and only if all vertices of V(P ) are non-degenerate. See Figure 6 for an
(solid dots).
edges and Delaunay edges are orthogonal. With a little thought, this
diagrams and their Delaunay tessellations are not only dual graphs,
diagrams, but have not yet discussed the meaning and results of this
(i) Qd = P
(ii) Cd = Q
14 ADAM DOBRIN
(iii)|Ed | = |E|
The points on a Delaunay circle, which are Delaunay vertices, are called
natural neighbors. Two vertices are natural neighbors if and only if they
have assumed that our generator points, besides their location, have
tor points can be more useful than having uniformly weighted points in
cable when looking at, for example, the population size of a settlement,
Recall from our definition of the Voronoi diagram from Section 3.4.
Let
\
Vw (pi ) = Domw (pi , pj ).
pj ∈P \{pi }
boat. Therefore, a faster boat will reach boats previously outside of its
region.
||~x − ~xi ||
dmw (p, pi ) = , where wi > 0.
wi
There are many names for the associated Voronoi diagram: Vmw , the
A REVIEW OF PROPERTIES AND VARIATIONS OF VORONOI DIAGRAMS17
not be using the last two, as they are highly specialized and this paper
diagrams.
Voronoi regions are not smaller than the weight of pi . Another way
may share disconnected edges. Edges in Vmw are circular arcs if and
if only if the weights of the two affective regions are not equal. Edges
in Vmw are straight lines if and only if the weights of the two affective
regions are equal. See Figure 9 for a diagram of the bisectors between
on ∂ CH(Pmax ).
islands. They all still travel at the same speed, but some boats begin
given by
given by
Domaw (pi , pj ) = {~x ||~x − ~xi || − ||~x − ~xj || ≤ wi − wj }, where i 6= j.
has positive area, then there exists at least one non-convex additively-
shaped with respect to its generator point. This means that from pi ,
we can draw a line to any point in the region Vaw (pi ), and the line will
section, we will look at two of these: the Plane-Sweep and the Tree
ing the Voronoi diagram, we draw all generator points, Voronoi vertices,
22 ADAM DOBRIN
given some initial assumptions. Therefore, we can copy our Voronoi di-
Let pi and pj be two generator points such that x(pi ) > x(pj ) where
ine that there is a current, moving directly to the right, that is faster
A REVIEW OF PROPERTIES AND VARIATIONS OF VORONOI DIAGRAMS23
point to the left of its island. Voronoi regions in φ(R2 ) are contained
For the plane-sweep method, we will need to make four initial as-
tor points align vertically. Lastly, any generator point and any of that
In this method, we take a vertical line and move it from left to right
over the plane. We update our data everytime there is an event, which
the line, from bottom to top. We will also use Q, a set of points in
the plane. It is a list of all possible points where events may occur.
Voronoi vertices φ(q) are added to Q as they are found by the sweepline.
24 ADAM DOBRIN
φ(V).
(1) Let Q = P .
4.b.iii.
(φ(V (pj )), h− (pi , pj ), φ(V (pi )), h+ (pi , pj ), φ(V (pj ))).
perbola on L.
and 4.c.iv.
(i) Replace the subsequence (h± (pi , pj ), φ(V (pj )), h± (pj , pk ))
h± (pj , pk ), and h.
(5) Report all half-hyperbolas that were ever listed on L, all the
From this algorithm, we have all edges, vertices, and generator points
d
(2) r(p) = .
2 cos(tan−1 (m))
Our point,
d
(4) p = (xφ − , yφ ).
2 cos(tan−1 (m))
task, involving much more work than is actually necessary. For simpler
information. Let qijk denote the vertex incident to V (pi ), V (pj ), and
V (pk ) in that order. For some point p = (x, y), we let H(pi , pj , pk , p) =
1 xi yi x2i+ yi2
2 2
1 xj yj xj + y j
(5) H(pi , pj , pk , pl ) = .
2 2
1 xk yk xk + y k
1 x y x2 + y 2
28 ADAM DOBRIN
(a) (b)
(c) (d)
also assume that Vl−1 divides the plane into i + 1 regions; call them
cells. Assume that every cell, besides infinite ones, is simply connected
We will let T be the set of vertices of Vl−1 that will be in V (pl ). As-
mized.
(2) Among the vertices qijk on ∂V (pi ), find the one that gives the
assumptions).
(4) For every generator pi associated with a vertex qijk that satis-
qikl . In the case that qijk lies on the outer circuit of Vl−1 , draw
30 ADAM DOBRIN
e(pi , pk ) be the vertex qikl . If the bisector does not intersect ei-
(5) Remove all vertices in T , and the edges incident to them, from
Vl−1 .
In Figure 14, we see the use of the Tree Expansion and Deletion
generator point affects only the central vertex. For the other three ver-
(a) (b)
(c)
are out. During the algorithm, vertices that are as of yet unchecked are
planning.
domain, will spit sand away from the center of its territory. For a high
of sand results in raised bars of sand, and, when viewed from above,
tessellation [3].
tessellation [3].
P and the Voronoi diagram V(P ). The same is not true for centroidal
population density, but did not explain why this was important. To
diagram has the property that each generator point is the centroid of
its Voronoi region. This implies that the problems that we will run
into when trying to find a CVT will lie solely in the generator set.
diagrams arises. The freedom that we have with regular Voronoi di-
distances given any set of data. CVT’s are limited in this sense, but, as
to find a generator set whose points are the mass centroids of their
respective Voronoi regions. There are many methods for doing so: the
A REVIEW OF PROPERTIES AND VARIATIONS OF VORONOI DIAGRAMS35
two that we will look at are Lloyd’s method and MacQueen’s method.
then move them incrementally such that the resulting Voronoi diagram
is also a CVT.
able random numbers and observing only the fraction of the numbers
use Lloyd or MacQueen’s method, and get a very similar, and for all
intents and purposes, the same result. The Monte Carlo method, in
effect, gives us a head start on creating the generator set for the CVT.
It gives us a set of n points that resembles our density function, ρ(x, y).
Because Lloyd and Macqueen’s methods take the n points and iterate
them to find a generator set for a CVT, which also resembles ρ, the
set found using a Monte Carlo method may be closer to our final set.
Carlo method.
36 ADAM DOBRIN
method as in Step 1.
(4) Set
ip∗ + y
p∗ =
i+1
and i = i + 1.
return to Step 2.
are very similar, but also very different. They differ just in their method
A REVIEW OF PROPERTIES AND VARIATIONS OF VORONOI DIAGRAMS37
are simply moving a point, then checking our results with the “prede-
one point at a time, instead of many points in the set. Lloyd’s method
of criteria, but only use one for any given construction of a CVT. One
criterion will be the total number of iterations, i. The higher the num-
(a) (b)
(c) (d)
is very rare that it would happen multiple times in a row. The lower
the error counter is, the closer V will be to an actual centroidal Voronoi
tessellation.
take note of the fact that we do not get CVT’s with either Lloyd or
(a) (b)
(c) (d)
using these methods close enough to true CVT’s because the differences
between them are minimal. The only way to guarantee that the Voronoi
This cannot be done for all practical purposes. Thus, we accept the
7.10. Figures 16 and 17. In both Figures 16 and 17, (a) and (b)
are computed using an iteration limit, while (c) and (d) are computed
using an error counter. In Figure 16, note the difference between the
two diagrams on the left and the two on the right. In (a) and (c), we
see Voronoi regions bunched together that are relatively far away from
the points (−1, −1) and (1, 1). Some of these problems are corrected in
(b) and (d), as we can observe a better order of smaller regions closer
to those points than in (a) and (c). In Figure 17, we can see the same
general trends as in Figure 16. The smaller regions are closer to the
origin in (b) and (d) than in (a) and (c). These results are expected
because (b) uses a higher iteration limit than does (a), and (d) uses a
lower error counter than does (c). These results follow directly from
8. Conclusion
grams. There are many topics on Voronoi diagrams that were not cov-
diagram, to find the locations of the generator set P [6]. Also fasci-
rithm for which mail centers deliver to which location, and in what
properties that are not included in this paper. For example, there
are many other metrics to use when computing weighted Voronoi di-
weighted diagrams are fairly complex, and are either part of a fourth-
the square of the Euclidean distance between the generator and the
lines. These two metrics of Voronoi diagrams are very useful in model-
ing real-world occurrences, as are the two covered in the paper, because
42 ADAM DOBRIN
system.
A REVIEW OF PROPERTIES AND VARIATIONS OF VORONOI DIAGRAMS43
References
[2] Geometry in Action. Drysdale, Scot. 14 Nov 1996. Donald Bren School of
//www.ics.uci.edu/~eppstein/gina/scot.drysdale.html.
[3] Du, Qiang, Vance Faber, and Max Gunzburger. “Centroidal Voronoi Tessel-
313-322.
[5] Ju, Lili, Qiang Du, and Max Gunzburger. “Probabilistic Methods for Cen-
Computer 28 (2002):1477-1500.
[6] Okabe, Atsuyuki, Barry Boots, Kokichi Sugihara, and Sung Nok Chiu. Spa-
1996.