Welchen EPSG-Code soll ich für das Koordinatensystem ETRS89 UTM Zone 33 nehmen?. Es gibt mehrere dafür. Es kursieren die Nummern 25833, 35833, 325833 und 5650 und es gibt unterschiedliche Angaben über die Verwendung von Zonenkennzeichen im Rechtswert der Koordinaten. Mal gar nicht, mal eine führende 3 und mal eine 33. In diesem Beitrag möchte ich ein wenig Klarheit in die Thematik bringen.
Ich will nicht die ganze Geschichte von Koordinatenreferenzsystemen hier beschreiben, das wäre vielleicht etwas für einen eigenen Blog-Beitrag, aber einige Begriffe, die ich hier verwenden möchte muss ich dennoch vorher erklären, denn: Begriffe ohne Erklärungen sind leer und Erklärungen ohne Begriffe sind schwer.
Mit EPSG-Codes meine ich die von der EPSG oder jetzt OGP zusammengestellten und bereitgestellten Nummern zur Identifizierung von Kartenprojektionen. Die Betonung liegt auf Projektion.
Die Projektion im Sinne vom EPSG definiert wie Punkte auf einem Ellipsoid in die Ebene abgebildet werden. Wer nicht genau weiß was ein Ellipsoid ist helfe sich vielleicht damit, dass er sich für die Erde eine Kugel vorstellt, die am Nord- und Südpol etwas zusammengedrückt ist. Von solchen Ellipsoiden gibt es allerdings einige, denn die sollen nicht nur am Nord- und Südpol möglichst gut passen, sondern auch jeweils da, wo man die Koordinaten verwendet, bzw. da wo man sie gemessen hat. Das ist historisch bedingt, weil man früher noch keine Satelliten hatte um ein globales Referenzsystem einfach verwenden zu können, wurden die Referenzellipsoide einfach dem jeweiligen Land bestmöglich angepasst. Die Festlegung eines solchen Referenzellipsoid wird geodätisches Datum genannt, und die Lagedifferenz zu dem mittleren Erdellipsoid, auf das sich mittlerweile alle geeinigt haben heißt Datumsshift.
So nun genug der Begriffe. Was genau wird denn nun in einem EPSG-Code beschrieben?
Den aktuellsten Stand was einem EPSG-Code zugeordnet wird findet man am besten in der EPSG Registry. Dort auf retrieve by code gehen und die Nummer eingeben und rechts auf view (Bild 1). Da steht nun eine ganze Menge und man kann noch zusätzliche Informationen einblenden und ausblenden. Was davon ist wichtig, bzw. was davon ist eigentlich das was zum EPSG-Code gehört und was sind Zusatzinformationen? Am besten kann man das sehen wenn man sich die Textrepräsentation der Definition ansieht. Dazu auf den Button WKT rechts klicken. Hier mal ein Beispiel für den EPSG-Code 25833:
PROJCRS["ETRS89 / UTM zone 33N",
BASEGEODCRS["ETRS89",
DATUM["European Terrestrial Reference System 1989",
ELLIPSOID["GRS 1980",6378137,298.257222101,LENGTHUNIT["metre",1.0]]]],
CONVERSION["UTM zone 33N",
METHOD["Transverse Mercator",ID["EPSG",9807]],
PARAMETER["Latitude of natural origin",0,ANGLEUNIT["degree",0.01745329252]],
PARAMETER["Longitude of natural origin",15,ANGLEUNIT["degree",0.01745329252]],
PARAMETER["Scale factor at natural origin",0.9996,SCALEUNIT["unity",1.0]],
PARAMETER["False easting",500000,LENGTHUNIT["metre",1.0]],
PARAMETER["False northing",0,LENGTHUNIT["metre",1.0]]],
CS[cartesian,2],
AXIS["easting (E)",east,ORDER[1]],
AXIS["northing (N)",north,ORDER[2]],
LENGTHUNIT["metre",1.0],
ID["EPSG",25833]]
Hier finden wir nun einiges von dem wieder was ich oben beschrieben habe:
- das Ellipsoid auf das sich die Koordinaten beziehen
- das Datum mit dem das Ellipsoid festgelegt ist
- und die Methode wie die Koordinaten von dem Ellipsoid in die Ebene abgebildet werden, also der Name und die Parameter der Projektion.
In diesem Beitrag will ich wie gesagt nicht genau darauf eingehen wie die Projektion funktioniert, nur soviel:
- der Datumsshift ist hier nicht angegeben
- die Koordinatenachsen sind beschrieben mit einer Verschiebung im Rechtswert (False easting)
Der Wert von False easting bedeutet in diesem Fall, dass die Rechtswerte im ebenen Koordinatensystem an dem genannten Meridian, hier Longitude of natural origin 15° östliche Länge, den Wert 500000 haben. Rechts davon steigen die Koordinaten an, links nehmen sie an. Das hat man gemacht, damit keine negativen Koordinaten vorkommen.
Also Punkt 1 hätten wir geklärt. Der EPSG-Code 25833 gilt nur für Koordinaten ohne Zonenkennzeichen im Rechtswert vor der 500000. Also 3500000 oder 33500000 wären falsch!!!
Schauen wir uns jetzt mal den EPSG-Code 5650 an, der von Landesamt für Innere Verwaltung Koordinierungsstelle für Geoinformationswesen im Dokument Empfohlene EPSG-Codes für Koordinatensysteme in Mecklenburg-Vorpommern für ETRS89 UTM Zone 33 empfohlen wird an.
PROJCRS["ETRS89 / UTM zone 33N (zE-N)",
BASEGEODCRS["ETRS89",
DATUM["European Terrestrial Reference System 1989",
ELLIPSOID["GRS 1980",6378137,298.257222101,LENGTHUNIT["metre",1.0]]]],
CONVERSION["UTM zone 33N with prefix",
METHOD["Transverse Mercator",ID["EPSG",9807]],
PARAMETER["Latitude of natural origin",0,ANGLEUNIT["degree",0.01745329252]],
PARAMETER["Longitude of natural origin",15,ANGLEUNIT["degree",0.01745329252]],
PARAMETER["Scale factor at natural origin",0.9996,SCALEUNIT["unity",1.0]],
PARAMETER["False easting",33500000,LENGTHUNIT["metre",1.0]],
PARAMETER["False northing",0,LENGTHUNIT["metre",1.0]]],
CS[cartesian,2],
AXIS["easting (E)",east,ORDER[1]],
AXIS["northing (N)",north,ORDER[2]],
LENGTHUNIT["metre",1.0],
ID["EPSG",5650]]
Die Definition ist fast identisch, bis auf den Namen, da ist der Zusatz (zE-N) drin und False easting ist 33500000. zE-N bedeutet mit Zonenkennziffer im Rechtswert und 33500000 dass im Rechtswert eine führende 33 zu stehen hat. Man braucht also um von 25833 zu 5650 zu gelangen nur bei allen Rechtswerten die Zahl 33000000 zu addieren. Erstaunlich ist, dass in dem Dokument vom 03.12.2012 nicht drin steht, dass das mit führenden 33 ist. Das müsste man im Erlass zu amtlichen geodätischen Bezugssystemen in Mecklenburg-Vorpommern (Landesbezugssystemerlass) vom 15.03.2005 nachlesen. In Absatz 3.2.2 steht, dass die auf den 9° Meridian bezogenen UTM Koordinaten die Kennziffer 32 haben und die auf 15° bezogenen die Kennziffer 33. Auch steht da, dass die Kennziffern dem Ostwert vorangestellt werden.
Und genau hier ist der Wiederspruch: Wenn man dem Ostwert die Kennziffer 33 voranstellt, sind es keine UTM Koordinaten mehr. Man muss also dringend berücksichtigen, dass wenn man von dem System ETRS89 UTM Zone 33 redet eigentlich ETRS89 UTM Zone 33 (zE-N) meint, wenn man über dass Bezugssystem vom M-V redet unbedingt den EPSG-Code 5650 verwenden, nicht 25833.
Hier noch mal dick hervorgehoben: EPSG:25833 = ETRS89 / UTM zone 33N
EPSG:5650 = ETRS89 / UTM zone 33N (zE-N)
Man hat hier also mal wieder eine Extrawurst gebraten.
--- Fortsetzung folgt --- |