Koordinatenumrechnung: Gauß Krüger => Dezimalgrad

OP #3318202
Lesenswert?

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
#3318319
Lesenswert?

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...
Gast #3318811
Lesenswert?

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.
OP #3318833
Lesenswert?

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.
#3318916
Lesenswert?

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.
OP #3319313
Lesenswert?

Wie versprochen hier mal die Berechnung als VBA-Code:
1
Private Sub GK_in_Dezimalgrad_umrechnen(rechts As Double, hoch As Double)
2

3
    Pi = 4 * Atn(1)
4
    rho = 180 / Pi
5
    
6
    e2 = 0.0067192188
7
    c = 6398786.849
8
    
9
    mKen = WorksheetFunction.RoundDown(rechts / 1000000, 0)
10
                
11
    rm = rechts - mKen * 1000000 - 500000
12
    
13
    bl = hoch / 10000855.7646
14
    
15
    bll = bl * bl
16
    
17
    bf_1 = 325632.08677 * bl * ((((((0.00000562025 * bll - 0.0000436398) * bll + 0.00022976983) * bll - 0.00113566119) * bll + 0.00424914906) * bll - 0.00831729565) * bll + 1)
18
                
19
    bf = bf_1 / 3600 / rho
20
                
21
    co = Cos(bf)
22
                
23
    g2 = e2 * co * co
24

25
    gl_1 = c / Sqr(1 + g2)
26
        
27
    t = Sin(bf) / Cos(bf)
28
                
29
    fa = rm / gl_1
30
                
31
    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)
32

33

34
'=================================
35
    lat = gb_1 * rho
36
'=================================
37
    
38
    
39
    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)
40
         
41
         
42
'=================================
43
    lon = (dl * rho / co) + (mKen * 3)
44
'=================================
45
    
46
    
47
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.
#3320133
Lesenswert?

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/

Antwort schreiben

Bitte melde dich an, um einen Beitrag zu schreiben.

oder

Mit Google-Account einloggen

Die Registrierung ist kostenlos und dauert nur eine Minute.

Jetzt registrieren