Excel-Tabelle GPS-WGS84 nach Gauss-Krüger-DHDN.

[moenk]

Das hier musste ich heute mal nachrechnen:

Stimmt so etwa. Die Tabelle gibts hier, falls jemand auch mal sowas umrechnen will:
http://www.geoclub.de/files/GPS_nach_GK.xls
Die letzten Stellen stimmen aber nicht. Liegt das wohl an anderen Translationsparametern oder an der Ungenauigkeit von Minuten/Sekunden?

 

[Prof. Y. Lupardi]

Die Hollandisch Coordinate Converter 4.1
(download http://www.rdnap.nl/download/setup%20cc41.exe )
macht dieses:
ETRS89 LatN: 53d 8m 25,2s LonE: 8d 13m 4,2s
Gauss-Krüger Zone 3: Hoch: 5890082,350 Rechts: 3447729,712

 

[geometer42]

Wenn man im Excel-Sheet den Parametersatz 4 (für Norddeutschland) einstellt, sollten die Spannungen im alten DHDN berücksichtigt sein. Leider weiß ich nicht wie das angegebene System LS100 (Lagestatus 100) in Niedersachsen definiert ist. Möglicherweise rechnet man da mit noch einem anderen Satz von Transformationsparametern.

Mich würde mal interessieren, wo das Excel-File herkommt. Wer hat es gemacht und welche Formeln wurden benutzt ?

Viele Grüße
Chris

 

[moenk]

Gauss-Krüger Zone 3: Hoch: 5890082,350 Rechts: 3447729,712

Das meinte ich. Den Rechtswert kann ich auf der Tafel dann ja noch als geodätisch gerundet akzeptieren. Der Hochwert liegt aber fast drei Meter daneben. So gross ist die Abweichung bei meiner Excel-Tabelle aber nicht.
Hab die Antwort aber nun auch gefunden: 0,1 Sekunden machen inder Länge etwa zwei Meter aus, in der Breite sogar drei. Dann ist die Angabe der Minuten/Sekunden ungenau und wohl aus den GK-Koordinaten gerechnet.

 

[moenk]

Wenn man im Excel-Sheet den Parametersatz 4 (für Norddeutschland) einstellt, sollten die Spannungen im alten DHDN berücksichtigt sein.

Mit dem Satz komme ich nun wieder auf ganz andere Werte. Auf dem Schild steht ja auch 2001, macht also Sinn. Die Tabelle habe ich aus ein paar Beispielen von Prof. Reinking aus der Vorlesung an der FH zusammengebastelt, die Formeln kann man ja recht gut erkennen.

 

[geometer42]

Vielen Dank für den Hinweis.

 

[moenk]

Wir hatten ja grad geklärt wie das gemeint war. Hab mich denn nun grad noch mal dran gemacht und die Gegenrichtung in eine Excel-Mappe verpackt:
http://www.geoclub.de/files/GK_nach_GPS.xls
So werden die Angaben vermutlich zustande gekommen sein: Der Messpunkt war vermutlich schon vom Kataster bekannt und wurde für die GPS-Empänger mit dem 2001er-Translationsdaten die damals aktuell waren umgerechnet.

 

[pfeffer]

in Deiner Excel-Mappe sind ja mehrere Parametersätze.
1. was bedeuten die Parametersätze darin?
sind das die Gauß-Krüger-Streifen?
Oder unterschiedliche Sätze wegen der Netzspannungen?
2. EDIT: Und bevor ich wieder stundenlag teste: wie ist die Beziehung zu dem online-Umrechnungstool vom geodatenzetrum http://www.geodatenzentrum.de/geodaten/gdz_home1.gdz_home_start?gdz_home_spr=deu&gdz_home_para0=0&gdz_home_para1=Koordinatentransformation&gdz_home_menu_nr=15&gdz_home_menu_nr2=1&gdz_home_para3=https://upd.geodatenzentrum.de/auftrupd/ktrans?sprache=deu kommt da genau das gleiche raus, wie bei Deiner Excelberechnung? ist Deine Berechnung evtl. sogar genauer?
3. sehe ich das richtig, dass es eine sog. "7-Parameter-Transformation" ist? - und sehe ich das weiter richtig, dass man für andere Koordinatensysteme (solange es natürlich keine andere Projektion ist) "einfach" andere Parameter einsetzen kann?

Gruß,
Pfeffer.

 

[pfeffer]

also: klar ist mir inzwischen: Gauß-Krüger-Streifen sind es nicht.
Wahrscheinlich beziehen sie sich alle auf GK 2?
und die unterschiedlichen Parameter für Nord-, Mittel- und Süddeutschland ergeben sich wegen der Netzspannungen?
Wo sollten da die Grenzen gezogen werden?
Gehe ich richtig in der Annahme, dass die Nord-, Mittel- und Süddeutschlandparameter die genauesten sind?

Beste Grüße,
Pfeffer.

 

[moenk]

Die speziellen Parameter sind besser an die tatsächliche Form der Erde in diesem Bereich angepasst, allgemein wird wohl der für Deutschland universell vereinbarte genommen. Die Vermesser in Niedersachsen haben wieder einen anderen Satz der sich LS100 nennt und der da noch nicht mal dabei ist (und den kenn ich auch nicht).

 

[Wallraff]

Hallo moenk,

(und den kenn ich auch nicht)

Du kannst ihn kennenlernen
http://cdl.niedersachsen.de/blob/images/C12283445_L20.pdf
wenn der noch gültig ist.

Jetzt aber schnell ein paar Proberechnungen ...

Grüße Wallraff

 

[pfeffer]

programmierst Du die Umrechnung dann?
Ich bräuchte sie in Java...

Schöne Grüße,
Pfeffer.

EDIT: Wen es interessiert: Die Transformationsparameter für Deutschland und NRW und Teilgebiete von NRW sind zu finden unter http://www.lverma.nrw.de/produkte/raumbezug/koordinatentransformation/Koordinatentransformation.htm

Dort steht, dass die Transformation für ganz Deutschland nur auf +/-1m genau ist. Die für ganz NRW immerhon auf +/- 34cm. Allerdings ist da was seltsam: die für ganz NRW ist angeblich genauer als die für die Teile II und III aus NRW - da scheint mir etwas inkonsistent zu sein.

moenk, wo hast Du denn die verschiedenen Transformationsparameter her?

Schöne Grüße,
Pfeffer.

 

[moenk]

moenk, wo hast Du denn die verschiedenen Transformationsparameter her?

Die kommen wie die ganzen Excel-Formeln aus der Vorlesung in Vermessungskunde bei Prof. Reinking an der FH Oldenburg/Ostfriesland.

 

[pfeffer]

So, habe sie nun in java nachprogrammiert (zunächst nur WGS84 -> Gauß-Krüger).

Dabei ist mir grad aufgefallen:
a) in GK_nach_GPS ist der Maßstab falsch herum: er muss einmal 1/ gerechnet werden.
b) man kann es leicht vereinfachen, wenn man die Translation vor der Drehung anwendet, dann muss man die Translationsparameter nicht auch drehen.

Ich habe eine Excel-Tabelle gemacht, in der ich diese Änderungen eingetragen habe. Ich könnte sie Dir mailen, wenn Du mir Deine Mail-Adresse sagt.

Vielen Dank für die Excel-Tabelle! - Sie ist eine hervorragende Vorlage zum Nachprogrammieren. Man kann schön jeden einzelnen Zwischenschritt nachprüfen, wenn man sich beim abschreiben der Formeln irgendwo vertan hat.

Der Code kann demnächst im CacheWolf-SVN abgerufen werden: Es ist dei Klasse navi.TransformCoordinates.java und Matrix.java die die Arbeit leisten, link: http://svn.berlios.de/wsvn/cachewolf/trunk/?rev=0&sc=0

Damit ist CacheWolf der Möglichkeit, georeferenzierte Karten von deutschen Landesvermessungsämtern online abzurufen einen großen Schritt näher gekommen.


Beste Grüße,
Pfeffer.

 

[pfeffer]

so, jetzt ist es fertig, in beide Richtungen

Allerdings: Genau jetzt hat natürlich das Bundesamt für Kartographie und Geodäsie eine völlig andere Methode und entsprechenden Datensätze veröffentlicht, die eine wesentlich exaktere Umrechnung ermöglichen ( siehe http://crs.bkg.bund.de/crs-eu/ click on "national CRS" -> germany -> DE_DHDN / GK_3 -> DE_DHDN (BeTA, 2007) to ETRS89 oder direkt: http://crs.bkg.bund.de/crseu/crs/descrtrans/BeTA/BETA2007dokumentation.pdf). hmpf. Naja, erstmal kann man meine aktuelle Umsetzung verwenden. Vielleicht setze ich später noch die für ganz Deutschland genauere Methode um. Mal sehen. Für NRW sind bereits mit der vorhanden Implementierungen Umrechnungen mit einer maximalen Abweichung von 34cm möglich. Vermutlich wäre ähnliches auch in anderen Bundesländern möglich. Dazu müsste man nur die entsprechenden Transformationsparameter erfragen. Das Landesvermessungsamt hat seine im Internet veröffentlicht (siehe http://www.lverma.nrw.de/produkte/druckschriften/verwaltungsvorschriften/images/gps/TrafopsNRW.pdf ). Bei anderen Landesvermessungsämtern habe ich noch nicht geguckt.

Schöne Grüße,
Pfeffer.

 

[pfeffer]

Habe inzwischen heruasgefunden, wo die Transformationsparameter für Nord/Mittel und Süddeutschland herkommen: http://crs.bkg.bund.de/crs-eu/ click on "national CRS" -> germany -> DE_DHDN / GK_3 -> DE_DHDN (Middle) to ETRS89

Die dortigen Infos habe ich umgesetzt:

public static GkPoint wgs84ToGkGermany(CWPoint ll) {
		if (FORMER_GDR.isInBound(ll)) wgs84ToGk(ll, GK_GERMANY_2001); // exlcude former GDR from the splitting germany in north/middel/south
		if (ll.latDec <= 55 && ll.latDec >= 52.33333334 ) return  wgs84ToGk(ll, GK_NORD_GERMANY);
		if (ll.latDec <= 52.33333334  && ll.latDec >= 50.33333334 ) return  wgs84ToGk(ll, GK_MID_GERMANY);
		if (ll.latDec <= 50.33333334  && ll.latDec >= 47) return  wgs84ToGk(ll, GK_MID_GERMANY);
		return  wgs84ToGk(ll, GK_GERMANY_2001); 	
	}

Damit bekommen wir eine Genauigkeit von weniger als einem Meter (vermtulcih so um 30cm) für Westdeutschland und gut einen Meter Genauigkeit für Ostdeutschland.

Gruß,
Pfeffer.[/code]

 

[hustelinchen]

Ich möchte gerne das hier umrechnen...

Geographische Koordinaten DHDN UTM - Kordinaten ETRS89
Länge: 7 ° 33,3' Ost (E): 32399592
Breite: 51° 28,7' Nord (N): 5703865

dazu finde ich Null im Netz und hoffte nun auf die Excel Tabelle. Aber dafür passen die Angaben nicht. Weiß jemand wie ich das umrechnen kann?

Mich interessiert vor allem die 32399592, wieso sind das 8 Stellen? So kann ich das in mein GPS nicht eingeben.

 

[pfeffer]

die 32 im Ost-Wert ist vermutlich die UTM-Zone,

Gruß,
Pfeffer.

PS: aber ich habe nicht verstanden, was Du genau hast und in was Du es umrechnen willst. Mit CacheWolf kann man ansonsten Gauß-Krüger <-> WGS84 (= ETRS89 bis auf weniger als 5cm) umrechnen.

 

[t31]

N 51° 28.700 E 007° 33.300
32U E 399653 N 5704010 - laut CW und damit sehr nahe an obiges:
32U E 399592 N 5703865

mit Excel http://www.geoclub.de/viewtopic.php?p=361641#p361641

 

[pfeffer]

ich weiß immernoch nicht genau, was Du willst: Sind die Grad-Angaben in Bessel, d.h. DHDN? und Du willst sie in UTM, basieren auf ETRS89 umrechnen? - Dazu hätte Cachewolf keine Funktion. DHDN-Koordinaten können in Cachewolf nur als Gauß-Krüger eingegeben werden.

UTM bezieht sich international fast immer auf WGS84 bzw. ETRS89, aber das ist nicht zwingend, UTM könnte sich theoretisch auch auf DHDN beziehen.

Oder sind die Gradangaben in WGS84/ETRS89 und Du willst sie in Gauß-Krüger DHDN umrechnen? Das macht cachewolf.

Und: was willst Du überhaupt umrechnen, Du scheinst ja beide Angaben bereits zu haben?

oder wie oder was - Fragen über Fragen.

Gruß,
Pfeffer.