sbornÃk
sbornÃk sbornÃk
Petr Kahánek, Alexej Kolcun 2 Triangulations Let S be a finite set of points in the Euclidean plane. A triangulation of S is a maximally connected straight line plane graph whose vertices are the points of S. By maximality, each face is a triangle except for the exterior face which is the complement of the convex hull of S [1]. 2.1 Delaunay triangulations Triangulation of a set of points S is a Delaunay triangulation, if and only if the circumcircle of any triangle does not contain any other points of S in its interior. It is Delaunay triangulation that is usually used when a triangulation is needed. It can be generated very effectively, in time O(n log n) for expected case, where n is the number of points to be triangulated. Moreover, a Delaunay triangulation maximizes the minimal angle. It even performs better than that: it lexicographically maximizes an angle vector V = (α 1 , α 2 , ..., α 3t ), α 1 ≤ α 2 ≤ ... ≤ α 3t , where t is the number of triangles. [4] A Delaunay triangulation is dual to another fundamental structure of computational geometry, Voronoi diagram. At this point, we will use Voronoi diagram to show the problem we are forced to deal with when using our mesh for interpolation purposes. When the values are interpolated from the vertices to the interior area, it is suitable to have triangles of aproximately equilateral shape.Triangulation of a set of points S is a Delaunay triangulation, if and only if the circumcircle of any triangle does not contain any other points of S in its interior. Figure 1: Voronoi diagrams for a) equilateral triangle, b) triangle with a "very sharp" angle, c) an obtuse-angle triangle See Figure 1. On the left (a), we can see a triangle that is just perfect for interpolation - all its vertices "have the same influence" over the triangle. 96
BRESENHAM'S REGULAR MESH DEFORMATION AND... For a triangle with one “very sharp” angle (b), the situation is not that good, but at least it is consistent in the following meaning: each vertex should affect primarily those points inside the triangle which are its close neighbours. Now, lets take a look at what happens when an obtuse angle appears (c). One can see that vertex A does not affect big part of line BC, although it is much closer to it than any of vertices B and C. This may cause big interpolation errors. So we will want to avoid "big" angles [2]. Unfortunately, although a Delaunay triangulation can help us a little in this aspect , it is generally not the one that minimizes the maximal angle. An example that illustrates this is on the Figure 2. Figure 2: Minmax triangulation that is not Delaunay. The diagonal divides the largest angle, producing a minmax triangulation. But there also is a triangle, which fails Delaunay criterion, because its circumcircle is not empty. 2.2 Minmax angle triangulations Minmax angle triangulation is the one that minimizes the minimal angle over all possible triangulations. It is therefore the one we would like to use for linear interpolation. The bad news is that the best known algorithm for minmax triangulation needs time O(n 2 log n) in the worst case [1]. However, for some special cases, a Delaunay triangulation satisfies minmax criterion. We will use this observation later. 3 Structured meshes Let us have a regular, orthogonal grid (its elements are identical squares). Now, lets deform the grid so that it respect the given geometry. The result 97
- Page 46 and 47: ÞÓÓÑÓÒÒ×ÓÙÖÒ×ÓÙÊ´
- Page 48 and 49: ÔÖØÓÑØÓÔÓÝÙ×ÓÙÖÓÚ
- Page 50 and 51: Bohumír Bastl The second reason fo
- Page 52 and 53: Bohumír Bastl 4 3 2 1 0 2 4 0 2 4
- Page 54 and 55: Bohumír Bastl the package that bot
- Page 56 and 57: Zuzana Benáková x B y z B B = u c
- Page 58 and 59: Zuzana Benáková Obrázek 4: třet
- Page 60 and 61: Michal Benes Since des t 2 = ds 2 +
- Page 62 and 63: Michal Benes The latter condition i
- Page 64 and 65: Michal Benes 4 Innitesimal deformat
- Page 66 and 67: Michal Benes where , 1, 2, ¢1, ¢2
- Page 68 and 69: Milan Bořík, Vojtěch Honzík sys
- Page 70 and 71: Milan Bořík, Vojtěch Honzík Dá
- Page 72 and 73: Milan Bořík, Vojtěch Honzík Obr
- Page 74 and 75: Jaromír Dobrý Example 2 Let’s c
- Page 76 and 77: Jaromír Dobrý M ϕ ϕ(M) ϑ V {o}
- Page 78 and 79: Jaromír Dobrý It is obvious from
- Page 80 and 81: Henryk Gliński distortion coeffici
- Page 82 and 83: Henryk Gliński polynomial. The coe
- Page 84 and 85: ÊÓÑÒÀõ ÎÖÚÑÓúÒÓÔÖÓ
- Page 86 and 87: ÊÓÑÒÀõ ÃõÒÐÓÝ´ÚÞÇÖ
- Page 88 and 89: ÊÓÑÒÀõ ÈÖÓÖÑÓÒ ÒÑÞ
- Page 90 and 91: Oldřich Hykš foundations and simp
- Page 92 and 93: Oldřich Hykš The choice of the tr
- Page 94 and 95: Oldřich Hykš discussed together w
- Page 98 and 99: Petr Kahánek, Alexej Kolcun of thi
- Page 100 and 101: Petr Kahánek, Alexej Kolcun Now we
- Page 102 and 103: Petr Kahánek, Alexej Kolcun In the
- Page 104 and 105: Mária Kmeťová Bernsteinove polyn
- Page 106 and 107: Mária Kmeťová H 0 O H 1 V 1 V 2
- Page 108 and 109: Mária Kmeťová Na obrázku 5 je k
- Page 110 and 111: Mária Kmeťová 6 Záver Program E
- Page 112 and 113: Milada Kočandrlová Pro body X eli
- Page 114 and 115: Milada Kočandrlová Podmínku cykl
- Page 116 and 117: Milada Kočandrlová 6 Homotetický
- Page 118 and 119: Jiří Kosinka medial axis transfor
- Page 120 and 121: Jiří Kosinka Figure 3: Asymptotic
- Page 122 and 123: Jiří Kosinka a) b) Figure 4: a) F
- Page 124 and 125: Iva Křivková • body (např. Bri
- Page 126 and 127: Iva Křivková V trojrozměrném eu
- Page 128 and 129: Iva Křivková asymptoty). V tomto
- Page 130 and 131: Karolína Kundrátová myslíme roz
- Page 132 and 133: Karolína Kundrátová Bázové fun
- Page 134 and 135: Karolína Kundrátová Obr. 2: Kři
- Page 136 and 137: Miroslav Lávička of this article
- Page 138 and 139: Miroslav Lávička if 〈x, x〉 M
- Page 140 and 141: Miroslav Lávička n: x - x = 0 0 4
- Page 142 and 143: Miroslav Lávička 5 Conclusion In
- Page 144 and 145: Pavel Leischner Věta 1 (ptolemaiov
Petr Kahánek, Alexej Kolcun<br />
2 Triangulations<br />
Let S be a finite set of points in the Euclidean plane. A triangulation of S is<br />
a maximally connected straight line plane graph whose vertices are the<br />
points of S. By maximality, each face is a triangle except for the exterior<br />
face which is the complement of the convex hull of S [1].<br />
2.1 Delaunay triangulations<br />
Triangulation of a set of points S is a Delaunay triangulation, if and only if<br />
the circumcircle of any triangle does not contain any other points of S in its<br />
interior.<br />
It is Delaunay triangulation that is usually used when a triangulation is<br />
needed. It can be generated very effectively, in time O(n log n) for expected<br />
case, where n is the number of points to be triangulated. Moreover, a<br />
Delaunay triangulation maximizes the minimal angle. It even performs<br />
better than that: it lexicographically maximizes an angle vector V = (α 1 ,<br />
α 2 , ..., α 3t ), α 1 ≤ α 2 ≤ ... ≤ α 3t , where t is the number of triangles. [4]<br />
A Delaunay triangulation is dual to another fundamental structure of<br />
computational geometry, Voronoi diagram. At this point, we will use<br />
Voronoi diagram to show the problem we are forced to deal with when<br />
using our mesh for interpolation purposes.<br />
When the values are interpolated from the vertices to the interior area, it<br />
is suitable to have triangles of aproximately equilateral shape.Triangulation<br />
of a set of points S is a Delaunay triangulation, if and only if the<br />
circumcircle of any triangle does not contain any other points of S in its<br />
interior.<br />
Figure 1: Voronoi diagrams for a) equilateral triangle, b) triangle with<br />
a "very sharp" angle, c) an obtuse-angle triangle<br />
See Figure 1. On the left (a), we can see a triangle that is just perfect for<br />
interpolation - all its vertices "have the same influence" over the triangle.<br />
96