Poprzeczny Mercator: seria Redfearn

Artykuł Projekcja poprzeczna Mercatora ogranicza się do ogólnych cech projekcji. Ten artykuł opisuje szczegółowo jedną z (dwóch) implementacji opracowanych przez Louisa Krügera w 1912 roku; wyrażone jako szereg potęgowy w różnicy długości geograficznej od południka środkowego. Szeregi te zostały ponownie obliczone przez Lee w 1946 r., przez Redfearna w 1948 r. i przez Thomasa w 1952 r. Są one często określane jako serie Redfearn lub serie Thomasa. Ta implementacja ma ogromne znaczenie, ponieważ jest szeroko stosowana w układzie współrzędnych płaszczyzny stanu USA, w krajowych (Wielka Brytania, Irlandia i wiele innych), a także międzynarodowych systemach mapowania, w tym w układzie współrzędnych Universal Transverse Mercator (UTM). Są one również włączone do konwertera współrzędnych Geotrans udostępnianego przez Narodową Agencję Wywiadu Geoprzestrzennego Stanów Zjednoczonych. W połączeniu z odpowiednim geodezyjnym układem odniesienia , seria zapewnia wysoką dokładność w strefach mniejszych niż kilka stopni w zakresie wschód-zachód.


Wstępy I: parametry odniesienia i elipsoidy

Serie muszą być używane z geodezyjnym układem odniesienia , który określa położenie, orientację i kształt elipsoidy odniesienia . Chociaż wzory projekcji zależą tylko od parametrów kształtu elipsoidy odniesienia, pełny zestaw parametrów odniesienia jest niezbędny do powiązania współrzędnych projekcji z rzeczywistymi pozycjami w przestrzeni trójwymiarowej. Układy odniesienia i elipsoidy odniesienia związane z poszczególnymi implementacjami formuł Redfearna są wymienione poniżej . Obszerna lista ważnych elipsoid jest podana w artykule na temat Figury Ziemi .

Przy określaniu normalne jest podawanie półosi wielkiej oś równikowa ) z spłaszczeniem lub ( oś biegunowa), a czasem jedno i serie wykorzystują zamiast Ponadto używają parametrów \ , zwany spłaszczeniem i mimośrodem . Istnieją tylko dwa niezależne parametry kształtu i istnieje między nimi wiele relacji: w szczególności

Wzory projekcji obejmują również promień południka (na szerokości i \ , promień krzywizny w głównym pionie . (Główna pionowa to pionowa płaszczyzna prostopadła do płaszczyzny południka w punkcie na elipsoidzie). Promienie krzywizny definiuje się w następujący sposób:

funkcje i są jako:

Dla zwartości normalne jest wprowadzenie następujących skrótów:


Eliminacje II: odległość południkowa

Odległość południkowa

Artykuł na temat południka kilka metod obliczania odległości od równika do punktu na szerokości geograficznej wyrażenia podane poniżej w faktycznej realizacji projekcji Transverse Mercator przez OSGB. Błąd obcięcia jest mniejszy niż 0,1 mm, więc seria jest z pewnością dokładna z dokładnością do 1 mm, czyli tolerancji projektowej implementacji OSGB.

gdzie współczynniki są podane na zamówienie (kolejność mi ) przez

Odległość południka od równika do bieguna wynosi

Postać szeregu określona dla UTM jest wariantem powyższych, wykazującym wyrazy wyższego rzędu z błędem obcięcia 0,03 mm.


Odwrotna odległość południka

Ani implementacje OSGB, ani UTM nie definiują szeregów odwrotnych dla odległości południka; zamiast tego używają schematu iteracyjnego. Dla danej odległości południka najpierw ustaw, następnie wykonaj iterację za

do mm.

Inwersję można przeprowadzić za pomocą szeregu, przedstawionego tutaj do późniejszego wykorzystania. Dla danej odległości południka określ szerokość geograficzną prostowania przez

Geodezyjna szerokość geograficzna odpowiadająca to (Snyder strona 17):

gdzie do ,


Zarys metody

Normalny aspekt rzutu Mercatora na kulę o promieniu jest opisany równaniami

gdzie izometryczna szerokość geograficzna jest dana przez

Na elipsoidzie staje się izometryczna szerokość geograficzna

Z rzut ze współrzędnych geodezyjnych ( współrzędne ( zgodny Jeśli współrzędne ( używane do zdefiniowania punktu na płaszczyźnie zespolonej, to każda funkcja analityczna zdefiniuje inną konforemną projekcję Metoda Krugera polega na poszukiwaniu określonego który generuje jednolitą skalę wzdłuż centralnego południka, . Osiągnął to, badając przybliżenie szeregu Taylora ze współrzędnymi projekcji podanymi przez:

część rzeczywista do funkcji odległości południka m (Złożone) współczynniki od pochodnych , można sprowadzić do pochodnych w odniesieniu do , (nie . Pochodne są zasadzie łatwe do oszacowania, ale wyrażenia stają się bardzo skomplikowane przy wyższych rzędach z powodu skomplikowanej relacji między i ψ . Rozdzielenie części rzeczywistych i urojonych daje szereg dla i , a dalsze pochodne dają skalę i współczynniki zbieżności.


Seria w szczegółach

Ta sekcja przedstawia serię ósmego rzędu opublikowaną przez Redfearn (ale z i i różnicą długości geograficznej od południka środkowego oznaczoną zamiast ). Równoważne serie ósmego rzędu, z różnymi zapisami, można znaleźć w Snyder (strony 60–64) oraz na wielu stronach internetowych, takich jak Ordnance Survey of Great Britain.

środkowego , wyrażonej w radianach: są opracowywane na podstawie Projekcja jest zwykle ograniczona do wąskich stref (w długości geograficznej), tak że oba parametry rozszerzania są zwykle mniejsze niż około 0,1, co gwarantuje szybką zbieżność . Na przykład w każdej UTM te parametry ekspansji są mniejsze niż 0,053, a dla brytyjskiej sieci krajowej ( NGGB ) są mniejsze niż 0,09. Wszystkie bezpośrednie serie dawania , , skala , zbieżność funkcje zarówno szerokości, jak i długości geograficznej oraz parametry elipsoidy: wszystkie szeregi odwrotne dają , , , są funkcjami zarówno jak i i parametry elipsoidy.

Seria bezpośrednia

W następującym szeregu różnica długości geograficznej dowolnego punktu i długości geograficznej wybranego południka środkowego: i jest dodatnia na wschód od południka środkowego. Współczynniki W funkcjami wymienionymi . Szereg dla się do wyskalowanej odległości południka, gdy .


Szeregi odwrotne

Szeregi odwrotne obejmują kolejną konstrukcję: szerokość geograficzną punktu stopy . Biorąc pod uwagę punkt na rzucie, punkt stopy , . Ponieważ skala na środkowym południku wynosi odległość południka od równika do stopy jest równa . Odpowiednia szerokość geograficzna punktu stopy lub odwrotny szereg odległości południka, jak opisano

Oznaczając funkcje oceniane w indeksie dolnym „1”, szeregi odwrotne to:


Skala punktowa i zbieżność

Skala punktowa dla transformacji konforemnej. Można go obliczyć pod względem współrzędnych geograficznych lub projekcji. Zauważ, że szereg dla zmniejsza się do , gdy albo x . Zbieżność można również obliczyć (w radianach) na podstawie współrzędnych geograficznych lub projekcji:

Współczynniki dla wszystkich szeregów


Dokładność serii

Dokładne rozwiązanie Lee-Thompsona, zaimplementowane przez Karneya (2011), ma dużą wartość w ocenie dokładności okrojonego szeregu Redfearna. Potwierdza to, że błąd obcięcia serii Redfearn (ósmego rzędu) jest mniejszy niż 1 mm na zewnątrz do różnicy długości geograficznej wynoszącej 3 stopnie, co odpowiada odległości 334 km od środkowego południka na równiku, ale zaledwie 35 km na północy granica strefy UTM.

Seria Redfearn staje się znacznie gorsza w miarę rozszerzania się strefy. Karney omawia Grenlandię jako pouczający przykład. Środek tego długiego i cienkiego lądu znajduje się na 42 W, aw najszerszym miejscu znajduje się nie dalej niż 750 km od tego południka, podczas gdy rozpiętość długości geograficznej sięga prawie 50 stopni. Seria Redfearn osiąga maksymalny błąd 1 kilometra .

Implementacje

Implementacje podane poniżej są przykładami wykorzystania serii Redfearn. Dokumenty definiujące w różnych krajach różnią się nieco zapisem i, co ważniejsze, zaniedbaniem niektórych małych terminów. Analiza małych terminów zależy od zakresów szerokości i długości geograficznej w różnych siatkach. Istnieją również niewielkie różnice we wzorach stosowanych na odległość południka: czasami do powyższego wzoru dodaje się jeden dodatkowy składnik, ale taki składnik jest mniejszy niż 0,1 mm.

OSGB

Implementacja poprzecznego odwzorowania Mercatora w Wielkiej Brytanii jest w pełni opisana w dokumencie OSGB Przewodnik po układach współrzędnych w Wielkiej Brytanii , Załączniki A.1, A.2 i C.

układ odniesienia: OSGB36
elipsoida: Airy 1830
oś główna: 6 377 563,396
oś pomocnicza: 6 356 256,909
długość geograficzna południka środkowego: 2°W
współczynnik skali południka środkowego: 0,9996012717
początek odwzorowania: 2°W i 0°N
początek siatki rzeczywistej: 2°W i 49°N
fałszywa współrzędna wschodnia rzeczywistego początku siatki, E0 (metry): 400 000
fałszywa współrzędna północna rzeczywistego początku siatki, N0 (metry): -100 000
E = E0 + x = 400000 + x
N = N0 + y -k0*m(49 °)= y - 5527063

Zasięg siatki wynosi 300 km na wschód i 400 km na zachód od południka środkowego i 1300 km na północ od fałszywego pochodzenia (OSGB sekcja 7.1), ale z wyłączeniem części Irlandii Północnej, Irlandii i Francji. Odniesienie do siatki jest oznaczone parą (E,N), gdzie E mieści się w zakresie od nieco powyżej zera do 800000m, a N w zakresie od zera do 1300000m. Aby zmniejszyć liczbę cyfr potrzebnych do podania odniesienia do siatki, siatka jest podzielona na 100-kilometrowe kwadraty, z których każdy ma dwuliterowy kod. Pozycje National Grid mogą być podane za pomocą tego kodu, po którym następuje współrzędna wschodnia i północna zarówno w zakresie 0, jak i 99999m.

Formuły projekcji różnią się nieco od przedstawionych tutaj formuł Redfearn. Zostały uproszczone przez zaniedbanie większości wyrazów siódmego i ósmego rzędu w : siódmego rzędu w szeregu dla pod względem . To uproszczenie opiera się na badaniu terminów Redfearn na podstawie rzeczywistych zasięg siatki. Jedyne inne to (a) absorpcja centralnego współczynnika skali w promieniach krzywizny i odległości południka (b) zastąpienie parametrem zdefiniowane powyżej ).

Podręcznik OSGB zawiera omówienie przekształceń Helmerta , które są wymagane do połączenia współrzędnych geodezyjnych na elipsoidzie Airy 1830 i na WGS84.

UTM

Artykuł na temat projekcji Universal Transverse Mercator zawiera ogólny przegląd, ale pełna specyfikacja jest zdefiniowana w instrukcjach technicznych US Defense Mapping Agency TM8358.1 i TM8358.2. Ta sekcja zawiera szczegółowe informacje dotyczące strefy 30 jako kolejnego przykładu formuł Redfearn (zwykle nazywanych formułami Thomasa w Stanach Zjednoczonych).

elipsoida: International 1924 (aka Hayford 1909)
oś główna: 6 378 388,000
oś podrzędna: 6 356 911,946
długość geograficzna południka środkowego: 3°W
początek rzutu: 3°W i 0°N
współczynnik skali południka środkowego: 0,9996
prawdziwy początek siatki: 3° W i 0°N
fałszywa współrzędna wschodnia rzeczywistego początku siatki, E0: 500 000
E = E0 + x = 500000 + x
półkula północna fałszywa współrzędna północna rzeczywistego początku siatki N0: 0
półkula północna: N = N0 + y = y
półkula południowa fałszywa północ prawdziwy początek siatki N0: 10 000 000
półkula południowa: N = N0 + y = 10 000 000 + y

Szeregi przyjęte dla odległości południka obejmują wyrażenia piątego rzędu w, podaje, że są one mniejsze niż 0,03 mm (TM8358.2 Rozdział 2). Wzory projekcji wykorzystują mimośrodowość (zdefiniowaną powyżej zamiast . Schematy odniesienia siatki są zdefiniowane w artykule Układ współrzędnych Universal Transverse Mercator . Deklarowana dokładność odwzorowań UTM wynosi 10 cm we współrzędnych siatki i 0,001 sekundy kątowej dla współrzędnych geodezyjnych.

Irlandia

Poprzeczna projekcja Mercatora w Irlandii i Irlandii Północnej (międzynarodowa implementacja obejmująca jeden kraj i część drugiego) jest obecnie wdrażana na dwa sposoby:

Irlandzki system odniesienia do siatki

punkt odniesienia: Irlandia 1965
elipsoida: zmodyfikowana Airy 1830
oś główna: 6 377 340,189
oś podrzędna: 6 356 034,447
współczynnik skali południka środkowego: 1,000035
początek rzeczywisty: 8° W i 53,5° N
fałszywa współrzędna wschodnia rzeczywistego początku siatki, E0: 200 000
fałszywa północ od prawdziwy początek siatki, N0: 250 000

Irlandzka siatka wykorzystuje formuły projekcji OSGB.

Irlandzki Merkator poprzeczny

punkt odniesienia: Irlandia 1965
elipsoida: GRS80
oś główna: 6 378 137
oś podrzędna: 6 356 752,314140
współczynnik skali południka środkowego: 0,999820
początek rzeczywisty: 8°W i 53,5°N
fałszywa współrzędna wschodnia rzeczywistego początku siatki, E0: 600 000
fałszywa północ względem prawdziwej siatki pochodzenie, N0: 750 000

Jest to interesujący przykład przejścia między użyciem tradycyjnej elipsoidy a nowoczesną globalną elipsoidą. Przyjęcie radykalnie różnych fałszywych początków pomaga uniknąć pomyłki między tymi dwoma systemami.

Zobacz też