Hallo zusammen,

ich bin gerade dabei mir via OpenOffice Calc ein paar Makros zurecht zuschneidern die mir ein paar Arbeitsschritte erleichtern sollen.
Darunter fallen auch einige Datumshifts. Soweit bekomme ich auch alles hin und, dank diesem Forum und anderen Seiten sind kleiner Probleme auch schnell beseitigt.
Allerdings hänge ich jetzt in den Seilen: die Berechnung von GK nach UTM läuft bis zu dem Punkt an dem die Transformation durchgeführt werden soll wie sie soll. Die kartesischen Koordinaten auf Bessel werden korrekt berechnet aber nach der Transformation sind die Koordinaten auf GRS80 um ca. 50m falsch.
Sowohl Formeln als auch Transformationsparameter habe ich von der CRSeu.
In Calc implementiert sieht das dann so aus:

Xgrs = Xbes + TX + ( (-1)*RZ*Ybes + RY*Zbes ) + SC*Xbes
Ygrs = Ybes + TY + ( RZ*Xbes + (-1)*RX*Zbes ) + SC*Ybes
Zgrs = Zbes + TZ + ( (-1)*RY*Xbes + RX*Ybes ) + SC*Zbes

Das eigenartige ist, dass die Kontrolldaten die auf der Seite der CRSeu korrekt umgewandelt werden. Andere Daten (z.B. habe ich noch welche von der AdV gefunden) werden eben mit dieser enormen Abweichen umgewandelt. Die Transformationsparameter habe ich mehrmals überprüft. Sie sind korrekt eingegeben und werden auch richtig in das Makro importiert.

Hat jemand ein Ahnung wo mein Fehler liegt?

Hallo Tom,

Allerdings hänge ich jetzt in den Seilen: die Berechnung von GK nach UTM läuft bis zu dem Punkt an dem die Transformation durchgeführt werden soll wie sie soll.

Kannst Du ein Zahlenbeispiel hier posten? Welchen Wert setzt Du für Z ein beim BKG-Beispieldatensatz ein?

Die Transformationsparameter habe ich mehrmals überprüft. Sie sind korrekt eingegeben und werden auch richtig in das Makro importiert.

Wenn ich einen Beispieldatensatz habe, würde ich die Formeln mal abtippen und Dir zumindest bestätigen, dass Du hier keinen Fehler gemacht hast.

Hat jemand ein Ahnung wo mein Fehler liegt?

Im Moment noch nicht.

Gruß Micha

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

Tags:
Gauß-Krüger, Transformation, UTM

Hi MichaeL,

für die ellipsoidische Höhe verwende ich immer 0.000 sofern ich keine Daten habe.
Zahlenbeispiele ellipsoidische Koordianten:

Bessel                      GRS80                       GRS80
B             L             B-Soll        L-Soll        B-Ist         L-Ist
48.18790000   8.27400000    48.18695833   8.27307777    48.18695793   8.27307781
53.50000000   7.48333333    53.49846114   7.48250602    53.49844824   7.48250197


Zahlenbeispiele kartesische Koordianten:

Bessel
Soll            Ist           
3794219.069     3794219.070
548024.691      548024.691
5079507.831     5079507.831

GRS80
Soll            Ist
3794799.376     3794850.222
548044.577      548051.928
5079891.705     5079959.773


Hallo Tom,

Zahlenbeispiele kartesische Koordianten:

Wo hast Du diese Zahlen her? Die Parameter, die auf der BKG-Seite angegeben werden, scheinen richtig zu sein. Auch bei Wikipedia sind diese so zu finden. Ich habe von Deinem Beispieldaten die Bessel (SOLL) Koordinaten nach GRS80 überführt und komme zu folgenden Zahlen:

3794850.22113067
548051.928821092
5079959.7726258

Das würde also Deine Berechnung bestätigen. Mehr kann ich dazu leider nicht sagen. Im Moment würde ich vermuten, dass der Beispieldatensatz nicht das ist, was Du vermutest.

Gruß Micha

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

Tags:
Gauß-Krüger, Transformation, UTM

Hallo,

weil ich grad hier bin.
Hab' reingeschaut und will morgen mal nachrechnen.

Roland

Den Datensatz habe ich aus einer - wie sich jetzt heraus stellt - nicht wirklich zuverlässigen Quelle.
Was ich jetzt bemerkt habe ist, dass auf der Seite der CRSeu es auch die Möglichkeit einer Online-Transformation gibt. Dort werde ich mir mal ein paar Referenzdaten besorgen. Damit denke ich bin ich auf der sicheren Seite.
Vielen Dank für deine Mühen.

Hallo,

mit 'Alltrans' kann ich die Bessel-GRS80 Transformation bestätigen und auch stichprobenweise die geogr. Koordinatenumwandlungen.

Mir kam aber eine Idee:
Wenn ich einen Punkt des Bessel-Ellipsoids mit der Höhe 0 aufs GRS80 mit Helmert transformiere, liegt der Punkt ca. 40-50 m über dem GRS80. Bei entsprechenden Rücktransformationen darf ich nicht 0 für die ell. Höhe ansetzen.

Der Punkt B= 48.18°, L= 8.27° hat eine ell. Höhe von 53.5 Metern, der Punkt bei 53.5° und 7.48° von 40,3 Metern ...

Grüße Roland