Jak dokładny jest przybliżenie Ziemi jako kuli?

63

Jaki poziom błędu napotkałem, przybliżając ziemię jako kulę? W szczególności, gdy mamy do czynienia z lokalizacją punktów i na przykład odległościami między dużymi okręgami.

Czy są jakieś badania dotyczące średniego i najgorszego błędu w porównaniu do elipsoidy? Zastanawiam się, ile dokładności poświęciłbym, gdybym poszedł z kulą dla łatwiejszych obliczeń.

Mój szczególny scenariusz dotyczy bezpośredniego mapowania współrzędnych WGS84, tak jakby były współrzędnymi na idealnej kuli (ze średnim promieniem określonym przez IUGG) bez żadnej transformacji.

Jeff Bridgman
źródło
Czy jesteś szczególnie zainteresowany modelem sferycznym, czy interesują Cię modele elipsoidalne? Wyobrażam sobie, że wielkość błędu będzie się znacznie różnić między kulą a elipsą.
Jay Laura
2
Powiązana analiza pojawia się w tej odpowiedzi . Aby jednak uzyskać odpowiedź na swoje pytanie, musisz określić, w jaki sposób przybliżona jest Ziemia jako kula. W użyciu jest wiele przybliżeń. Wszystkie są równoznaczne z nadaniem funkcji f '= u (f, l) i l' = v (f, l) gdzie (f, l) są współrzędnymi geograficznymi kuli, a (f ', l') są współrzędnymi geograficznymi elipsoida. Patrz sekcja 1.7 („Transformacja ... elipsoidy obrotu na powierzchnię kuli”) w Bugayevskiy & Snyder, Map Projections, A Reference Manual . Taylor i Francis [1995].
whuber
Jest to podobne do wczesnej debaty nad projekcją Google / Bing EPSG 900913 (która wykorzystuje współrzędne WGS84, ale projektuje tak, jakby znajdowały się na kuli), a błędy prawdopodobnie odpowiadają początkowej odrzuceniu projekcji przez EPSG, dopóki nie ulegną presji ze strony programistów. Nie chcąc nadmiernie cię rozpraszać, śledzenie niektórych z tej debaty może zwiększyć zakres informacji w doskonałym linku dostarczonym przez whuber.
MappaGnosis
@ Jzl5325: Tak, miałem na myśli ścisłą sferę, a nie elipsoidę, zredagowałem pytanie, aby zapewnić nieco więcej kontekstu.
Jeff Bridgman
1
Myślę, że powinieneś przeczytać to: en.wikipedia.org/wiki/Haversine_formula
longtsing

Odpowiedzi:

83

Krótko mówiąc, odległość może być błędna do około 22 km lub 0,3%, w zależności od danych punktów. To jest:

  • Błąd można wyrazić na kilka naturalnych, użytecznych sposobów , takich jak (i) błąd resztkowy, równy różnicy między dwoma obliczonymi odległościami (w kilometrach), i (ii) błąd względny, równy różnicy podzielonej przez „poprawna” (elipsoidalna) wartość. Aby uzyskać liczby wygodne w użyciu, mnożę te współczynniki przez 1000, aby wyrazić względny błąd w częściach na tysiąc .

  • Błędy zależą od punktów końcowych. Ze względu na symetrię obrotową elipsoidy i kuli oraz ich dwustronne (północ-południe i wschód-zachód) symetria, możemy umieścić jeden z punktów końcowych gdzieś wzdłuż pierwszej południka (długość 0) na półkuli północnej (szerokość geograficzna między 0 a 90 ) i drugim punktem końcowym na półkuli wschodniej (długość geograficzna między 0 a 180).

Aby zbadać te zależności, narysowałem błędy między punktami końcowymi przy (lat, lon) = (mu, 0) i (x, lambda) w funkcji szerokości x między -90 a 90 stopni. (Wszystkie punkty są nominalnie na elipsoidalnej wysokości zero.) Na rysunkach rzędy odpowiadają wartościom mu w {0, 22,5, 45, 67,5} stopniach, a kolumny wartościom lambda w {0, 45, 90, 180} stopnie. To daje nam dobry widok na spektrum możliwości. Zgodnie z oczekiwaniami, ich maksymalne rozmiary to w przybliżeniu spłaszczenie (około 1/300) razy większa od osi głównej (około 6700 km) lub około 22 km.

Błędy

Błędy resztkowe

Błędy względne

Błędy względne

Działka konturowa

Innym sposobem na zwizualizowanie błędów jest naprawienie jednego punktu końcowego i pozwolenie, aby drugi się zmieniał, konturując pojawiające się błędy. Tutaj na przykład jest wykres konturowy, w którym pierwszy punkt końcowy znajduje się na 45 stopniach szerokości geograficznej północnej i 0 stopniach długości geograficznej. Tak jak poprzednio, wartości błędów podano w kilometrach, a błędy dodatnie oznaczają, że obliczenia sferyczne są zbyt duże:

Działka konturowa

Może być łatwiej czytać, gdy jest zawinięty na całym świecie:

Działka globu

Czerwona kropka na południu Francji pokazuje lokalizację pierwszego punktu końcowego.

Dla przypomnienia, do obliczeń wykorzystano kod Mathematica 8 :

WGS84[x_, y_] := GeoDistance @@ (GeoPosition[Append[#, 0], "WGS84"] & /@ {x, y});
sphere[x_, y_] := GeoDistance @@
   (GeoPosition[{GeodesyData["WGS84", {"ReducedLatitude", #[[1]]}], #[[2]], 0}, "WGS84"] & /@ {x, y});

I jedno z poleceń drukowania:

With[{mu = 45}, ContourPlot[(sphere[{mu, 0}, {x, y}] - WGS84[{mu, 0}, {x, y}]) / 1000, 
                   {y, 0, 180}, {x, -90, 90}, ContourLabels -> True]]
Whuber
źródło
23
Co za kiepska odpowiedź @whuber
Ragi Yaser Burhum,
21

Ostatnio badałem to pytanie. Myślę, że ludzie chcą wiedzieć

  1. jakiego promienia sferycznego powinienem użyć?
  2. jaki jest wynikowy błąd?

Rozsądną miarą jakości aproksymacji jest maksymalny bezwzględny błąd względny w odległości wielkiego koła

err = |s_sphere - s_ellipsoid| / s_ellipsoid

z maksimum ocenianym dla wszystkich możliwych par punktów.

Jeśli spłaszczenie f jest małe, promień sferyczny, który minimalizuje błąd, jest bardzo bliski (a + b) / 2, a wynikowy błąd wynosi około

err = 3*f/2 = 0.5% (for WGS84)

(oceniane z losowo wybranymi parami punktów 10 ^ 6). Czasami sugeruje się użycie (2 * a + b) / 3 jako promienia sferycznego. Powoduje to nieco większy błąd, err = 5 * f / 3 = 0,56% (dla WGS84).

Geodezja, której długość jest najbardziej niedoceniana przez przybliżenie sferyczne, leży w pobliżu bieguna, np. (89,1) do (89,180). Geodezja, której długość jest najbardziej zawyżona w przybliżeniu sferycznym, jest południkowa w pobliżu równika, np. (-0,1,0) do (0,1,0).

DODATEK : Oto inny sposób rozwiązania tego problemu.

Wybierz pary równomiernie rozmieszczonych punktów na elipsoidzie. Zmierz elipsoidalną odległość s oraz odległość na kuli jednostkowej t . Dla dowolnej pary punktów s / t daje równoważny promień sferyczny. Uśrednij tę ilość we wszystkich parach punktów, co daje średni równoważny promień sferyczny. Jest pytanie, jak dokładnie należy zrobić średnią. Jednak wszystkie wybory próbowałem

1. <s>/<t>
2. <s/t>
3. sqrt(<s^2>/<t^2>)
4. <s^2>/<s*t>
5. <s^2/t>/<s>

wszystkie wyszły w odległości kilku metrów od zalecanego średniego promienia IUGG, R 1 = (2 a + b ) / 3. Zatem wartość ta minimalizuje błąd RMS w obliczeniach odległości sferycznej. (Jednak powoduje to nieco większy maksymalny błąd względny w porównaniu do ( a + b ) / 2; patrz powyżej.) Biorąc pod uwagę, że R 1 może być wykorzystywany do innych celów (obliczenia powierzchni i tym podobne), istnieje dobry powód, aby trzymaj się tego wyboru do obliczania odległości.

Dolna linia :

  • Do każdego rodzaju systematycznej pracy, w której tolerujesz 1% błąd w obliczeniach odległości, użyj kuli o promieniu R 1 . Maksymalny błąd względny wynosi 0,56%. Używaj tej wartości konsekwentnie, gdy zbliżasz ziemię do kuli.
  • Jeśli potrzebujesz dodatkowej dokładności, rozwiąż elipsoidalny problem geodezyjny.
  • W przypadku obliczeń z tyłu obwiedni użyj R 1 lub 6400 km lub 20000 / pi km lub a . Powoduje to maksymalny błąd względny około 1%.

KOLEJNY DODATEK : Możesz wycisnąć nieco większą dokładność z odległości dużego koła, używając μ = tan- 1 ((1 - f ) 3/2 tanφ) (szerokość korygująca biedaka) jako szerokość geograficzną w obliczeniu wielkiego koła. Ogranicza to błąd względny od 0,56% do 0,11% (przy R 1 jako promień kuli). (Nie jest jasne, czy naprawdę warto zastosować to podejście, a nie bezpośrednio obliczać elipsoidalną odległość geodezyjną).

cffk
źródło