Hallo,

so, der Hinweis mit Deinem c=0 war der entscheidene, ich konnte diese Angabe sonst nirgens finden.

Der Hinweg geht nun bei mir. Ich habe dazu zwei kleine Klassen geschrieben, wobei der Rückweg - wie man sich denken kann - noch fehlt.

Die Handhabung sollte recht einfach sein, für Deine Beispiele:

    CoordsConverter CC = null;

    Ellipsoid Bessel = new Ellipsoid("Bessel", 3, false);
    CC = new CoordsConverter(4581441.0, 5936991.0, Bessel);
    System.out.println(CC.getLatitude()+"  "+CC.getLongitude());
    
    Ellipsoid Krassowski = new Ellipsoid("Krassowski", 6, false);
    CC = new CoordsConverter(3382537.0, 5938294.0, Krassowski);
    System.out.println(CC.getLatitude()+"  "+CC.getLongitude());
    
    Ellipsoid GRS80 = new Ellipsoid("GRS80", 6, true);
    CC = new CoordsConverter(33382586.0, 5935821.0, GRS80);
    System.out.println(CC.getLatitude()+"  "+CC.getLongitude());

3° Potsdam-Datum (= Bessel-Ellipsoid):
GK ; geographisch (Soll laut Top50) ; geographisch (Ist laut meinem TR)
4581441 ; 13,22919° ; 13,22918°
5936991 ; 53,55962° ; 53,55962°

6° 42/83 (= Krassowski-Ellipsoid):
GK ; geographisch (Soll laut Top50) ; geographisch (Ist laut meinem TR)
3382537 ; 13,22741° ; 13,22740°
5938294 ; 53,55806° ; 53,55806°

WGS84/ETRS89-Ellipsoid:
UTM ; geographisch (Soll laut Top50) ; geographisch (Ist laut meinem TR)
33382586 ; 13,22741° ; 13,22741°
5935821 ; 53,55812° ; 53,55812°

Dabei erhalte ich:

53.559616969942574  13.229180899939035
53.55806194137487   13.227411172382293
53.55811998388578   13.227409699981028

Sieht so schlecht nicht aus

Formeln habe ich aus dem [FRÖHLICH/KÖRNER, 2001] genommen, da deutlich Übersichtlicher, wie Du Dir denken kannst.

Micha

--
applied-geodesy.org - OpenSource Least-Squares Adjustment Software for Geodetic Sciences