Koordinatenumrechnung: Gauß Krüger => Dezimalgrad

[Matt B. (mattb)]

Hallo,

ich brauche eine Formel um Gauß Krüger Koordinaten in Dezimalgrad 
umzurechnen. Da ich 3000 Koordinatensätze umrechnen muss geht es nicht 
manuell mit einem Onlinerechner. Wahrscheinlich hab ich mich bei der 
Suche im Netz zu dumm angestellt und habe deswegen nichts gefunden... 
:-/

Danke schonmal
Gruss
Matt

[Purzel H. (hacky)]

Was spricht denn gegen ein kleines Prograemmchen? Oder eine Formel im 
Excel ?

[Bernd S. (bernds1)]

Siebzehn mal Fuenfzehn schrieb:
> Was spricht denn gegen ein kleines Prograemmchen? Oder eine Formel
> im Excel ?

Gar nichts, deswegen sucht er ja auch genau diese Formel.

Matt B. schrieb:
> ich brauche eine Formel um Gauß Krüger Koordinaten in Dezimalgrad
> umzurechnen.

Erna sagt zu Max: "Ich brauche ein neues Kleid!"
Sagt Max zu Erna: "Was spricht gegen ein neues Kleid?"

Erst lesen, dann posten! Steht sogar als Überschrift beim Posting 
verfassen oben drüber...

[Bernd S. (bernds1)]

Das einzige, was ich anbieten kann ist dieses PDF. Falls dir das was 
nützt. Ich habe es vor längerer Zeit mal runtergeladen, aber selbst noch 
nichts damit gemacht.

Gruß Bernd

[Purzel H. (hacky)]

Passt das etwa?

>http://en.wikipedia.org/wiki/Transverse_Mercator

[Timm T. (Gast)]

Bernd S. schrieb:
> Gar nichts, deswegen sucht er ja auch genau diese Formel.

Gibt es nicht.

Da Gauss-Krüger auf jeweils 3° breiten Meridanstreifen entlang von 
Mittelmeridianen beruht, muss man erstmal schauen, auf welchen 
Mittelmeridian sich die Koordianten beziehen. Dann kann man die 
Koordinaten zerlegen und von diesem Mittelmeridian aus umrechnen. Das 
kann freilich ein Programm machen, aber eine Formel, die man einfach in 
Excel klopft, wird das nicht.

[Matt B. (mattb)]

Ok, dann hab ich mich bei der Suche nicht zu dumm angestellt ! ;-)

Ich habe aber trotzdem was gefunden:
www.kram.de/pca/koordinatenumrechnung.xls

Volker Kram hat hier eine Excel-Datei veröffentlicht, in der er mittels 
Funktionen im Arbeitsblatt Gauß Krüger Koordinaten umgerechnet hat.

Anhand dieser Datei habe ich mir eine Funktion zusammengebastelt. Damit 
waren die 3000 Datensätze in 2 Sekunden umgerechnet. Hat sich also 
gelohnt!

Nachdem ich die Koordinaten in Google Maps eingetragen hatte, bemerkte 
ich, dass die Punkte um einen konstanten Offset zur Karte verschoben 
sind. Dem muss ich noch nachgehen. Vielleicht stimmt etwas an meiner 
Formel nicht.

Die Funktion selbst stelle ich hier in den nächsten Tagen rein. Erklären 
werde ich sie aber nicht.

[Bernd W. (berndwiebus)]

Hallo Bernd.

Bernd S. schrieb:
> Das einzige, was ich anbieten kann ist dieses PDF. Falls dir das was
> nützt. Ich habe es vor längerer Zeit mal runtergeladen, aber selbst noch
> nichts damit gemacht.

Herzlichen Dank! Sowas habe ich selber auch gesucht.

Mit freundlichem Gruß: Bernd Wiebus alias dl1eic
http://www.dl0dg.de

[Uhu U. (uhu)]

Matt B. schrieb:
> Nachdem ich die Koordinaten in Google Maps eingetragen hatte, bemerkte
> ich, dass die Punkte um einen konstanten Offset zur Karte verschoben
> sind. Dem muss ich noch nachgehen. Vielleicht stimmt etwas an meiner
> Formel nicht.

Die Umrechnung von Gauss-Krüger in irgend was "normales" ist offenbar 
auch in ausgewachsenen GIS-Programmen ziemlich ungenau. Ich habe mal ca. 
100 umgerechnete Koordinaten "repariert", indem ich die Punkte in Google 
Maps mit verschieblichen Markern angezeigt, zurechtgeschoben und die 
Koordinaten dann exportiert habe.

[Arno H. (arno_h)]

Die Formeln im EXCEL-Sheet beziehen sich auf ein Buch von 1964, GPS und 
neuere Karten  benutzen das WGS84-System.
Check doch mal mit was Neuerem wie z.B. dem Konverter auf 
http://sky-maps.de/wgs84-to-gk-gk-to-wgs84/ gegen.

Arno

[Matt B. (mattb)]

Wie versprochen hier mal die Berechnung als VBA-Code:

Private Sub GK_in_Dezimalgrad_umrechnen(rechts As Double, hoch As Double)

    Pi = 4 * Atn(1)
    rho = 180 / Pi
    
    e2 = 0.0067192188
    c = 6398786.849
    
    mKen = WorksheetFunction.RoundDown(rechts / 1000000, 0)
                
    rm = rechts - mKen * 1000000 - 500000
    
    bl = hoch / 10000855.7646
    
    bll = bl * bl
    
    bf_1 = 325632.08677 * bl * ((((((0.00000562025 * bll - 0.0000436398) * bll + 0.00022976983) * bll - 0.00113566119) * bll + 0.00424914906) * bll - 0.00831729565) * bll + 1)
                
    bf = bf_1 / 3600 / rho
                
    co = Cos(bf)
                
    g2 = e2 * co * co

    gl_1 = c / Sqr(1 + g2)
        
    t = Sin(bf) / Cos(bf)
                
    fa = rm / gl_1
                
    gb_1 = bf - (fa * fa * t * (1 + g2) / 2) + (fa * fa * fa * fa * t * (5 + (3 * t * t) + (6 * g2) - (6 * g2 * t * t)) / 24)


'=================================
    lat = gb_1 * rho
'=================================
    
    
    dl = fa - (fa * fa * fa * (1 + (2 * t * t) + g2) / 6) + (fa * fa * fa * fa * fa * (1 + (28 * t * t) + (24 * t * t * t * t)) / 120)
         
         
'=================================
    lon = (dl * rho / co) + (mKen * 3)
'=================================
    
    
End Sub

In dem Code werden die Lat & Lon Werte nicht weiter verarbeitet. Aber 
das sollte jeder selbst hinbekommen.


Ich habe allerdings bemerkt, dass alle damit berechneten Punkt um einen 
konstanten Offset verschoben sind. Den Code habe ich noch nicht 
kontrolliert. Vielleicht habe ich darin einen Fehler.

[Ronald R. (hidden_impact)]

Bei der Umrechnung musst du noch folgendes beachten. Die GK-Koordinaten 
basieren (auf deutschen Karten) auf dem Bessel-Ellipsoid. Zusätzlich 
wird dieser noch gekippt und verschoben. Das nennt sich dann 
Potsdam-Datum. Um also in Google-Maps (WGS84) etwas sinnvolles 
anzuzeigen, musst du diese Transformationen ausführen. Infos gibt z.B. 
es auf:

http://www.ottmarlabonde.de/