institut für kartographie und geoinformation prof. dr. lutz plümer geoinformation ii vorlesung 8...
Post on 05-Apr-2015
109 Views
Preview:
TRANSCRIPT
Institut für Kartographie und GeoinformationProf. Dr. Lutz Plümer
Geoinformation IIVorlesung 8
08.06.00
Foliendesign: cand. geod. Jörg Steinrücken
Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00
22
Zur Erinnerung: Punktsuche in Landkarten
In welcher Masche liegt q?
Außen
x
y
Landkarte S mitn Kanten
q
Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00
33
Außen
q
Lösung 1: Zerlegung in vertikale Streifen
x
y
Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00
44
Nachteil: O(n²) Speicherplatz im Worst Case
4
n
4
n
Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00
55
Lösung 2: Trapezkarte
R
Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00
66
Trapezkarte: Konstruktionsprinzip
• Die Trapezkarte T(S) einer Landkarte S wird wie folgt konstruiert:– Bilde ein umschließendes Rechteck R
• Vermeidung der unbeschränkten Masche „Außen“– Konstruiere für jeden Endpunkt P eines Segments aus S
eine obere und eine untere vertikale Extension (Linie); diese Linien enden am Schnittpunkt mit dem nächsten Segment aus S oder an R
• Diese Konstruktion zerlegt R in disjunkte Trapeze, die von (höchstens) 4 Seiten begrenzt werden: – Ein oder zwei vertikale Seiten, gebildet aus den Extensionen– Eine obere und eine untere Seite, gebildet aus den
Segmenten
Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00
77
Trapezkarte: vereinfachende Annahme
• Es gibt keine zwei Knoten mit gleicher x-Koordinate, d.h. es gibt keine vertikalen Kanten
• Falls diese Annahme nicht erfüllt ist, kann sie durch Rotation um den Ursprung + Scherung mit einem sehr kleinen Winkel , der topologisch „nichts kaputt macht“, hergestellt werden
• Später werden wir sehen, dass diese Transformation rein virtuell ist
Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00
88
5 Fälle für die vertikalen (linken) Kanten
1. Kante entartet zu einem Punkt
2. Die untere vertikale Er-weiterung trifft auf eine Kante von S
3. Die obere vertikale Er-weiterung trifft auf eine Kante von S
4. Kante besteht aus oberer und unterer Erweiterung
5. Kante besteht aus einer Kante von R
Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00
99
Bezeichnungen
leftp
rightp
top
bottom
Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00
1010
Trapezkarte - Eigenschaften
Satz: Eine Trapezkarte einer Landkarte S mit n Kanten enthält höchstens (a) 6n + 4 vertikale Linien und (b) 3n + 1 Trapeze.
Beweis: (a) Ein Knoten der Trapezkarte ist entweder– ein Eckpunkt von R 4– ein Knoten der Karte S 2n– ein Punkt, bei dem eine Erweiterung auf eine
Kante von S trifft 4n
Insgesamt: 6n + 4
Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00
1111
Trapezkarte - Eigenschaften
Satz: Eine Trapezkarte einer Landkarte S mit n Kanten enthält höchstens 6n + 4 vertikale Linien und sowie 3n + 1 Trapeze.
Beweis: (b): aus a) mit Eulers Formel
Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00
1212
Trapezkarte - Eigenschaften
Zwei Trapeze T und T‘ heißen adjazent, wenn sie sich entlang einer vertikalen Linie berühren
T‘T
Es gilt entweder top(T) = top(T‘) oder
bottom(T) = bottom(T‘)
T T‘
T‘T
Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00
1313
Probleme und weiteres Vorgehen
Probleme:• Konstruktion der Trapezkarte T(S)• Unterstützung der Suche in einer Trapezkarte
Idee für das weitere Vorgehen:• Unterstützung der Suche durch eine Art „binärer
Suchbaum“ D(S) mit 2 Sorten von Knoten–X-Knoten für Endpunkte / Extensionen
• Links oder Rechts?–Y-Knoten für Segmente
• Oben oder Unten?
• Trapezkarte und „Baum“ werden simultan konstruiert
Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00
1414
T(S) und D(S)
p1
A q1
s1
B
C
p2
q2
s2
s2
D
E
F
G
D
E
F
G
A
B
C
p1
q1s1
p2
q2
s2
Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00
1515
Schrittweise Konstruktion von T(S) und D(S)
Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00
1616
Trapezkarte und Suchstruktur zweier Kanten
X
p1
Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00
1717
Trapezkarte und Suchstruktur zweier Kanten
p1p1
p1
Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00
1818
Trapezkarte und Suchstruktur zweier Kanten
X
A
X
p1
q1
A
p1
Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00
1919
Trapezkarte und Suchstruktur zweier Kanten
A
q1
p1
A q1
s1s1
s1 q1
p1
Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00
2020
X
Trapezkarte und Suchstruktur zweier Kanten
A
p1
A q1
s1
B
B X
Y
Yp2
s1
p1
q1
Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00
2121
Trapezkarte und Suchstruktur zweier Kanten
A
B
p1
A q1
s1
B
Y
Yp2p2
s1 q1
p2
p1
Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00
2222
Trapezkarte und Suchstruktur zweier Kanten
A
B
p1
A q1
s1
B
Y
Yp2
CC
XX
q2
s1 q1
p2
p1
Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00
2323
Trapezkarte und Suchstruktur zweier Kanten
A
B
C X
p1
A q1
s1
B
C
p2
X
q2
q2q2
s2
s1
p2
q1
p1
Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00
2424
A
B
C
Trapezkarte und Suchstruktur zweier Kanten
p1
A q1
s1
B
C
p2
q2
p2
s2
s2
p1
q1s1
q2
s2
s2
s2
Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00
2525
D
E
F
G
A
B
C
Trapezkarte und Suchstruktur zweier Kanten
p1
A q1
s1
B
C
p2
q2
s2
s2
p1
q1s1
D
E
F
Gp2
q2
s2
Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00
2626
Beachte:
• T(S) und D(S) werden simultan konstruiert• D(S) ist kein Baum, sondern ein „DAG“, ein „directed
acyclic graph“, ein gerichteter azyklischer Graph• Dieser DAG ist zusammenhängend und hat genau
eine Wurzel• Unterschied zum Baum: innere Knoten können
mehrere Vorgänger haben• Die Blätter von D(S) und die Trapeze von T(S)
referenzieren sich gegenseitig• Wie stets hängt die Tiefe (=Güte) des Baumes von
der Reihenfolge der Bearbeitung der Segmente ab• Idee: Zufällige Permutation der Segmente von S
Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00
2727
Datenstruktur für T(S)
• Möglich wäre eine doppelt verkettete Kantenliste (s. Overlay)• Wegen der einfachen Struktur der Trapeze bietet sich folgende
Alternative an:• Knoten (mit Koordinaten)• Segmente (mit Referenzen auf Knoten)• Trapeze mit Referenzen auf:
– Top
– Bottom
– Leftp
– Rightp
– alle (maximal 4) Nachbarn
• Beachte: Die Geometrie der Trapeze ist nur implizit, kann aber in konstanter Zeit rekonstruiert werden
Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00
2828
Skizze: Konstruktion von T(S) und D(S)
Input: Eine Menge S von n Segmenten
1. Konstruiere ein umschließendes Rechteck R
2. Berechne eine Permutation s1, s2,...,sn von S
3. for i= 1 to n doKonstruiere T(Si) und D(Si) mit Si = {s1, s2,...,si} unter Verwendung von T(Si-1) und D(Si-1)
„Schleifeninvariante“: T(Si-1) ist eine Trapezkarte und D(Si-1) ist eine Suchstruktur für diese Trapezkarte
Der Unterschied zwischen „i-1“ und „i“ betrifft genau die Trapeze in T(Si-1) , die von si geschnitten werden
Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00
2929
Von T(Si-1) und D(Si-1) zu T(Si) und D(Si)
Fall 1: si liegt in genau einem Trapez
Input: Das Segment si, T(Si-1) und D(Si-1)
Idee: Nutze die Suchstruktur D(Si-1), um schnell zu finden
Vorgehen: Modifiziere T(Si-1) und D(Si-1) in der auf den folgenden Folien beschriebenen Weise:
Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00
3030
Einfügen einer Kante (I)
D(Si-1)
pi
T(Si-1)
Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00
3131
Einfügen einer Kante (I)
D(Si-1)
pipi
sisi
si
qiqi
pi
qi
Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00
3232
Einfügen einer Kante (I)
C D
B
A
D(Si-1)
pi
qi
si
AC
BD
pi
qisi
T(Si) D(Si)
Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00
3333
Von T(Si-1) und D(Si-1) zu T(Si) und D(Si) Fall2
Fall 2: si liegt in mehrere Trapezen 1, 2,... n
Input: Das Segment si, T(Si-1) und D(Si-1)
Teilziel: Bestimmung von 1, 2,... n Beobachtungen:• i ist rechter Nachbar von i-1
• Jedes Trapez hat maximal 2 rechte Nachbarn• Der richtige kann mit rightp() rasch identifiziert werden:
– wenn rightp() oberhalb von si, dann wähle den unteren rechten Nachbarn, sonst den oberen
– Der Übergang von i zu i+1 ist also einfach
• Nutze D(Si-1) um 1 schnell zu finden
Modifiziere T(Si-1) und D(Si-1) nun auf die folgende Weise:
Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00
3434
Einfügen einer Kante (II)
3
D(Si-1)
1 20 3
0
1
2
qi
T(Si-1)
pi
Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00
3535
Einfügen einer Kante (II)
0
1
2 qi
si
qi
qi
pi
D(Si-1)
1 20
Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00
3636
Einfügen einer Kante (II)
D(Si-1)
qi
si
si sisi
qi
pi
si si si
sisi
Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00Lutz Plümer - Geoinformation II - SS 2000 - Vorlesung 8 - 08.06.00
3737
Einfügen einer Kante (II)
D(Si-1)
qisi sisiA
B
DE
FB D si
A
C
F
E
Csipi
qi
T(Si)D(Si)
Schönen Dank für Ihre Aufmerksamkeit und
Auf Wiedersehen
top related