Hallo,

ich damit schon zurechtkommen. Es wurden ja auch schon alle Variablen
erklärt, wobei mir nicht klar ist, was "E,F = 1./2. Exzentrizität" ist.

            ---------
           /  2     2 |
          /  a  -  b
e = --   /  ---------
       /        2
               a 

                        ---------
                      /  2     2 |
                     /  a  -  b
e' = f(Flip) = --   /  ---------
                  /        2
                          b 

Wo wir schon bei Excel sind auch gleich etwas zu MichaeL: Schade, das du
kein Excel hast. Meiner Meinung nach das mächtigste Office Programm von
Microsoft.

Stimmt, das beste Office von _Microsoft_. Aber ich setze da auf eine Lösung von der Softwareschmiede SUN.

So langsam glaube ich, das ich mir das Buch von dem Herrn Fröhlich doch
kaufen sollte. Mal sehen ob mein Chef mich sponsoren will. Ich werd mal
nett lächeln ...

Probier es mit Alkohol

Könnt Ihr mir vielleicht noch was zum Ellipsoid sagen?

Sie sind durch ihre Halbachse a und b definiert - siehe Häufig verwendete Referenzellipsoide.

Gauß-Krüger/Krassowski-Ellipsoid umschalte stimmen die angegebenen
Koordinaten. Beim Umschalten auf Gauß-Krüger/Bessel-Ellipsoid jedoch kommen
ganz andere Werte heraus. Komisch.

Gut, das ist logisch, also das verschiedene Werte raus kommen; Wann welche stimmen, da warte ich Flips nachträgliche Antwort mal ab, er ist ja nun wieder Motiviert
Was bekommst Du denn raus?

Bei mir sieht es nun so aus (hin und zurück):

    PointRHBL P = null;
    CoordsConverter CC = null;
    
    Ellipsoid Bessel = new Ellipsoid("Bessel", 3, false);
    P = new CoordsConverter(Bessel).RH2BL(4581441.0, 5936991.0);
    System.out.println(P.toString());
    P = new CoordsConverter(Bessel).BL2RH(P.getLatitude(), P.getLongitude());
    System.out.println(P.toString());
    System.out.println("----------n" );
    
    Ellipsoid Krassowski = new Ellipsoid("Krassowski", 6, false);
    P = new CoordsConverter(Krassowski).RH2BL(3382537.0, 5938294.0);
    System.out.println(P.toString());
    P = new CoordsConverter(Krassowski).BL2RH(P.getLatitude(), P.getLongitude());
    System.out.println(P.toString());
    System.out.println("----------n" );
    
    Ellipsoid GRS80 = new Ellipsoid("GRS80", 6, true);
    P = new CoordsConverter(GRS80).RH2BL(33382586.0, 5935821.0);
    System.out.println(P.toString());
    P = new CoordsConverter(GRS80).BL2RH(P.getLatitude(), P.getLongitude());
    System.out.println(P.toString());
    System.out.println("----------n" );

liefert:

East:   4581441.0
North:  5936991.0
Lat:    53.55961697043168
Lon:    13.229180899953215

East:   4581441.000003739
North:  5936990.999897883
Lat:    53.55961697043168
Lon:    13.229180899953215
----------

East:   3382537.0
North:  5938294.0
Lat:    53.55806194187098
Lon:    13.227411172361553

East:   3382536.9999718675
North:  5938293.999899098
Lat:    53.55806194187098
Lon:    13.227411172361553
----------

East:   3.3382586E7
North:  5935821.0
Lat:    53.558119984382245
Lon:    13.227409699960269

East:   3382585.9999718764
North:  5935820.999899081
Lat:    53.558119984382245
Lon:    13.227409699960269
----------

By The Way: Wie macht Ihr das eigentlich immer mit der grauen,
eingerückten Schrift, wenn Ihr irgendworauf bezug nehmt?

Wir lassen beim Antworten dieses Zeichen > am Anfang der Zeile stehen, den Rest macht die tolle Forensoftware

Aber was heißt bei Deinem TR-Programm bitte "Abs"?

Absolut/absoluter Betrag

a = |-3| = |3| = 3;

Micha

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

Hallo,

So, ich bin einen riesigen Schritt weiter.

Das freut mich!

war der entscheidende Tipp, da hier schon Umrechnungen in einer
Basic-Sprache vorliegen.

Ich wusst, Du ergötzt Dich dran. Bin durch Zufall drauf gestoßen als ich auch was suchte.

Y = Rechtswert

Ok

X = Hochwert

Ok

D = Y (wofür auch immer das gut sein mag)

Hat keinen Sinn, wie ich bemerkt habe, aber ok

A = 6377397.155 (große Halbachse des Bessel-Ellipsoid)

Sieht gut aus

B = 0.003342773154 (kleine Halbachse oder Abplattung des
Bessel-Ellipsoid)

Das ist schon ein Unterschied! Schau noch einmal die kleine Halbachse ist b = 6.356.078,965.

Stell Dir nur die Ellipse (Ellpsoid) vor mit einer Achse unter einem Meter - verdammt kleine Welt, in der wir leben
Die Abplattung kommt hin 1/299,1528128; Flip meint aber b in seinem Tool.

Pi = 3.14159265358979 (das ist wohl klar)

Ich denke, geprüft habe ichs nicht

1. wofür stehen die Fälle W=1 und W=2

W1 = 3° Streifen; W2 = 6° Streifen; U = 2 ist wohl UTM - das kann ich zumindets aus dem Code so entziffern.

2. was ist die Variable F

Siehe mein letztes Posting, 2. num. Exzentrizität e'. Das f(Flip) sollte F bei Philipp heißen - kam wohl nicht so gut an...

3. In der Zeile "A/(1-E²(sin Q)²)^(1/2)->N" wird da nach dem "1-E²"
multipliziert?

Ich zeichen mal wieder mit dem analogen Formeleditor:

             A
N = --------------------
       ---------------
   -  /     2     2
    / 1 - E * sin Q

Es gilt, Punktrechnung vor Strichrechnung - ganz klar.

4. Mit den letzten drei Formeln ab "lambda=":L kann ich in Visual Basic
nichts anfangen. Was wird da gemacht?

Dort wird zum einen der Wert auf dem Display ausgegeben und gleichzeitig in der Variabel L gespeichert, um, so ist zu vermuten, nachfolgenden Berechnungen zur Verfügung zu stehen - das hätte ich zumindest so gemacht.

5. Wie erfolgt die Ausgabe bzw. in welchen Variablen stehen die
geographischen Koordinaten?

siehe 4. in L, wie Lamda und P wie phi vgl. Griechisches Alphabet

Erstmal war ich einfach nur begeistert.

Freude und....

Aber dann kam die Ernüchterung.

Leid liegen Dich beisammen.

Overlay-Format nur geographische Koordinaten akzeptiert,
die sich auf das WGS84-Ellipsoid beziehen

Von uns (beiden) kannte keiner das Format, sodass wir Dich hier nicht wissentlich ins Leere liefen lassen.

Folgerichtig funktioniert die Umrechnung meiner UTM-Koordinaten exakt
richtig.

Die beziehen sich auf GRS80 und das entspricht WGS84

Die Fett markierte Zahl ist wohl die große Halbachse im Bessel-Ellipsoid,
allerdings kann man die nicht einfach durch die des WGS84-Ellipsoids
ersetzen.

Nein, Du erkennst langsam den Unterschied zwischen einer einfachen Umrechnung und einer Transformation. Für die Umrechnung musst Du auf einem Ellispoid bleiben. Deshalb gibt es hier auch die geschlossenen Formeln, die Flip oder das Basic-Tool verbaut haben. Beachte noch einmal den ersten Absatz aus seinem ersten Posting.

Irgendwas muß noch geändert werden,

Nein, Du bist hier Chancenlos, mit dem derzeitigen Ansatz wirst Du keinen Datumsübergang hinbekommen.

Schönes Wochenende

Dir auch!

Micha

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

Ich will ja hier nicht schon wieder den Lustlosen spielen - aber wenn die Umsetzung des TR-Codes wegen einer besseren Umsetzung anderen Codes nicht weiterverfolgt werden soll, würde ich mich mit einer kurzen Analyse begnügen wollen, da auch ich mich andernfalls erstmal ausführlichst in meinen eigenen TR-Code zurück-einarbeiten müsste (die Programmierung ist jetzt so ziemlich genau 3 Jahre und 2 Monate her!).

Hier also mal auf die Schnelle...

Zuallererst sei nochmals darauf hingewiesen, daß das beileibe nicht der ganze Code war. Es war ungefähr ein Viertel davon. Das Programm kann neben GK und UTM auch mit Soldnerkoordinaten umgehen. Während GK und UTM ähnliche Vorgehensweisen benötigen, ist Soldner schon wieder ein anders gelagertere Fall - somit sind einige Sachen speziell für diesen Fall programmiert worden und somit hier unwichtig (wie geklärt z.B. alles, was mit der Variablen C zusammen hängt)!
Zum vollständigen Verständnis wäre die Gesamt-Veröffentlichung sicherlich hilfreich (wenn gewünscht, mach ich das hier - s. dazu meine Ausführungen zu Beginn dieses Beitrags).

Y und X sind klar, D ist nochmals Y, weil das Y kurz danach "bearbeitet" wird, aber auch der Original-Y-Wert (nun als D) noch gebraucht wird.
A stimmt (genauer: 6377397.15507605), B nicht (B=6356078.96289778). W müsste auf den ersten schnellen Blick die Streifenfrage klären (W=1=3°, W=2=6°, W=3=Soldner) - bei UTM also immer W=2, bei GK im Normalfall W=1 (aber auch =2 möglich). F ist E' (=2. Exzentrizität) - das hatte "mein" Micha schon in seiner letzten Antwort vor deinem Beitrag (gleich zu Beginn) geklärt.

Wenn kein Zeichen steht, wird multipliziert, ja. Man beachte aber Punkt- vor Strichrechnung! Und deine beiden letzten Fragen erfordern eine gemeinsame Antwort - diese Formeln bedeuten die Ausgabe der Werte. Was in Anführungszeichen steht, wird im TR als Text ausgegeben, gespeichert wird der angezeigte Wert dann in der Variablen am Ende der Zeilen (also L (lambda) und P (phi)).

Und um noch kurz was zu deiner Basic-Geschichte zu sagen (den Hauptteil überlasse ich da großzügig deinem Namensvetter ):
Das UTM sich auf WGS84 bezieht und GK auf Bessel (oder eigenlich ja eher Krassowski in MV) hast du richtig erkannt - folgerichtig passen die UTM-Sachen gut und die GK-Sachen nicht (rein rechnerisch stimmt die GK-Sache zwar auch, aber der Übergang von lamda/phi in Bessel oder Krassowski zu lamda/phi in WGS84 fehlt dann halt noch (womit wir dann tatsächlich bei der Transformation angekommen wären)).
53,1234° in Bessel sind dann halt 53,1987° in WGS84 und somit folgerichtig leicht woanders in deiner TOP50-Karte (ist jetzt natürlich nur ein "Phantasie-Unterschied"!!!).

Wie dem auch sei: Diese Transformation ist nicht so einfach - finde ich jedenfalls. Micha sah das hier ja schon anders, womit er da vll. weiter ansetzen sollte. Die ganze 7-Parameter-Helmert-Transformation mit irgendwie herbeigezauberten 7 Parameterwerten war mir für eine Umsetzung immer zu hoch...