[C# => Delphi] Geokoordinaten nach Mercator

[Mithrandir]

Hi ihr,

Eigentlich sollte es nicht zu schwer sein, aber ich bekomme es trotzdem nicht hin...

Ich habe mir von hier den Code geschnappt, um von Geokoordinaten auf die Mercator-Projektion zu kommen. Dieses Listing habe ich nach Delphi umgeschrieben.

Folgendes ist gegeben:

Latitude: 53.5582°
Longitude: 7.9479°

Das soll in die Mercator-Projektion umgewandelt werden. Hier gibt es einen Konverter: http://mal2.ch/stuff/latlontomercator.html

Heraus kommt:

Latitude(Y): 7086925;
Longitude(X): 884756;

Ich tippe einfach mal, dass das stimmt. Und jetzt komm ich:

Latitude(Y): -10266881
Longitude(X): 884756

Wie man unschwer erkennen kann, geht bei der Berechnung der Latitude was in die Hose.

Aber wo?

[Mithrandir]

Oookeee...

Also, ich hatte einen Zeilendreher drin. In der orpMercator.pas finden sich in der Funktion LatToY die folgenden Zeilen:

Delphi-Quellcode:

  ts := Tan(0.5 * ((PI * 0.5) - phi)) / con;
con := Power(((1.0 - con) / (1.0 + con)), fCom);

diese müssen vertauscht werden:

Delphi-Quellcode:

con := Power(((1.0 - con) / (1.0 + con)), fCom);
ts := Tan(0.5 * ((PI * 0.5) - phi)) / con;

Dann erhält man folgendes:

Latitude(Y): 7052526
Longitude(X): 884756

Wie man erkennt, passt der Y-Wert immer noch nicht. Ich habe deshalb testweise mal den PHP-Code laufen lassen:

Code:

Array
(
    [x] => 884756.180876
    [y] => 7052526.40477
)

Wie zu erkennen ist, liefert diese ähnliche Werte. Es ist also davon auszugehen, dass die Codes stimmen, und das es einfach daran liegt, dass die Algorithmen genauer sind.