Hallo Dieter,
nochmal Danke für die Anregung 
Der Teil Loxodrom hat gepasst (nachdem ich Wikipedia ignoriert habe - ist OMO mathematisch nicht ganz sauber. Mathematische Herleitung bei rainerstumpe s.u. gut erläutert, Aufbau jedoch für Schiffsnavigation mit Hilfsmittel Rechenschieber)
Bei der idealem Kugel (Erdabplattung noch ignoriert) lässt sich die Berechnung relativ einfach automatisieren. Daher habe ich ausschließlich das Bogenmaß genutzt; für Strecken km.
a) aus der Mercatorprojektion ergibt sich für die einzelnen Koordinaten (WP)
- WP (lat;lon) 'Wegpunkt in Bogenmaß angeben
- M_WPx_lat = LN(TAN(lat/2+PI/4) ' Projektion der Breite
- d_lat = M_WPziel_lat - M_WPstart_lat
- d_lon = lon_ziel - lon_start
b) Berechnung des Kurswinkels KW und der Strecke
- KW = PI/2 - ARCTAN2(d_lon;d_lat) ' Kurswinkel ARCTAN2 hat den Ergebnisbereich von -PI bis +PI, keine weitere Winkelbetrachtung notwendig
- KW = ARCTAN(d_lon/d_lat) ' Bewertung der Winkels durchführen, da der Ergebnisbereich zwischen -PI/2 bis +PI/2 liegt, ggf muss zum Winkel PI addiert werden
- für den Erdradius gibt es verschiedene Betrachtungsmöglichkeiten.
Nach Erd-Ellipsoid WGS 84
er_a = 6378,137 km ' Äquator
er_b = 6356,752315 km ' Pol
er_m = 6367,4446575 km ' gemittelter Wert aus er_a + er_b
- Strecke = (WPziel_lat - WPstart_lat) / COS(KW) * er_m
c) Berechnung einer Zielkoordinate von WP mit KW + Strecke
- Nach Formelumstellung ergibt sich (o.a. gilt noch
- WPziel_lat = Strecke/er_m * COS(KW) + WPstart_lat
- WPziel_lon = TAN(KW) * (M_WPziel_lat - M_WPstart_lat) + WPstart_lon
Überlegung: Noch genauere Strecke: mit er_a arbeiten und eine Funktion mitlaufen lassen, die die Länge des Erdradius in Abhängigkeit des Breitengrades berücksichtigt.
a) f(x) * cos(?)
b) würde ein Mittelwert aus den beiden Breitengraden ein hinreichend genauer Faktor sein
Quellen:
http://www.rainerstumpe.de/HTML/navigation.html
http://wikipedia.de/
ETH Zürich, Hans Walser, Mecatorkarte Lernumgebung (gefunden per Google Suche)
Hoffe, es sind keine Tippfehler drin 