Beobachtete Koordinaten einführen

by PeterW, Friday, March 22, 2019, 08:17 (32 days ago)

Hallo Micha,

ich hätte da mal ein Sonderproblem...
Ich möchte zwei freie Netze (Epo0 und Epo1) mittels einer zweistufigen Defo-Analyse prüfen, wobei ich entsprechend Stütz- und Objektpunkte habe. Nun sind einige meiner Stützpunkte auch nicht wirklich stabil, ich kenne aber das Bewegungsverhalten einzelner Punkte, d.h. ich kann vorausberechnen, wo sich manche Stützpunkte zum Zeitpunkt der Epo1 befinden werden. Dieses Wissen würde ich gerne in der Defo-Analyse nutzen, um die Punkte weiterhin als Stützpunkte nutzen zu können (d.h. mein Hypothesentest sollte die Abweichung der Epo1 von der prognostizierten Bewegung zwischen Epo0 und Epo1 testen).
Da ich kein Ausgleichungsprogramm kenne, welches zeitabhängige Punktbewegungen mit modellieren kann, hatte ich an folgende Vorgehensweise gedacht:

  • Ich berechne die Koordinaten der Punkte mit bekannter Bewegungscharakteristik extern und füge sie dann der Epo0 als „beobachtete Koordinaten“ mit entsprechender Standardabweichung und neuer Punktnummer hinzu. Diese Punkte sind zwar mangels verknüpfter Beobachtungen nicht kontrollierbar, aber das stört mich ja nicht - Hauptsache, sie verfügen über eine stochastische Beschreibung.
  • In der Epo1 ersetze ich die Punktnummern der tatsächlichen Beobachtungen durch diejenigen meiner fiktiven Punkte, so dass ein ganz normales freies Netz berechnet werden kann.
  • Jetzt sollte, da für beide Epochen ein vollständiges stochastisches Modell vorliegt und meine bewegten Punkte auch in Epo0 enthalten sind, eine Defo mit dem Standard-Algorithmus möglich sein.

Beobachtete Koordinaten in das AGL-Modell einzuführen, ist ja eigentlich trivial - trotzdem kenne ich kein Ausgleichungsprogramm, welches das unterstützt...
Ist dieses Vorgehen mit JAG3D realisierbar, bzw. hast Du eine bessere Idee für eine Lösung?

Schöne Grüße
Peter

Beobachtete Koordinaten einführen

by Micha ⌂, Friday, March 22, 2019, 11:30 (32 days ago) @ PeterW

Hallo Peter,

Da ich kein Ausgleichungsprogramm kenne, welches zeitabhängige Punktbewegungen mit modellieren kann, hatte ich an folgende Vorgehensweise gedacht:

Du kennst ja die absolute Position nicht, sondern modellierst nur die Verschiebung zwischen den Epochen. Du bestimmst aus Deinem Bewegungsmodell also nur δX und δY. Diese Verschiebungsvektoren würde ich entweder in eine polare Darstellung konvertieren (Richtungen und Strecken) oder direkt als relative GNSS-Basislinie in JAG3D importieren. Letzteres erscheint mir am einfachsten. Das stochastische Modell kannst Du dann vermutlich direkt Deinem externen Bewegungsmodell entnehmen.

Beobachtete Koordinaten in das AGL-Modell einzuführen, ist ja eigentlich trivial - trotzdem kenne ich kein Ausgleichungsprogramm, welches das unterstützt...

In JAG3D heißen diese Punkte Anschlußpunkte. Im Kontext der Deformationsanalyse werden diese Punktbeobachtungen aber nicht unterstützt, da ein Netz mit Anschlußpunkten kein freies Netz mehr darstellt. Anschlußpunkte werden in JAG3D aber im Rahmen einer klassischen Ausgleichung auf Veränderung hin getestet.

Ist dieses Vorgehen mit JAG3D realisierbar, bzw. hast Du eine bessere Idee für eine Lösung?

Ich würde vermutlich die erste Variante versuchen, da Du ein Bewegungsmodell hast. Unter einem solchen Modell stelle ich mir eine relative Information vor, sodass sie von der Charakteristik der relativen GNSS-Beobachtung entspricht. Ich hoffe, das hilft Dir bei Deiner Überlegung weiter.

Viele Grüße
Micha

--
applied-geodesy.org - OpenSource Least-Squares Adjustment Software for Geodetic Sciences

Tags:
JAG3D, datum, deformation analysis, time-dependent deformation

RSS Feed of thread