UNIVERZITA OBRANY Fakulta vojenských technologií Základy matematické kartografie (Skripta) Autor: plukovník doc. Ing. Václav TALHOFER, CSc. BRNO 2007 Talhofer, V.: Základy matematické kartografie Oponenti: prof. Ing. Bohuslav VEVERKA, DrSc, katedra mapování a kartografie, Fakulta stavební, ČVUT Praha Ing. Petr BUCHAR, CSc, katedra mapování a kartografie, Fakulta stavební, ČVUT Praha Skripta byla schválena na zasedání katedry dne 12. října 2007 ISBN: 978-80-7231-297-9 2 Talhofer, V.: Základy matematické kartografie Obsah 1. Referenční plochy a souřadnicové soustavy..............................................................................................8 LI Referenční plochy.....................................................................................................................................8 1.1.1 Referenční elipsoid..........................................................................................................................8 1.1.2 Referenční koule..............................................................................................................................9 1.1.3 Referenční rovina.......................................................................................................................... 10 1.2 Souřadnicové soustavy............................................................................................................................10 1.2.1 Souřadnicové soustavy na referenčním elipsoidu......................................................................... 11 1.2.1.a Výpočet délky poledníkového a rovnoběžkového oblouku......................................................................12 1.2.1. b Izometrické souřadnice.............................................................................................................................14 1.2.2 Souřadnicové soustavy na referenční kouli................................................................................... 15 1.2.2. a Určeni polohy kartografického pólu........................................................................................................18 1.2.3 Souřadnicové soustavy v zobrazovací rovině................................................................................20 2. Dělení a klasifikace zobrazení..................................................................................................................21 2.1 Základní transformace mezi referenčními plochami a rovinnými souřadnicovými systémy...................21 2.2 Základní vlastnosti jednoduchých zobrazení..........................................................................................22 2.3 Základní vlastnosti nepravých zobrazení................................................................................................24 2.4 Základní charakteristiky obecných zobrazení.........................................................................................25 2.5 Klasifikace zobrazení podle zkreslení.....................................................................................................25 3. Zákony zkreslení........................................................................................................................................27 3.1 Délkové zkreslení....................................................................................................................................27 3.1.1 Délkové zkreslení na referenční kouli...........................................................................................29 3.1.2 Extrémní délkové zkreslení............................................................................................................30 3.1.3 Extrémní délkové zkreslení na referenční kouli.............................................................................32 3.2 Uhlové zkreslení......................................................................................................................................32 3.2.1 Úhlové zkreslení na referenční kouli.............................................................................................35 3.2.2 Extrémní úhlové zkreslení.............................................................................................................35 3.3 Plošné zkreslení......................................................................................................................................37 3.4 Zákony zkreslení při užití polárních rovinných souřadnic......................................................................38 3.5 Vizualizace průběhu zkreslení................................................................................................................40 4. Teorie zobrazení........................................................................................................................................43 4.1 Ekvidistantní zobrazení...........................................................................................................................43 4.2 Ekvivalentní zobrazení............................................................................................................................44 4.3 Konformní zobrazení..............................................................................................................................45 4.4 Konformní zobrazení geodetické čáry....................................................................................................47 4.4.1 Průběh obrazu geodetické čáry v rovině konformního zobrazení.................................................47 4.4.2 Směrová korekce geodetické čáry.................................................................................................49 4.4.1 Délková korekce geodetické čáry..................................................................................................55 5. Zobrazení referenčního elipsoidu na kouli..............................................................................................57 .5.1 Základní vztahy a vzorce.........................................................................................................................57 5.2 Zobrazení se zachovanými zeměpisnými souřadnicemi..........................................................................58 5.3 Konformní zobrazení elipsoidu na kouli.................................................................................................59 6. Jednoduchá válcová zobrazení.................................................................................................................61 6.1 Základní vztahy a vzorce.........................................................................................................................61 6.2 Ekvidistantní válcové zobrazení..............................................................................................................63 6.3 Ekvivalentní válcově zobrazení...............................................................................................................65 6.4 Konformní válcové zobrazeni.................................................................................................................66 6.5 Šikmá poloha válcového zobrazení.........................................................................................................68 7. Jednoduchá kuželová zobrazení...............................................................................................................68 7.1 Základní vztahy a vzorce.........................................................................................................................69 7.2 Ekvidistantní kuželové zobrazení............................................................................................................ 71 7.2.1 Ekvidistantní kuželové zobrazení s jednou nezkreslenou rovnoběžkou......................................... 72 7.2.2 Ekvidistantní kuželové zobrazení se dvěmi nezkreslenými rovnoběžkami..................................... 75 7.3 Ekvivalentní kuželové zobrazení............................................................................................................. 78 7.3.1 Ekvivalentní kuželové zobrazení s jednou nezkreslenou rovnoběžkou.......................................... 79 7.3.2 Ekvivalentní kuželové zobrazení se dvěmi nezkreslenými rovnoběžkami......................................80 7.4 Konformní kuželové zobrazení................................................................................................................81 3 Talhofer, V.: Základy matematické kartografie 7.4.1 Konformní kuželové zobrazení s jednou nezkreslenou rovnoběžkou.............................................82 7.4.2 Konformní kuželové zobrazení se dvěmi nezkreslenými rovnoběžkami.........................................82 7.5 Šikmá poloha kuželového zobrazení.......................................................................................................84 8. Jednoduchá azimutální zobrazení............................................................................................................84 8.1 Základní vztahy a vzorce.........................................................................................................................85 8.2 Ekvidistantní azimutální zobrazení.........................................................................................................86 8.3 Ekvivalentní azimutální zobrazení..........................................................................................................88 8.4 Konformní azimutální zobrazení.............................................................................................................90 8.5 Azimutální projekce................................................................................................................................92 8.5.1 Gnomonická projekce....................................................................................................................93 8.5.2 Stereografická projekce.................................................................................................................94 8.5.3 Ortografická projekce...................................................................................................................95 9. Nepravá zobrazení.....................................................................................................................................96 9.1 Nepravá válcová zobrazení.....................................................................................................................96 9.1.1 Nepravá válcová sinusoidální zobrazení.......................................................................................97 9.1.1. a Mercator-Sansonovo (Flamsteedovo) zobrazení.......................................................................................97 9.1.1 .b Eckertovo sinusoidální zobrazení...........................................................................................................100 9.1.2 Nepravá válcová eliptická zobrazení........................................................................................... 101 9.1.2. a Mollweidovo zobrazení..........................................................................................................................101 9.2 Nepravá kuželová zobrazení................................................................................................................. 103 9.2.1 Bonneovo nepravé kuželové zobrazení........................................................................................ 105 9.3 Nepravá azimutální zobrazení.............................................................................................................. 106 9.3.1 Werner-Stabovo nepravě azimutální zobrazení........................................................................... 107 9.3.2 Ginzburgovo zobrazení............................................................................................................... 108 9.3.3 Modifikovaná azimutální zobrazení............................................................................................ 109 9.3.3. a Aitovovo nepravé azimutální zobrazení.................................................................................................109 9.3.3.b Hammerovo zobrazení............................................................................................................................110 9.3.3.c Wagnerovo zobrazení.............................................................................................................................111 9.4 Polykáni cká zobrazení.......................................................................................................................... 112 10. Gaussovo zobrazení.................................................................................................................................115 10.1 Základní charakteristiky zobrazení.................................................................................................. 116 10.2 Zobrazovací rovnice......................................................................................................................... 119 10.2.1 Zobrazovací rovnice UTM...................................................................................................... 122 10.3 Inverzní funkce k zobrazovacím rovnicím........................................................................................ 122 10.3.1 Inverzní funkce k zobrazovacím rovnicím UTM.................................................................... 126 10.4 Meridiánová konvergence................................................................................................................ 126 10.4.1 Meridiánová konvergence v UTM.......................................................................................... 129 10.5 Zákony zkreslení............................................................................................................................... 129 10.5.1 Zákony zkreslení v UTM......................................................................................................... 132 10.6 Směrová a délková korekce geodetické čáry.................................................................................... 133 10.6.1 Směrová korekce v zobrazení UTM........................................................................................ 136 10.7 Délková korekce geodetické čáry..................................................................................................... 136 10.8 Mezipásmové transformace.............................................................................................................. 137 11. Křovákovo zobrazení..............................................................................................................................137 11.1 Základní charakteristiky zobrazení.................................................................................................. 138 11.2 Zobrazovací rovnice......................................................................................................................... 138 11.2.1 Zobrazení referenčního elipsoidu na referenční kouli............................................................ 138 11.2.2 Transformace zeměpisných souřadnic na referenční kouli na kartografické souřadnice...... 138 11.2.3 Transformace do zobrazovací roviny..................................................................................... 140 11.2.4 Převod rovinných polárních souřadnic na pravoúhlé............................................................ 141 11.3 Inverzní funkce k zobrazovacím rovnicím........................................................................................ 141 11.4 Meridiánová konvergence................................................................................................................ 142 11.5 Zákony zkreslení............................................................................................................................... 143 12. Používaná zobrazení v Armádě České republiky a v NATO...............................................................144 12.1 Zobrazení UTM................................................................................................................................ 144 12.2 Zobrazení UPS.................................................................................................................................144 12.2.1 Zobrazovací rovnice zobrazení UPS...................................................................................... 145 12.2.2 Inverzní funkce k zobrazovacím rovnicím............................................................................... 146 12.3 Lambertovo konformní kuželové zobrazení...................................................................................... 148 13. Transformace zobrazení.........................................................................................................................149 13.1 Prostorové transformace.................................................................................................................. 151 4 Talhofer, V.: Základy matematické kartografie 13.1.1 Prostorové pravoúhlé souřadnice.......................................................................................... 151 13.1.2 Tříprvková prostorová transformace...................................................................................... 153 13.1.3 Sedmiprvková prostorová transformace................................................................................. 154 13.1.4 Moloděnského transformace................................................................................................... 154 13.1.5 Zjednodušená Moloděnského transformace........................................................................... 155 13.2 Rovinné transformace......................................................................................................................155 13.2.1 Shodnostní transformace........................................................................................................155 13.2.2 Podobnostní transformace...................................................................................................... 156 13.2.3 Afinní transformace................................................................................................................ 156 13.2.4 Interpolační metody................................................................................................................ 157 14. Aplikace zobrazení v nástrojích GIS.....................................................................................................158 14.1 Volba geodetického referenčního systému....................................................................................... 158 14.2 Transformace mezi geodetickými referenčními systémy.................................................................. 159 14.3 Volba zobrazení............................................................................................................................... 160 14.4 Vizualizace matematických prvků.................................................................................................... 161 5 Úvod Základním úkolem geoinformatiky j e vytváření a správa modelů krajinné sféry, jejích objektů a jevů. Zabývá se vztahem skutečnosti a jejím modelem ve všech aspektech, které s touto činností souvisí a je zároveň chápána jako vědní obor i praktická činnost. Geoinformatika je široký obor, který vychází z řady vědních oborů a má na ně silné vazby. Z hlediska klasické kartografie se geoinformatiky zabývá: • naukou o mapách, která zahrnuje historii kartografie, tříděním a dokumentací map a atlasů, jejich povšechným studiem; • kartografickou vizualizací, která řeší otázky kartografického jazyka a generalizace obsahu; • kartografickou tvorbou - vlastním zpracováváním obsahu map; • kartografickou polygrafií a reprografií, tedy způsoby rozmnožování map; • kartometrií a kartografickým výzkumem - způsoby analýzy obsahu map a syntézy zjištěných výsledků. Digitální geoinformatika vychází z obecné informatiky a zabývá se zejména: • definováním objektů a jevů a jejich vztahů v geografické realitě; • aplikacemi databázových přístupů k tvorbě digitálních geodatabází; • datovou analýzou; • prezentací dat a způsoby j ej ich zobrazování; Klasická i digitální geoinformatika se zabývá i řízením celého procesu modelováním včetně zjišťování a objektivizací uživatelských potřeb na vytvářené modely. Všechny modelované objekty a jevy je nutné mít lokalizovány na povrchu Země či v jejím blízkém okolí. Základní lokalizace je především otázkou topografického nebo tematického mapování zpravidla ve výchozím referenčním rámci, který je dán zvolených geodetickým referenčním systémem. Při jejich vizualizaci (zpravidla grafické trvalé nebo virtuální) je však nutné zvolit jeho rovinné zobrazení. Metodami zobrazování geodetických systémů do roviny se zabývá matematická kartografie. Matematická kartografie je tedy částí kartografie a obecně geoinformatiky zabývající se matematickými a geometrickými základy kartografických děl v obecném slova smyslu. Matematická kartografie studuje proces transformace prostorových souřadnic objektů a jevů na referenčních plochách do roviny. Zkoumá jeho zákonitosti, zkreslení, která při transformacích vznikají, jejich prostorové závislosti a poskytuje i metodiku výběru vhodných transformací pro modelovaná území. Matematická kartografie se zabývá i speciálními úkoly, jako je rovinné zobrazování bodů, čar a ploch, které se uplatňují například při zobrazování stran trigonometrických sítí, drah letadel, raket a kosmických těles, drah šíření elektromagnetických signálů radiotechnických prostředků apod. Výsledkem matematické kartografie jsou kartografická zobrazení (krátce zobrazeni) jako matematický aparát pro výše uvedené transformace. Součástí kartografických zobrazení j sou i charakteristiky zkreslení, které při transformaci prostorových souřadnic do roviny vznikají. Tyto studijní texty jsou určeny ke studiu základů matematické kartografie studované v rámci předmětu kartografie v bakalářském studijním programu vojenské technologie v oboru vojenská geografie a meteorolog. Mohou být však využity i pro jiné obory, které se zabývají teorií a praxí kartografických zobrazení. Talhofer, V.: Základy matematické kartografie Texty jsou členěny do 14 kapitol včetně úvodní kapitoly. Ve první kapitole jsou souhrnně definovány používané referenční plochy a jsou zde definovány základní souřadnicové soustavy na těchto referenčních plochách a v zobrazovacích rovinách. Druhá kapitola je věnována základních vlastnostem jednotlivých zobrazení a klasifikaci těchto zobrazení. Tato kapitola je zde zařazena i z terminologických důvodů, protože v cizojazyčné literatuře se lze setkat i s jinými názvy uváděných zobrazení a projekcí. Stěžejní kapitolou pro pochopení celé matematické kartografie je třetí kapitola, věnovaná zákonům zkreslení. Jsou zde vysvětleny příčiny zkreslení daných transformacemi prostorových těles (elipsoidu, koule) do roviny. Na tuto kapitolu navazuje kapitola vysvětlující princip odvozování zobrazovacích rovnic jednotlivých druhů a typů zobrazení. V páté až deváté kapitole jsou uvedeny jednotlivé druhy zobrazení, které jsou používány především v praxi při tvorbě map menších měřítek, zpravidla nástěnných a atlasových, kdy se jako výchozí referenční plocha většinou používá koule. Desátá a jedenáctá kapitola jsou věnovány zobrazením používaným při tvorbě státního mapového díla v České republice, závazných geoinformačích systémů (GIS) a v geodetické praxi. Jsou uváděny jak celosvětový systém WGS84 a jeho zobrazení UTM, tak i systém S-JTSK a Křovákovo zobrazení. Tyto kapitoly navazují na předmět geodézie. Dvanáctá kapitola je věnována používaným zobrazením v Armádě České republiky a v NATO. Předposlední kapitola je zaměřena na transformaci zobrazení mezi sebou. Poslední kapitola se zabývá některými aplikace matematické kartografie v programových prostředcích geografických informačních systémech se zaměřením na systém ArcGIS® firmy ESRI. Ve studijních textech nejsou vzhledem k jej ich zaměření uvedeny podrobnější informace týkající se zejména zobrazování křivek a čar v konformních zobrazeních s aplikací na Gaussovo zobarzení a zobrazení UTM. Stejně tak řada použitých vzorců není plně odvozena. K jejich bližšímu studiu je možné využít některé materiály uvedené v seznamu literatury. V textu jsou některé vybrané termíny uváděny i v anglické verzi (kurzívou v závorce za českým termínem). Důvodem bylo jak obecná znalost anglické terminologie z oblasti matematické kartografie, tak i jejich používání v programových nástrojích GIS. plk. doc. Ing. Václav Talhofer, CSc. 7 Talhofer, V.: Základy matematické kartografie 1. Referenční plochy a souřadnicové soustavy Fyzický povrch zemského tělesa je velice složitý a členitý a v modelech krajinné sféry je těžko zobrazitelný. Proto je pro vytváření těchto modelů nahrazován topografickou plochou, která je spojitou plochou vyhlazující mikrostrukturu a ty terénní tvary, které jsou z hlediska rozlišovací úrovně modelu bezvýznamné. Topografická plocha je však stále poměrně složitá pro přímé zobrazování do map nebo pro definování digitálních modelů. Pro účely mapování a tvorby modelů terénu se tato plocha nahrazuje referenčními plochami, které jsou jednodušší a jsou matematicky nebo fyzikálně přesně definované. Tyto referenční plochy jsou potom součástí definovaného geodetického referenčního systému (Datum, Geographic Coordinate System). 1.1 Referenční plochy Referenční plochou pro výšková měření je geoid. Geoid je definován jako plocha, na které všechny body mají stejný geopotenciál a která nejlépe odpovídá nerušené střední hladině světových moří, protažené i pod kontinenty. Tato plocha je ve všech bodech kolmá na směr tíže. Protože geoid je definován jako fyzikální těleso, jeho matematické vyjádření je značně složité. Pro potřeby praktické geodézie, mapování, kartografie i celé geoinformatiky je proto nahrazován referenčním elipsoidem (spheroid), referenční koulí (sphere) nebo i referenční rovinou. Vztahy mezi fyzickým povrchem Země, geoidem, resp. kvazigeoidem a elipsoidem jsou znázorněny na následujícím obrázku (Obr. 1-1). geoid, resp. kvazigeoid—-—^^ H/y e/ipso/ciN. i/\ 'tížnice\\ \\ fyzický povrch Země Obr. 1-1 Vztahy mezi fyzickým povrchem Země, geoidem, resp. kvazigeoidem a elipsoidem 1.1.1 Referenční elipsoid Výchozí referenční plochou v matematické kartografii je rotační elipsoid. Parametry rotačního elipsoidu jsou voleny tak, aby v maximální míře nahrazoval geoid v zájmové části Země nebo aby nahrazoval celý geoid. Elipsoid je plně definován dvěma parametry, kterými mohou být: • a, b - velikost hlavní a vedlejší poloosy (semimajor axis, semiminor axis), • a, e - velikost hlavní poloosy a numerická výstřednost (excentricita, eccentricity), • a, e' - velikost hlavní poloosy a druhá excentricita, • a, f- velikost hlavní poloosy a zploštění (flattening). Mezi jednotlivými parametry platí vztahy ( 1-1 ): 8 Talhofer, V.: Základy matematické kartografie a2-b2 a2 a2-b2 b2 a - b a (1-1) Referenční elipsoidy jsou jako výchozí referenční plocha používány zejména tehdy, pokud je nutné definovat zobrazení s minimálními hodnotami zkreslení rovinného obrazu. Tento způsob se volí u kartografických zobrazení používaných při definici státních souřadnicových systémů nebo mezinárodních systémů. Současně se používá i při tvorbě státních mapových děl. Do současné doby byla odvozena řada referenčních elipsoidů. Na území České republiky se používá pro civilní státní mapová díla Besselův elipsoid, pro bývalé vojenské topografické mapy v souřadnicovém systému S-1942/83 (používané do roku 2005) elipsoid Krasovského a pro současné vojenské mapové dílo a pro celosvětový systém WGS84 elipsoid WGS84.. Parametry uvedených elipsoidů jsou uvedeny v následující tabulce (Tabulka 1-1): Tabulka 1-1 Parametry referenčních elipsoidů používaných na území České republiky Elipsoid Besselův Krasovského WGS84 (GRS80) Velká poloosa a [m] 6 377 397,1550 6 378 245 6 378 137 Malá poloosa b [m] 6 356 078,9629 6 356 863,0188 6 356 752,3142 Druhá mocnina excentricity -e2 0,006 674 372 2 0,006 693 421 6 0,006 694 380 Druhá mocnina druhé excentricity -e'2 0,006 719 218 8 0,006 738 525 4 0,006 739 496 7 Reciproká hodnota zploštění M f 299,152 812 853 298,300 003 2 298,257 223 6 Poznámka: Elipsoid GRS80 je součástí geodetického referenčního systému ETRS-89, který se též používá v rezortu Českého úřadu zeměměřického a katastrálního. Jeho parametry jsou v rámci v tabulce uváděné přesnosti prakticky shodné s elipsoidem WGS84. 112 Referenční koule Není-li vyžadována vysoká přesnost prostorové lokalizace modelovaných objektů a jevů, je často používána jako referenční plocha koule. Uplatňuje se zejména při tvorbě map malých měřítek, při vizualizaci digitálních dat s menšími nároky na minimalizaci zkreslení a při řešení jednodušších navigačních úloh. Zvláštním případem je použití referenční koule při tzv. dvojitém zobrazení, kdy je referenční elipsoid nejprve zobrazen na kouli, která se poté zobrazuje do roviny. Tento postup je používán zejména při obecné poloze konstrukční osy zobrazení. Poloměr referenční koule je možné volit na základě různých hledisek. Je-li zobrazované území rozloženo podél rovnoběžky o zeměpisné šířce (po, je vhodné zvolit poloměr koule rovný příčnému poloměru křivosti elipsoidu (1-2 ): R = Nn (1-2 ) Při tomto řešení zůstává zachována původní délka rovnoběžky g^na elipsoidu (Obr. 1-2). 9 Talhofer, V.: Základy matematické kartografie Pro území kruhového tvaru se volí poloměr koule rovný střednímu poloměru křivosti rovnoběžky procházející jeho těžištěm ( 1-3 ): R = ylM0No (l-3^ Obě tělesa se poté v okolí těžiště velmi těsně ve všech směrech přimykají (Obr. 1-3). F 3S Ps Po Po /^"** / \ R=No /[
\ \\ \\ \\ rovník \ \ / \ R=MoNo S \ \ \ \ \ \ \ \^ \
i \ ' \ \ i \ N i \ 1 \ f t / i / i / -D* sin6^9 + iľ* sin8#>-
(1-13)
V následující tabulce (Tabulka 1-2) jsou uvedeny hodnoty koeficientů pro používané elipsoidy v ČR:
Tabulka 1-2 Koeficienty pro výpočet délky poledníkových oblouků referenčních elipsoidů používaných v ČR
Elipsoid A*[m] B* [m] C* [m] D*[m] £*[m]
Besselův 111120,61960 15988,63853 16,72995 0,02178 3,07731.10"5
Krasovského 111134,86108 16036,48027 16,82807 0,02198 3,11311.10 5
WGS84 111132,95255 16038,50866 16,83261 0,02198 3,11485.10"5
Dosadí-li se do vzorce ( 1-13 ) zeměpisná šířka pólu (cp = 90°), vypočítá se délka zemského kvadrantu. Ta bude pro:
• Besselův elipsoid 10 000 855,764 metrů,
• Krasovského elipsoid 10 002 137,497 metrů,
• elipsoid WGS84 10 001 965,729 metrů.
Poznámka: Pro určení délky metru jako desetimiliónté části zemského kvadrantu stanovil Delambre koncem 18. století rozměry elipsoidu, jehož délka kvadrantu byla 10 000 000 metrů.
Méně často je nutné stanovit i délku oblouku rovnoběžky. Poloměr rovnoběžky v zeměpisné šířce (p]t r = Ncosp. Délka jejího oblouku sr mezi dvěma body o zeměpisné délce Ai a Az (vyjádřené v radiánech) je:
13
Talhofer, V.: Základy matematické kartografie
sr = JVcos
N2 cos2 q>
(1-16)
Je možné zavést novou souřadnici jako funkci zeměpisné šířky q> . Tato souřadnice se nazývá izometrická šířka q. Její diferenciál bude:
dq
Mdcp
(1-17)
(1-18)
Ncosp
Rovnice (1-16 ) potom nabude tvaru
ds2 = N2 cos2 (p(dq2 + dA2)
a souřadnice q, A vytvoří na referenčním elipsoidu soustavu izometrických souřadnic. Bude-li dq = d A, potom na povrchu referenčního elipsoidu vznikne síť diferenciálních čtverců, jejichž velikost se bude s rostoucí zeměpisnou šířkou zmenšovat v závislosti na výrazu N2cos2(p.
Vzorec pro výpočet izometrické šířky se odvodí integrací výrazu ( 1-17 ):
Mdcp Ncosp
a(l-e2)
(l-e2sin2 , někdy i v rozsahu <-180°; 180°>, tedy obdobně jako u zeměpisných délek. Poloha počátku Vmůže být pevná nebo se může měnit v závislost na hodnotě zeměpisné šířky.
Pokud se v některých zobrazeních ztotožňují počátky obou soustav (Obr. 1-16), potom je výhodnější měřit polární úhel s od kladného směru osy X Pro transformaci mezi soustavami poté platí vtahy ( 1-29 )
x = p cos s y = psins ( 1-29)
i x
v \
X, YP P'
Yp 1 1 Xp 1 ' s / 1 Xp
0 0
Obr. 1-15 Polární souřadnicová soustava s různým Obr. 1-16 Polární souřadnicová soustava s totožným počátkem než pravoúhlá soustava počátkem jako pravoúhlá soustava
Počátek rovinných souřadnicových soustav se zpravidla volí uprostřed zobrazovaného území. Z hlediska konstrukce map, jejich používání nebo používání prostorových geoinformací je však výhodné, aby celé území leželo pouze v 1. kvadrantu. Proto se často k vypočteným souřadnicím přičítají vhodné konstanty Ax (falše northing) a Ay (falše easting) (Obr. 1-17).
20
Talhofer, V.: Základy matematické kartografie
(X) f X1 i
Y
i l o f/Y -^
AX v
' 1 (Y)
(0) ^- -^ AY
Obr. 1-17 Posun počátku pravoúhlé souřadnicové soustavy mimo zobrazované území
Poznámka: Orientace os X, Y nemusí být vždy stejná jako na předchozích obrázcích. Některé systémy, používané zejména pro státní mapy, mohou mít orientaci otočenou například o 180° (v ČR).
2. Dělení a klasifikace zobrazení
Kartografické zobrazení (map projection, projection) je dáno matematickým vyjádřením závislostí mezi zeměpisnými souřadnicemi na referenční ploše a souřadnicemi v zobrazovací rovině. Při definici uvedené závislosti je možné využít několika způsobů.
2.1 Základní transformace mezi referenčními plochami a rovinnými souřadnicovými systémy
Obr. 2-1 ukazuje možné způsoby transformace zeměpisných souřadnic na referenčních plochách do rovinných souřadnic. Výchozími souřadnicemi jsou zpravidla zeměpisné souřadnice na referenčním elipsoidu % Ä, v některých případech, zejména u maloměřítkových map, i zeměpisné souřadnice na referenční kouli U, V. Konečné souřadnice jsou vždy rovinné pravoúhlé souřadnice x, y.
REFERENČNÍ ZOBRAZOVACÍ PLOCHA ROVINA
21
Talhofer, V.: Základy matematické kartografie
Obr. 2-1 Způsoby transformace souřadnic mezi referenčními plochami a zobrazovací rovinou
V praxi se lze setkat se všemi kombinacemi transformace. Například zobrazení vojenských topografických map je přímou transformací mezi zeměpisnými souřadnicemi q>, X na rovinnými pravoúhlými souřadnicemi x, y (resp. N, E). Zobrazení základních map České republiky je naopak postupnou transformací od zeměpisných souřadnic na referenčním elipsoidu, přes zeměpisné souřadnice na referenční kouli, kartografické souřadnice, polární souřadnice k výsledným rovinným pravoúhlým souřadnicím.
Výchozí referenční plochou při kartografickém zobrazování je referenční elipsoid nebo referenční koule. Referenční elipsoid je zpravidla používán tehdy, pokud je požadavek na minimalizaci zkreslení rovinného obrazu. Využívá se zejména při zobrazení státních mapových děl, vizualizaci objektů a jevů databází státních informačních systémů apod. Referenční koule se využívá jako výchozí plocha zejména při tvorbě map menších měřítek (v atlasech, nástěnných map apod.) či při vizualizaci digitálních dat s menší rozlišovací úrovní. Referenční koule se používá též při řešení jednodušších navigačních úloh. Je ji však možné využívat při zobrazení státních mapových děl s vysokými požadavky na minimalizaci zkreslení rovinného obrazu, potom ovšem ve variantě dvojitého zobrazení (například u Křovákova zobrazení, které je popsáno v kapitole 11).
Kartografické zobrazení může být definováno geometrickou nebo matematickou cestou.
Zobrazení definovaná geometrickou cestou se odvozují z matematického popisu perspektivní projekce referenčních těles (v podstatě však výhradně koule) na plochy rozvinutelné do roviny. Tato zobrazení jsou označována jako projekce a jsou v současné době používány poměrně zřídka. V podstatě všechna dnes používaná zobrazení j sou definována matematickou cestou.
Poznámka: V anglické terminologii jsou však pod pojmem projekce (projection) uvažovány jak projekce ve významu uvedeném v předchozím textu, tak i všechna ostatní zobrazení.
Zobrazení se třídí podle různých hledisek, z nichž nej významnější jsou vlastnosti zkreslení obrazu a tvar zeměpisné sítě v rovině. Dalšími hledisky je i tvar zobrazovacích rovnic, poloha konstrukční osy, počet na sebe navazujících částí, na které je povrch zobrazován apod.
2.2 Základní vlastnosti jednoduchých zobrazení
Významnou třídou jsou zobrazení jednoduchá. Jejich charakter je možné přibližně vyjádřit pomocí geometrické představy promítání referenční plochy na plochy rozvinutelné do roviny, což se používá při rámcovém popisu zobrazení.
Poznámka: Projekce i jednoduchá zobrazení mají totožné obecné tvary zobrazovacích rovnic, proto jsou projekce často zahrnovány do třídy jednoduchých zobrazení.
Pod pojmem plochy rozvinutelné do roviny se rozumí plášť válce, kužele nebo rovina sama. Jednoduchá zobrazení se podle druhu zobrazovací roviny dělí na válcová (cylindrical), kuželová (conic) a azimutální (planar, azimuthal).
Charakter zobrazení je výrazně ovlivněn vzájemnou polohou referenční plochou a konstrukční osou zobrazovací plochy. Konstrukční osa je u válcových zobrazení osou válce, u kuželových zobrazení osou kužele a u azimutálních zobrazení normálou k tečné rovině v tečném bodě (nebo ve středu zobrazovaného území). Je-li konstrukční osa totožná s osou rotace Země, je zobrazení označováno jako pólové (normální, polar), leží-li konstrukční osa v rovině rovníku, potom je zobrazení nazýváno příčné (rovníkové, transverzální, transversal), při obecné poloze konstrukční osy se zobrazení nazývá obecné (šikmé, oblique).
22
Talhofer, V.: Základy matematické kartografie
Obrazem zeměpisné sítě jednoduchých válcových zobrazení v pólové poloze soustava vzájemně ortogonálních přímek (Obr. 2-2). Souřadnice bodů na referenční ploše se přímo transformují na rovinné pravoúhlé souřadnice. Obecné rovnice pro toto zobrazení lze vyjádřit pro referenční elipsoid vztahy (2-1 ), pro referenční kouli potom vztahy ( 2-2 ):
(2-1)
(2-2)
Obr. 2-2 Princip jednoduchého válcového zobrazení (převzato z [23])
Obraz zeměpisné sítě je u jednoduchým kuželových a azimutálních zobrazení v pólové poloze tvořen soustavou polopřímek vycházejících z jednoho bodu (poledníky) a oblouků soustředných kružnic (u kuželových zobrazení) nebo celých soustředných kružnic (u zobrazení azimutálních) s totožným středem (Obr. 2-3, Obr. 2-4). Obrazy poledníků a rovnoběžek jsou na sebe vzájemně kolmé.
Obr. 2-3 Princip jednoduchého kuželového zobrazení (převzato z [23])
Zobrazovací rovnice u obou typů zobrazení mají podobný tvar a vycházejí nejprve z transformace souřadnic na referenční ploše na rovinné polární souřadnice podle obecných vztahů ( 2-3 ) pro referenční elipsoid a ( 2-4 ) pro referenční kouli:
P = f( ,X) y = f((p,A)
budou mít tyto transformované body rovinné souřadnice P '[x, y] a Q '[x+dx, y+dyj (Obr. 3-2). Jejich vzdálenost v rovině bude opět délkovým elementem dS, jehož zeměpisný azimut v rovině bude A 'a jeho směrník ď.
Délkový element na referenční ploše lze vyjádřit podle obrázku (Obr. 3-1) rovnicí
ds2 = M2d(p2 + N2 cos2 (pdÁ2 (3-2}
Tomuto elementu bude podle (Obr. 3-2) odpovídat délkový element v rovině:
dS2 =dx2+dy2
Hodnoty dx a dy se určí jako totální diferenciály z obecných zobrazovacích rovnic ( 2-9 ):
(3-3)
dx ôx
dx = —d(p-\--d A
8
Tangenta úhlu je v intervalu 360° dvojznačná (v I. a III. kvadrantu kladná a ve II. a IV. kvadrantu záporná). Proto rovnice ( 3-24 ) určuje dva azimuty:
• 2Aa,
• 2Ab = \m°-2Aa
tedy:
• Aa,
• Ab = 90° -Aa
Azimuty Aa, A b jsou měřeny na referenčním elipsoidu. Po jejich zobrazení do zobrazovací roviny budou označeny A'a, A 'b, přitom obecně platí:
Aa* A'a,Ah*A'h
Dosadí-li se hodnoty Aa, A b do rovnice ( 3-14 ), získají se dvě rovnice pro extrémní délková zkreslení ma a mb, která jsou ve vzájemně kolmých směrech:
p
m2 =m2cos2Aa+-sin24, +m2 sin2 An (3-25)
MNcoscp
a r
30
Talhofer, V.: Základy matematické kartografie
ml = m2 cos2 Ab H--sin2^ +m2 sin2 Ab
MNcoscp
S uvážením, že Ab = 90° - Aa může mít druhá rovnice i tvar:
p
m2h = m2n sin2 An +-ún2An +m2r cos2 An ( 3"26 )
MNcoscp
Směry a, b, ve kterých jsou extrémní délková zkreslení, se nazývají hlavní paprsky zkreslení. Hlavní paprsky zkreslení jako jediné ortogonální směry na referenční ploše zůstávají vzájemně kolmé i po zobrazení do roviny.
Na obrázku (Obr. 3-3) je kružnice na referenčním elipsoidu o poloměru ds. Po jejím zobrazení do rovina se kružnice změní na elipsu, jejíž poloosy leží ve směrech hlavních paprsků zkreslení (Obr. 3-4). Velikosti jejich poloos budou a = ma ds a b = mb ds. Tato elipsa se nazývá elipsa zkreslení, též Tissotova elipsa, Tissotova indikatrix. Tvar elipsy a orientace jejích os umožňuje posuzovat hodnoty zkreslení a orientaci jeho extrémů v různých částech zobrazovaného území.
Obr. 3-3 Diferenciální kružnice na referenčním Obr. 3-4 Obraz diferenciální kružnice v zobrazovací elipsoidu rovině
Uvažuj e-li se elipsa zkreslení, lze obecný vzorec pro výpočet délkového zkreslení definovat i pomocí hodnot extrémů zkreslení a souřadnic definovaných v ortogonální soustavě hlavních paprsků zkreslení (Obr. 3-4):
dS2 =m2ada2 +m2bdb2
kde podle (Obr. 3-3) je:
da = ds cos pi db = ds sin pi
Po dosazení do ( 3-1) se získá jednodušší vzorec pro výpočet délkového zkreslení:
2 2 2 2*2 ( '-1-97 ^
m =ma cos pi + mb sin pi y J L' >
kde //je směrník uvažovaný od hlavního paprsku zkreslení a počítaný podle vztahu:
31
Talhofer, V.: Základy matematické kartografie
ju = A-Aa
3.1.3 Extrémní délkové zkreslení na referenční kouli
Pokud se jako referenční těleso uvažuje referenční koule o poloměru R, potom rovnice ( 3-24 ), ( 3-25 ) a ( 3-26 ) budou mít tvar:
tg2Aa
2F
[m2p -m2)R2cosU
m2 = m2 cos2 A„ +
F
a R2cosU
sin 24, +m2sm2 A„
F
2 2 • 2 a
m, = wi sin Aa+ —
p R2cosU
3.2 Úhlové zkreslení
-sin2A„ +m2 cos2 A„
( 3-28 ) ( 3-29 ) ( 3-30 )
Úhlově zkreslení Ac^je definováno vztahem ( 3-31 ):
Aco = co' - co (3-31)
kde co] e úhel na referenční ploše mezi dvěma směry PQ\ a PQi, co')q odpovídající úhel po zobrazení do zobrazovací roviny.
Uhel co lze však vyjádřit i jako rozdíl dvou azimutu jak na referenční ploše (co=A2 -A{), tak i v zobrazovací rovině (co'=A i-A j), jak je ukázáno na obrázcích (Obr. 3-5 a Obr. 3-6).
X Qi X' yQ'1
\__A2 \^_A'2 /
P \ \ 00 ' Q2 \ co' _ Q'2
p\
Obr. 3-5 Úhel jako rozdíl dvou azimutu na referenční Obr. 3-6 Úhel jako rozdíl dvou azimutu v zobrazovací ploše rovině
Zkreslení úhluje potom možné vyjádřit jako:
Aco = (Ä2-A\)-(A2 -Al) = (Ä2-A2)-(A\-A,) = AA2 -AAl kde AA je zkreslení azimutu vyjádřené obecným vzorcem:
AA = A-A (3-32)
32
Talhofer, V.: Základy matematické kartografie
Zkreslení azimutu je možné odvodit z obrázku (Obr. 3-7) zobrazující rovinný obraz zeměpisného poledník A a libovolného směru s, jehož azimut v zobrazovací rovině je A' a na němž diferenciálně blízku od výchozího bodu /"leží bod Q'. Podle obrázku platí:
Obr. 3-7 Azimut obecného směru v zobrazovací rovině Ä=\S0°-(ďp-ď)
Lze psát:
tgA = -tg(ďp-ď)
Podle vzorce pro tangentu rozdílu dvou úhlů lze předchozí vzorec upravit:
tgďp-ígď _ tgď-tgďp
tgA
\ + tgv'ptgď l + tgďptgď
( 3-33 )
Pro určení azimutu A 'je tedy nutné stanovit tangenty směrníků směru s a poledníku A. Podle obrázku (Obr. 3-7) lze směrník c'vyjádřit jako
dy
tgrj
dx
kde dx a dy je možné vyjádřit vztahy (3-3 ) a ( 3-4 ). Potom bude:
tgď
—dcp +—dA dq> dA
dx dx
—dcp-\--dA
dq> dA
(3-34)
Dosadí-li se za diferenciály zeměpisných souřadnic tvary ( 3-12 ), rovnice ( 3-34 ) bude mít tvar:
tgď
dy cos A . dy sin A , dy .T . dy , . . .
—--ds + —--ds —Ncos(pcosA + —MsmA
ô d A dA dq>
tgA
( 3-37 )
f dx'] ( dy^ — + — d(p \dq>
M . dx dx dy dy
—cos q> cot gA -\----1---
N dep dA dep d A
Využijí-li se Gaussovy koeficienty ( 3-7 ) ( 3-8 ):
E
fdx\2
F
dep dxdx
dy d(P.
dydy
--1--
dcpdA dcpdA
a zavede-li se čtvrtý Gaussův koeficient H
^ _ dx dy dx dy d(p dA dA d(p
potom lze výraz ( 3-37 ) psát ve tvaru:
H
( 3-38 )
tgA'--
h — cos cp cot gA + r
N
( 3-39 )
Pomocí rovnice ( 3-39 ) je tedy možné vypočítat k azimutu A obecného směru s na referenční ploše azimut A' obrazu tohoto směru v zobrazovací rovině a tím lze vypočítat i zkreslení azimutu AA. Ze znalosti zkreslení azimutu lze vypočítat i zkreslení obecného úhlu A ca
Zkreslení azimutu je možné vypočítat i z extrémních hodnot délkového zkreslení. Protože podle obrázků Obr. 3-3 a Obr. 3-4 platí:
db
tg V
da
mbdb mda
potom lze psát:
jfi
tgjU '=—^tgjU
( 3-40 )
Pokud je /u = A - Aa, lze tedy z něho vypočítat //' = A '- A 'a a tím i zkreslení směrníku /u i azimutu A podle vztahů (3-41 ) a ( 3-42 ):
AjU = JLÍ-JU
(3-41)
34
Talhofer, V.: Základy matematické kartografie
AA = Ä-A = Aju + Aa-Aa (3-42)
3.2.1 Úhlové zkreslení na referenční kouli
Rovnice (3-38)a(3-39) budou mít na referenční kouli tvar:
dx dy dx dy (3-43)
H
dU dV dV dU
tgA'=-—- ( 3"44 )
EcoslI cotgA + F Ostatní odvozené vztahy jsou beze změny.
3.2.2 Extrémní úhlové zkreslení
Z hlediska tvarů rovnic ( 3-39 ) a ( 3-40 ) v elipse zkreslení existují symetrické směry, ve kterých úhlové zkreslení dosahuje extrémů. Označí-li se tyto symetrické směry symbolem s, potom jejich směrníky v ortogonální soustavě hlavních paprsků zkreslení budou /4 a ju's a jim odpovídající azimuty AsaA 'e. Velikost extrémního zkreslení směrníku bude A/4 =Ju,s-Jus a velikost extrémního zkreslení azimutu AA =A 'e-Ae,.
Podle pravidla pro hledání extrémů funkce platí:
d(ju'-ju) = 0
z toho
dju'=dju
Diferencováním rovnice ( 3-40 ) se získá:
1 _mb 1
cos2/// macos2jUE
( 3-45 )
Protože obecně platí, že
l + tg2x
cos2 x
lze rovnici ( 3-45 ) upravit na tvar:
i+fcV=—(i+£2/0 <3-46)
a po dosazení za tg2ju'z rovnice ( 3-40 ) pro hodnotu ju'e se obdrží: , m} 2 m, i 2 \
a odtud
2,. Mb
tg M.
m
a \ma j
™L_i
35
Talhofer, V.: Základy matematické kartografie
Lze tedy vypočítat:
Pokud by se do výrazu ( 3-46 ) dosadila hodnota z výrazu ( 3-40 ) tgju, potom se získá obdobný výraz jako ( 3-47 ):
(3"48)
Vzorce ( 3-47 ) a ( 3-48 ) určují jak na referenčním elipsoidu, tak i v zobrazovací rovině čtyři symetrické směry s dané dvojznačností tangenty v intervalu 0°- 360° a kladnou a zápornou hodnotou výrazu ///„ ///,,, resp. ///,, ///„ po odmocnění. V těchto směrech dosahuje úhlové zkreslení svých extrémů.
Ze známých hodnot jUs a ju's se vypočítá na základě vzorce ( 3-41 ) velikost extrémního zkreslení směrníku
AME =Ms'-Ms
Uvedený vzorec je možné dále upravovat
a po dosazení za tgjus a tg/u'e ze vzorců ( 3-47 ) a ( 3-48 ) se získá tvar:
mh — wi
j^A^e=—5=2r (3-49)
který je ještě možné upravit podle obecného vztahu
tgx
sin x ■
2x
fi+tg na tvar
m — m
sinA//E= b a (3-50) mb +ma
Ze znalosti hodnoty AjUs je možné vypočítat hodnotu extrémního zkreslení azimutu podle vztahu ( 3-42 ):
ME=AuE+Äa-Aa (3-51)
Pomocí vyjádření extrémního zkreslení směrníku uvažovaného v ortogonální soustavě hlavních paprsků zkreslení je možné odvodit vzorec pro výpočet extrémního zkreslení obecného úhlu. K jeho odvození se využije elipsa zkreslení (Obr. 3-8), na níž jsou vyznačeny jak směry s vymezené v zobrazovací rovině směrníky ju'e, tak i odpovídající směry na referenčním elipsoidu (zde označené vymezené s), které mají směrníky jue.. Všechny čtyři dvojice svírají stejný úhel AjUs. Z obrázku je patrné, že právě úhly sevřené jednotlivými směry s budou nejvíce zkreslené a platí pro ně vztah:
36
a lze tedy vypočítat extrémní úhlové zkreslení
Talhofer, V.: Základy matematické kartografie
(3-52)
x84 a El .
S3'/
Obr. 3-8 Směry extrémů úhlového zkreslení
Při vyhodnocování extrémního úhlového zkreslení jeho znaménko nemá praktický význam a proto se většinou uvádí v absolutní hodnotě.
Dosadí-li se výraz ( 3-52 ) do výrazu ( 3-50 ), získá se po úpravě nejvíce používaný vztah pro výpočet extrémního zkreslení:
A&> m, -m„
sin-
mb +ma
(3-53)
Poznámka: Při rozborech kartografických zobrazení se úhlové zkreslení vyhodnocuje pouze výjimečně a často se uvažují pouze jeho extrémní hodnoty. Proto se zpravidla vynechává znak extrému s a symbolem Aco se rozumí přímo extrémní úhlové zkreslení.
3.3 Plošné zkreslení
Plošné zkreslení je definováno výrazem:
dP dp
m
Pi
(3-54 )
kde dp je diferenciální plocha na referenční ploše,
dP]e odpovídající diferenciální plocha v zobrazovací rovině.
Diferenciální plochu lze na referenční ploše vymezit diferenciálně blízkými poledníky a rovnoběžkami jako čtyřúhelník (Obr. 3-9). Tato plocha bude mít velikost:
dp = MN cos(pd(pdA
Po zobrazení uvedené diferenciální plochy do roviny bude mít obecně tvar rovnoběžníka, ve kterém bude úhel rovnoběžky a poledníku^ ŕ. Z obrázku (Obr. 3-10) plyne
dP = mpmrMNr cos (pd(pdA sin Ar'
Dosadí-li se výše uvedené hodnoty do výrazu ( 3-54 ), obdrží se vzorec pro výpočet plošného zkreslení:
37
Talhofer, V.: Základy matematické kartografie
m
pi
mp mrMN cosqd (pdX sin Ar MNcosgxIgxIA
m
Pi
■ mpmr sin Ar
(3-55)
Obr. 3-9 Diferenciální plocha na referenčním elipsoidu Obr. 3-10 Obraz diferenciální plochy v zobrazovací
rovině
Plošné zkreslení lze vypočítat i z poměru plochy diferenciální kružnice na referenční ploše a jejího obrazu v zobrazovací rovině - elipsy zkreslení. Plocha kruhu na obrázku (Obr. 3-3) je
dp = Tľds1
Plocha elipsy na obrázku (Obr. 3-4) potom je
dP = 7unadsmbds
Po dosazení uvedených výrazů do vzorce ( 3-54 ) lze vypočítat plošné zkreslení i ze známých hodnot extrémů délkového zkreslení
mPi=mamb
( 3-56 )
3.4 Zákony zkreslení při užití polárních rovinných souřadnic
Ke zobrazení referenční plochy do roviny a k vyjádření zákonů zkreslení se v případě kuželových a azimutálních zobrazení využívají polární rovinné souřadnice p, s, do nichž jsou transformovány výchozí zeměpisné souřadnice q>, A. Polární rovinné souřadnice jsou teprve následně transformovány do pravoúhlé souřadnicové soustavy - do souřadnic x, y. Mezi jednotlivými souřadnicemi platí vztahy (viz odstavec 2.4 , vztahy (2-11 ) a odstavec 1.2.3 , vztahy ( 1-28 )):
P = f( p -> p -> X
-> X -> s -> -> X
Schéma je možné využít jako pomůcku pro parciální derivace rovnic ( 3-58 ):
dx _ dx dxv dx dp dx ds d dp dq> ds dq>
dx _ dx dp dx ds dk dp dk ds dk
dy _dy dp ^dy ds d G = jV2cos> (4-8)
Ncoscp
Pro referenční kouli vztah ( 4-8 ) platí obdobně
4g
RcosU
l^g = R2cos2U (4"9)
Hodnoty Gaussových koeficientů E, g budou mít tvary dané parametry konkrétní referenční plochy a konkrétními zobrazovacími rovnicemi použitého zobrazení.
4.2 Ekvivalentní zobrazení
Základní vlastností ekvivalentního zobrazení je, že se při jeho použití nezkreslují plochy zobrazovaných objektů a jevů, případně jsou tyto plochy konstantně zkreslené v celém zobrazovaném území. Tato varianta je však pouze modifikací (měřítkovou změnou) základního ekvivalentního zobrazení a proto ji není nutno uvažovat jako zvláštní případ.
Podmínku zachování velikosti ploch je možné vyjádřit z obecné rovnice plošného zkreslení (viz ( 3-55 )):
mpl=\ (4-10)
kterou lze psát ve tvaru:
44
Talhofer, V.: Základy matematické kartografie
m„mr smA =1 (4-11)
p r
Vzorec ( 4-11 ) se výhodně používá i pro rozpoznávání ekvivalentních zobrazení, a to zejména pro jednoduchá. U těchto zobrazení platí A'r = 90°; výše uvedený vzorec potom nabývá tvaru:
mpmr=l (4-12)
Jiné vyjádření podmínky ekvivalentního zobrazení vychází ze vztahů pro plošné zkreslení využívající Gaussovy koeficienty. Pro ekvivalentní zobrazení z referenčního elipsoidu lze psát:
H \^>H =MNcoscp (4-13)
MN cos cp
Pro ekvivalentní zobrazení z referenční koule bude mít podmínka tvar:
H o (4-14)
:1=>H = R2cosU
R2cosU
Hodnota koeficientu H je vyjádřena podle druhu referenční plochy a typu zobrazení. Uvedená podmínka se využívá zejména u nepravých nebo obecných zobrazení.
4.3 Konformní zobrazení
Vlastností konformního zobrazení je, že nezkresluje úhly. Tuto vlastnost je možné vyjádřit vztahem:
Acv = 0 (4-15) Uhlové zkreslení je obecně dáno vzorcem (( 3-53 ), viz. odstavec 3.2 ):
. Acos mb-ma (A^f,^ sin-= -2- (4-16)
2 mb+ma
Vzhledem ke vzorci ( 3-53 ) bude podmínka (4-15 ) splněna pouze za předpokladu, že:
ma =mb (4-17 )
Hodnoty ma a ntb jsou extremními hodnotami délkového zkreslení. Pokud se mají tyto dvě hodnoty rovnat, potom je délkové zkreslení konstantní a nezávislé na směru azimutu délkového elementu. Elipsa zkreslení se tudíž zobrazuje jako kružnice.
Z obecné rovnice délkového zkreslení dané vztahem ( 3-23 ) (viz odstavec 3.1):
2 2 2 ^T * 2*2
m =m„cos A-\--sin2A + mr sin A
MNcos cp
je zřejmé, že konstantní délkové zkreslení bude pouze v případě platnosti podmínek:
1. mp = mr
2. F = 0
První podmínku lze vyjádřit i pomocí Gaussových symbolů. Pro referenční elipsoid lze psát:
45
Talhofer, V.: Základy matematické kartografie
4Ě= 4g _e= M2 (4.18)
M Ncosy G N2 cos2 q) Pro referenční kouli lze psát obdobně:
■Je _ Vg e _ i (4.19)
"^"iŕcosf/^G'cos2^
Hodnoty Gaussových koeficientů budou vyjádřeny opět podle druhu a charakteru daného zobrazení.
Konformní zobrazení lze definovat i pomocí izometrických souřadnic. V kapitole Referenční plochy a souřadnicové soustavy byly odvozeny vztahy pro výpočet izometrické šířky na referenčním elipsoidu a na referenční kouli.
Vyjde-li se z obecné rovnice délkového zkreslení ve tvaru:
m2
dS2
ds
2
a dosadí-li se za diferenciály délek jejich tvary vyjádřené v diferenciálech izometrických souřadnic (viz kapitola 1), lze psát:
m1
dx2+dy2 (4"2°)
N2 cos2 cp{dq2 +dA2)
V konformním zobrazení nesmí být zkreslení délkového elementu závislé na jeho azimutu. Bude uváženo, kdy bude rovnice ( 4-20 ) vyhovovat této podmínce. Azimut délkového elementu na referenční ploše lze vyjádřit:
N cos cpdA
tgA =-
Mdcp
neboli:
. dA dq
Směrník v zobrazovací rovině se vyjádří vztahem:
, dy dx
Je tedy zřejmé, že závislost délkového zkreslení na směru délkového elementu vyjadřují poměry diferenciálů zeměpisných (izometrických) a rovinných souřadnic. Aby zobrazení bylo konformní, nesmí se ve výrazu ( 4-20 ) uvedené diferenciály vyskytovat. To bude splněné pouze za předpokladu využití obecných zobrazovacích rovnic ve tvaru:
x + iy = f(q + iA) (4-21)
x-iy = f(q-iA) (4-22)
Přitom pro praktické použití stačí uvažovat pouze jednu z uvedených funkcí.
46
Talhofer, V.: Základy matematické kartografie
4.4 Konformní zobrazení geodetické čáry
Většina úloh praktické geodézie je řešena v rovině konformního zobrazení. Při jejich řešení je nutné znát průběh rovinného obrazu geodetických čar, kterými jsou poledníky, rovník a všechny ostatní nej kratší spojnice dvou bodů na referenčním elipsoidu, tedy i strany trigonometrických sítí. Změny průběhu uvedených geodetických čar, které jsou výsledkem jejich zobrazení do roviny konformního zobrazení, se projevují v tzv. korekcích ze zobrazení, tedy jako směrová a délková korekce geodetické čáry [14].
4.4.1 Průběh obrazu geodetické čáry v rovině konformního zobrazení
Podle [22] lze prostorové čáře definovat geodetickou křivost y jako křivost pravoúhlého průmětu této čáry do tečné roviny vdaném bodě, pro který je křivost počítána. Dále lze definovat geodetickou čáru (geodetickou křivku) jako čáru, jejíž geodetická křivost y je v každém jejím bodě rovna 0.
Protože konformní zobrazení není ortogonální, geodetické čáry se v rovině tohoto zobrazení obecně zobrazují jako křivky, jejichž křivost se v každém bodě mění. Definici změny křivosti je nutné vyjít ze situace ilustrované na obrázku (viz Obr. 4-1)
dS' _B'
Obr. 4-lPrůběh geodetické čáry v rovině konformního zobrazení
Mějme na referenčním elipsoidu diferenciální pravoúhlý čtyřúhelník PQBA omezený stranami ds a dt, které jsou částí geodetických čar. Po jeho zobrazení do roviny konformního zobrazení se sice zachovají pravé úhly u vrcholů čtyřúhelníka, avšak každá strana bude jinak zakřivena a bude délkově zkreslena. Po zobrazení budou tyto strany označeny dS, dS', dT a dT'. Na základě obrázku (Obr. 4-1) lze zkoumat křivost Y elementu dS v bodě P'. Z obrázku je zřejmé, že:
dS = -da. (4-23)
r
Pokud se zanedbá minimální rozdíl mezi délkami křivek A 'B' a A 'C, potom lze psát:
dS'=dS + dodT. (4"24) Po dosazení za duz (4-23) do (4-24) se získá vztah:
47
Talhofer, V.: Základy matematické kartografie
dS'= dS(l + TdT). Pokud se vydělí obě strany rovnice ( 4-33 ) hodnotou ds, potom bude:
*l = *(i+rv/r)
ds ds
(4-25)
(4-26)
Výraz dS/ds je výrazem pro délkové zkreslení m v bodě P' a výraz dSVds vyjadřuje délkové zkreslení v bodě A', jehož hodnota je proti hodnotě délkového zkreslení v bodě i5' změněna o dm. Tuto změnu je možné vyjádřit jako:
ds
= m + dm
(4-27)
m + dm = m(í + rdT) (4-28)
Po dosazení výrazů pro délkové zkreslení a jeho změnu do (4-26) se získá výraz: a z něho po úpravě:
1 dm
r =
(4-29 )
m dT
kde m je tedy délkové zkreslení a dm/dT]e změna zkreslení ve směru kolmém ke geodetické čáře.
Poznámky:
1. Délkové zkreslení m je počítáno pro délková element ds mezi body P a Q, který má azimute. Změnu tohoto zkreslení ve směru kolmém na element ds lze vypočítat pomocí stejného výrazu, pokud se změní hodnota azimutu a 90° v příslušném směru.
2. Je-li geodetická čára vedena kolmo k ekvideformátám, potom dm/dT = 0 a proto její obraz bude přímka. Pokud je geodetická čára ve směru ekvideformát, potom změna zkreslení v kolmém směru bude maximální a tedy bude i její obraz maximálně zakřiven (viz obrázek Obr. 4-2).
Obr. 4-2 Průběh geodetické čáry vzhledem k ekvideformátám
Délkové zkreslení se mění v závislosti na poloze konkrétního bodu na dané geodetické čáře. Pro výpočty zejména v geodetické praxi je ale nutné znát především tvar geodetické čáry na jejím počátečním a koncovém bodě a délku jejího obrazu. Tyto vlastnosti lze určit pomocí výpočtů tzv. směrové a délkové korekce geodetické čáry.
48
Talhofer, V.: Základy matematické kartografie
4.4.2 Směrová korekce geodetické čáry
Strany trojúhelníků trigonometrických sítí jsou geodetickými čarami. Na kouli jsou tyto čáry součástí hlavních kružnic, na referenčním elipsoidu jsou částmi zmíněných geodetických čar. K pochopení principu výpočtu směrové korekce v jakémkoliv zobrazení (tedy nejen v konformním) je možné využít jeden sférický trojúhelník ABC na kouli, jehož vrcholové úhly jsou a, /?, y (Obr. 4-3) (podle [15]). Při zobrazení tohoto trojúhelníka do roviny některou ze zobrazovacích metod, budou všechny úhly zkresleny a nabydou hodnot a', /?', y' (Obr. 4-4). Zkreslení těchto úhluje možné vyjádřit rovnicemi:
Aa = ď-a, A/3 = j3'-j3, Ay = y-y (4"30)
S výjimkou gnóm oni cké projekce se části oblouků hlavních kružnic zobrazí jako křivky.
Obr. 4-3 Sférický trojúhelník ABC na kouli
Obr. 4-4 Obraz sférického trojúhelníka A'B'C'v zobrazovací rovině
Ve sférickém trojúhelníku platí pro součet jeho vrcholových úhlů vztah:
a + j3 + y = lS0° + s,
kde e je sférický exces.
Pokud se řeší vztahy v trojúhelníku v rovině, což je typické pro všechny úlohy praktické geodézie, nepočítá se s křivkami, ale s přímými spojnicemi vrcholů A'B', A 'C, B'C. Součet vrcholových úhlů, které jsou mezi uvedenými přímými spojnicemi, potom bude:
a0+j80 + y0=\S0°
Přímé spojnice vrcholů trojúhelníka A 'B 'C svírají s tečnami ke křivkám (obrazům hlavních kružnic) v bodech A', B', C malé úhly ô, které se nazývají směrové korekce. Aby bylo možné trojúhelník řešit v rovině, je nutné směrníky příslušných trigonometrických stran opravit o příslušnou směrovou korekci, pomocí které se převede křivka na přímku. Podle obrázku (viz Obr. 4-4) budou využity následující směrové korekce:
• wboděA':SAB;SAC;
• v bodě 5': SBA-, SBC;
• v bodě C: óCA; óCB.
Potom lze vypočítat vrcholové úhly v rovinném trojúhelníku A'B'C pomocí následujících vztahů:
a0 = a±SA.B. ± 8A.C.
(4-31)
49
Talhofer, V.: Základy matematické kartografie
P0=P'±ÔKA,±ÔB,C,
r0=r'±ôc,A,±sc,B,
Vztahy ( 4-31 ) jsou platné pro všechny zobrazení s výjimkou gnómonické projekce.
U konformních zobrazení se úhly nezkreslují, proto všechny rozdíly úhlů počítané podle ( 4-30 ) jsou nulové a potom lze tedy psát:
a0=a±8A.B. ±8A.C.
P0=P±ÔB,A,±ÔB,C, (4-32)
n=r±8c.A.±sc.B.
Vrcholové úhly se počítají z naměřených hodnot úhlů v terénu, s nichž se získávají směrníky geodetických čar o. Příslušný vrcholový úhel je tedy získán z rozdílů směrníků jednotlivých stran daného trojúhelníku.
V další části je tedy uvedený pouze postup odvození výpočtu směrové korekce jedné strany trojúhelníka. Vzhledem k obvyklé délce stran trojúhelníků v trigonometrické síti nebo velikosti měřených délek je možné použít výše uvedené závěry i pro plochu referenčního elipsoidu. Výchozí situace je na obrázku (viz Obr. 4-5),
1 k —~iV- r.dX-~ zde bude nahrazena šířkou izometrickou q. Obecné zobrazovací rovnice potom budou mít tvar:
X = f^V (10-4)
y = f (q, X)
Jelikož zobrazení je konformní, je možné je definovat pomocí obecných rovnic konformního zobrazení odvozených v kapitole Teorie zobrazení:
x + iy=f(q + iA) (10-5)
X-jy=f(q-jÄ) (10-6)
Vzhledem k tomu, že poledníkový pás je velice úzký, je hodnota A vzhledem ke q diferenciálně malá, je možné pravou stranu rovnice ( 10-5 ) rozvést v Taylorovu řadu:
119
Talhofer, V.: Základy matematické kartografie
i2 A2 ,,„,, x ŕ A3 m) , s ŕ ?i m , s ŕ'A5
xXy = f(q) + fXq)iA + fXqý-^ + rXqy-X- + f> (q)^ + fJ (<ľ)^p +
ŕ (9)
6!
Rovnice ( 10-6 ) by se upravila stejně. Pro odvození zobrazovacích rovnic však stačí uvažovat pouze první z nich.
Uváží-li se mocniny imaginárního čísla i (i2 = -1, ŕ = -i, atd.) a oddělí-li se reálná a imaginární část, potom lze psát obecné rovnice:
x = fti)-f\q)^+r iq)yA-ŕ (i)4s+- (1°"7)
y = f\q)A-r\q)^f (q)~... ^ 6 120
Při původním odvození zobrazení byla stanovena podmínka, že osový poledník zůstane nezkreslený. Jelikož osový poledník má X = 0°, rovnice ( 10-7 ) a ( 10-8 ) pro bod P'0 ležící v zeměpisné šířce (p na tomto poledníku potom nabývají tvaru:
*o = f (a)
Jo=0
Z podmínky jeho nezkreslení vyplývá, že:
xo =SP = $ p
kde sp a Sp je délka oblouku osového poledníku od rovníku k bodu P'o na referenčním elipsoidu a v zobrazovací rovině. S uvážením výše uvedeného lze vyjádřit funkci f(q):
f(q) = Sp (1 f2
ae2 sin(pcos2 (p-aún(p{^-e2 sin2 (p) _ aún(p(\-e2) _ . (10-14)
L 2-2 V/2 =~Z 2 • 2 V/2 = ~M sin^
\l- e sin (p) \l- e sin q>)
Po dosazení do (10-13 ) se získá druhý člen obecné zobrazovací rovnice pro souřadnici x: f"(q) = -N cos (p sirup (10-15)
Obdobně se získávají další členy obecných zobrazovacích rovnic derivacemi rovnice ( 10-15 ). Jednotlivé členy mají tvary:
fm(q) = -Ncos3 (ól-58ŕ2+ŕ4 +27072 -33072ŕ2) (10-19) v nichž jsou použity symboly n = e cos^a t = \%q>.
V rovnicích ( 10-18 ) a ( 10-19 ) jsou zanedbány členy s rf a n6, protože jejich hodnoty jsou již zanedbatelné.
Po dosazení derivací funkce f(q) do obecných zobrazovacích rovnic ( 10-7 ) a ( 10-8 ) se získají základní tvary zobrazovacích rovnic Gaussovo zobrazení:
A2 i \ A4
x = Sp + jVcos^sin^ — + iVsin^cos3
která se dosadí do druhého členu rovnice (10-21). Po úpravě se získá rovnice:
y = Ncos(pA + ^—(l-t2 + n2)
ze které se vypočítá hodnota A ve druhém přiblížení:
Ä = ^---^(l-ŕ+n2)
N cos (p 6N cos (p
a znovu se vypočte A3 a A5. V uvedených mocninách stačí uvažovat hodnoty y pouze do páté mocniny. Tyto mocniny budou mít tvar:
A3 =--£3---/-^{\-t2+n2)
jV3cos> 2jV5cos>v ;
A>- y5
N5 cos5 cp
Hodnoty se opět dosadí do ( 10-21 ) a vypočte se hodnota A ve třetím, posledním přiblížení, přičemž se u páté mocniny y zanedbávají členy s rf a n4, jejichž hodnoty jsou zanedbatelné:
A = ---4-(l-t2 + n2)+-4-(5 " 2*2 9ŕ ) (10-24}
Ncoscp 6N3cos(py ' \2W5cos(px '
K tomu, aby bylo možné hodnotu A vypočítat, je nutné znát hodnotu zeměpisné šířky q>. Ta se vypočítá pomocí pomocné hodnoty q>f (viz. Obr. 10-6). Z bodu i5' se spustí kolmice k ose X, která tuto osu protne v bodě P'f. Protože se osový poledník délkově nezkresluje, je možné k hodnotě souřadnice x hned vypočítat (nebo nalézt v tabulkách) hodnotu zeměpisné šířky (pf (viz kapitola 1).
Poznámka: Pro výpočet hodnoty f- = sin[(pf - ((pf - cp)\
Hodnota ((pf - cp) je obecně ve srovnání s hodnotou (p/ diferenciálně malá, proto je možné pravou stranu rovnice rozvinout v Taylořovu řadu, přičemž se vzhledem k jejich velikosti uváží pouze první dva členy této řady:
sin cp = sin cpf - cos cpf (pf - (p)
Pokud se za (sf>f~ .r ' ''^V' Obdobně se upraví i zbylé goniometrické funkce. Po úpravě bude:
2
cos (o = cos (o f + sin (o f—--tf (10-33)
Vf ^2MfNf f
i. ^ v2
t = t
2
ÍX-t2f)^--tf (10-34)
f v J/2MfNf
Dosadí-li se výrazy ( 10-32 ), ( 10-33 ) a ( 10-34 ) do ( 10-30 ) a ( 10-24 ), získají se po úpravách konečné vztahy:
.2 ,.4
je konstantní, derivace budou pouze pro A a vzhledem k velikosti jednotlivých členů zobrazovacích rovnic stačí uvažovat výsledné členy nejvýše s A3. Derivace mají tvar:
dx í
— = 7Vcos^sin^/l + 7Vsin^cos3^(5-ř2 +9^2
dA
( 10-41)
127
Talhofer, V.: Základy matematické kartografie
^ = iVcos^ + iVcos>(l-ř2 + ^2)— (10-42) dA 2
Po dosazení výrazů (10-41) a (10-42 ) do (10-40 ) se získá vztah:
5 -12 + 9r/2 + Arf )— sin cpA + sin $>cos2 q>\5 -11 + 9r/2 + Arf )— 6 _ 6
tgy.
(\A? í \?i
l-t2+r/2)— \ + cos2 q)\l-t2+r/2j—
Po další úpravě se získá rovnice:
/ \A3 tgy = sirupA + sin cpcos2 (l + 3?/2 + 2n4)— + sin(pcos4 (p\2 -12)— ^10-44 ^
Pro praktické výpočty, pokud není požadována vysoká přesnost, je možné použít pouze první člen rovnice ( 10-44 ).
Meridiánovou konvergenci lze počítat i z rovinných pravoúhlých souřadnic X, Y . Pokud se dosadí do výrazu ( 10-44 ) vztahy ( 10-36 ) a ( 10-37 ), získá se po úpravě rovnice:
/ = — tf-^tÁl + t2-n2-2n4f)+^^tÁ2 + 5ň+3řf) (10-45)
Nf f 3N} /V f f fJ \5N5f /V f fJ
Uvedený vzorec zabezpečí přesnost výpočtu konvergence v prostoru České republiky (na úrovni zeměpisné rovnoběžky 50°) v jednom poledníkovém pásu 0,0005" .V případě, že je požadována přesnost výpočtu do 3", lze v rozmezí jednoho poledníkového pásu použít i zjednodušený vzorec
y = 2-tgq> (10"46)
N
Poznámka: Meridiánová konvergence v rámci jednoho poledníkového pásu nabývá jak kladných (na východní části pásu), tak i záporných hodnot (na západní části pásu). Při cp = 50° je její absolutní hodnota na okrajích pásu přibližně 2°18'.
128
Talhofer, V.: Základy matematické kartografie
10.4.1
Meridiánová konvergence v UTM
Výpočet meridiánové konvergence ze zeměpisných souřadnic je totožný jako v původním Gaussově zobrazení. Pokus se k výpočtům použijí rovinné pravoúhlé souřadnice N, E, je nutné uvážit i měřítkový faktor m0. Rovnice ( 10-45 ) potom bude mít tvar:
y.
E
E3
NeIfm0 f 3/V>03 A / //
E5
f> 15/V>05 /
^(2 + 54+3^) (10-47)
a obdobně zjednodušená rovnice (10-46 ) ze změní na:
E
ľ
NelmQ
tg(P
( 10-48)
10.5 Zákony zkreslení
V Gaussově zobrazení stačí vypočítat pouze délkové zkreslení m. Plošně zkreslení bude jeho kvadrátem a úhlově zkresleniJe zde nulové.
Délkové zkreslení je možné vypočítat ze zeměpisných nebo rovinných pravoúhlých souřadnic. Pokud se počítá ze zeměpisných souřadnic, lze použít obecné výrazy (viz kapitola Zákony zkreslení):
mr
VĚJ
yd(p) [d = 50° umožní přesnost výpočtu zkreslení na setiny milimetru. Pro tuto hodnotu zeměpisné šířky člen s A4 dosahuje hodnoty 4.10"8, tedy 0,04 mm.km"1. Proto pokud není požadována taková přesnost výpočtu, je možné tento člen zanedbat. Stejně tak, pokud se pro tuto zeměpisnou šířku zanedbá hodnota výrazu (1 + n2), vzniklá chyba nepřesáhne hodnotu 15.10"7 (tedy 1,5 mm.km"1). Proto se často v praxi používá pouze zjednodušený vzorec:
m = l + cos (p
2 %
( 10-50)
129
Talhofer, V.: Základy matematické kartografie
Délkové zkreslení je možné vypočítat i z rovinných pravoúhlých souřadnic. V tomto případě se volí jednodušší cesta, při níž se referenční elipsoid v daném bodě nahradí koulí o poloměru
r = 4mň~
a místo Gaussova zobrazení se uvažuje jednoduché válcové konformní zobrazení v příčné poloze (Mercatorovo) (viz Obr. 10-10).
Obr. 10-10 Náhrada referenčního elipsoidu koulí pro výpočet délkového zkreslení Gaussova zobrazení
S ohledem na tvar rovnice délkového zkreslení je možné psát:
1
m =--
cos S
a tento vzorec upravit rozvojem kosinu v řadu:
1
1--+--...
2 24
V rámci šestistupňového pásu není příliš velký rozdíl v délce zobrazené části oblouku původní rovnoběžky, části oblouku kartografického poledníku a souřadnice y (viz Obr. 10-11).
X RŠ \
Pj
Obr. 10-11 Ilustrace výpočtu délkového zkreslení Gaussova zobrazení
Proto lze psát:
y = rš
a tento výraz se dosadí do ( 10-51 ), přičemž stačí uvažovat pouze do mocniny/*. Výsledný vzorec bude:
130
Talhofer, V.: Základy matematické kartografie
m
y
y
( 10-52)
Jelikož výraz
y
2RL 24i?4
y
2R1 24R4
dosahuje velmi malých hodnot, lze vzorec dále upravit:
m = 1 +
y
y
( 10-53 )
2R1 24R4
Pokud se opět v rámci jednoho šestistupňového pásu zanedbá ve vzorci ( 10-53 ) třetí člen (sy4), maximální rozdíl na území ČR dosáhne hodnoty 5.10"8 (tedy 0,05 mm.km"1). Proto se v praxi často používá pouze zjednodušený výraz:
.2
m = 1 +
y
2R2
( 10-54)
Průběh délkového zkreslení v rámci jednoho šestistupňového pásu je zřejmý z obrázku (viz. Obr. 10-12).
Graf na následujícím obrázku (viz Obr. 10-13) znázorňuje závislost délkového zkreslení na zeměpisné šířce a hodnotě A. Další obrázek (viz Obr. 10-14) ilustruje průběh ekvideformát délkového zkreslení Gaussova zobrazení na území střední Evropy. Na území ČR dosahuje délkové zkreslení maximálních hodnot na okrajích 3. poledníkového pásu kolem 0,58 m.km"1.
F >' 1 \ 1 l\
/I 1 1 /III MM 1 1 ! 1 1 1 1 1* 1 1 !\ 1 1 l\ 1 1 1 \ Ml 1 M 11 r1111
1 1 M E 1 1 1 I =i 1111 1 M 1 1
1 1 ! 1 m>i I M I i i t I 1 i i i i m>i MM MM
t ! 1 I MM I MM I I M I
Délkové zkreslení Gaussova zobrazení
10 20 30 40 50 60 70 80 84
Obr. 10-12 Zobrazení ekvideformát délkového zkreslení Gaussova zobrazení (převzato z [23])
Obr. 10-13 Graf délkového zkreslení v Gaussově zobrazení
131
Talhofer, V.: Základy matematické kartografie
Obr. 10-14 Průběh ekvideformát délkového zkreslení Gaussova zobrazení na území střední Evropy
(převzato z [23])
10.5.1
Zákony zkreslení v UTM
V zobrazení UTM se délkové zkreslení ze zeměpisných souřadnic počítá podle vzorce:
\a2 í í~ . ?\^4
m = m,
1 + cos2 < Cr'l2>°'l2
1 + ^12 -#12
Cr'l2<°'l2
Pokud se zobrazí informace z tabulky graficky, je zřejmé, že se obraz geodetické čáry vždy konkávne zakřivuje k osovému poledníku daného pásu (Obr. 10-18).
135
Talhofer, V.: Základy matematické kartografie
10.6.1 Směrová korekce v zobrazení UTM
Směrová korekce v zobrazení UTM se počítá v podstatě podle stejných vzorců jako v Gaussově zobrazení, avšak s uvážením měřítkového faktoru m0. Rovnice postačující pro běžnou geodetickou praxi mají tvar:
Sl2=—^r(Nl-N2)(2El+E2)
6mLR
^i=-4^(^2-^i)(2£2+£1)
( 10-68)
( 10-69)
Pokud by byla nutná vyšší přesnost, je možné využít přesnější vztahy uvedené například v [16]. Rovněž zakřivení obrazu geodetické čáry je stejné jako na obrázku (viz Obr. 10-18).
10.7 Délková korekce geodetické čáry
Hodnoty délkového zkreslení se počítají podle jednoho ze vzorců ( 10-52 ) až ( 10-54 ) u Gaussova zobrazení, resp. ( 10-55 ) až ( 10-58 ) u zobrazení UTM. Pro praktické výpočty za použití rovinných pravoúhlých souřadnic lze výrazy pro výpočet délkové korekce dále upravovat.
Vyjádří-li se m pomocí těchto souřadnic, lze pro Gaussovo zobrazení psát: m,=l + yi
2R2
yi + y2 '2
i v - y i y\ +2 y\y 2 + y\
m i / = 1 + ---^— = l + +±-—^-
2R2 SR2
m1=\+ yi
2" 2R2
a po dosazení například do (4-64 ) a úpravě se získá rovnice:
Sl2=sl2+^{yl + yiy2 + yl) oo-™) 6R
kde hodnotu R lze počítat ke středu čáry a zaokrouhlovat ji na celé kilometry. Výsledná rovnice pro výpočet délkové korekce potom získá tvar:
SSii=Sl2-sl2=^{yl+yiy2+yl) (10-71)
Pokud se v zobrazení UTM použije obdobný postup, potom se vypočítá délka obrazu geodetické čáry podle upraveného vztahu ( 10-72 ) a z takto určené hodnoty délky obrazu geodetické čáry se vypočítá výsledná délková korekce.
Sl2=m0sl2+-^{E2+ElE2+E22) ( 10-72)
6m0R
136
Talhofer, V.: Základy matematické kartografie
Uvedené vztahy pro výpočet délkové korekce jsou použitelné pro běžné geodetické výpočty do vzdálenosti 20 km. Pro delší vzdálenosti nebo k dosažení vyšší přesnosti výpočtu je nutné použít přesnější vzorce uvedené například v [16] a [23].
10.8 Mezipásmové transformace
Protože v Gaussově zobrazení (i v jeho variantě UTM) má každý pás svoji souřadnicovou soustavu, je nutné v praxi poměrně často řešit transformaci souřadnic bodů ze souřadnicové soustavy jednoho pásu do souřadnicové soustavy druhého pásu. Tuto transformaci je možné řešit několika způsoby. V dřívějším období se nejčastěji používaly různé varianty rovinné transformace (viz. například kapitola Transformace zobrazení v [9]a [23]). Pro tyto transformace byly zpracovány i výpočetní tabulky.
K rovinným transformacím lze počítat i grafickou transformaci použitou na vojenských topografických mapách, kde v rámu těchto map (v tzv. překrytovém pásmu) jsou vykresleny rysky kilometrových čar sousedního poledníkového pásu. Jejich spojením se na mapu vykreslí celá kilometrová síť a pomocí ní je potom možné odečítat souřadnice i v souřadnicovém systému tohoto pásu.
V současné době je nejběžnější univerzální metoda transformace podle schématu: x\ý -+q>,A^>xn ,y", resp.
N1,E1 -^(p,Ä^Nn\En
Pro jednotlivé kroky se použijí vztahy ( 10-20 ), ( 10-21 ) a ( 10-35 ), ( 10-36 ) pro Gaussovo zobrazení a pro zobrazení UTM potom vztahy ( 10-23 ), ( 10-22 ) a ( 10-37 ),( 10-38 ).
Poznámka: Při transformaci jednoho bodu do souřadnicového systému jiného pásuje nutné uvážit rychlý nárůst délkového zkreslení, což v důsledku může ovlivnit i přesnost výpočtů v rovinných souřadnicích.
11. Křovákovo zobrazení
Po vzniku Československé republiky v roce 1918 byly budovány i nové geodetické a kartografické základy nového státního mapového díla, které se měly použít i pro katastrální účely. V roce 1922 navrhl Křovák (Josef Křovák 1884 až 1951) konformní kuželové zobrazení v obecné poloze jako součást geodetického referenčního systému jednotné trigonometrické sítě katastrální (S-JTSK). Toto zobrazení se využívalo do roku 1938 a bylo znovu zavedeno po druhé světové válce. S výjimkou padesátých a šedesátých let 20. století se používá dodnes. Zobrazení bylo definováno s ohledem na protáhlý a mírně vůči zeměpisným rovnoběžkám stočený tvar území bývalé Československé republiky (včetně tzv. Zakarpatské Ukrajiny) tak, aby minimalizovalo na tomto území délkové zkreslení. Dnes je používáno pouze v České a Slovenské republice.
V současné době jsou v tomto zobrazení vydávána státní mapová díla určená pro státní správu a samosprávu (viz Nařízení vlády ČR č.430/2006 - [17]). Jedná se zejména o Státní mapu v měřítku 1 : 5 000, Základní mapy ČR v měřítkách 1:10 000, 1 : 50 000, 1 : 100 000 nebo 1 : 200 000 a Mapu ČR v měřítku 1 : 500 000. V tomto zobrazení j sou také poskytována digitální data z databáze ZABAGED.
137
Talhofer, V.: Základy matematické kartografie
11.1 Základní charakteristiky zobrazení
V dále uvedených vzorcích jsou použité původní symboly, které zavedl Křovák. Zejména pro rovinné polární souřadnice se používají symboly R, D' namísto p, s a pro poloměr referenční koule r namísto původního R.
Křovákovo zobrazení je dvojité zobrazení, které je možné vyjádřit schématickým zápisem:
U,V-> Š,D ->R, D'->x,y
Výchozí referenční plochou je Besselův elipsoid, který je nejprve konformně zobrazen na referenční kouli. Na ní jsou definovány kartografické souřadnice, pomocí kterých je povrch koule transformován do zobrazovací roviny konformním kuželovým zobrazením. Poslední fází je transformace z polárních rovinných souřadnic na pravoúhlé. V následujících odstavcích jsou popsány jednotlivé fáze.
11.2 Zobrazovací rovnice
11.2.1
Zobrazení referenčního elipsoidu na referenční kouli
V první fázi je Besselův elipsoid konformně zobrazen na referenční kouli s jednou nezkreslenou rovnoběžkou (po = 49°30', která probíhá přibližně středem území původní Československé republiky. Poloměr referenční kouleje
r = jM0N0
K transformaci souřadnic jsou použity rovnice odvozené v kapitole 5:
SI f+ 45°
tg
a| ^- + 45°
f l-eúrup^ 2 l + eúxup
V = aA
Konstanty zobrazení jsou:
r = 6380 703,6105 m k= 1,00341 91640 cc = 1,00059 7498372
Po transformaci odpovídá původní hodnotě (po nezkreslené rovnoběžky hodnota na referenční kouli Uo = 49°27'35",84625.
11.2.2 Transformace zeměpisných souřadnic na referenční kouli na kartografické souřadnice
Na referenční kouli je definována souřadnicová soustava kartografických souřadnic Š, D. Tato soustava vyhovuje protáhlému a mírně stočenému tvaru původní republiky. Osu území tvoří základní kartografická rovnoběžka So, z jejíhož tvaru byla vypočítána poloha kartografického pólu K podle postupu uvedeném v odstavci 1.2.2.a . Na této kartografické rovnoběžce byl za nejvýchodnějším cípem republiky, který tvořil okraj tehdejší speciální mapy 1:75 000, zvolen bod A, jehož zeměpisné souřadnice jsou:
138
Talhofer, V.: Základy matematické kartografie
01 = 48° 15'
Aa = 42°30' východně Ferra (24°50' východně Greenwich).
Tento bod má na referenční kuli souřadnice:
f/y4 = 48°12'42",69689 J^ = 42°31'31",41725
Z polohy základní kartografické rovnoběžky byla vypočítána poloha kartografického pólu K. Pól leží na stejném poledníku jako bod A je od něho na sever o 11°30\ Jeho zeměpisné souřadnice na kouli jsou:
£4 = 59°42'42",69689 ř/ff = 42°31'31",41725
Základní kartografická rovnoběžka má hodnotu Šo = 78°30'. Celé území bývalého Československa leželo potom v úzkém pásu vymezeném dvěma kartografickými rovnoběžkami v relativně malé vzdálenosti AS = 2°3ľ, což je asi 280 km. Uvedené hlavní prvky zobrazení dokumentuje obrázek (viz. Obr. 11-1).
Zeměpisné souřadnice U, V jsou transformovány na kartografické souřadnice Š, D pomocí rovnic ( 1-24 ) a ( 1-25 ). Křovák jej pouze upravil zavedením zenitové vzdálenosti a (a = 90°- UK) kartografického pólu (viz Obr. 11-2).
Obr. 11-1 Základní prvky Křovákova zobrazení
139
Talhofer, V.: Základy matematické kartografie
Ps
a=90°-Ut
90'-U
í 180°-D K
I '°
Šo í /
A
Obr. 11-2 Transformace zeměpisných souřadnic na kartografické
Upravené rovnice potom budou:
sin Š = smll cos a + cos U sin a cos(F — Vk)
. cosíl .
smD =-¥sm(V-Vk)
cosi
(11-1)
(11-2)
11.2.3
Transformace do zobrazovací roviny
Pro zobrazení referenční kouleje použito jednoduché konformní kuželové zobrazení s jednou nezkreslenou rovnoběžkou, která se z důvodů zmenšení absolutní hodnoty zkreslení dodatečně zkresluje pomocí měřítkového faktoru m0 = 0,9999. Pro výpočty se používají vztahy ( 7-31 ), ( 7-5 ), ( 7-34 ) a ( 7-15 ). Protože se však jedná o obecnou polohu zobrazení, mají zobrazovací rovnice a další vztahy následující tvary:
R = R0
tg
^ + 45 2
j
tg
2
D' = nD
+ 45°
(11-3)
kde:
Rq = m0rcotgŠ0 n = sin Šn
(11-4)
(11-5) (11-6)
Hodnota R0 = 1 298 039,0046 m a n = 0,97992 47046.
Poznámka: Použitím měřítkového se zobrazení mění v zobrazení se dvěmi nezkreslenými kartografickými rovnoběžkami o hodnotách & = 79°18'03" aŠ2 = 77°40'50".
140
Talhofer, V.: Základy matematické kartografie
11.2.4
Převod rovinných polárních souřadnic na pravoúhlé
Polární souřadnice R, D' jsou transformovány na rovinné pravoúhlé x, y v souřadné soustavě, kde osa X]e umístěna v obraze poledníku Ak a její počátek je v obraze kartografického pólu K. Kladný směr osy je na jih. Osa Y je na ní kolmá a její kladná orientace je na západ. Polární souřadnice j sou transformovány podle vzorců:
x = RcosD' (H-7)
y = RsinD'
Celé území republiky potom leží v prvním kvadrantu (viz Obr. 11-3).
(11-8)
Obr. 11-3 Poloha rovinného pravoúhlého systému v Křovákově zobrazení
Obrazem kartografických poledníků jsou polopřímky vycházející z obrazu kartografického pólu, obrazem kartografických rovnoběžek jsou soustředné kružnice se středem opět v obraze kartografického pólu. Obrazem zeměpisných poledníků a rovnoběžek jsou složité křivky, které však na zobrazovaném území České a Slovenské republiky mohou být na mapách středních měřítek nahrazeny přímkami (poledníky) nebo soustřednými kružnicemi (rovnoběžky), jejichž zakřivení je téměř totožné se zakřivením obrazu rovnoběžek u Gaussova zobrazení.
11.3 Inverzní funkce k zobrazovacím rovnicím
Inverzní funkce k zobrazovacím rovnicím se řeší postupně podle schématu (Kratochvíl: Polohové geodetické sítě, 2000):
x, y —> R, D' —> Š, D —> U, V-> (p, A
Nejprve se ze vztahů (11-7 ) a (11-8 ) vypočítají polární souřadnice:
(11-9)
R = i
x2+y2
D' = arctg
y
(11-10)
141
Talhofer, V.: Základy matematické kartografie
Poté se s využitím vztahů ( 7-31 ) a ( 7-5 ) vypočítají kartografické souřadnice na referenční kouli:
Š = 2 jarctan
f č
tg
So
+ 45°
j
-45c
D
D' sin Šn
(11-11)
(11-12)
Transformaci kartografických souřadnic na zeměpisné lze řešit pomocí vztahů ze sférické trigonometrie (viz. Obr. 11-2):
U = arcsin (cos a sin Š - sin a cos Š cos D)
f
V = Vh - arcsin
cos Š cosll
sin D
(11-13) (11-14)
Výpočet zeměpisné šířky je nutné provést v několika iteracích, protože argument m je na obou stranách rovnice. Jeden z možných postupuje následující:
. Zeměpisná délka se potom určí podle schématu:
A = A\ pokud jmenovatel ve zlomku je (12-10)
kladný
(12-11) ( 12-12)
á = x + A', pokud A£>0 a/l'<0
A = -7i + A\ pokud AE <0 a A' > 0
V případě, že se pro výpočet použije funkce arctg se dvěma argumenty (čitatel i jmenovatel vstupují do výpočtu samostatně), potom A je přímo rovna A' v rozsahu <-n, it>.
Zeměpisná šířka q> se počítá postupně s výjimkou hodnoty AN = 0, kdy q> = 90°. Nejprve se vypočítá p podle j ednoho za vztahů:
p = |A£|, jestliže AjV=0 (12-13)
p = \AN\, jestliže AE=0 (12-14)
P
AE
ve všech ostatních případech
(12-15)
sin/t
Dále se vypočítá hodnota z a izometrické šířky q:
z = 2arctg—^— (12-16)
m0C0
q = --z (12-17)
2
Výsledná zeměpisná šířka se potom vypočítá podle vztahu:
l + eúxup
( 12-23 )
Konstanty zobrazení jsou počítány pro každou vrstvu samostatně, a to jak pro vrstvy s intervalem 8°, tak pro vrstvy s intervalem 4°. Pro vrstvy, ve kterých leží Česká republika, parametry a konstanty zobrazení mají následující hodnoty:
Tabulka 12-2 Hodnoty základních parametrů pro Lambertovo konformní kuželové zobrazení pro území ČR
Interval vrstvy 0
y
cos (p0 cos A cos (p0 sin A
Hodnoty q>, HeiaNse dále počítají pro i-tou iteraci následovně:
q>{ = arctg N
z ^-1+^-1
a
(13-4)
■e2 sin2 q>{
H.
cos cpi cos A
y
cos sin A
Při prvním výpočtu zeměpisné šířky je možné zanedbat elipsoidickou výšku bodu P. Pomocí (po se vypočítají první aproximace iVo a Hei0. Jejich dosazením do rovnic ( 13-4 ) se získají zpřesněné hodnoty q>,, N, a Hei,. Iterační výpočet se ukončí, pokud je rozdíl mezi předcházející a počítanou hodnotou menší než požadovaná přesnost výpočtu.
Obr. 13-1 Prosotorové pravoúhlé souřadnice
Obdobné vztahy jako ( 13-2 ) a ( 13-3 ) platí i pro transformaci zeměpisných souřadnic U, V a výšky H na referenční kouli o poloměru R na prostorové pravoúhlé souřadnice. Výšku HJe možné uvažovat pouze jako nadmořskou, protože tato transformace se používá především pro méně přesné úlohy. Pokud se uvažuje stejná poloha a orientace souřadnicových os jako v případě referenčního elipsoidu, lze transformaci vyjádřit vzorci:
152
Talhofer, V.: Základy matematické kartografie
X = (R + H)cos U COS V y = (R + H)cosUsmV z = (R + H)smU
Zpětný převod je potom možný podle vztahů:
13.1.2
V = arctg
(f)
f
U = arctg
H
x2+/
(13-5)
(13-6)
R
1
4
x2+y2
R
sin U cos U Tříprvková prostorová transformace
Nejj ednodušší transformací mezi referenčními systémy je tříprvková prostorová transformace. Rozdíl mezi původním a novým referenčním systémem geocentrických souřadnic je pouze v lineárním posunu počátků obou systémů. Počátek nového souřadnicového systému je pouze vzhledem k původnímu systému posunut o hodnoty dx, dy a dz. Vlastní transformaci lze vyjádřit vztahem ( 13-7 ):
Xn dx X
yn = dy + y
dz z
(13-7)
Graficky lze podstatu transformace vyjádřit následujícím obrázkem (Obr. 13-2):
Zn
Původní souřadnicový systém
On
Nový souřadnicový systém
Yn
Xn
dz
i-' dx
dy
Obr. 13-2 Tříprvková prostorová transformace Hodnoty souřadnic i lineárních posunů se vyjadřují v metrech.
153
Talhofer, V.: Základy matematické kartografie
13.1.3 Sedmiprvková prostorová transformace
Přesnější a komplexnější transformace využívající prostorové pravoúhlé souřadnice je sedmiprvková prostorová transformace, někdy nazývaná i jako prostorová podobnostní transformace. Vedle lineárních posunuje zde uvažováno i se třemi rotacemi kolem původních os (rx, ry a rz) a se změnou měřítka - měřítkovým faktorem m=l+ju.
Transformaci je možné vyjádřit následující rovnicí ( 13-8 ):
Xn dx m 0 0" 1 X
yn = dy + 0 m 0 1 y (13-8)
Zn dz 0 0 m ry 1 z
Graficky lze podstatu transformace opět vyjádřit následujícím obrázkem (Obr. 13-3):
Obr. 13-3 Sedmiprvková prostorová transformace
Hodnoty souřadnic a lineárních posunů se opět uvádějí v metrech, hodnoty rotací v desetinách vteřin a měřítkový faktor bývá v řádech 10"6 až 10"5. V anglické literatuře bývá uváděn v jednotkáchppm (partsper milion).
Hodnoty rotací jsou definovány dvojím způsobem. Pokud se díváme směrem k původním osám X, Y, Z, je rotace jsou buďto ve směru nebo proti směru hodinových ručiček. Přitom je možné přiřadit kladné nebo záporné znaménko oběma směrům. Například v Austrálii byla pro rotace při definici referenčního souřadnicového systému použita kladná znaménka pro směry otáčení hodinových ručiček, v Evropě tomu bylo naopak [10]. Před použitím popsané transformace je nutné zjistit, jaký systém rotací byl použit. Pokud by byl použit nesprávný, výsledné transformované hodnoty jsou chybné.
13.1.4 Moloděnského transformace
Moloděnského transformace umožňuje přímou transformaci zeměpisných souřadnic definovaných v souřadnicových systémech, aniž je nutný jejich převod do prostorových pravoúhlých souřadnic. K této transformaci je nutná znalost parametrů původního elipsoidu
154
Talhofer, V.: Základy matematické kartografie
(velikost poloos a, b), lineárních posunů dx, dy a dz a rozdílů parametrů použitých referenčních elipsoidů (původního a nového) - velké poloosy Aa a zploštění Af.
Transformace je dána vztahy ( 13-9 ):
A(p-
, e smcpcoscp . sin ^>cos Adx - sin