Orthodrome

Die Orthodrome (von Vorlage:GrcS) ist die kürzeste Verbindung zweier Punkte auf einer Kugeloberfläche (d. h. nicht die gerade Strecke durch die Kugel hindurch).

Die Orthodrome ist eine Geodäte für den speziellen Fall einer Kugeloberfläche. Die Orthodrome ist immer ein Teilstück eines Großkreises. In der Luftfahrt fliegt man meist entlang dieser Orthodrome, um die geringste Flugstrecke zurücklegen zu können. Die umgangssprachlich häufiger gebrauchte synonyme Bezeichnung ist daher Luftlinie.

Orthodrome auf der Erdkugel zwischen Los Angeles und London

Der kürzeste Weg auf der Kugeloberfläche zwischen Punkt A und B ist eine Orthodrome.

Definitionsbeschränkung

Die Orthodrome ist nicht definiert, wenn die gegebenen Punkte übereinstimmen oder auf der Kugeloberfläche gegenüberliegen (auf der Erdoberfläche also für antipode Punkte).

Berechnung

Grundlage für die folgenden Berechnungen sind die Formeln aus der sphärischen Trigonometrie.

Verwendete Variablen	Bedeutung
$\varphi$	Geographische Breite
$\lambda$	Geographische Länge
$A({\varphi }_A,{\lambda }_A)$	Anfangspunkt
$B({\varphi }_B,{\lambda }_B)$	Endpunkt
$N$	Nordpol
$\zeta$	Zentriwinkel (Strecke AB, ausgedrückt als Winkel)
$\alpha$	Innenwinkel von Dreieck ABN bei A
$\beta$	Innenwinkel von Dreieck ABN bei B
${\omega }_A$	Kurswinkel bei A
${\omega }_B$	Kurswinkel bei B
$P_N({\varphi }_N,{\lambda }_N)$ $P_S({\varphi }_S,{\lambda }_S)$	Nördlichster oder südlichster Punkt der Orthodrome

Dabei ist $\lambda$ in Richtung Westen negativ, Richtung Osten positiv; $\varphi$ ist positiv für Breiten der Nordhemisphäre und negativ auf der Südhalbkugel.

Strecke

Als Maß für die Entfernung zweier Punkte dient in der sphärischen Geometrie der zugehörige Mittelpunktswinkel (mit dem Kugelmittelpunkt als Scheitel). Dieser Winkel, hier mit $\zeta$ bezeichnet, lässt sich berechnen durch

$\cos \zeta =\sin ({\varphi }_A)\cdot \sin ({\varphi }_B)+\cos ({\varphi }_A)\cdot \cos ({\varphi }_B)\cdot \cos ({\lambda }_B-{\lambda }_A).$

Sind die gegebenen Punkte gleich, so erhält man durch diese Formel $\zeta =0^{\circ }$. Liegen die beiden Punkte auf der Kugeloberfläche gegenüber, so ergibt sich $\zeta =180^{\circ }$. In diesen Fällen sind Kurswinkel und die Lage der Orthodrome nicht definiert.

Um die Distanz zwischen den zwei Punkten zu berechnen, muss $\zeta$ noch mit dem mittleren Erdradius (rund 6.371 km) multipliziert werden (für $\zeta$ im Bogenmaß; falls $\zeta$ in Grad angegeben ist, muss noch zusätzlich mit $\pi /180$° multipliziert werden).

Der Winkel $\zeta$ kann über das Skalarprodukt der Ortsvektoren von $A$ und $B$ berechnet werden. Die obige Formel ergibt sich dann durch Umformungen mit Hilfe geometrischer Additionstheoreme für Sinus und Kosinus. Alternativ kann die Formel hergeleitet werden, indem der Seiten-Kosinussatz der sphärischen Trigonometrie auf das aus den Punkten $A$ und $B$ und dem Nordpol gebildete Dreieck angewendet wird.

Kurswinkel und rechtweisende Kurse

Der Kurswinkel (das Azimut) gibt die Bewegungsrichtung an und wird ab Nordrichtung im Uhrzeigersinn gezählt. Zur Berechnung kann man für die nicht-trivialen Fälle das oben erwähnte Dreieck verwenden. $\alpha$ und $\beta$ seien die Innenwinkel dieses Dreiecks bei A bzw. B (mit $0^{\circ }<\alpha ,\beta <180^{\circ }$). Sie können z. B. mithilfe des Seiten-Kosinussatzes berechnet werden.

$\cos \alpha =\frac{\sin ({\varphi }_B)-\sin ({\varphi }_A)\cdot \cos \zeta }{\cos ({\varphi }_A)\cdot \sin \zeta }$

$\cos \beta =\frac{\sin ({\varphi }_A)-\sin ({\varphi }_B)\cdot \cos \zeta }{\cos ({\varphi }_B)\cdot \sin \zeta }$

Fällt einer der gegebenen Punkte mit dem Nord- oder Südpol zusammen, so ist die entsprechende Gleichung wegen Division durch 0 unbrauchbar.

Die beiden Parameter $\alpha$ und $\beta$ lassen sich auch direkt aus den Breiten- und Längengraden ${\varphi }_A$ bzw. ${\varphi }_B$ und ${\lambda }_A$ bzw. ${\lambda }_B$ bestimmen:

$\cos \alpha =\frac{\cos ({\varphi }_A)\cdot \sin ({\varphi }_B)-\cos ({\lambda }_A-{\lambda }_B)\cdot \cos ({\varphi }_B)\cdot \sin ({\varphi }_A)}{\sqrt{1-(\cos ({\lambda }_A-{\lambda }_B)\cdot \cos ({\varphi }_A)\cdot \cos ({\varphi }_B)+\sin ({\varphi }_A)\cdot \sin ({\varphi }_B){)}^2}}$

$\cos \beta =\frac{\cos ({\varphi }_B)\cdot \sin ({\varphi }_A)-\cos ({\lambda }_A-{\lambda }_B)\cdot \cos ({\varphi }_A)\cdot \sin ({\varphi }_B)}{\sqrt{1-(\cos ({\lambda }_A-{\lambda }_B)\cdot \cos ({\varphi }_A)\cdot \cos ({\varphi }_B)+\sin ({\varphi }_A)\cdot \sin ({\varphi }_B){)}^2}}$

Kurswinkel

Für die Kurswinkel ${\omega }_A$ (am Anfangspunkt) und ${\omega }_B$ (am Endpunkt) gilt $0^{\circ }\le {\omega }_A,{\omega }_B<360^{\circ }$. Es ergibt sich:

Voraussetzung	Kurswinkel bei A	Kurswinkel bei B
B östlich von A	${\omega }_A=\alpha$	${\omega }_B=180^{\circ }-\beta$
B westlich von A	${\omega }_A=360^{\circ }-\alpha$	${\omega }_B=180^{\circ }+\beta$
Gleicher Meridian, B nördlich von A	${\omega }_A=0^{\circ }$	${\omega }_B=0^{\circ }$
Gleicher Meridian, B südlich von A	${\omega }_A=180^{\circ }$	${\omega }_B=180^{\circ }$

Rechtweisende Kurse A → B (Hinweg)

$rw{K}_A={\omega }_A$

$rw{K}_B={\omega }_B$

Rechtweisende Kurse B → A (Rückweg)

$rw{K_B}^{\prime }={\omega }_B\pm 180^{\circ }$

$rw{K_A}^{\prime }={\omega }_A\pm 180^{\circ }$

Das Rechenzeichen ($+$ oder $-$) ist hier so zu wählen, dass $0^{\circ }\le rw{K_B}^{\prime },rw{K_A}^{\prime }<360^{\circ }$ erfüllt ist.

Nördlichster oder südlichster Punkt

In einer gnomonischen Projektion werden Orthodromen stets als gerade Strecke abgebildet

Ob eine Orthodrome zwischen den Grenzpunkten einen nördlichsten oder südlichsten Punkt (Scheitelpunkt) hat, hängt vom Kurswinkel am Anfangspunkt (${\omega }_A$) und am Endpunkt (${\omega }_B$) ab.

Ein nördlicher Scheitelpunkt existiert, wenn entweder zugleich $0^{\circ }<{\omega }_A<90^{\circ }$ und $90^{\circ }<{\omega }_B<180^{\circ }$ oder zugleich $270^{\circ }<{\omega }_A<360^{\circ }$ und $180^{\circ }<{\omega }_B<270^{\circ }$ gilt. Man erhält unter dieser Voraussetzung:

${\varphi }_N=\arccos (|\sin ({\omega }_A)|\cdot \cos ({\varphi }_A))$

${\lambda }_N=\{\begin{array}{cc}{\lambda }_A+\arccos \left(\frac{\tan ({\varphi }_A)}{\tan ({\varphi }_N)}\right),& \text{falls}{\omega }_A<180^{\circ }\\ {\lambda }_A-\arccos \left(\frac{\tan ({\varphi }_A)}{\tan ({\varphi }_N)}\right),& \text{falls}{\omega }_A>180^{\circ }\end{array}$

Ein südlicher Scheitelpunkt existiert, wenn entweder zugleich $90^{\circ }<{\omega }_A<180^{\circ }$ und $0^{\circ }<{\omega }_B<90^{\circ }$ oder zugleich $180^{\circ }<{\omega }_A<270^{\circ }$ und $270^{\circ }<{\omega }_B<360^{\circ }$ gilt. In diesem Fall ergibt sich:

${\varphi }_S=-\arccos (|\sin ({\omega }_A)|\cdot \cos ({\varphi }_A))$

${\lambda }_S=\{\begin{array}{cc}{\lambda }_A+\arccos \left(\frac{\tan ({\varphi }_A)}{\tan ({\varphi }_S)}\right),& \text{falls}{\omega }_A<180^{\circ }\\ {\lambda }_A-\arccos \left(\frac{\tan ({\varphi }_A)}{\tan ({\varphi }_S)}\right),& \text{falls}{\omega }_A>180^{\circ }\end{array}$

Wegen der Festlegung $-180^{\circ }\le \lambda \le +180^{\circ }$ müssen die Ergebnisse für ${\lambda }_N$ und ${\lambda }_S$ unter Umständen korrigiert werden: Bei Werten über $180^{\circ }$ muss $360^{\circ }$ subtrahiert werden, bei Werten unter $-180^{\circ }$ muss $360^{\circ }$ addiert werden.

Die Formeln lassen sich begründen durch Anwendung der neperschen Regel auf das rechtwinklige Kugeldreieck, das durch Punkt A, den Scheitelpunkt und den Nordpol bestimmt ist.

Beispiel Berechnung der Entfernung Berlin–Tokio

Geographische Koordinaten der Anfangs- und Endpunkte:

• Berlin
  • 52° 31′ 0″ N = 52,517°
  • 13° 24′ 0″ E = 13,40°
• Tokio
  • 35° 42′ 0″ N = 35,70°
  • 139° 46′ 0″ E = 139,767°

Winkelberechnung

${\varphi }_A=52,517^{\circ }$

${\lambda }_A=13,40^{\circ }$

${\varphi }_B=35,70^{\circ }$

${\lambda }_B=139,767^{\circ }$

$\begin{array}{cc}\zeta & =\arccos (\sin ({\varphi }_A)\sin ({\varphi }_B)+\cos ({\varphi }_A)\cos ({\varphi }_B)\cos ({\lambda }_B-{\lambda }_A))\\ & =\arccos (\sin (52,517^{\circ })\sin (35,70^{\circ })+\cos (52,517^{\circ })\cos (35,70^{\circ })\cos (139,767^{\circ }-13,40^{\circ }))\\ & =\arccos (0,79353\cdot 0,58354+0,60853\cdot 0,81208\cdot (-0,59296))\\ & =\arccos (0,1700)\\ & =80,212^{\circ }\end{array}$

bzw. $\zeta =1,400$ im Bogenmaß

Streckenberechnung

Zur Vereinfachung wird von einer Erdkugel mit dem Umfang 40.000 km bzw. dem Radius 6366 km ausgegangen.

$\begin{array}{cc}L& =\frac{\zeta }{360^{\circ }}\cdot 40000\mathrm{k}\mathrm{m}\\ & =\frac{80,212^{\circ }}{360^{\circ }}\cdot 40000\mathrm{k}\mathrm{m}\\ & =8912\mathrm{k}\mathrm{m}\end{array}$

Oder für $\zeta$ im Bogenmaß:

$\begin{array}{cc}L& =\zeta \cdot 6366\mathrm{k}\mathrm{m}\\ & =8912\mathrm{k}\mathrm{m}\end{array}$

Das sind aufgrund der idealisierten Geodaten selbstverständlich nur zwei Näherungen. Die tatsächliche Entfernung zwischen den beiden angegebenen Punkten in Berlin und Tokio kann bei Verwendung des WGS84-Referenzellipsoids zu 8941,2 km genauer berechnet werden, also mit einer Abweichung von etwa 23 km oder 0,26 % im Vergleich zur zweiten Näherung.

Genauere Formel zur Abstandsberechnung auf der Erde

Mit folgenden Formeln kann der Abstand zwischen zwei Standorten auf der Erde auf 50 Meter genau berechnet werden, siehe dazu auch Thaddeus Vincenty. Dabei wird keine Kugel, sondern das WGS84-Ellipsoid zugrunde gelegt. Sollten Koordinaten eines anderen Referenzellipsoids verwendet werden, müssen die Parameter $a$ (Radius) und $f$ (Abplattung) angepasst werden.

Seien ${\varphi }_A$ und ${\lambda }_A$ die geografische Breite und Länge von Standort A, ${\varphi }_B$ und ${\lambda }_B$ die geografische Breite und Länge von Standort B im Gradmaß. Der Abstand zwischen beiden Standorten berechnet sich wie folgt:

Abplattung der Erde: $f=\frac{1}{298,257223563}$

Äquatorradius der Erde: $a=6378,137\mathrm{k}\mathrm{m}$

$F=\frac{{\varphi }_A+{\varphi }_B}{2}$, $G=\frac{{\varphi }_A-{\varphi }_B}{2}$, $l=\frac{{\lambda }_A-{\lambda }_B}{2}$

Zunächst wird der grobe Abstand D ermittelt:

$S=(\sin G{)}^2\cdot (\cos l{)}^2+(\cos F{)}^2\cdot (\sin l{)}^2$

$C=(\cos G{)}^2\cdot (\cos l{)}^2+(\sin F{)}^2\cdot (\sin l{)}^2$

$w=\arctan \sqrt{\frac{S}{C}}$

$D=2\cdot w\cdot a$

Dabei ist $w$ im Bogenmaß einzusetzen.

Der Abstand $D$ wird durch die Faktoren $H_1$ und $H_2$ korrigiert:

$T=\frac{\sqrt{S\cdot C}}{w}$

$H_1=\frac{3\cdot T-1}{2\cdot C}$

$H_2=\frac{3\cdot T+1}{2\cdot S}$

Der Abstand $s$ in Kilometern berechnet sich abschließend wie folgt:

$s=D\cdot (1+f\cdot H_1\cdot (\sin F{)}^2\cdot (\cos G{)}^2-f\cdot H_2\cdot (\cos F{)}^2\cdot (\sin G{)}^2)$

Berechnungsbeispiel Berlin – Tokio

$\begin{array}{ccc}{\varphi }_A& =& 52,516666666666667^{\circ }\\ {\lambda }_A& =& 13,400^{\circ }\\ {\varphi }_B& =& 35,700^{\circ }\\ {\lambda }_B& =& 139,766666666666667^{\circ }\\ \\ f& =& 0,00335281066474748\\ a& =& 6378,137\mathrm{k}\mathrm{m}\\ F& =& 44,108333333333333^{\circ }\\ G& =& 8,408333333333333^{\circ }\\ l& =& -63,183333333333333^{\circ }\\ S& =& 0,41498261872684\\ C& =& 0,58501738127316\\ w& =& 0,699965690768276\\ D& =& 8928,9541420394\mathrm{k}\mathrm{m}\\ T& =& 0,70391883329502\\ H_1& =& 0,95019099899696\\ H_2& =& 3,74926124548527\\ \\ s& =& 8941,20250458698\mathrm{k}\mathrm{m}\end{array}$

Der Abstand $s$ ist also auf etwa 50 m genau zu 8.941,2 km bestimmt worden.

Loxodrome

Gegenüberstellung von Loxodromen (rot) und Orthodromen (blau) auf einer Mercatorkarte, mit Weglängen in Kilometern.

Weg	Lox.	Orth.	Diff.
NY-MO	8359 km	7511 km	10,1 %
NY-DA	6207 km	6150 km	Vorlage:00,9 %
DA-MO	6596 km	6509 km	Vorlage:01,3 %

Loxodromeverlängerung relativ zur Orthodrome entlang des 50. Breitengrades in Prozent.

Bei der Navigation von Punkt A nach B mit einem Kompass eignet sich die Loxodrome besser, da sie die Meridiane immer im gleichen Winkel kreuzt, man also den einmal eingestellten (Kompass-)Kurs einfach beibehalten kann.

Bei kurzen Strecken ist eine Loxodrome nur unwesentlich länger als eine Orthodrome. Bei hoher Breite und bei Entfernungen unterhalb von 30 Längengraden liegt der relative Längenunterschied bei weniger als 1 %. Danach steigt er deutlich an. Eine Reise entlang des 50. Breitengrades über 180 Längengrade ist 45 % länger als der Weg über einen Großkreis, der dann über den Pol verläuft.

Siehe auch

• Geodätische Linie (Geodäte)
• Abweitung

Weblinks

• Berechnung von Entfernung und Anfangskurs für beliebige Rotationsellipsoide (englisch)
• Great Circle Mapper – Great Circle mapper including ETOPS ranges (englisch)
• Erweitertes Messwerkzeug – mit Vergleich von Orthodrome und Loxodrome zwischen zwei Punkten auf OSM-Weltkarte
• Entfernungsberechnung zwischen zwei Orten der Erde – auf Basis einer idealen Kugel mit einem Radius von 6378,388 km.
• Ausführliche Großkreisberechnung für die Praxis (englisch) – Orthodrome und Loxodrome Distanz, Anfangs und Endkurs, Mischsegeln, Etmal, Meridianschnittpunkte

Literatur

• J. Meeus: Astronomical Algorithms. , Astronomische Algorithmen (Deutsch) 2. Auflage. Willmann-Bell, Richmond 2000, ISBN 0-943396-61-1, S. 85.