Bis jetzt nur http://www.delphi-treff.de/tipps/mathematik/wiki/Geographische in Gauß-Krüger-Koordinaten umrechnen/
Da steht
c := 6398786.849;
das ist aber weder der terrestrische
Polradius von 6735 km = 6735000 m,
noch der Äquatorradius von 6378 km = 6378000 m
Der Wert in der Delphi-Rechnung weicht also mindestens um ca 0.33% vom Realwert ab.
---
Wenn ich das Verfahren richtig deute, wird der Kugelstreifen auf einen Zylinder, und dieser auf eine Ebene (linkshändig x-y) projeziert.
Ich frag mich grade wo eigentlich der Unterschied zu einer direkten Projektion des Kugelstreifens auf die x-y Ebene ist und seh grade keinen.
Bei direkter Projektion wären die Formeln einfach:
Höhenwert y = (sin(Breitengrad))*Radius;
Rechtswert x = (sin(Längengrad-Längengrad(Mittelmeridian))*Radius
Vielleicht mal mit Radius = 6378000 m für einen Punkt rund um Greenwich testen.
---
Es mag hilfreich sein, für Überschlagsrechnungen in Radians statt Grad zu rechnen
rad(degree) = degree/(360°/(2*PI)) = (degree/180°) * PI
da bei kleinen Winkeln (hier bei den Längengraden L mit |L| <= 1.5°) gilt
sin(rad) "ungefähr =" rad.
Grüsse
Solkar