Referenzellipsoid
Bei gegebenem <math>\varphi</math> ergibt sich die Höhe h als
- <math>h=\frac{\sqrt{X^2+Y^2}}{\cos\varphi} - N_\varphi</math>
Obwohl diese Beziehung exakt ist, bietet sich die Formel
- <math>h\approx\sqrt{X^2+Y^2}\cdot\cos\varphi + Z \sin\varphi - a^2 / N_\varphi</math>
eher für praktische Berechnungen an, da der Fehler <math>\Delta h \approx \frac{1}{2} (h+\sqrt{X^2+Y^2})(\Delta\varphi)^2</math> nur quadratisch vom Fehler in <math>\varphi</math> abhängt.<ref>Bowring: The accuracy of geodetic latitude and height equations (Survey Review, Vol. 28)</ref> Das Ergebnis ist somit um einige Größenordnungen genauer.
Für die Berechnung von <math>\varphi</math> muss auf Näherungsverfahren zurückgegriffen werden. Aufgrund der Rotationssymmetrie wird das Problem in die X-Z-Ebene verlegt (<math>\lambda=0</math>). Für den allgemeinen Fall wird dann X durch <math>\sqrt{X^2+Y^2}</math> ersetzt.
Das Lot des gesuchten Punktes <math>(X,Z)</math> auf die Ellipse hat den Anstieg <math>\tan\varphi</math>. Das verlängerte Lot geht durch den Mittelpunkt M des Krümmungskreises, welcher die Ellipse im Lotfußpunkt berührt. Die Koordinaten des Mittelpunktes lauten
- <math>\begin{pmatrix}X_M\\ Z_M\end{pmatrix}=\begin{pmatrix}\varepsilon^2 a \cos^3 t\\ - {\tilde \varepsilon}^2 b \sin^3 t\end{pmatrix}</math>
mit
- <math>t</math> - parametrische Breite, d. h., Punkte auf der Ellipse sind durch <math>\left(\begin{smallmatrix}a\cos t\\b\sin t\end{smallmatrix}\right)</math> beschrieben
- <math>\tilde \varepsilon = \frac{\sqrt{a^2 - b^2}}{b}</math>
Damit gilt
- <math>\tan\varphi = \frac{Z+{\tilde\varepsilon}^2 \sin^3 t}{X - \varepsilon^2 \cos^3 t}</math>
Dies ist eine Iterationslösung, da <math>\varphi</math> und <math>t</math> über <math>\tan \varphi = \frac{a}{b} \tan t</math> in Beziehung stehen. Ein naheliegender Anfangswert wäre
- <math>\tan t_0=\frac{a}{b} \frac{Z}{X}</math>
Mit dieser Wahl erreicht man nach einem Iterationsschritt eine Genauigkeit von <math>\Delta\varphi \approx \frac{3}{2} \varepsilon^3 \frac{a h^2}{(a+h)^3} \sin^3\varphi \cos^3\varphi</math>.<ref>Bowring: Transformation from Spatial to Geographical coordinates (Survey Review, Vol. 23)</ref> Das heißt, auf der Erdoberfläche ergibt sich für <math>\varphi</math> ein maximaler Fehler von 0,00000003″ und das globale Maximum des Fehlers (bei <math>h=2a</math>) beträgt 0,0018″.
Bei günstiger Wahl von <math>t_0</math> kann auch der maximale Fehler für Punkte im Weltraum noch weiter reduziert werden. Mit
- <math>\tan t_0=\frac{b}{a}\frac{Z}{X}\left(1 + \tilde{\varepsilon} \frac{b}{\sqrt{X^2+Z^2}}\right)</math>
ist durch einmaliges Einsetzen in die Iterationsformel der Winkel <math>\varphi</math> (für die Parameter der Erde) auf 0,0000001″ genau bestimmt (unabhängig vom Wert von h).
Wichtige Referenzellipsoide
Die Form und Größe der in verschiedenen Regionen verwendeten Ellipsoide werden im Allgemeinen durch ihre große Halbachse <math>a</math> und die Abplattung <math>f</math> (engl. flattening) festgelegt. Ferner ist noch jener zentral gelegene „Fundamentalpunkt“ zu definieren, auf dem das Referenzellipsoid das Geoid berührt und ihm damit eine unzweideutige Höhenlage gibt. Beide Festlegungen zusammen werden „geodätisches Datum“ genannt.
Auch wenn zwei Länder dasselbe Ellipsoid verwenden (z. B. Deutschland und Österreich das Bessel-Ellipsoid), unterscheiden sie sich doch in diesem Zentralpunkt bzw. Fundamentalpunkt. Daher können sich die Koordinaten der gemeinsamen Grenzpunkte um bis zu einem Kilometer unterscheiden.
Die Achsen der Ellipsoide sind je nach der Region, aus deren Messungen sie bestimmt wurden, um bis zu 0,01 % verschieden. Die Genauigkeitssteigerung bei der Bestimmung der Abplattung <math> f=(a-b)/a </math> (Differenz der Ellipsoid-Achsen rund 21 km) hängt mit dem Start der ersten künstlichen Satelliten zusammen. Diese zeigten sehr deutliche Bahnstörungen bzgl. der Bahnen, die man vorausberechnet hatte. Anhand der Fehler konnte man zurückrechnen und die Abplattung genauer bestimmen.
Die Tabelle zeigt regionale Ellipsoide 1810 bis 1906 und global bestimmte Erdellipsoide von 1924 bis 1984 und die Entwicklung der Kenntnis vom mittleren Äquatorradius und der Erdabplattung.
Ellipsoid | Jahr | große Halbachse a in Meter | kleine Halbachse b in Meter | Numerus = 1/Abplattung (n = 1/f = a/(a−b)) | Anmerkungen | EPSG-Code |
---|---|---|---|---|---|---|
Delambre, Frankr. | 1810 | 6.376.985 | 308,6465 | Pionierarbeit | ||
Schmidt | 1828 | 6.376.804,37 | 302,02 | Pionierarbeit | ||
G.B. Airy | 1830 | 6.377.563,4 | 6.356.256,91 | 299,3249646 | ||
Airy 1830 modifiziert | 1830 | 6.377.340,189 | 6.356.034,447 | 299,3249514 | EPSG::7002 | |
Everest (Indien) | 1830 | 6.377.276,345 | 300,8017 | EPSG::7015 | ||
Bessel 1841<ref>[2], Constants for Reference Ellipsoids used for Datum Transformations From Bessel 1841 to WGS 1984 und Built-in Ellipsoids geben den auf 1 mm gerundeten Wert für b an, ausgehend von den in EPSG:7004 definierten Parametern a und f</ref> | 1841 | 6.377.397,155 | 6.356.078,963<ref><math>b=-(\tfrac{a}{n}-a)=-(\tfrac{6.377.397,155\,\text{m}}{{\color{Green}299,1528128}}-6.377.397,155\,\text{m})=6.356.078,96282m\simeq6.356.078,{\color{Red}963}\,\text{m}</math></ref> | 299,1528128<ref>Fälschlicherweise auch <math>n=\tfrac{1}{f}=\tfrac{a}{a-b}=\tfrac{6.377.397,155\,\text{m}}{6.377.397,155\,\text{m}-6.356.078,{\color{Red}963}\,\text{m}}={\color{Green}299,15281}{\color{Red}53513205998}</math></ref> | ideal angepasst in Eurasien oft benutzt in Mitteleuropa |
EPSG::7004<ref>[3] [4] EPSG:7004 nutzt <math>n=\tfrac{1}{f}={\color{Green}299,1528128}</math> Dieser Wert stammt aus dem US Army Map Service Technical Manual; 1943. “Remarks: Original Bessel definition is a=3272077.14 and b=3261139.33 toise. This used a weighted mean of values from several authors but did not account for differences in the length of the various toise: the "Bessel toise" is therefore of uncertain length.”</ref> |
Clarke | 1866 | 6.378.206,400 | 294,9786982 | ideal angepasst in Asien | EPSG::7008 | |
Clarke 1880 /IGN | 1880 | 6.378.249,17 | 6.356.514,99 | 293,4663 | EPSG::7011 | |
Friedrich Robert Helmert | 1906 | 6.378.200,000 | 298,3 | EPSG::7020 | ||
Australian Nat. | 6.378.160,000 | 298,25 | EPSG::7003 | |||
Modif. Fischer | 1960 | 6.378.155,000 | 298,3 | |||
Internat. 1924 Hayford | 1924 | 6.378.388,000 | 6.356.911,946 | 297,0 | ideal angepasst in Amerika bereits 1909 publiziert |
EPSG::7022 |
Krassowski | 1940 | 6.378.245,000 | 6.356.863,019 | 298,3 | EPSG::7024 | |
Internat. 1967 Luzern | 1967 | 6.378.165,000 | 298,25 | |||
SAD69 (South America) | 1969 | 6.378.160,000 | 298,25 | |||
WGS72 (World Geodetic System 1972) | 1972 | 6.378.135,000 | ≈ 6.356.750,52 | 298,26 | EPSG::7043 | |
GRS 80 (Geodätisches Referenzsystem 1980) | 1980 | 6.378.137,000 | ≈ 6.356.752,3141 | 298,257222101 | EPSG::7019 | |
WGS84 (World Geodetic System 1984) | 1984 | 6.378.137,000 | ≈ 6.356.752,3142 | 298,257223563 | für GPS-Vermessungen | EPSG::7030 |
Das Bessel-Ellipsoid ist für Eurasien ideal angepasst, sodass sein „800-m-Fehler“ für die Geodäsie Europas günstig ist – ähnlich wie die gegenteiligen 200 m des Hayford-Ellipsoids (nach John Fillmore Hayford) für Amerika.
Für viele Staaten Mitteleuropas ist das Bessel-Ellipsoid wichtig, ferner die Ellipsoide von Hayford und Krassowski (Schreibweise uneinheitlich), und für GPS-Vermessungen das WGS84.
Die Resultate von Delambre und von Schmidt sind Pionierarbeiten und beruhen auf nur begrenzten Messungen. Hingegen entsteht der große Unterschied zwischen Everest (Asien) und Hayford (Amerika) durch die geologisch bedingte Geoid-Krümmung verschiedener Kontinente. Einen Teil dieses Effekts konnte Hayford durch mathematische Reduktion der Isostasie eliminieren, sodass man dessen Werte damals für besser hielt als die europäischen Vergleichswerte.
Literatur
- Wolfgang Torge: Geodesy. 3. completely revised and extended edition. De Gruyter-Verlag, Berlin u. a. 2001, ISBN 3-11-017072-8.
- J. Ihde et al.: European Spatial Reference Systems – Frames for Geoinformation Systems (PDF).
Weblinks
- Kartenbezugssysteme, Ellipsoide, Geoide und topografische Oberflächen
- MapRef – Europäische Referenzsysteme und Kartenprojektionen
- CRS-EU - Information and Service System for European Coordinate Reference Systems.
- http://www.euref-iag.net/euref_links.html - EUREF-Links zu weiteren geodätischen Informationen
- http://www.epsg-registry.org/ - Datenbank zu Referenzsystemen und Koordinatentransformationsparametern
Einzelnachweise
<references />