Lösen von Problemen der Sphärischen Trigonometrie unter Zuhilfenahme von Computeralgebra

1. Lösen von Problemen der Sphärischen Trigonometrieunter Zuhilfenahme von Computeralgebra Parametrisierung und Visualisierung Walter Wegscheider, PI-Hollabrunn (Niederösterreich)

2. Suche nach einem Koordinatensystem • Verschiedene Weltbilder (geozentrisch, heliozentrisch) • Orientierung auf der Erdkugel – keine 1:1 Übertragung von einfachen rechtwinkeligen Systemen • schon frühe (Griechen, Chinesen u.a.) Bestrebungen für Koordinaten auf der Kugel • Schließlich durchgesetzt hat sich • Breitengrad (Latitude) und • Längengrad (Longitude)(und Höhe)

3. Längengrad - Breitengrad Bezugsgrößen: • Geogr. Breite • Pole, Äquator • Sonnenhöchststand • Datum / Schiefe der Ekliptik • Geogr. Länge • Nullmeridian (Greenwich – London), willkürlich • Sonnenhöchststand • Genaue Uhrzeit

4. Geogr. Breite - Breitenkreise • Seit dem späten 15. Jhdt. berechenbar • Astrolabium, Jakobsstab, Sextant • Auch Parallelkreise genannt • Parallel zum Äquator • Kleinkreise (außer Äquator) • Eckdaten • 90° S(üdpol) • 0° - Äquator • 90° N(ordpol) • Abstand eines Punktes vom Äquator Sextant

5. Geogr. Länge - Meridiane • Erst seit 1774 berechenbar (der Chronometer – eine genügend robuste und genau gehende Uhr war erfunden worden) • Meridiane sind Halbkreise (Großkreise), die die beiden Pole miteinander verbinden • Willkürliche Festlegung von Greenwich (London) als Nullmeridian. • 180° W (-180°) – 0° (Greenwich) – 180° O • Übergang bei Datumsgrenze

6. Probemessung • Am 21. Juni messen wir zu Mittag einen Sonnenhöchststand von 53.5° • Geogr. Breite 90° - 53.5° + 23.5° (Ekliptik) = 60° • Voraussetzung: Nordhalbkugel! • Der Chronometer zeigt an, dass dieser Sonnenhöchststand um 14.00 Uhr Nachmittag erreicht wurde. • Geogr. Länge: 30° Ost • 360° / 24 Stunden  15° pro Stunde • Wir sind in St. Petersburg!

8. Visualisierung – 3D Plots • Wenn wir versuchen, die Kugel (Modell der Erde) und Großkreise bzw. Meridiane zu visualisieren, bietet sich ein CAS wie DERIVE an. • Recht bald zeigen sich aber Schwierigkeiten bei gewohnten Schreibweisen • Beispiel 1: Kugel mit Radius 1 • Beispiel 2: Kreis als Schnitt einer Ebene mit einer Kugel! DERIVE v01

9. Koordinatentransformation • Statt mit kartesischen Koordinaten rechnen wir mit Kugelkoordinaten! • Ein beliebiger Punkt P wird im sphärischen Koordinatensystem durch den Radius r und zwei Winkel (r, a, d) festgelegt. Dabei steht d für den Winkel zur xy-Ebene (geogr. Breite), a für den Winkel innerhalb der xy-Ebene (geogr. Länge).

10. Umrechnungsformeln 1 Der Ortsvektor zu P hat die Länge r. Somit gilt: Der Winkel a (analog dazu bei d) kann in der xy-Ebene über trigonometrische Grundbeziehungen beschrieben werden:

11. Umrechnungsformeln 2 Die Projektion von P in die xy-Ebene erhält man über . - die weitere Projektion auf die x-Achse durch Analog dazu ergibt die Projektion auf die y-Achse:

12. Beliebiger Punkt auf der Kugel • Eine Möglichkeit, einen beliebigen Punkt auf der Kugel zu definieren – alternativer Zugang: • Wir betrachten den vertikalen Großkreis der Kugel mit Radius r, der beim Schnitt mit y=0 entsteht (entspricht Nullmeridian – wobei ein Meridian eigentlich nur ein Halbkreis ist!):

13. Nun lassen wir den entstanden Kreis um einen beliebigen Winkel s rotieren. Dazu multiplizieren wir q(t) mit der Rotationsmatrix Q. Ein Punkt auf der Kugel mit den Koordinaten r, t (geogr. Breite, 0 bis 2p bzw. –p bis p) und s (geogr. Länge, 0 bis 2p bzw. –p bis p) genügt daher folgender Parameterdarstellung. DERIVE v02

14. Visualisierung von Breiten- und Längenkreisen bzw. Raumpunkten • Wenn man in der Parameterdarstellung der Kugel einen der beiden Winkel durch einen fixen Wert ersetzt, erhält man Breiten- bzw. Längenkreise! • Wenn beide Parameter durch Werte ersetzt werden, erhält man die kartesischen Koordinaten eines Raumpunkts! • Beispiele: • t ersetzen durch d = 48° (N) – Plot des Breitenkreises • s ersetzen durch a = 16° (O) – Plot des Längenkreises • Plot des Raumpunkts Wien(6371,16°,48°)! DERIVE v03

15. Großkreise - Orthodrome • Ein Großkreis entsteht durch Schnitt der Kugel mit einer Ebene durch den Kugelmittelpunkt. Ein Großkreis wird daher über die beiden folgenden Gleichungen definiert: Eine Parameterdarstellung geht von der Darstellung eines Kreises in der Ebene aus: Im Raum in der xy-Ebene!

16. Für die Definition eines beliebigen Großkreises ist also eine entsprechende Orthonormalbasis erforderlich: In der sphärischen Geometrie interessieren wir uns in erster Linie für Großkreise, die durch zwei Punkte am Kreis definiert sind. Wir müssen also, um eine entsprechende Orthonormalbasis für eine Parameterdarstellung zu finden, die beiden Ortsvektoren zu den Punkten P und Q einbinden! … steht normal auf die Ortsvektoren der Punkte P und Q und damit normal auf die Ebene des Großkreises … steht normal auf die Normale des Großkreises, befindet sich also wieder in der Ebene des Großkreises und steht normal auf den Ortsvektor von P

17. Normiert auf Länge 1 erhalten wir Und damit den Großkreis mit: DERIVE v04

18. Länge der Orthodrome Die kürzeste Verbindung zwischen zwei Punkten, als Teil eines Großkreises, nennt man Orthodrome (Geradlaufende). Für die sphärische Entfernung zwischen zwei Punkten P und Q gilt: Wenn man nun mit Hilfe der Meridianbögen durch P und Q ein spärisches Dreieck erzeugt, kann man mit Hilfe der Sätze im Dreieck die Länge von c berechnen! Seiten-Cosinussatz:

19. Loxodrome - Kursgleiche • Sphärische Kurve, die entsteht, wenn auf einer Kugel mit Radius r ein fixer Winkele zu den Meridianen eingehalten wird! • Man spricht auch von einem konstanten Kurswinkel oder Azimut. Die Bedeutung der Loxodrome liegt in der praktikableren Navigation (mit Kompass) – bei speziellen Kartenentwürfen (Mercator) werden Loxodrome als Gerade abgebildet!

20. Wir betrachten zwei knapp nebeneinander liegende Punkte P und P1 auf einer Loxodrome und konstruieren mit Hilfe der Meridiane durch P und P1 ein spezielles Dreieck mit einem Hilfspunkt Q am Breitenkreis von P.

21. Wenn die beiden Punkte P und P1 unendlich knapp nebeneinander gewählt werden, entartet das sphärische Dreieck PQP1 zu einem rechtwinkeligen ebenen Dreieck mit der Hypothenus Ds! Die Beziehungen in diesem Dreieck lauten (a und d im Bogenmaß, r Radius der Kugel): Die Umformung dieser Beziehung liefert folgende Differentialgleichung: DERIVE v05

22. Wir erhalten damit die Gleichung einer Loxodrome(nschar) mit Kurswinkel e: Für die Berechnung einer bestimmten Loxodrome durch zwei Punkte P(aP, dP) und Q(aQ, dQ) kann man daraus ableiten – die Steigung tan e der Loxodrome beträgt:

23. Loxodrome - Parameterdarstellung • Die Parameterdarstellung der Loxodrome mit Steigung lautet: Die Loxodrome ist ein Spezialfall einer logarithmischen Spirale!

24. Loxodrome - Bogenlänge Für die Berechnung der Länge zwischen P und Q betrachten wir wieder das Dreieck PQP1 und berechnen die Länge von ds. Die Bogenlänge der Loxodrome zwischen den Punkten P und Q beträgt: DERIVE v06 (Vergleich der Längen: Orthodrome, Loxodrome)

25. Kartenentwürfe • Bildet man die Erdoberfläche oder einen Teil von ihr auf eine Ebene ab, so entsteht eine geographische Karte! • Problem: Längentreue in alle Richtungen ist unmöglich (Beweis durch Euler 1777) • Möglich sind Längentreue in bestimmten Richtungen (z.B. entlang der Meridiane, Flächentreue und Winkeltreue) • Verschiedenste Entwürfe (z.B. perspektivische, Zylinder- und Kegelentwürfe)

26. Mercator-Entwurf, Winkeltreue Der Mercator-Entwurf (Gerhard Mercator, 1512 – 1594) ist eine Zylinderprojektion. Die Erdoberfläche wird auf einen Zylindermantel projiziert, welcher die Erde am Äquator berührt. Der Zylinder wird nach der Projektion längs einer Mantellinie aufgeschnitten und abgerollt. Die Besonderheit gegenüber der klassischen Zylinderprojektion ist die rechnerische Korrektur, um Winkeltreue zu erzielen.

27. Verzerrung bei Zylinderprojektion abhängig von der geogr. Breite d, mit Hilfe der ersten Fundamentalform von Flächen erhält man die dazugehörigen Verzerrungsfaktoren: Verzerrung bei Mercatorprojektion Idee: Streckung in senkrechter Richtung – für winkeltreue Abbildung muss die Verzerrung gleich sein! Kompensation um:

28. Die neue senkrechte Komponente sei also v = f(d) mit: Über mehrere Umformungen kommt man zu folgender Differentialgleichung: Die Lösung dieser Gleichung führt zu einer schon bekannten Lösung (vgl. Loxodrome): Eine Loxodrome wird damit in der Mercator-Projektion als Gerade abgebildet! – Bedeutung für Navigation immens!