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

Wieso willst du überhaupt Transformieren? Nur wegen dem Maßstabsfaktor? Wenn sowohl Autocad davon ausgeht, dass der Maßstab 1 ist als auch das Tachymeter dann sind nur die Fehlervektoren vielleicht was größer.
Wie sieht es mit Open Surce Software aus? http://derletztekick.com/software/koordinatenumformung
Hab ich aber auch noch nicht getestet! Falls du positive Erfahrungen sammeln solltest kannst du ja nochmal bescheid sagen.

Gruß Markus

Hallo zusammen
ich bin neu und entschuld mich wenn ich noch Fragen stelle die für viele selbstverständlich sind.

Ich benötige ein kleines Programm das mir GK-Koordinaten in Länge und Breite umrechnet.
habe in meinen alten Unterlagen noch ein paar Formeln gefunden und diese auch umgesetzt jedoch komme ich nicht auf die richtige Ergebnisse. Vielleicht kann mir jemand da weiter helfen. Ich bin auch für andere Lösungen offen. es ist schon ein bisschen her als ich damit aktiv zu tun hatte.

Mein Code bis jetzt

private MyPunkt2D berechnung(double x, double y)
        {
            string breite= "";
            string länge= "";
            double roh = 180 / Math.PI;
        // Besselellipsoid
            double halbachse = 6377397.155;     // große Halbachse in Meter
            double f = 1/299.152813;            // Abplattung
            double lh = 9;                  // Hauptmeridian
            double mh = 1;
   // start Berechnungen
        // Vorberechnungen
            double g = x / mh;
            double b = halbachse - (f - halbachse);
            double e2 = (Math.Pow  (halbachse,2) - Math.Pow(b,2)) / (Math.Pow(halbachse,2));
            double e_strich = Math.Sqrt((Math.Pow(halbachse,2) - Math.Pow(b,2)) / (Math.Pow(b,2)));
            double c = Math.Pow(halbachse,2) / b;
 
            double e0 = c * (1 - 3 / 4 * Math.Pow(e_strich, 2) + 45 / 64 * Math.Pow(e_strich, 4) - 175 / 256 * Math.Pow(e_strich, 6) + 11025 / 16384 * Math.Pow(e_strich, 8));
            double f2 = 3 / 8 * Math.Pow(e_strich, 2) - 3 / 16 * Math.Pow(e_strich, 4) + 213 / 2048 * Math.Pow(e_strich, 6);
            double f4 = 21 / 256 * Math.Pow(e_strich, 4) - 21 / 256 * Math.Pow(e_strich, 6);
            double f6 = 151 / 6144 * Math.Pow(e_strich, 6);
 
            double a = g / e0;
            double bf = a + f2 * Math.Sin(2 * a) + f4 * Math.Sin(4 * a) + f6 * Math.Sin(6 * a);
            double tf = Math.Tan(bf);
            double nf2 = Math.Pow(e_strich, 2) * Math.Pow(Math.Cos(bf), 2);
            double rnf = a / (Math.Sqrt(1 - e2 * Math.Pow(Math.Sin(bf), 2)));
 
            double b2 = -(tf*(1+nf2))/(2*Math.Pow(mh,2)*Math.Pow(rnf,2));
            double b4 = (tf * (5 + 3 * Math.Pow(tf, 2) + 6 * nf2 * (1 - Math.Pow(tf, 2)))) / (24 * Math.Pow(mh, 4) * Math.Pow(rnf, 4));
            double b6 = -(tf * (61+ 90 * Math.Pow(tf, 2) + 45* Math.Pow(tf, 4))) / (720 * Math.Pow(mh, 6) * Math.Pow(rnf, 6));
 
            double l1 = 1 / (rnf * mh * Math.Cos(bf));
            double l3 = -(1 + 2 * Math.Pow(tf, 2) + nf2) / (6 * Math.Pow(mh, 3) * Math.Pow(rnf, 3) * Math.Cos(bf));
            double l5 = (5 + 28 * Math.Pow(tf, 2) + 24 * Math.Pow(tf, 4)) / (120 * Math.Pow(mh, 5) * Math.Pow(rnf, 5) * Math.Cos(bf));
 
     // Berechnung länge und Breite
            double breite_dezimal = bf + b2 * Math.Pow(y, 2) + b4 * Math.Pow(y, 4) + b6 * Math.Pow(y, 6);
            double länge_dezimal = lh + l1 * Math.Pow(y, 1) + l3 * Math.Pow(y, 3) + l5 * Math.Pow(y, 5);
 
            int breite_grad = Convert.ToInt16(breite_dezimal);
            int breite_minute = Convert.ToInt16(((breite_dezimal-Convert.ToDouble(breite_grad))*100*0.6));
            int breite_sekunde = Convert.ToInt16(((((breite_dezimal-Convert.ToDouble(breite_grad))*100*0.6)-Convert.ToDouble(breite_minute))*100*0.6));
            breite += Convert.ToString(breite_grad) + "°" + Convert.ToString(breite_minute) + "'" + Convert.ToString(breite_sekunde) + "''";
 
            int länge_grad = Convert.ToInt16(länge_dezimal);
            int länge_minute = Convert.ToInt16(((länge_dezimal - Convert.ToDouble(länge_grad)) * 100 * 0.6));
            int länge_sekunde = Convert.ToInt16(((((länge_dezimal - Convert.ToDouble(länge_grad)) * 100 * 0.6) - Convert.ToDouble(länge_minute)) * 100 * 0.6));
            länge += Convert.ToString(länge_grad) + "°" + Convert.ToString(länge_minute) + "'" + Convert.ToString(länge_sekunde) + "''";
 
            MyPunkt2D punkt_LB = new MyPunkt2D (länge,breite);
            return punkt_LB;

vielen Dank schon mal