Identifikation identischer Punkte in Höhennetzen mehrerer ...–Rückführung auf Ausgleichung nach...

29
Titelmaster Geodätische Woche 2010, Köln Identifikation identischer Punkte in Höhennetzen mehrerer Epochen Sebastian Halsig Sebastian Halsig, Andreas Ernst, Wolf-Dieter Schuh 07. Oktober 2010

Transcript of Identifikation identischer Punkte in Höhennetzen mehrerer ...–Rückführung auf Ausgleichung nach...

Page 1: Identifikation identischer Punkte in Höhennetzen mehrerer ...–Rückführung auf Ausgleichung nach Bedingungen 14 Modell: Ausgleich nach Bedingungen mit folgender Bedingungsmatrix

Titelmaster

Geodätische Woche 2010, Köln

Identifikation identischer Punktein Höhennetzen mehrerer Epochen

Sebastian Halsig

Sebastian Halsig, Andreas Ernst, Wolf-Dieter Schuh

07. Oktober 2010

Page 2: Identifikation identischer Punkte in Höhennetzen mehrerer ...–Rückführung auf Ausgleichung nach Bedingungen 14 Modell: Ausgleich nach Bedingungen mit folgender Bedingungsmatrix

Motivation

• Häufiges Problem in der Geodäsie:

– Einhängen eines geodätischen Netzes in ein Netz höherer Ordnung

2

Sebastian Halsig 07.10.2010 Identifikation identischer Punkte in Höhennetzen mehrerer Epochen

Page 3: Identifikation identischer Punkte in Höhennetzen mehrerer ...–Rückführung auf Ausgleichung nach Bedingungen 14 Modell: Ausgleich nach Bedingungen mit folgender Bedingungsmatrix

Motivation

• Häufiges Problem in der Geodäsie:

– Einhängen eines geodätischen Netzes in ein Netz höherer Ordnung

• Problem: Deformation an der Erdoberfläche

– häufig in Bergbauregionen

3

Sebastian Halsig

– häufig in Bergbauregionen

– oftmals ruckartig und räumlich eingegrenzt

07.10.2010 Identifikation identischer Punkte in Höhennetzen mehrerer Epochen

Page 4: Identifikation identischer Punkte in Höhennetzen mehrerer ...–Rückführung auf Ausgleichung nach Bedingungen 14 Modell: Ausgleich nach Bedingungen mit folgender Bedingungsmatrix

Motivation

• Häufiges Problem in der Geodäsie:

– Einhängen eines geodätischen Netzes in ein Netz höherer Ordnung

• Problem: Deformation an der Erdoberfläche

– häufig in Bergbauregionen

4

Sebastian Halsig

– häufig in Bergbauregionen

– oftmals ruckartig und räumlich eingegrenzt

• Punktauswahl:

– Punktidentitäten

– Welche Punkte sind über einen längeren Betrachtungszeitraum unbewegt?

07.10.2010 Identifikation identischer Punkte in Höhennetzen mehrerer Epochen

Page 5: Identifikation identischer Punkte in Höhennetzen mehrerer ...–Rückführung auf Ausgleichung nach Bedingungen 14 Modell: Ausgleich nach Bedingungen mit folgender Bedingungsmatrix

Ziel und Problemstellung

• Ziel: Entwicklung eines Detektionsverfahrens

– Identitätsmodell

– Suche nach Punktidentitäten5

Sebastian Halsig

• Problem: Lagerung des Höhennetzes

– Detektionsverfahren muss unabhängig von Datumspunkt sein!

07.10.2010 Identifikation identischer Punkte in Höhennetzen mehrerer Epochen

Page 6: Identifikation identischer Punkte in Höhennetzen mehrerer ...–Rückführung auf Ausgleichung nach Bedingungen 14 Modell: Ausgleich nach Bedingungen mit folgender Bedingungsmatrix

Umsetzung

6

• Datengrundlage:

• Überwachungsmessungen im Rahmen des Leitnivellements der GEOBasis NRW

• Standard-Datenverarbeitung

Sebastian Halsig 07.10.2010 Identifikation identischer Punkte in Höhennetzen mehrerer Epochen

• Leitnivellement Ruhrgebiet

• Datumspunkt

in Essen-Stoppenberg

Page 7: Identifikation identischer Punkte in Höhennetzen mehrerer ...–Rückführung auf Ausgleichung nach Bedingungen 14 Modell: Ausgleich nach Bedingungen mit folgender Bedingungsmatrix

Detektionsverfahren

• Aufbau des Detektionsverfahrens

– Zweistufiges Verfahren

– Gauß-Markov-Modell mit Restriktionen

7

1. Auswertungsstufe 2. Auswertungsstufe

Sebastian Halsig 07.10.2010 Identifikation identischer Punkte in Höhennetzen mehrerer Epochen

Einzelepochenauswertung Ausgleich nach Bedingungen

Freie Netzausgleichung

� Eindeutige LagerungHypothesentests zur Detektion von Deformationen

Netzanalyse Algorithmen zur Punktauswahl

Page 8: Identifikation identischer Punkte in Höhennetzen mehrerer ...–Rückführung auf Ausgleichung nach Bedingungen 14 Modell: Ausgleich nach Bedingungen mit folgender Bedingungsmatrix

Detektionsverfahren

• 1. Auswertungsstufe

– Einzelepochenausgleich

8

Modell:

Reduzierte Beobachtungen u. Verbesserungen

Parameterzuschläge

Designmatrix

Sebastian Halsig 07.10.2010 Identifikation identischer Punkte in Höhennetzen mehrerer Epochen

Page 9: Identifikation identischer Punkte in Höhennetzen mehrerer ...–Rückführung auf Ausgleichung nach Bedingungen 14 Modell: Ausgleich nach Bedingungen mit folgender Bedingungsmatrix

Detektionsverfahren

• 1. Auswertungsstufe

– Einzelepochenausgleich

9

Modell: Normalgleichungssystem:Reduzierte Beobachtungen u. Verbesserungen

Parameterzuschläge

Designmatrix

Sebastian Halsig 07.10.2010 Identifikation identischer Punkte in Höhennetzen mehrerer Epochen

Page 10: Identifikation identischer Punkte in Höhennetzen mehrerer ...–Rückführung auf Ausgleichung nach Bedingungen 14 Modell: Ausgleich nach Bedingungen mit folgender Bedingungsmatrix

Detektionsverfahren

• 1. Auswertungsstufe

– Einzelepochenausgleich

10

Modell: Normalgleichungssystem:Reduzierte Beobachtungen u. Verbesserungen

Parameterzuschläge

Designmatrix

Sebastian Halsig 07.10.2010 Identifikation identischer Punkte in Höhennetzen mehrerer Epochen

Page 11: Identifikation identischer Punkte in Höhennetzen mehrerer ...–Rückführung auf Ausgleichung nach Bedingungen 14 Modell: Ausgleich nach Bedingungen mit folgender Bedingungsmatrix

Detektionsverfahren

• 1. Auswertungsstufe

– Einzelepochenausgleich

11

Modell: Normalgleichungssystem:Reduzierte Beobachtungen u. Verbesserungen

Parameterzuschläge

Designmatrix

Sebastian Halsig 07.10.2010 Identifikation identischer Punkte in Höhennetzen mehrerer Epochen

Verbesserung und empirischer Varianzfaktor der ersten Stufe:

Page 12: Identifikation identischer Punkte in Höhennetzen mehrerer ...–Rückführung auf Ausgleichung nach Bedingungen 14 Modell: Ausgleich nach Bedingungen mit folgender Bedingungsmatrix

Detektionsverfahren

• 2. Auswertungsstufe

– Rückführung auf Ausgleichung nach Bedingungen

12

Modell:

Sebastian Halsig 07.10.2010 Identifikation identischer Punkte in Höhennetzen mehrerer Epochen

Page 13: Identifikation identischer Punkte in Höhennetzen mehrerer ...–Rückführung auf Ausgleichung nach Bedingungen 14 Modell: Ausgleich nach Bedingungen mit folgender Bedingungsmatrix

Detektionsverfahren

• 2. Auswertungsstufe

– Rückführung auf Ausgleichung nach Bedingungen

13

Modell: Ausgleich nach Bedingungen mit folgender Bedingungsmatrix

Sebastian Halsig 07.10.2010 Identifikation identischer Punkte in Höhennetzen mehrerer Epochen

Page 14: Identifikation identischer Punkte in Höhennetzen mehrerer ...–Rückführung auf Ausgleichung nach Bedingungen 14 Modell: Ausgleich nach Bedingungen mit folgender Bedingungsmatrix

Detektionsverfahren

• 2. Auswertungsstufe

– Rückführung auf Ausgleichung nach Bedingungen

14

Modell: Ausgleich nach Bedingungen mit folgender Bedingungsmatrix

Sebastian Halsig 07.10.2010 Identifikation identischer Punkte in Höhennetzen mehrerer Epochen

Auftransformation auf alle als identisch angenommenen Punkte

Page 15: Identifikation identischer Punkte in Höhennetzen mehrerer ...–Rückführung auf Ausgleichung nach Bedingungen 14 Modell: Ausgleich nach Bedingungen mit folgender Bedingungsmatrix

Detektionsverfahren

• 2. Auswertungsstufe

– Rückführung auf Ausgleichung nach Bedingungen

15

Modell: Ausgleich nach Bedingungen mit folgender Bedingungsmatrix

Sebastian Halsig 07.10.2010 Identifikation identischer Punkte in Höhennetzen mehrerer Epochen

Auftransformation auf alle als identisch angenommenen Punkte

Residuen:

ausgeglichene Parameter in zweiter Stufe:

empirischer Varianzfaktor der zweiten Stufe:

Page 16: Identifikation identischer Punkte in Höhennetzen mehrerer ...–Rückführung auf Ausgleichung nach Bedingungen 14 Modell: Ausgleich nach Bedingungen mit folgender Bedingungsmatrix

Detektionsverfahren

• Hypothesentest

– Wiederholte Anwendung unter Verwendung eines speziellen Algorithmus

16Nullhypothese und Alternativhypothese:

Sebastian Halsig 07.10.2010 Identifikation identischer Punkte in Höhennetzen mehrerer Epochen

Page 17: Identifikation identischer Punkte in Höhennetzen mehrerer ...–Rückführung auf Ausgleichung nach Bedingungen 14 Modell: Ausgleich nach Bedingungen mit folgender Bedingungsmatrix

Detektionsverfahren

• Hypothesentest

– Wiederholte Anwendung unter Verwendung eines speziellen Algorithmus

17Nullhypothese und Alternativhypothese:

Realisierung der Testgröße

Sebastian Halsig 07.10.2010 Identifikation identischer Punkte in Höhennetzen mehrerer Epochen

Realisierung der Testgröße

über die empirischen Varianzfaktoren von erster und zweiter Stufe

Fisher-Verteilung n … Anzahl Beobachtungen

m … Anzahl Parameter

b … Anzahl Datumsrestriktionen

p … Anzahl Restriktionsgleichungen

Page 18: Identifikation identischer Punkte in Höhennetzen mehrerer ...–Rückführung auf Ausgleichung nach Bedingungen 14 Modell: Ausgleich nach Bedingungen mit folgender Bedingungsmatrix

Detektionsverfahren

• Hypothesentest

– Wiederholte Anwendung unter Verwendung eines speziellen Algorithmus

18Nullhypothese und Alternativhypothese:

Realisierung der Testgröße

Sebastian Halsig 07.10.2010 Identifikation identischer Punkte in Höhennetzen mehrerer Epochen

Realisierung der Testgröße

über die empirischen Varianzfaktoren von erster und zweiter Stufe

Hypothesentest:

Nullhypothese wird widerlegt, wenn

Fisher-Verteilung

Quantil der Fisher-Verteilung mit entsprechenden Freiheitsgraden und festgelegtem Signifikanzniveau

n … Anzahl Beobachtungen

m … Anzahl Parameter

b … Anzahl Datumsrestriktionen

p … Anzahl Restriktionsgleichungen

Page 19: Identifikation identischer Punkte in Höhennetzen mehrerer ...–Rückführung auf Ausgleichung nach Bedingungen 14 Modell: Ausgleich nach Bedingungen mit folgender Bedingungsmatrix

Algorithmus

19

Ausgangslage:

Punktgruppe mit fester Punktanzahl n

Hinzunahme weiterer Punktgruppe mit fester

Punktanzahl n

Entferne Punkt aus „Testmenge“

Gesamtheit aller Punkte

Abbruch, wenn alle Punkte durchlaufen

Strenges Kriterium zur Auswahl notwendig!

Einzelepochen-ausgleich

1. Auswertungsstufe

Sebastian Halsig 07.10.2010 Identifikation identischer Punkte in Höhennetzen mehrerer Epochen

Hypothesen-test

Nullhypothese kann nicht widerlegt werden!

Nullhypothese abgelehnt –

Alternativhypothese angenommen

Testergebnis

Ausgleich nach Bedingungen

2. Auswertungsstufe

Page 20: Identifikation identischer Punkte in Höhennetzen mehrerer ...–Rückführung auf Ausgleichung nach Bedingungen 14 Modell: Ausgleich nach Bedingungen mit folgender Bedingungsmatrix

Ergebnisdiskussion

20

Darstellung der

• Punktidentitäten (grün)

•Bewegten Punkte (blau)

•Betrag der Deformation

Sebastian Halsig 07.10.2010 Identifikation identischer Punkte in Höhennetzen mehrerer Epochen

•Betrag der Deformation

Epochen: 2002-2008

Datumspunkt: 4508900382

Signifikanzniveau: 95 %

���� Anzahl der Punktidentitäten: 238

(von 1060 gemeinsamen Punkten)

Page 21: Identifikation identischer Punkte in Höhennetzen mehrerer ...–Rückführung auf Ausgleichung nach Bedingungen 14 Modell: Ausgleich nach Bedingungen mit folgender Bedingungsmatrix

Ergebnisdiskussion

21

Darstellung der

• Punktidentitäten (grün)

•Bewegten Punkte (blau)

•Betrag der Deformation

Sebastian Halsig 07.10.2010 Identifikation identischer Punkte in Höhennetzen mehrerer Epochen

•Betrag der Deformation

Epochen: 2004-2008

Datumspunkt: 4508900382

Signifikanzniveau: 95 %

���� Anzahl der Punktidentitäten: 532

(von 1205 gemeinsamen Punkten)

Page 22: Identifikation identischer Punkte in Höhennetzen mehrerer ...–Rückführung auf Ausgleichung nach Bedingungen 14 Modell: Ausgleich nach Bedingungen mit folgender Bedingungsmatrix

Ergebnisdiskussion

22

Signifikanzniveau: 99 %

���� Anzahl der Punktidentitäten: 587

(von 1205 gemeinsamen Punkten)

Sebastian Halsig 07.10.2010 Identifikation identischer Punkte in Höhennetzen mehrerer Epochen

Epochen: 2004-2008

Datumspunkt: 4508900382

Signifikanzniveau: 95 %

���� Anzahl der Punktidentitäten: 532

(von 1205 gemeinsamen Punkten)

Page 23: Identifikation identischer Punkte in Höhennetzen mehrerer ...–Rückführung auf Ausgleichung nach Bedingungen 14 Modell: Ausgleich nach Bedingungen mit folgender Bedingungsmatrix

Ergebnisdiskussion

23

Signifikanzniveau: 99 %

���� Anzahl der Punktidentitäten: 587

(von 1205 gemeinsamen Punkten)

• Erhöhung des Signifikanzniveaus

– Anzahl der gefundenen Punktidentitäten steigt

Sebastian Halsig 07.10.2010 Identifikation identischer Punkte in Höhennetzen mehrerer Epochen

Epochen: 2004-2008

Datumspunkt: 4508900382

Signifikanzniveau: 95 %

���� Anzahl der Punktidentitäten: 532

(von 1205 gemeinsamen Punkten)

Page 24: Identifikation identischer Punkte in Höhennetzen mehrerer ...–Rückführung auf Ausgleichung nach Bedingungen 14 Modell: Ausgleich nach Bedingungen mit folgender Bedingungsmatrix

Ergebnisdiskussion

24

Signifikanzniveau: 99 %

���� Anzahl der Punktidentitäten: 587

(von 1205 gemeinsamen Punkten)

• Erhöhung des Signifikanzniveaus

– Anzahl der gefundenen Punktidentitäten steigt

• Hypothesentest verliert an

Sebastian Halsig 07.10.2010 Identifikation identischer Punkte in Höhennetzen mehrerer Epochen

Epochen: 2004-2008

Datumspunkt: 4508900382

Signifikanzniveau: 95 %

���� Anzahl der Punktidentitäten: 532

(von 1205 gemeinsamen Punkten)

• Hypothesentest verliert an Sensitivität!

– Annahmebereich vergrößert

– Kleinere Abweichungen spielen weniger eine Rolle

– Wahrscheinlichkeit für Fehler 2. Art wächst

Page 25: Identifikation identischer Punkte in Höhennetzen mehrerer ...–Rückführung auf Ausgleichung nach Bedingungen 14 Modell: Ausgleich nach Bedingungen mit folgender Bedingungsmatrix

Ergebnisdiskussion

25

Darstellung der

• Punktidentitäten (grün)

•Bewegten Punkte (blau)

•Betrag der Deformation

Sebastian Halsig 07.10.2010 Identifikation identischer Punkte in Höhennetzen mehrerer Epochen

•Betrag der Deformation

Epochen: 2006-2008

Datumspunkt: 4508900382

Signifikanzniveau: 95 %

���� Anzahl der Punktidentitäten: 599

(von 1307 gemeinsamen Punkten)

Page 26: Identifikation identischer Punkte in Höhennetzen mehrerer ...–Rückführung auf Ausgleichung nach Bedingungen 14 Modell: Ausgleich nach Bedingungen mit folgender Bedingungsmatrix

Ergebnisdiskussion

26

Darstellung der

• Punktidentitäten (grün)

•Bewegten Punkte (blau)

•Betrag der Deformation

Aktuelle Deformation, aufgetreten zwischen 2006 und 2008

Sebastian Halsig 07.10.2010 Identifikation identischer Punkte in Höhennetzen mehrerer Epochen

•Betrag der Deformation

Epochen: 2006-2008

Datumspunkt: 4508900382

Signifikanzniveau: 95 %

���� Anzahl der Punktidentitäten: 599

(von 1307 gemeinsamen Punkten)

Page 27: Identifikation identischer Punkte in Höhennetzen mehrerer ...–Rückführung auf Ausgleichung nach Bedingungen 14 Modell: Ausgleich nach Bedingungen mit folgender Bedingungsmatrix

Fazit und Ausblick

Fazit:

– Detektion von Punktidentitäten• durch strengen statistischen Test

• Steuerung der Punktauswahl durch Algorithmus

• Lokalisierung möglich!

• Bedingte Aussage über Zeitpunkt der Deformation

27

Sebastian Halsig

• Bedingte Aussage über Zeitpunkt der Deformation

– Ergebnis Ruhrgebiet: • flächenmäßig größte Deformation: 2002-2004

• mehr Punktidentitäten im Westen des Ruhrgebietes

• teilweise sprunghafte Senkungen im dm-Bereich

07.10.2010 Identifikation identischer Punkte in Höhennetzen mehrerer Epochen

Page 28: Identifikation identischer Punkte in Höhennetzen mehrerer ...–Rückführung auf Ausgleichung nach Bedingungen 14 Modell: Ausgleich nach Bedingungen mit folgender Bedingungsmatrix

Fazit und Ausblick

Ausblick:

– bisher noch keine Validierung • Möglichkeit über GPS-Permanentstationen

– Realisierung des bestmöglichen Festpunktnetzes

28

Sebastian Halsig

– Realisierung des bestmöglichen Festpunktnetzes über alle Punkte mit allen Epochen• unter Verwendung der Punktidentitäten

– Punktidentitäten • Einhängen in das Deutsche Gesamtnetz

07.10.2010 Identifikation identischer Punkte in Höhennetzen mehrerer Epochen

Page 29: Identifikation identischer Punkte in Höhennetzen mehrerer ...–Rückführung auf Ausgleichung nach Bedingungen 14 Modell: Ausgleich nach Bedingungen mit folgender Bedingungsmatrix

Vielen Dank29

Sebastian Halsig

für Ihre Aufmerksamkeit!

07.10.2010 Identifikation identischer Punkte in Höhennetzen mehrerer Epochen