Hallo Thomas,

Zunächst konvertiere ich die UTM-Koordinaten zu geographischen Koordinaten und diese dann ins Gauß-Krüger-System.

Wie realisierst Du den Ellipoidübergang?

Ich habe mit Deinen Werten mal schnell die Umrechnung in GeoTra durchgeführt. Ich komme hier zumindest auf Deine Soll-Werte. Die Methode zum Umformen sieht so aus:

 
 private Point2D BL2XY(double phi, double dLambda) {
  double a  = this.ellipsoid.getMajorAxis();
  double b  = this.ellipsoid.getMinorAxis();
  double e2 = this.ellipsoid.getE2Square();
 
  // Meridianbogelnlaenge
  double Bphi = this.ellipsoid.getMeridianArcLengthFromLatitude(phi);
 
  double eta2 = e2*Math.cos(phi)*Math.cos(phi);
  double t    = Math.tan(phi);
  double t2   = t*t;
 
  // Normalkruemmungsradius fuer phi
  double N = a*a/(b*Math.sqrt(1.0+eta2));
 
  double x = Bphi + t/2.0*N*Math.pow(Math.cos(phi),2)*Math.pow(dLambda,2) +
             t/   24.0*N*Math.pow(Math.cos(phi),4)*(   5.0 -        t2 + 9.0*eta2+4.0*eta2*eta2)*Math.pow(dLambda,4) +
             t/  720.0*N*Math.pow(Math.cos(phi),6)*(  61.0 -   58.0*t2 + t2*t2 + 270.0*eta2 - 330*t2*eta2)*Math.pow(dLambda,6) +
             t/40320.0*N*Math.pow(Math.cos(phi),8)*(1385.0 - 3111.0*t2 + 543.0*t2*t2 - Math.pow(t,6))*Math.pow(dLambda,8);
 
  double y = N*Math.cos(phi)*dLambda +
             1.0/   6.0*N*Math.pow(Math.cos(phi),3)*( 1.0 -       t2 + eta2)*Math.pow(dLambda,3) +
             1.0/ 120.0*N*Math.pow(Math.cos(phi),5)*( 5.0 -  18.0*t2 + t2*t2 + 14.0*eta2 - 58.0*t2*eta2)*Math.pow(dLambda,5) +
             1.0/5040.0*N*Math.pow(Math.cos(phi),7)*(61.0 - 479.0*t2 + 179.0*t2*t2 - Math.pow(t,6))*Math.pow(dLambda,7);
 
  return new Point2D("BL2XY", x, y);
 }

Kennziffer und die 500.000 sind noch zu berücksichtigen. Die Formeln stammen aus Bernhard Hofmann-Wellenhof, Gerhard Kienast, Herbert Lichtenegger: GPS in der Praxis. Ein Buch, das ich hiermit gleich mal mit empfehlen möchte.

Schöne Grüße
Micha

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

Tags:
Gauß-Krüger, UTM, Java, GeoTra, Koordinaten, Bessel