wie GPS-Koordinaten (WGS84) nach Gauß-Krüger umrechnen?

[Husi]

Moin moin,

im Augenblick komme ich mal wieder nicht weiter und bräuchte eure Unterstützung.

Ich habe viele GPS-Logs, die sind im WGS84 Format (Beispiel: 50.10173°N ; 8.65043°E) gespeichert. Da ich die Logs nicht im Maps via *.kmz Datei betrachten möchte, sondern die Daten im CAD weitervarbeiten möchte, brauche ich sie in einem winkeltreuen Koordinatensystem. Meine erste Wahl wäre das Gauß-Krüger Koordinatensystem. Dort bzw. im UTM-System sähe die Koordinate (Beispiel von oben) so aus: 32U 475000 5550000
Also der Streifen Nr 32, ein Rechtswert von 475000 und ein Hochwert von 5550000

Bei Wiki findet man z.B. hier einiges: https://de.wikipedia.org/wiki/Referenze ... ellipsoide

@Tracer, das Bild ist nicht von Wiki, ich habe es selbst gepinselt...

Skizze_Koordinaten_umrechnen.gif (19.61 KiB) 4539 mal betrachtet

Ich schaffe es aber nicht von der Koordinate ( 50.10173°N ; 8.65043°E) auf (32U 475000 5550000) im z.B. Excel selbst zu kommen. Online-Rechen Tools helfen mir nicht, da ich wirklich viele Koordinaten habe, die ich umrechnen muß (und nicht noch einen weiteren Zwischenschritt beim Import der Logdaten haben möchte).

Im Anhang ist meine Excel-Datei (ich mußte die Dateiendung umbenennen, damit ich sie hochladen durfte) mit den entsprechenden Formeln. Wer findet den Fehler, bzw. kann mir erklären, was ich falsch mache?


Vielen Dank
Mirko

[dilg]

Hallo;


Eine schöne Kopfnuss für die Nacht.....


Ich hab da ein paar Fragen:
-von den Formeln sieht alles mathematisch korrekt aus, nur wie ist die Eingabe?
Das check ich noch nicht ganz.... für was ist Lambda? gibst du bei phi X oder Y ein? Sicher nur Y oder..?

Wie ist die Eingabe aufgebaut? kann es sein das du da einen Fehler drin hast weil die Grad nicht dezimal sind?

Die Formeln sind zwar schön und gut, aber ich sehe im Wiki nicht in welchen Format die Daten eingegeben werden sollen.....


Und dann hab ich noch drei Prinzipfragen:

-wie hast du im Excel die Variable phi global definiert?

-Warum das überhaupt?
ist eine Kugelform nicht genau genung, ich mein die Abplattung ist ja im Vergleich zum Erdraius lächerlich klein....

-Was hast du vor?

[Husi]

Moin michl,

danke für deine Gegenfragen.
Soweit ich das kann versuche ich mich mal an Antworten.

-wie hast du im Excel die Variable phi global definiert?

Das geht recht einfach, du markierst im Excel eine Zelle und trägst oben links den gewünschten Namen ein.

Excel_Zelle_benennen.jpg (68.45 KiB) 4458 mal betrachtet

-Warum das überhaupt?

Damit die Formeln "etwas" leichter zu lesen sind. Ich weiß, das hätte ich mit den anderen Zellen auch machen sollen...

ist eine Kugelform nicht genau genung, ich mein die Abplattung ist ja im Vergleich zum Erdraius lächerlich klein....

Das kann ich z.Z. noch nicht eindeutig beantworten. Solange ich aber Mist raus bekomme, sollte ich mich aber schon an die "Vorgaben" halten.

-Was hast du vor?

Das ist die schwierigste Frage. Ich möchte meine Log-Flüge im CAD als 2D und auch als 3D Flüge importieren und mich an den Buschreihen, Hochspannungsmasten etc. in der Flugplatznähe orientieren. Nebenbei könnte ich mit der selben Programm-Routine, Mittags bei Termikflügen einfach die Positionen der guten Thermik (mit ihrer Uhrzeit) in eine "Thermik-Karte" einzeichnen lassen.
Bzw. du schaust mal hier: https://www.rchelifan.org/viewtopic.p ... 2&t=101576

warum_das_alles_02.jpg (115.95 KiB) 4457 mal betrachtet

Blöder weise wird die Geometrie verzerrt, wenn man die WGS84 Koordinten als in ein rechwinkliges Koordinatensystem importiert.
Unser Modellflugplatz ist in Wirklichkeit nicht so ein starkes Parallelogramm, wie es oben im AutoCAD zu sehen ist, sondern eher ein Rechteck...
Wenn ich zusätzlich noch ein Screenshot von Maps unten drunter lege (aber ich hierfür mal aus geschaltet, ich weiß nicht wie das mit den Bildrechten aussieht) sieht man ganz deutlich die Verzerrung.

nur wie ist die Eingabe?

Jetzt verstehe ich deine Frage nicht. Meinst du die Eingabei in meiner Excel Datei, oder allgemein, bei der "Wiki-Anleitung"?
In der Excel-Datei sind es die drei grauen Zellen B49(=Phi) + B50(=Lamda) + B51(=h)

für was ist Lambda?

Ich habe da mal wieder ein Bildchen gemalt, ich denke daraus müßte es deutlich werden.

Was_ist_Lambda.jpg (135.73 KiB) 4460 mal betrachtet

Es ist der Winkel, der die Schnittposition in Bezug zu der Drehung um die Z-Achse festlegt.

gibst du bei phi X oder Y ein? Sicher nur Y oder..?

Siehe oben graue Zellen (nein nicht deine, oder meine grauen Zellen im Kopf)

Wie ist die Eingabe aufgebaut? kann es sein das du da einen Fehler drin hast weil die Grad nicht dezimal sind?

Das glaube ich nicht. Lägen die Winkel in diesem Grad Minuten Sekunden Format vor, dann würde ich sie auch erst einmal dezimal-Werte umrechnen, weil ich nicht weiß, wie ich den cos bzw. sin von diesen Minuten+Sekunden Winkeln ausrechnen soll.
dezimal Winkel =(Grad*3600+Minuten*60+Sekunden)/3600

Dies ist immer ein und die selbe Koordinate (siehe auch: https://tools.retorte.ch/map/?wgs84=50. ... map=hybrid):
Rechtwinklig:
Gauss-Krüger (DHDN/Bessel) => 3475065 5551781
UTM (WGS84) => 32U 475000 5550000

nicht Rechtwinklig:
Geogr. Koordinaten (WGS84)
Grad, Minuten, Sekunden (sexagesimal) =>50° 6´6.22´´ / 8° 39´1.53´´
Grad, Minuten Dezimal (nautisch) => 50°6.104´ / 8° 39.026´

Viele Grüße
Mirko

[CMaxX]

Moin ....!!!

Die Abplattung der Erde würde ich nicht gerade vernachlässigen ...!!! Sind immerhin ja rund 21km. Das wirkt sich doch recht stark auf die Projektion aus.

Mit was hast du denn die Daten aufgenommen? Und wieviele Koordinaten must du umrechnen?

Eigentlich bräuchtest du für die Transformation vom WGS ins RD83 Koordinatensystem Identische Koordinaten. Sprich, du brauchst Passpunkte (Punkte die in beiden Koordinatensystemen bekannt sind) am besten mit einer Überbestimmung (um eine Aussage über die Restklaffen (Genauigkeit) zu machen).

[dilg]

Hallo;


Also, die Frage nach dem Warum, ja da meinte ich warum die Umrechnung; du hast sowieso bei der besten Rechnung +/- 1m. Ob dann das Rechteck deines Flugplatzes stimmt wage ich zu bezweifeln.
Und selbst dann, könnstest mit dem Befehl Varia auch was gerade rücken.


Und die Frage nach Lamdba (oder eigentlich gamma) war, das ganze ist ja rotationssmetrisch...
Ich hab grad keinen Überblick zum Nachdenken, aber was gefunden.
Schau mal da rein
http://moenk.de/plugin/dlfile_13
is ein Ecxel, auf Tabelle 6 ist alles drin. ich frag mich warum und wieso da eine Bogenlänge drin is, vielleicht komm ich noch drauf. (bin halt doch kein Vermessungstechniker.....)

[Husi]

Hey Danke,

auf den ersten Blick ist die Excel Tabelle das, was ich gesucht habe. Jetzt muß ich es nur noch verstehen und "Nachprogrammieren".

Mit was hast du denn die Daten aufgenommen? Und wieviele Koordinaten must du umrechnen?

Die GPS-Koordinaten sind mit dem GPS-Logger 2 von SM-Modellbau aufgenommen, meine restlichen Daten mit allerlei dazu gestricktem.
Pro Sek. Flugzeit werden 10 Koordinaten gespeichert. Bei den normalen Logflügen sind das ca. 10 000 Koordinaten pro Flug. Davon mache ich in der Regel pro Meßtag vier Stück.

Danke
Mirko

[CMaxX]

GPS-Koordinaten sind mit dem GPS-Logger 2

Nur so als kleine Anmerkung: Du schaffst damit in der Regel ne Genauigkeit von 5 bis 10m. Das könnte auch zur Verzerrung der Daten beitragen.

[tracer]

@Tracer, das Bild ist nicht von Wiki, ich habe es selbst gepinselt...

Wäre in dem Fall auch OK gewesen, ist noch liberaler als eine CC-Lizenz:
"I grant anyone the right to use this work for any purpose, without any conditions, unless such conditions are required by law."

[Husi]

Sooo, fertig bin ich noch lange nicht, aber wieder mal einen Schritt weiter.

Mittlerweile habe ich durch die oben verlinkte Excel-Datei verstanden, das man zuerst die Geographischen Koordinaten in Geozentrische Koordinaten umrechtet. Das habe ich in meiner (ersten) Excel-Datei (im ersten Post) gemacht. Damit bin ich leider mit der Aufgabe noch nicht fertig. Denn die Ausgangskoordinaten sind ja in dem WGS84 Koordinatensystem gespeichert. Da ich aber ein anderes (Gaus-Krüger = Bessel-Koordinatensystem) haben möchte, müssen die ermittelsten X, Y, Z Koordinaten (Geozentrische Koordinaten) in das Ziel-Koordianntensystem (Bessel) transformiert werden.
Diese Transformation ist hier beschrieben:
https://de.wikipedia.org/wiki/Helmert-Transformation
Kurze Erklärung, die zwei Ellipsoide sind leicht gegeneinander verdreht, verschoben und skaliert. Das WGS84-Ellipsoid soll möglichst wenig Abweichung auf der ganzen Erde haben, während das Bessel-Ellipsoid möglichst wenig Abweichung in Deutschland haben soll. Ist halt schon doof, wenn die Erde eine so verbeulte Eierform hat.

Nachdem die Transformation gelaufen ist, können aus den neuen X, Y, Z-Koordinaten neue Geographische Koordinaten (Phi, Lambda, Höhe) berechnet werden. Diese werden dann im letzten Schritt zu den eingeebenten Koordinaten y,x,h (Rechts- und Hochwert und Höhe) gemacht.

Viele Grüße
Mirko

[echo.zulu]

Hallo Mirko.

Ist halt schon doof, wenn die Erde eine so verbeulte Eierform hat.

Dann nimm dir nen Spaten und fang schon mal an.

Ganz ehrlich bin ich immer wieder erstaunt was Du Dir so alles an Aufgaben auferlegst. Für mich ist dann immer die Frage, ob das Ergebnis den Aufwand rechtfertigt. Ich bin auf jeden Fall gespannt wohin Dich der Weg noch so führt.

[CMaxX]

Nachdem die Transformation gelaufen ist, können aus den neuen X, Y, Z-Koordinaten neue Geographische Koordinaten (Phi, Lambda, Höhe) berechnet werden. Diese werden dann im letzten Schritt zu den eingeebenten Koordinaten y,x,h (Rechts- und Hochwert und Höhe) gemacht.

Falls du nicht weiter kommst, sag bescheid. Da schaue ich mal in meinen Unterlagen nach.

Für was brauchst du sowas eigentlich?

[Husi]

Falls du nicht weiter kommst, sag bescheid. Da schaue ich mal in meinen Unterlagen nach.

Danke, aber jetzt bin ich damit schon fertig. Im Augenblick sitze ich noch daran, die Rechnung für große Koordinaten Arrays zu optimieren.

Für was brauchst du sowas eigentlich?

Ich möchte einen Screenshot von Google Maps hinter unseren Modellflugplatz legen. Wenn ich ihn in eine Richtung anpasse, dann stimmt die andere Richtung nicht. Siehe rote Pfeile.

AutoCAD_Koordinaten_umrechnen.jpg (444.98 KiB) 4233 mal betrachtet

Darum wollte ich die GPS-Koordinaten in ein Rechtwinkliges Koordinatensystem umrechnen, damit das zueinander passt.
Die Feldwege etc. passen exakt zu meinen Flügen. Denn ich bin sie mit den beiden GPS Loggern von SM-Modellbau abgefahren, da ist eine Abweichung von nur wenigen cm zu sehen.

Viele Grüße
Mirko

[Husi]

Für mich ist dann immer die Frage, ob das Ergebnis den Aufwand rechtfertigt.

Nö, aber ist das mit dem kompletten Hobby nicht so? Mit möglichst viel Aufwand und Geld, möglichst wenig erreichen...

Ich bin auf jeden Fall gespannt wohin Dich der Weg noch so führt.

Seit vielen Jahren habe ich zwei fest gesteckte Ziele. Alles was ich "zwischendurch" gemacht habe sind Zwischenstritte um die Ziele verwirklichen zu können.
Aber danke für die netten Worte!

Viele Grüße
Mirko

[echo.zulu]

Ziele sind im Leben immer wichtig. Oft sind die Wege dort hin verschlungen und man darf sie auch nicht aus den Augen verlieren. Aber die kleinen Abstecher machen es oft um so interessanter.

[Husi]

Es klappt

Ich habe mir gerade eine Kataster Import Routine geschrieben, die automatisch die Koordinaten auf Gauß-Krüger umrechnet.
Und siehe da, schon passt der Google-Maps Screnshot exakt zu meinen Katasterdaten.

es_klappt.jpg (453.64 KiB) 4161 mal betrachtet

Jetzt muß ich nur noch die Transformations-Routine an meinen Logdaten-Import ankoppeln...

Viele Grüße
Mirko