Transformace geodetických souřadnicových systémů

Aktuální verze stránky ještě nebyla zkontrolována zkušenými přispěvateli a může se výrazně lišit od verze recenzované 27. září 2020; kontroly vyžadují 8 úprav .

V geodézii , úloha přechodu mezi různými souřadnicovými systémy vyvstává z existence rozmanitých souřadnicových systémů , které se vynořily po celém světě v průběhu doby. Využití různých souřadnicových systémů při řešení praktických problémů geodézie , kartografie , navigace a v geografických informačních systémech je nevyhnutelné. Existuje několik typů transformací souřadnic: přechod mezi různými formáty souřadnic , přechod mezi různými souřadnicovými systémy a mapovými projekcemi a transformace dat . Všechny tyto typy transformace budou diskutovány v tomto článku.[jeden]

Změna formátu a jednotek

Označit zeměpisné místo obvykle znamená sdělit zeměpisnou šířku a délku místa. Číselné hodnoty pro zeměpisnou šířku a délku mohou být reprezentovány v několika různých typech jednotek a formátů: [2]

sexagesimální : stupně, minuty a sekundy: 40° 26′ 46″ N 79° 58′ 56″ W

stupně a desetinné minuty: 40° 26,767′ N 79° 58,933′ zd

desetinné stupně : 40,446° N 79,982° W

Stupeň má 60 minut a minuta 60 sekund. Proto pro převod ze stupňů/minut/sekund na desetinné stupně můžete použít vzorec:

desetinné stupně=stupně+minuty/60+sekundy/3600.

Chcete-li převést zpět z desetinného formátu stupňů na formát stupňů/minut/sekund, můžete použít vzorce:

stupně = [desetinné stupně]

minuty =[60*(desetinné stupně-stupně)]

sekund = 3600*(desetinné stupně-stupně)-60*minut

kde zápis [ x ] znamená, že musíte vzít celou část x a odkázat na „ funkci police “.

Přechod mezi různými souřadnicovými systémy

Transformace souřadnicového systému  je přechod z jednoho souřadnicového systému do druhého, přičemž oba souřadnicové systémy jsou založeny na stejném geodetickém základu. Úkolem transformace je často změna z geodetického souřadnicového systému na pravoúhlé souřadnice nebo změna z jedné projekce mapy na jinou.

Od geodetického souřadnicového systému k pravoúhlému

Pravoúhlé souřadnice bodů v prostoru lze vypočítat ze známých geodetických souřadnic těchto bodů (zeměpisná šířka B, zeměpisná délka L, výška H) pomocí vzorců: [3]

kde

kde a  jsou rovníkové (hlavní poloosa) a polární poloměry (vedlejší poloosa).  je druhá mocnina první excentricity elipsoidu. poloměr křivosti první svislice je vzdálenost podél normály k elipsoidu od průsečíku plochy elipsoidu s normálou k ose oZ (obr. 1).

Od kartézského ke geodetickému

Při přechodu z pravoúhlých prostorových souřadnic do geodetického souřadnicového systému (jako je WGS84 ) je často nutné geodetické zeměpisné šířky B a výšky H vypočítat iterativně, to znamená prováděním postupných aproximací. Pokud jde o zeměpisné délky L, počítají se obvyklým způsobem.

Existuje několik metod pro výpočet geodetických zeměpisných šířek a výšek, my zvážíme dvě z nich.

Newton-Raphsonova metoda

Následující iracionální Bowringova rovnice [4] pro geodetickou šířku je řešena Newton-Raphsonovou iterační metodou : [5] [6]

kde ,

Zeměpisnou šířku B lze zjistit z rovnice .

Výška H se vypočítá takto:

Iteraci lze převést do následující podoby:

kde

Konstanta je dobrou výchozí hodnotou pro iteraci, když . Bowring ukázal, že v takových případech již první iterace poskytuje dostatečně přesné řešení. Ve své původní formulaci použil další goniometrické funkce .

Ferrari rozhodnutí

Výše uvedená rovnice pro může být vyřešena Ferrari metodou : [7] [8]

Použití Ferrariho rozhodnutí

Existuje celá řada metod a algoritmů, ale nejpřesnější je podle Zhu [9] následující sekvence stanovená Heikkinenem [10] . Předpokládá se, že geodetické parametry jsou známy.

Poznámka: arctan2 [Y, X] je zadní tečna ke čtyřem kvadrantům.

Mocninná řada

Pro malé e 2 začíná výkonová řada od

Přechod z geodetického souřadnicového systému do ENU a naopak

Převod z geodetických souřadnic na topocentrické souřadnice ENU se skládá ze dvou kroků:

  1. Převod souřadnic z geodetického systému na pravoúhlý.
  2. Převod souřadnic z pravoúhlého na topocentrický souřadnicový systém ENU.
Převod souřadnic z pravoúhlých na topocentrické souřadnice ENU

Pro převod pravoúhlých souřadnic na topocentrické souřadnice potřebujete znát počáteční bod topocentrického souřadnicového systému, obvykle se nachází v nějakém pozorovacím bodě. Pokud je pozorování provedeno v bodě a pozorovaný objekt je v, pak má vektor poloměru tohoto směru v souřadnicovém systému ENU tvar:

Transformace souřadnic z topocentrického souřadnicového systému ENU na pravoúhlý.

Inverzní transformací souřadnic z pravoúhlého systému získáme topocentrický souřadný systém:

Přepnutí na jinou projekci mapy

Převod souřadnic a poloh na mapě mezi různými projekcemi mapy , vázanými na stejný geodetický povrch , lze provést buď pomocí vzorců pro přímý přechod z jedné projekce do druhé, nebo se nejprve projekce převede na mezilehlý souřadnicový systém, například obdélníkový, a již z ní do projekce . Použité vzorce mohou být složité, v některých případech transformace nemá uzavřené tvarové řešení a je třeba použít přibližné metody. Obvykle se k provádění úloh transformace souřadnic používají počítačové programy, například s programem GEOTRANS podporovaným DoD a NGA. [jedenáct]

Počáteční transformace

Transformace mezi vztažnými body lze provádět různými způsoby. Existují transformace, které vám umožňují provést přímý přechod z geodetických souřadnic jednoho základu na geodetické souřadnice jiného základu. Existuje méně přímých přechodů, které převádějí geodetické souřadnice na geocentrické (ECEF), převádějí geocentrické souřadnice z jednoho základu na jiný a poté převádějí geocentrické souřadnice jiného údaje zpět na geodetické. Existují také promítací transformace, které umožňují provést přímý přechod z jednoho (datum, projekce) páru do druhého (datum, projekce) páru.

Promítací transformace

Transformace projekce umožňují přímý přechod ze souřadnic na mapě pro jeden pár (projekce mapy, datum) k souřadnicím na mapě pro jiný pár (projekce mapy, datum). Příkladem je metoda NADCON pro převod z 1927 North American Date (NAD) na 1983 NAD datum [12] . High Accuracy Reference Network (HARN), vysoce přesná verze transformací NADCON, má přesnost přibližně 5 centimetrů. Národní transformace verze 2 ( NTv2 ) je kanadská verze NADCON pro přechod mezi NAD 1927 a NAD 1983 . Metody HARN jsou také známé jako NAD 83/91 a High Precision Grid Networks (HPGN) [13] . Následně Austrálie a Nový Zéland pro sebe přijaly formát NTv2, aby vytvořily metody transformace projekce pro přechody mezi jejich vlastními místními daty.

Stejně jako transformace využívající více regresních rovnic, projekční metody používají k transformaci mapových souřadnic interpolaci nízkého řádu, ale ve dvou prostorech místo tří. NOAA poskytuje software (jako součást NGS Geodetic Toolkit) pro vytváření NADCON transformací. [14] [15]

Molodenského transformace

Molodenského transformace umožňuje provést přímý přechod mezi geodetickými souřadnicemi různých datových bodů bez nutnosti přechodu na geocentrické souřadnice. [16] Vyžaduje tři odsazení mezi středy souřadnicových systémů a rozdíly mezi hlavními poloosami a kompresními parametry referenčních elipsoidů.

Molodenského transformaci používá National Geospatial-Intelligence Agency (NGA) ve své bílé knize TR8350.2 a také v programu GEOTRANS podporovaném NGA. [17] Molodenského transformace byla populární před nástupem moderních počítačů a metoda je součástí mnoha geodetických programů.

Vícenásobné regresní rovnice

Transformace dat pomocí empirických metod vícenásobné regrese byly navrženy tak, aby bylo dosaženo větší přesnosti pro malé geografické oblasti než standardní Molodenského transformace. Transformační data se používají k převodu místních dat, které jsou generovány pro kontinenty nebo menší regiony, na globální data, jako je WGS 84 . [18] NIMA TM 8350.2, příloha D [19] uvádí transformace pomocí vícenásobných regresních rovnic z několika lokálních dat do WGS 84 s přesností asi 2 metry. [dvacet]

Metoda vícenásobných regresních rovnic umožňuje přímou transformaci geodetických souřadnic bez mezilehlého převodu na geocentrické souřadnice. Geodetické souřadnice v novém základu B jsou modelovány jako polynomy až do devátého stupně v geodetických souřadnicích původního základu A. Přírůstek lze například rozložit jako (zobrazeno je pouze kvadratické rozšíření):

kde

for a podobné rovnice jsou sestaveny. Při dostatečném počtu dvojic souřadnic (A, B) pro body v obou vztažných bodech se pro dobrou statistiku používají metody vícenásobné regrese k přizpůsobení parametrů těchto polynomů. Polynomy spolu s proloženými koeficienty tvoří vícenásobné regresní rovnice.

Helmertova transformace

Použití Helmertovy transformace při přechodu z geodetických souřadnic referenčního bodu ke geodetickým souřadnicím referenčního bodu probíhá ve třech krocích:

1 Převeďte geodetické souřadnice počátku na geocentrické;

2 Použití Helmertovy transformace s vhodnými transformačními parametry pro přechod z geocentrických výchozích souřadnic na geocentrické výchozí souřadnice ;

3 Převod geocentrických souřadnic na geodetické souřadnice pro základ .

Pro geocentrické souřadnice XYZ má Helmertova transformace tvar: [21]

Helmertova transformace je sedmiprvková transformace se třemi offsetovými parametry , třemi rotačními parametry a jedním měřítkovým parametrem . Helmertova transformace je přibližná metoda, kterou lze považovat za přesnou pouze tehdy, když jsou parametry transformace malé ve srovnání s hodnotami vektorů geocentrického souřadnicového systému. Za těchto podmínek lze transformaci považovat za vratnou. [22]

Čtrnáctiparametrová Helmertova transformace s lineární časovou závislostí každého parametru může být použita k pozorování časových změn zeměpisných souřadnic v důsledku geomorfologických procesů , jako je kontinentální drift [23] a zemětřesení . [24] Byl převeden na software, jako je nástroj Horizontal Time Dependent Positioning (HTDP) v softwaru US NGS. [25]

Molodensky-Badekasova transformace

K oddělení offsetů a rotací Helmertovy transformace lze použít tři další parametry, aby se nový střed rotace XYZ dostal blíže k transformovaným souřadnicím. Tato desetiparametrová transformace se nazývá Molodensky-Badekasova transformace a neměla by být zaměňována s jednodušší transformací Molodenského .

Stejně jako při použití Helmertovy transformace se i použití Molodenského-Badekasovy transformace skládá ze tří kroků:

  1. Převod geodetických souřadnic referenčního bodu na geocentrické.
  2. Použití Molodensky-Badekasovy transformace s vhodnými transformačními parametry pro přechod z geocentrických výchozích souřadnic na geocentrické výchozí souřadnice .
  3. Převeďte geocentrické souřadnice na geodetické souřadnice pro základ .

Transformace má tvar [26] :

kde  je počátek obrácení a transformace měřítka a  je faktor měřítka.

Transformace Molodensky-Badekas se používá k převodu místních geodetických údajů na globální údaje, jako je WGS 84 . Na rozdíl od Helmertovy transformace je Molodensky-Badekasova transformace nevratná, protože počátek obrácení se vztahuje k původnímu datu.

Viz také

Literatura

  1. Roger Foster Dan Mullaney Článek 018 základní geodézie: Konverze a transformace (4. března 2014). Získáno 9. prosince 2019. Archivováno z originálu dne 27. listopadu 2020.
  2. Velká Británie Ordnance Survey. souřadnicový transformátor . Získáno 9. prosince 2019. Archivováno z originálu 12. srpna 2013.
  3. B. Hofmann-Wellenhof H. Lichtenegger J. Collins. GPS - teorie a praxe. — 282 s. — ISBN 3-211-82839-7 .
  4. Bowring BR Transformace z prostorových na zeměpisné souřadnice // Surv. Rev.. - 1976. - V. 23 , č. 181 . - S. 323-327 . - doi : 10.1179/003962676791280626 .
  5. Fukušima, T. Rychlá transformace z geocentrických na geodetické souřadnice  //  J. Geod. : deník. - 1999. - Sv. 73 , č. 11 . - S. 603-610 . - doi : 10.1007/s001900050271 . (Příloha B)
  6. Sudano, JJ (1997). „Přesný převod ze souřadnicového systému zaměřeného na Zemi na zeměpisnou šířku, délku a výšku“. doi:10.1109/NAECON.1997.622711
  7. Přímá transformace z geocentrických na geodetické souřadnice // Vermeille, HH J. Geod .. - 2002. - T. 76 . - S. 451-454 . - doi : 10.1007/s00190-002-0273-6 .
  8. Irene PoloBlanco Gonzalez-Vega. Symbolická analýza Vermeilleho a Borkowského polynomů pro transformaci 3D kartézských souřadnic na geodetické souřadnice // J. Geod.. - 2009. - V. 83 . - S. 1071-1081 . - doi : 10.1007/s00190-009-0325-2 .
  9. J.Zhu. Převod Země-středných Země-pevných souřadnic na geodetické souřadnice // IEEE Transactions on Aerospace and Electronic Systems. - 1994. - T. 30 . - S. 957-961 . - doi : 10.1109/7.303772 .
  10. M.Heikkinen. Geschlossene formeln zur berechnung räumlicher geodätischer koordinaten aus rechtwinkligen koordinaten // Z. Vermess .. - 1982. - T. 107 . - S. 207-211 .
  11. MSP GEOTRANS 3.3 (Geographic Translator) (downlink) . NGA: Pobočka analýzy souřadnicových systémů. Získáno 9. prosince 2019. Archivováno z originálu 15. března 2014. 
  12. Nápověda ArcGIS 10.1: Metody založené na mřížce . ESRI. Získáno 9. prosince 2019. Archivováno z originálu dne 4. prosince 2019.
  13. Metoda posunu data NADCON/HARN . bluemarblegeo.com. Datum přístupu: 9. prosince 2019. Archivováno z originálu 6. března 2014.
  14. NADCON - Verze 4.2 . NOAA. Získáno 9. prosince 2019. Archivováno z originálu dne 6. května 2021.
  15. Donald M. Mulcare. NGS Toolkit, Část 8: Národní geodetický průzkum NADCON Tool (nedostupný odkaz) . Profesionální časopis zeměměřičů. Archivováno z originálu 6. března 2014. 
  16. Nápověda ArcGIS 10.1: Metody založené na rovnicích . ESRI. Získáno 9. prosince 2019. Archivováno z originálu dne 4. prosince 2019.
  17. Počáteční transformace . National Geospatial-Intelligence Agency. Získáno 9. prosince 2019. Archivováno z originálu dne 9. října 2014.
  18. Uživatelská příručka o transformacích vztažných bodů zahrnujících WGS 84(3. vyd.), Zvláštní publikace č. 60, Monaco: International Hydrographic Bureau, srpen 2008 , < https://web.archive.org/web/20160412230130/http://www.iho.int/iho_pubs/standard/S60_Ed3Eng.pdf > . Staženo 10. ledna 2017. . 
  19. ODDĚLENÍ GEODETICKÉHO SYSTÉMU OBRANY SVĚTA 1984 Jeho definice a vztahy s místními geodetickými systémy . Národní agentura pro snímkování a mapování (NIMA). Získáno 9. prosince 2019. Archivováno z originálu 11. dubna 2014.
  20. Taylor Chuck. Vysoce přesné referenční transformace . Datum přístupu: 9. prosince 2019. Archivováno z originálu 4. ledna 2013.
  21. Rovnice používané pro transformace vztažných bodů . Informace o zemi Nový Zéland (LINZ). Datum přístupu: 9. prosince 2019. Archivováno z originálu 6. března 2014.
  22. Geomatics Guidance Note číslo 7, část 2 Převody a transformace souřadnic včetně vzorců (odkaz není k dispozici) . Mezinárodní asociace producentů ropy a zemního plynu (OGP). Archivováno z originálu 6. března 2014. 
  23. Paul Bolstad. Základy GIS, 4. vydání . — Atlasové knihy. — 93 str. - ISBN 978-0-9717647-3-6 .
  24. Dodatek k NIMA TR 8350.2: Implementace světového geodetického systému 1984 (WGS 84) Referenční rám G1150 . National Geospatial-Intelligence Agency. Získáno 9. prosince 2019. Archivováno z originálu 11. května 2012.
  25. HDDP – horizontální časově závislé polohování . US National Geodetic Survey (NGS). Získáno 9. prosince 2019. Archivováno z originálu dne 25. listopadu 2019.
  26. Molodensky-Badekas (7+3) Transformace . National Geospatial Intelligence Agency (NGA). Získáno 9. prosince 2019. Archivováno z originálu 19. července 2013.

Poznámky