155KAR1 Kartografie 1

Z GeoWikiCZ

Anotace

Význam matematické kartografie. Referenční plochy a souřadnicové soustavy. Kartografická zkreslení. Klasifikace kartografických zobrazení. Zobrazení elipsoidu na kouli. Jednoduchá zobrazení kuželová, válcová a azimutální. Nepravá, polykónická, polyedrická a obecná zobrazení. Přehled zobrazení užitých na území ČR a ve světě. Volba, identifikace a hodnocení zobrazení. Referenční souřadnicové systémy v GIS.

Harmonogram přednášek

Přednáší Doc. Ing. Jiří Cajthaml, Ph.D.

pondělí 10:00 C208

  1. Význam matematické kartografie. Referenční plochy. Souřadnicové soustavy. Důležité křivky. (PDF)
  2. Kartografická zobrazení. Délkové zkreslení. Podmínky konformity. Zkreslení azimutu a úhlu. Hlavní paprsky. Plošné zkreslení. Zkreslení při známých hlavních paprscích. (PDF)
  3. Klasifikace zobrazení. Jednoduchá zobrazení. Zobrazení elipsoidu na kouli. (PDF)
  4. Kuželová zobrazení ekvidistantní, ekvivalentní, konformní. (PDF)
  5. Válcová zobrazení ekvidistantní, ekvivalentní, konformní. Válcové projekce. Azimutální zobrazení. Azimutální projekce. (PDF)
  6. Zobrazení použitá na území ČR. Cassiniho zobrazení. Křovákovo zobrazení. (PDF)
  7. Zobrazení použitá na území ČR. Gaussovo zobrazení elipsoidu v poledníkových pásech. Zobrazení UTM. (PDF)
  8. Nepravá zobrazení kuželová, azimutální a válcová. Jejich modifikace. (PDF)
  9. Polykónická zobrazení. Neklasifikovaná zobrazení. Polyedrické aplikace. (PDF)
  10. Volba, užití a identifikace zobrazení. Hodnocení zobrazení. Kritéria pro optimalizaci zobrazení. (PDF)
  11. Referenční souřadnicové systémy v GIS. Systém JTSK/05. Přehled v praxi užívaných zobrazení. Transformace rovinných souřadnic. (PDF)



Literatura :

  • Buchar P.: Matematická kartografie (skriptum)
  • Fiala F.: Matematická kartografie (učebnice)
  • Hojovec a kol.: Kartografie (učebnice)

Program cvičení

datum náplň cvičení
24.9. zadání 1. úlohy
1.10. zadání 2. úlohy
8.10. odpadá
15.10. konzultace
22.10. konzultace
29.10. odpadá
5.11. zadání 3. úlohy, test úloh 1+2
12.11. konzultace
19.11. konzultace
26.11. zadání 4. úlohy
3.12. konzultace
10.12. konzultace
17.12. test úloh 3+4


Každá úloha bude obsahovat :

  • zadání
  • matematickou formulaci
  • výpočty včetně potřebných mezivýpočtů
  • zdrojové kódy skriptů
  • obrazové přílohy
  • přehledný souhrn výsledků v tabulkách
  • zhodnocení výsledků a závěr

Úlohy nesprávně numericky vyřešené či nedbale vyhotovené nebudou přijaty. Úlohy 1+2, 3+4 se po ukončení testují. Pro udělení zápočtu je nezbytně nutné odevzdání všech úloh a úspěšné vykonání požadovaných testů.

Literatura je uvedená v harmonogramu přednášek.

Zadání cvičení

ÚLOHA 1 - Loxodroma a ortodroma

Vypočítejte zeměpisné souřadnice mezilehlých bodů loxodromy a ortodromy dané koncovými body a na referenční kouli. Jako interval mezilehlých bodů použijte zeměpisné délky.

Proveďte zákres obou křivek v následujících kartografických zobrazeních. Jako podklad použijte síť poledníků a rovnoběžek vygenerovanou skriptem v prostředí MATLAB (ukázka).

* Válcové konformní zobrazení s jednou nezkreslenou rovnoběžkou (Mercatorovo)


* Azimutální gnómonická projekce


* Válcová stereografická projekce

* Mercator-Sansonovo nepravé válcové sinusoidální zobrazení


Dále vypočtěte délky obou křivek. Výpočty provádějte na referenční kouli s poloměrem s přesností na úhlové vteřiny a .

Číselné zadání:

číslo
1 10 10 50 80
2 10 10 50 90
3 10 10 50 100
4 10 10 50 110
5 10 10 50 120
6 10 10 50 130
7 10 10 60 80
8 10 10 60 90
9 10 10 60 100
10 10 10 60 110
11 10 10 60 120
12 10 10 60 130
13 10 10 60 140
14 10 10 70 80
15 10 10 70 90
16 10 10 70 100
17 10 10 70 110
18 10 10 70 120
19 10 10 70 130
20 10 10 70 140
21 10 20 50 80

ÚLOHA 2 - Výpočet zkreslení

Je dáno zobrazení koule o poloměru R=6380 km do roviny mapy zobrazovacímí rovnicemi:

Určete souřadnice obrazu daného bodu geografické sítě. Zjistěte vlastnosti daného zobrazení. V daném bodě vypočtěte veškeré charakteristické hodnoty:

délkové zkreslení v poledníku ,

délkové zkreslení v rovnoběžce ,

úhel mezi obrazem poledníku a rovnoběžky ,

azimuty extrémního délkového zkreslení v originále , ,

azimuty extrémního délkového zkreslení v obraze , ,

hodnoty extrémních délkových zkreslení (poloosy Tissotovy indikatrix) , ,

maximální úhlové zkreslení ,

plošné zkreslení .

Vypočtené hodnoty zakreslete do Tissotovy indikatrix, která bude zobrazena ve vhodném měřítku nad geografickou sítí celého světa v daném zobrazení (interval 10 stupňů). Požadovaná přesnost výpočtu je u měřítek 6 desetinných míst, a pro úhly 1 úhlová vteřina. Číselné výsledky porovnejte s výsledky určenými v programu PROJ.4 (parametry)

Číselné zadání:

číslo
1 10 70
2 10 80
3 10 90
4 10 100
5 10 110
6 10 120
7 10 130
8 10 140
9 10 150
10 10 160
11 20 70
12 20 80
13 20 90
14 20 100
15 20 110
16 20 120
17 20 130
18 20 140
19 20 150
20 20 160
21 30 70

ÚLOHA 3 - Srovnání zobrazení

Pro zadané státní území vzájemně porovnejte z hlediska hodnot délkového zkreslení:

  • válcové konformní zobrazení o dvou nezkreslených rovnoběžkách (Mercatorovo)
  • stereografickou azimutální projekci s nezkreslenou rovnoběžkou

Navrhujte obecnou polohu zvolených zobrazení. Nezkreslené rovnoběžky volte tak, aby vliv délkového zkreslení ve středu území a na okraji byl v absolutní hodnotě stejný.

Vyhotovte obrázky se zákresem státních hranic a zeměpisné sítě s intervalem 1 stupně.

Rozhodněte, kolik bude zapotřebí souřadnicových soustav za podmínky, že vliv délkového zkreslení nepřekročí 20 cm/km.

V uvedených řešeních uvažujte referenční kouli. Výpočet zkreslení proveďte s přesností 6 desetinných míst, úhlové výpočty na minuty.

Shapefile se zeměpisnou sítí po 1°

Geodatabáze hranic států

Doporučený postup v ArcGIS pro Mercatorovo zobrazení:

  1. otevřít shapefile s hranicemi států, vybrat stát a exportovat do samostatného souboru
  2. nastavit souřadnicový systém okna ArcGIS na Mercator
    1. Data Frame Properties - Coordinate system - Predefined - Projected CS -- Mercator (world)
  3. vytvořit minimální ohraničující obdélník s nejmenší šířkou
    1. Data Management Tools - Features - Minimum Bounding Geometry
    2. nastavit vstup, výstup, RECTANGLE BY WIDTH
    3. nastavit okno aplikace, aby tam byl celý stát a v Environments nastavit u Output Coordinates, u Processing Extent a u Cartography souřadnicový systém jako "Same as Display"
    4. spustit funkci
  4. získat souřadnice bodů uprostřed kratších stran (krajní body ortodromy, která prochází středem státu)
    1. Feature To Line
    2. Split Line at Vertices
    3. spustit Editor, Split (Percentage 50)
    4. vytvořit novou bodovou vrstvu, nastavit Snapping, vytvořit 2 body
    5. v atributové tabulce přidat sloupce fi, lambda a vypočítat souřadnice pomocí Calculate geometry nad sloupcem
  5. odečíst souřadnice bodu na jedné z krajních rovnoběžek (roh obdélníku)
  6. pomocí sférické trigonometrie vypočítat Š tohoto bodu
  7. vypočítat délkové zkreslení na krajní rovnoběžce (závisí pouze na Š)
  8. vypočítat kolik je potřeba souřadnicových soustav, aby zkreslení nepřesáhlo požadovanou hodnotu
  9. vytvořit mapu
    1. nastavit souřadnicový systém okna
    2. Data Frame Properties - Coordinate system - New - Projected CS
    3. do Geographic CS nastavit WGS84
    4. Projection: Hotine Oblique Mercator Two Point Center
    5. nastavit dva krajní body hlavní ortodromy (šířka a délka) a dále šířku bodu zhruba uprostřed území
    6. v layout view exportovat mapu s hranicemi státu, zeměpisnou sítí, minimálním obdélníkem

Konkrétní zadání:

číslo stát
1 NORSKO
2 CHILE
3 FINSKO
4 ŠVÉDSKO
5 VIETNAM
6 LITVA
7 LOTYŠSKO
8 ESTONSKO
9 ŘECKO
10 THAJSKO
11 LAOS
12 KAMBODŽA
13 DÁNSKO
14 EKVÁDOR
15 BOLÍVIE
16 NOVÝ ZÉLAND
17 PORTUGALSKO
18 BULHARSKO
19 PARAGUAY
20 TUNISKO
21 ŠPANĚLSKO

ÚLOHA 4 - Výpočet souřadnic bodů v kartografických zobrazeních

Je dán bod P svými zeměpisnými souřadnicemi na elipsoidu WGS84. Vypočtěte pravoúhlé souřadnice jeho rovinného obrazu a uveďte hodnotu délkového zkreslení v systému JTSK. Použijte vlastní skript v Matlabu.

Pomocí software PROJ.4 vypočtěte rovinné pravoúhlé souřadnice bodu P a hodnotu délkového zkreslení:

  • v systému JTSK
  • v systému UTM (33. šestistupňový pás, počátek posunut od průsečíku základního poledníku a rovníku o 500 km západně)
  • v systému S-42 (3. šestistupňový pás, počátek posunut od průsečíku základního poledníku a rovníku o 500 km západně)

Použití PROJ.4:

  • (cs2cs +proj=latlong +datum=WGS84 +to +proj=krovak +ellps=bessel +towgs84=570.8,85.7,462.8,4.998,1.587,5.261,3.56)
  • (proj +proj=utm +datum=WGS84 +zone=33)
  • (cs2cs +proj=latlong +datum=WGS84 +to +proj=tmerc +ellps=krass +lon_0=15 +x_0=500000 +towgs84=26.0,-121.0,-78.0,0,0,0,0)


Pomocí software ArcGIS vypočtěte rovinné pravoúhlé souřadnice bodu P:

  • v systému JTSK (transformace S_JTSK_To_WGS_1984_1)
  • v systému UTM (33. šestistupňový pás, počátek posunut od průsečíku základního poledníku a rovníku o 500 km západně)
  • v systému S-42 (3. pás šestistupňový pás, počátek posunut od průsečíku základního poledníku a rovníku o 500 km západně)(transformace Pulkovo_1942_To_WGS_1984_5)


Pro kartografická zobrazení používejte příslušné geodetické datumy (elipsoidy). Parametry transformace mezi nimi:

  • Bessel2WGS: 570.8,85.7,462.8,4.998,1.587,5.261,3.56
  • Krasovsky2WGS: 26.0,-121.0,-78.0,0,0,0,0

Požadovaná přesnost výsledků : (souřadnice na cm), (zkreslení na 6 desetinných míst).

Číselné zadání :

číslo
1 50°10' 15°10'
2 50°10' 15°20'
3 50°10' 15°30'
4 50°10' 15°40'
5 50°10' 15°50'
6 50°20' 15°10'
7 50°20' 15°20'
8 50°20' 15°30'
9 50°20' 15°40'
10 50°20' 15°50'
11 50°30' 15°10'
12 50°30' 15°20'
13 50°30' 15°30'
14 50°30' 15°40'
15 50°30' 15°50'
16 50°40' 15°10'
17 50°40' 15°20'
18 50°40' 15°30'
19 50°40' 15°40'
20 50°40' 15°50'
21 50°50' 15°10'


Hodnocení odevzdaných úloh

úloha č. 1

jméno číslo zadání stav komentář
Brhlíková 1
ok
stejný závěr jako Nápravník
Formánek 2
opraveno, ok
odevzdáno pozdě, bez obálky, nepopsané obrázky a křivky, špatně Mercator-Sansonovo z., délky nejsou na mm, úhly nejsou na úhlové vteřiny
Hrdina 3
opraveno, ok
délky nejsou na mm
Lauwereys 4
opraveno, ok
odevzdáno pozdě, špatně zaokrouhlování 60 vteřin
Masák 5
ok
odevzdáno pozdě
Nápravník 6
opraveno, ok
špatně minuty a vteřiny, špatně délka ortodromy, stejný závěr jako Brhlíková
Pudil 7
ok
skript stejný jako Nedbalová
Tomášková 8
ok
Zbíral 9
opraveno, stále chybí délky křivek, po opravě stále nejsou délky na mm, již ok
chybí celá technická zpráva, výsledky nejsou na úhlové vteřiny
Brezničanová 10
opraveno, ok
délky nejsou na mm, špatně obrázek Mercatora, chybí popis křivek v obrázcích
Beliančinová 11
opraveno, délky stále nejsou na mm
délky nejsou na mm, špatná kvalita obrázků (jpg)
Děkanová 12
ok
Kouba 13
ok
délka loxodromy 10x větší
Merdymshaeva 14
ok
Nedbalová 15
opraveno, ok
špatně průsečíky loxodromy, špatně zaokrouhlování 60 vteřin, skript stejný jako Pudil
Novotný 16
opraveno, ale stále špatná kvalita obrázků
prohozené křivky v obrázcích, špatná kvalita (jpg), Mercator není celá Země
Talapková 17
opraveno ok
délky nejsou na mm
Usik 18
opraveno, ok
špatně průsečíky a délka ortodromy, špatná kvalita obrázků (jpg)
Vynikal 19
opraveno, ok
špatně Mercator-Sansonovo z.
Zeman 20
ok
Zíma 21
opraveno, ještě opravit délky, ok
průsečíky loxodromy mírně jinak, délky nejsou na mm, používán poloměr 6378km

úloha č. 2

jméno číslo zadání stav komentář
Brhlíková 1
opraveno, ok
odevzdáno pozdě, chybí skript pro zakreslení elipsy
Formánek 2
opraveno, ok
odevzdáno pozdě, souřadnice nejsou na cm, zkreslení nejsou na 6 desetinných míst, chybí jednotky u souřadnic, chybí skripty
Hrdina 3
ok
Lauwereys 4
opraveno, ok
úhly převést z rad na stupně
Masák 5
ok
odevzdáno pozdě
Nápravník 6
opraveno, ok
odevzdáno pozdě, elipsa ve špatné poloze, úhly nejsou na vteřiny, špatné jednotky u a,b
Pudil 7
ok
Tomášková 8
ok
odevzdáno pozdě
Zbíral 9
opraveno, ok
špatně vypočítané Ae2 v obraze
Brezničanová 10
opraveno, ok
velmi špatná kvalita obrázku (nedělejte ho v JPG, nebo použijte lepší kvalitu)
Beliančinová 11
ok
Děkanová 12
opraveno, ok
špatně azimuty extrémních zkreslení v obraze, prohozené značky pro minuty a vteřiny (' ")
Kouba 13
opraveno, ok
chybí výsledky X,Y bodu
Merdymshaeva 14
opraveno, stále nejsou souřadnice na mm, již ok
souřadnice nejsou na mm, špatné jednotky u a,b, úhly nejsou na stupně-minuty-vteřiny
Nedbalová 15
opraveno, ok
chybí výsledky X,Y bodu
Novotný 16
ok
Talapková 17
ok
Usik 18
opraveno, souřadnice stále nejsou na mm
souřadnice nejsou na mm
Vynikal 19
ok
Zeman 20
ok
Zíma 21
ok

úloha č. 3

jméno číslo zadání stav komentář
Brhlíková 1
ok
odevzdáno pozdě
Formánek 2
opraveno, ok
pokud je měřítko 0.9984, znamenalo by to 160cm/km a nevyhovovalo by zadání
Hrdina 3
ok
Lauwereys 4
ok
Masák 5
NEODEVZDÁNO
Nápravník 6
opraveno, ok
opsáno, nesedí závěr výpočet, špatně počet pásů u Mercatora
Pudil 7
ok
Tomášková 8
ok
odevzdáno pozdě
Zbíral 9
opraveno, ok
chybí závěr kolik souřadnicových systémů je potřeba v obou zobrazeních; nejdená se o gnónomickou ale stereografickou projekci
Brezničanová 10
opraveno, ok
špatně počty soustav, u válce stačí 3 pásy, u azimutálního zobrazení chybí obrázek se zákresem kružnic, protože to nelze jen tak spočítat
Beliančinová 11
ok
Děkanová 12
opraveno, ok
špatně počet soustav, u válce stačí 2 pásy (pás je celá soustava, ne jen její polovina), u kružnic chybí výpočet poloměru - soustav stačí určitě méně než 5 (možná průměr x poloměr)
Kouba 13
ok
Merdymshaeva 14
ok
Nedbalová 15
ok
Novotný 16
opraveno, ok
odevzdáno pozdě, chybí obrázek s kružnicemi, které splňují podmínku zkreslení
Talapková 17
opraveno, ok
špatně počet soustav u azimutálního zobrazení, chybí přehledně mezivýsledky a tak není možné zkontrolovat, kde je chyba
Usik 18
ok
Vynikal 19
ok
Zeman 20
ok
Zíma 21
ok

úloha č. 4

jméno číslo zadání stav komentář
Brhlíková 1
opraveno, ok
špatně výsledky z ArcGIS
Formánek 2
opraveno, ok
úplně chybí zpracování v Matlab, opsaná 2x hohnota zkreslení UTM a S42
Hrdina 3
opraveno, ok
špatně délková zkreslení u S-42 a UTM
Lauwereys 4
opraveno, ok
chybí některé výsledky (mUTM, třetí souřadnice JTSK)
Masák 5
NEODEVZDÁNO
Nápravník 6
ok
vymyšlené hodnoty zkreslení z ArcMap
Pudil 7
ok
Tomášková 8
ok
odevzdáno pozdě
Zbíral 9
ok
Brezničanová 10
ok
Beliančinová 11
opraveno, ok
špatně délkové zkreslení u S-42
Děkanová 12
opraveno, ok
prohozené výsledky UTM a S42, špatně délkové zkreslení u S-42, překlep o 100000 u Y v S42
Kouba 13
ok
Merdymshaeva 14
ok
vymyšlené hodnoty zkreslení z ArcMap
Nedbalová 15
ok
Novotný 16
ok
Talapková 17
ok
Usik 18
opraveno, ok
špatně délkové zkreslení u UTM
Vynikal 19
ok
špatně opsané hodnoty zkreslení z obrazovek do tabulek
Zeman 20
ok
Zíma 21
opraveno, ok
chybí přehledná tabulka výsledků, chybí délkové zkreslení z vlastního skriptu u JTSK, špatně výsledky S42 (PROJ i ArcGIS)