hybride Ausgleichung mit GNSS und terrestrischen Messungen

by htw9056, (8 days ago)
edited by htw9056,

Hallo Micha,

wir hatten schon einmal das Thema in 2023 und ich bin nach Jahren der Untätigkeit hier wieder dran.

Kannst du mir das Beispiel von Wetzel zukommen lassen, dass ich mir das noch einmal näher betrachten kann?

Gern als fertiges Jag3D Projekt und auch die Rohmessdaten und GNSS Fest-/Neupunkte.

Noch eine weitere Frage:
Wie werden die Basislinien verarbeitet? In WGS84 DeltaY,DeltaX,DeltaZ oder als lokal kartesische Koordinaten Nord,Ost,Höhe?

Bsp. WGS84 DeltaY,DeltaX,DeltaZ:
0082 0101 -20206.6154 18466.4881 -10908.0587
0082 0101 -20206.6140 18466.4930 -10908.0494
0082 0101 -20206.6143 18466.4906 -10908.0561
0082 0101 -20206.6165 18466.4898 -10908.0568

Bsp. lokal kartesische Koordinaten:
0082 0101 -17534.0441 -23682.5053 -31.1746
0082 0101 -17534.0385 -23682.5009 -31.1831
0082 0101 -17534.0328 -23682.5009 -31.1806
0082 0101 -17534.0249 -23682.5110 -31.1401

LG
Stefan

hybride Ausgleichung mit GNSS und terrestrischen Messungen

by Micha ⌂, Bad Vilbel, (8 days ago) @ htw9056

Hallo,

Kannst du mir das Beispiel von Wetzel zukommen lassen, dass ich mir das noch einmal näher betrachten kann?

Die Daten, die wir für unsere Untersuchung und die Veröffentlichung genutzt haben, sind online verfügbar. Du kannst diese direkt herunterladen.

Noch eine weitere Frage:
Wie werden die Basislinien verarbeitet? In WGS84 DeltaY,DeltaX,DeltaZ oder als lokal kartesische Koordinaten Nord,Ost,Höhe?

Das ist fast egal: Es werden Koordinatenunterschiede verarbeitet. Ich persönlich würde eine Darstellung wählen, die dicht am gewählten Datum liegt. Das erleichtert meiner Meinung nach die Interpretation erheblich. Wenn es starke Unterschiede gibt, müsstest Du u.U. geeignete Drehwinkel vorgeben als Startwerte. JAG3D wird als Startwert per Default von Null ausgehen.

Viele Grüße
Micha

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

hybride Ausgleichung mit GNSS und terrestrischen Messungen

by htw9056, (5 days ago) @ Micha

Hallo Micha,

vielen Dank für die schnelle Antwort und den Hinweis auf die Daten. Ich habe inzwischen selbst nachgerechnet und möchte mein Anliegen genauer erklären, weil die Festpunktanalyse bei mir nicht plausibel wirkt.

Was ich gemacht habe:
Die WGS84-Basislinien (Reihenfolge dY, dX, dZ) rechne ich mit derselben Drehmatrix um, mit der die Festpunkte ins lokal ellipsoidische System gebracht wurden (Fundamentalpunkt B = 51,3539773354°, L = 12,3740940991°), nur ohne Verschiebung. Ergebnis sind dE, dN, dh. Kontrolle: Diese Basislinien stimmen mit den Koordinatendifferenzen der lokalen Festpunkte auf unter 0,1 mm überein. Die Höhen der Festpunkte unterscheiden sich um höchstens 203,86 m, und die lokalen dh der Basislinien liegen zwischen -203,9 m und +184,7 m. Das passt also zusammen.

Zusätzlich habe ich die Mittelwerte der Basislinien zwischen den 8 Festpunkten (14 Paare, 24 bis 56 km) mit den Festpunktkoordinaten verglichen. Maßstab und Drehungen sind praktisch null (Maßstab 0,01 ppm, Drehungen unter 0,1 Bogensekunden). Die Restabweichungen liegen bei etwa 6 mm (Ost), 10 mm (Nord) und 17 mm (Höhe). Die einzelnen Epochen streuen in der Höhe mit etwa 25 mm, in der Lage mit 8 und 11 mm.

Warum ich deine Aussage "fast egal" nicht ganz teile:
Für die Koordinatendifferenzen stimme ich dir zu. Aber nach dem Wiki geht in das stochastische Modell der Wert der jeweiligen Vektorkomponente als Abstand d ein, und die Kovarianzmatrix der Basislinie ist diagonal. Dann ist die Darstellung für die Gewichte nicht egal:

1. Im lokalen System ist dh nur Meter bis Dezimeter groß, dE und dN sind dagegen mehrere Kilometer. Bei einem entfernungsabhängigen Anteil bekommt die Höhe damit eine deutlich kleinere Standardabweichung als Ost und Nord, obwohl GNSS die Höhe in Wirklichkeit schlechter bestimmt. In WGS84 (dX, dY, dZ) ist die Verteilung eine ganz andere. Verstehe ich das richtig?

2. Die Korrelationen zwischen dE, dN und dh aus der GNSS-Auswertung gehen dabei verloren. Gibt es eine Möglichkeit, die volle 3x3-Kovarianzmatrix je Basislinie zu verwenden, oder sollte man die Genauigkeiten besser einzeln vorgeben?

3. Meine Datei enthält pro Basislinie etwa 1080 Zeilen (Epochen). Werden diese als unabhängige Beobachtungen behandelt, wird die Netzgenauigkeit meiner Ansicht nach viel zu optimistisch. Würdest du die Epochen vorher zu einem Mittelwert mit realistischer Standardabweichung zusammenfassen?

4. Mir ist aufgefallen, dass in meinem ersten Beitrag das Beispiel "lokal kartesisch" zu einem ganz anderen Fundamentalpunkt gehört (rechnerisch etwa B = 51,63°, L = 12,32°, außerdem in der Reihenfolge Nord, Ost, Höhe). Dadurch unterscheiden sich dE, dN und dh um 12 m, 16 m und 72 m. Die Drehung muss also zum selben Fundamentalpunkt wie die Koordinaten passen. Das ist bei meinen aktuellen Daten der Fall.

Was die Festpunktanalyse bei mir zeigt:
[HIER ERGEBNIS EINFÜGEN: welche Festpunkte werden verworfen, in Lage oder Höhe, mit welchen Testwerten und welche Einstellungen (Maßstab/Drehwinkel mitschätzen ja/nein, Genauigkeiten der Basislinien).]

Kannst du mir sagen, wie du die Genauigkeiten der Basislinien (Punkte 1 bis 3) in so einem Fall am besten ansetzen würdest? Ich vermute, dass hier die Ursache für die auffälligen Ergebnisse liegt.

Viele Grüße
Stefan

hybride Ausgleichung mit GNSS und terrestrischen Messungen

by Micha ⌂, Bad Vilbel, (5 days ago) @ htw9056

Hallo,

Die WGS84-Basislinien (Reihenfolge dY, dX, dZ) rechne ich mit derselben Drehmatrix um, mit der die Festpunkte ins lokal ellipsoidische System gebracht wurden (Fundamentalpunkt B = 51,3539773354°, L = 12,3740940991°), nur ohne Verschiebung. Ergebnis sind dE, dN, dh. Kontrolle: Diese Basislinien stimmen mit den Koordinatendifferenzen der lokalen Festpunkte auf unter 0,1 mm überein. Die Höhen der Festpunkte unterscheiden sich um höchstens 203,86 m, und die lokalen dh der Basislinien liegen zwischen -203,9 m und +184,7 m. Das passt also zusammen.

Bitte beachte, dass der Fundamentalpunkt - sofern Du nichts umstellst - in Neugrad (gon) einzugeben ist.

1. Im lokalen System ist dh nur Meter bis Dezimeter groß, dE und dN sind dagegen mehrere Kilometer. Bei einem entfernungsabhängigen Anteil bekommt die Höhe damit eine deutlich kleinere Standardabweichung als Ost und Nord, obwohl GNSS die Höhe in Wirklichkeit schlechter bestimmt. In WGS84 (dX, dY, dZ) ist die Verteilung eine ganz andere. Verstehe ich das richtig?

Bei Basislinien wird die Komponentenlänge als Abstand verwendet und nicht die Norm des Vektors selbst, um diese Missverhältnisse in den Längen entgegenzuwirken. Im Allgemeinen steht es Dir aber frei, individuelle Unsicherheiten zu spezifizieren, wenn Du mit einem implementierten Ansatz unzufrieden bist.
Im Wiki heißt es: "Der Abstand d zur Bildung des entfernungsabhängigen Anteils im stochastischen Modell entspricht hierbei der zugehörigen Basislinenvektorkomponente."

2. Die Korrelationen zwischen dE, dN und dh aus der GNSS-Auswertung gehen dabei verloren. Gibt es eine Möglichkeit, die volle 3x3-Kovarianzmatrix je Basislinie zu verwenden [...]?

Nein, gegenwärtig nicht. Man könnte das System aber so drehen, dass keine Korrelationen mehr auftreten. Den Aufwand würde ich aber nicht treiben.

3. Meine Datei enthält pro Basislinie etwa 1080 Zeilen (Epochen). Werden diese als unabhängige Beobachtungen behandelt, wird die Netzgenauigkeit meiner Ansicht nach viel zu optimistisch. Würdest du die Epochen vorher zu einem Mittelwert mit realistischer Standardabweichung zusammenfassen?

GNSS kommt in meinem beruflichen Alltag praktisch nicht vor, aber prinzipiell handhabe ich das bei allen Beobachtungen so, ja.

Ich hoffe, dass hilft Dir weiter.

Viele Grüße
Micha

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

Tags:
JAG3D, WGS84, GNSS, Ausgleichung, Fundamentalpunkt

RSS Feed of thread