Geoinformation II

53
Institut für Kartographie und Geoinformation Prof. Dr. Lutz Plümer Geoinformation II Vorlesung 6 SS 2000 Konstruktion der Voronoi-Diagramme II

description

Geoinformation II. Vorlesung 6 SS 2000. Konstruktion der Voronoi-Diagramme II. Übersicht I. Divide and Conquer: Merge Konstruktion des Voronoi-Diagramms „Merge“ Extrempunkte von CH(P 1 )  CH(P 2 ) Tangente von CH(P 1 )  CH(P 2 ) Nochmals zur konvexen Hülle CH Tangente - PowerPoint PPT Presentation

Transcript of Geoinformation II

Page 1: Geoinformation II

Institut für Kartographie und GeoinformationProf. Dr. Lutz Plümer

Geoinformation IIVorlesung 6

SS 2000

Konstruktion der Voronoi-Diagramme II

Page 2: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

2 2

Übersicht I

• Divide and Conquer: Merge• Konstruktion des Voronoi-Diagramms• „Merge“

• Extrempunkte von CH(P1) CH(P2)

• Tangente von CH(P1) CH(P2)

• Nochmals zur konvexen Hülle CH• Tangente• Nachfolger - Bestimmung• Nachfolger• Bestimmung der (oberen) Tangenten der konvexen

Hüllen

Page 3: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

3 3

Übersicht II

• Extrempunkte

• 2 vertikal monotone Kantenzüge

• Tangente

• Bestimmung des Nachfolgers

• Bestimmung des Nachfolgers

• Bestimmung des Nachfolgers

• Bestimmung des Nachfolgers

• Konvexe Hülle

• Bestimmung des Nachfolgers

• Konvexe Hülle

Page 4: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

4 4

Übersicht III

• Konstruiere den trennenden Kantenzug als Verbindung der beiden Mittelsenkrechten

• Vereinigung

• Konstruiere den trennenden Kantenzug als Verbindung der beiden Mittelsenkrechten

• Länge des Kantenzuges im Worst Case

• Größenordnung des Kanten-Umrings im worst case

• O(n) * O(n) = O(n2) ?

• „Investitionen müssen sich amortisieren“

Page 5: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

5 5

Divide and Conquer: Merge

Page 6: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

6 6

Konstruktion des Voronoi-Diagramms

„Divide and Conquer“• Input: Gegeben ist eine Menge P von mindestens 2 Punkten

• Split: Zerlege P in zwei etwa gleich große Teilmengen P1 und P2

• Rekursiv: Berechne Voronoi-Diagramme von VD(P1) und VD(P2)

• Merge: Verknüpfe VD(P1) und VD(P2)

Page 7: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

7 7

„Merge“

• Die Voronoi-Diagramme VD(P1) und VD(P2) sind bereits berechnet.

• Die konvexen Hüllen CH(P1) und CH(P2) seien ebenfalls an dieser Stelle bekannt.

1. Bestimme die oberen und unteren Extrempunkte und die beiden oberen und unteren Tangenten von CH(P1) CH(P2)

2. Konstruiere CH(P1 P2)

3. Bilde die Mittelsenkrechten zu den beiden neu eingeführten Kanten

4. Konstruiere den trennenden Kantenzug als Verbindung der beiden Mittelsenkrechten

5. Entferne die überstehenden Kanten

6. Bilde die neu entstandenen Voronoi-Regionen (Maschen)

Page 8: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

8 8

max y

min ymin y

max y

Extrempunkte von CH(P1) CH(P2)

Page 9: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

9 9

Tangente von CH(P1) CH(P2)

Page 10: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

10 10

Nochmals zur konvexen Hülle CH

Was wissen wir über die „konvexe Hülle“ CH(P) einer Punktmenge P?

Die Extrempunkte sind die Knoten auf der Grenze von CH.

• Zu je zwei Punkten P1 und P2 ist die verbindende Kante ganz in CH enthalten.

• Der obere und der untere Extrempunkt zerlegen die Grenze von CH in zwei vertikal monotone Kantenzüge.

• Die Verbindungskante k zweier Punkte P1 und P2 aus P definiert eine Randkante von CH genau dann, wenn alle übrigen Punkte von P auf der gleichen Seite von k liegen.

• P2 ist genau dann Nachfolger von P1 auf dem Rand von CH, wenn der zugehörige polare Winkel von P2 minimal ist.

Page 11: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

11 11

Tangente

Page 12: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

12 12

Nachfolger - Bestimmung

Winkel minimal

P1

P2

Page 13: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

13 13

Nachfolger

Winkel minimal

P2

P1

Page 14: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

14 14

Bestimmung der (oberen) Tangenten der konvexen Hüllen

• Bestimme die oberen und unteren Extrempunkte von CH(P1), CH(P2) und CH(P1) CH(P2)

• Betrachte die oberen Extrempunkte P1 und Q1 und die Nachfolger P2 und Q2 im Uhrzeigersinn, und sei P1 höher als Q1

• Bestimme das Minimum der mit P1P2, P1Q1 und P1Q2 assoziierten Winkel

• Fälle:– P1 Q1 ist minimal: Tangente gefunden, fertig

– P1 P2 minimal: ersetze P1 durch P2 und P2 durch P3 (wandere auf der linken konvexen Hülle im Uhrzeigersinn)

– P1 Q2 minimal: ersetze Q1 durch Q2 und Q2 durch Q3 (wandere auf der rechten konvexen Hülle im Uhrzeigersinn)

• Der Fall der unteren Tangente ist symmetrisch

Page 15: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

15 15

Extrempunkte

Page 16: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

16 16

2 vertikal monotone Kantenzüge

Page 17: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

17 17

Tangente

Page 18: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

18 18

Bestimmung des Nachfolgers

Winkel nicht minimal

Page 19: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

19 19

Bestimmung des Nachfolgers

Winkel minimal

Page 20: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

20 20

Bestimmung des Nachfolgers

Page 21: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

21 21

Bestimmung des Nachfolgers

Page 22: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

22 22

Konvexe Hülle

Page 23: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

23 23

Bestimmung des Nachfolgers

Page 24: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

24 24

Konvexe Hülle

Page 25: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

25 25

Konstruiere den trennenden Kantenzug als Verbindung der beiden Mittelsenkrechten

Page 26: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

26 26

Vereinigung

Mittelsenkrechte bilden

Page 27: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

27 27

Vereinigung

Page 28: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

28 28

Vereinigung

Aktive Voronoi-Diagramme

Schnittpunkte mit Seg-menten suchen

Page 29: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

29 29

Vereinigung

Aktive Voronoi-Diagramme

Schnittpunkte mit Seg-menten suchen

Neues aktives VD

Page 30: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

30 30

Vereinigung

Aktive Voronoi-Diagramme

Schnittpunkte mit Seg-menten suchen

Neues aktives VD

Mittelsenkrechte zuwischenden aktiven VD

Page 31: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

31 31

Vereinigung

Schnittpunkte suchen

Page 32: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

32 32

Vereinigung

Schnittpunkte suchen

Neues aktives VD suchen

Page 33: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

33 33

Vereinigung

Schnittpunkte suchen

Neues aktives VD suchen

Page 34: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

34 34

Vereinigung

Schnittpunkte suchen

Neues aktives VD suchen

Mittelsenkrechte deraktiven VD

Page 35: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

35 35

Vereinigung

Schnittpunkte suchen

Page 36: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

36 36

Vereinigung

Schnittpunkte suchen

Neues aktives VD suchen

Page 37: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

37 37

Vereinigung

Schnittpunkte suchen

Neues aktives VD suchen

Mittelsenkrechte deraktiven VD

Page 38: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

38 38

Vereinigung

Nächsten relevanten Schnittpunkte suchen

Neues aktives VD suchen

Page 39: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

39 39

Vereinigung

Nächsten relevanten Schnittpunkte suchen

Neues aktives VD suchen

Mittelsenkrechte deraktiven VD

Page 40: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

40 40

Vereinigung

Nächsten relevanten Schnittpunkte suchen

Neues aktives VD suchen

Page 41: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

41 41

Vereinigung

Nächsten relevanten Schnittpunkte suchen

Neues aktives VD suchen

Mittelsenkrechte deraktiven VD

Page 42: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

42 42

Vereinigung

Nächsten relevanten Schnittpunkte suchen

Neues aktives VD suchen

Page 43: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

43 43

Vereinigung

Nächsten relevanten Schnittpunkte suchen

Neues aktives VD suchen

Mittelsenkrechte deraktiven VD

Page 44: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

44 44

Vereinigung

Nächsten relevanten Schnittpunkte suchen

Neues aktives VD suchen

Page 45: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

45 45

Vereinigung

Nächsten relevanten Schnittpunkte suchen

Neues aktives VD suchen

Mittelsenkrechte deraktiven VD

Page 46: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

46 46

Vereinigung

Nächsten relevanten Schnittpunkte suchen

Neues aktives VD suchen

Page 47: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

47 47

Vereinigung

Nächsten relevanten Schnittpunkte suchen

Neues aktives VD suchen

Verknüpfung mit der Mittel-senkrechten vom Anfang

Page 48: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

48 48

Konstruiere den trennenden Kantenzug als Verbindung der beiden Mittelsenkrechten

• gegeben: die beiden oberen und unteren Mittelsenkrechten g und g*• die zugehörigen oberen Voronoi-Regionen seien P und Q• Solange die untere Mittelsenkrechte noch nicht erreicht ist

– Bestimme für die aktuelle Mittelsenkrechte

• die Austrittspunkte p und q aus den aktuellen Voronoi-Regionen,

• die zugehörigen Kanten

• die zugehörigen Nachbarn P‘ und Q‘

– wenn p höher ist als q

• ersetze P durch P‘ und schneide g an der Stelle p ab

– wenn q höher als p

• ersetze Q durch Q‘ und schneide g an der Stelle q ab

– bestimme die aktuelle Mittelsenkrechte g des neuen Paares P, Q

Page 49: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

49 49

O(n)

Länge des Kantenzuges im Worst Case

Page 50: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

50 50

Größenordnung des Kanten-Umrings im worst case

O(n)

Page 51: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

51 51

O(n) * O(n) = O(n2) ?

Voronoi-Regionen sind konvex

Kantenzug ist monoton

war jetzt alles umsonst?

Page 52: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

52 52

O(n) * O(n) = O(n2) ?

Voronoi-Regionen sind konvex

Kantenzug ist monoton

Keine Kante öfter als zwei mal anfassen!

Page 53: Geoinformation II

Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6 Lutz Plümer - Geoinformation - 6. Semester - SS 2001 - Vorlesung 6

53 53

„Investitionen müssen sich amortisieren“

• Ziel: keine Kante mehr als zwei mal „anfassen“• Es gibt insgesamt höchstens 3* n – 6 Kanten O(n)• Konvexität der Voronoi-Regionen höchstens

zwei Schnittpunkte mit der aktiven Halbgeraden• Es genügt, die linken (grünen) Kantenumringe im

Uhrzeigersinn und die rechten (roten) Kantenumringe gegen den Uhrzeigersinn zu durchlaufen und den zuletzt gefundenen und verworfenen Schnittpunkt als Haltepunkt zu merken!