155GIS2 / 2. cvičení: Porovnání verzí
Řádek 62: | Řádek 62: | ||
{| | {| | ||
|{{fig|grass-gis2-cv4|Příklad určení výměry území s nadmořskou výšku mezi 500 a 700 metry}} | |{{fig|grass-gis2-cv4|Příklad určení výměry území s nadmořskou výšku mezi 500 a 700 metry}} | ||
|{{fig|grass-gis2- | |{{fig|grass-gis2-cv2|Reklasifikace rastrové mapy sklonu svahu v systému GRASS|size=700}} | ||
|} | |} | ||
</center> | </center> |
Verze z 17. 2. 2015, 15:32
Základy mapové algebry, reklasifikace rastrových dat
Základní pojmy
- Mapová algebra
- Mapová algebra
- Raster Operations
- Rastrový datový model a mapová algebra
- Velikostní stupnice (poměrová - ratio, intervalová - interval, pořadová, kategoriální)
- Mapová algebra
- Lokální funkce (Cell Statistics)
- jedna rastrová vrstva (goniometrické funkce, reklasifikace, ...)
- více rastrových vrstev (min, max, mean, sum, product, majority, minority, variety, ...)
- Fokální funkce (Neighborhood Statistics)
- Fokalní min, max, mean, sum,...
- Kernelové operace (density, slope, aspect)
- Inkrementální funkce (globalní funkce)
- Např. vzdálenostní funkce, obalová zóna, cost distance, cost path
- Zonální funkce (Zonal Statistics)
- Zonalní sum, mean, min, max, ..
- Lokální funkce (Cell Statistics)
- Reklasifikace
Poznámky k systému ArcGIS
- Mapová algebra
- Zapnout nadstavbu - Customize • Extensions • Spatial Analyst
- Základní mapová algebra (Math)
- Pokročilá mapová algebra, rastrový kalkulátor (Map Algebra • Raster Calculator)
- Spatial Analyst • Map Algebra • Raster Calculator
- Spatial Analyst • Extraction • Exctract By Mask
- Spatial Analyst • Raster Creation
- Spatial Analyst • Zonal • Zonal Geometry
- 3D Analyst • Functional Surface • Add Surface Information
Odkazy
- Mapová algebra
- 3D Analyst. Vizualizace dat v ArcScene. Mapová algebra
- Zonal Geometry
- How Cell Statistics works
- How Focal Statistics works
- How Zonal Statistics works
Poznámky k systému GRASS
Návod na řešení úloh z tohoto cvičení zde.
Úlohy
V následujících úlohách budeme, pokud nebude uvedeno jinak, používat data z datové sady ArcCR500. Tato datová sada je poskytovaná firmou ArcData zdarma.
V učebně B-870 jdou data nainstalována v adresáři
S:\K155\Public\data\ArcGIS\ArcCR500 3.3
Přístupová práva k tomu adresáři jsou omezena pouze na čtení! Nelze tedy datové vrstvy z ArcCR500 upravovat (včetně např. vytváření pyramid).
Data meteorologických stanic byla převzata z http://www.in-pocasi.cz/pocasi-na-web/meteostanice.php. Na základě ukázkového XML souboru byl vytvořen soubor ve formátu MS Excel, který je ke stažení zde.
- Souřadnice meteorologických stanic jsou uvedny ve WGS-84 (zeměpisná délka a šířka).
- Rastr interpolovaný z naměřených hodnot meteorologických stanic počítejte v prostorovém rozlišení 1km (pokud není uvedeno jinak) a v souřadnicovém systému S-JTSK. Rastr bude vždy ohraničen státní hranici ČR.
1.
- Jaká je plocha území v ha s nadmořskou výškou mezi 500 a 700m?
- Datové vrstvy: DigitalniModelReliefu
- Výsledek: 1 970 760ha
("DigitalniModelReliefu" > 500) & ("DigitalniModelReliefu" < 700)
2.
- Jaká je výměra území v ha pro kterou platí, že leží v nadmořské výšce nad 700m a má sklon svahu větší než 25 gonů?
- Datové vrstvy: DigitalniModelReliefu
- Výsledek: 12 401ha
("DigitalniModelReliefu" > 700) & (Slope("DigitalniModelReliefu") * 10 / 9 > 25)
3.
- Jaký průměrný sklon mají svahy, které jsou vzdáleny do 10km od státní hranice. Jak velký rozdíl to je oproti průměrné hodnotě počítané pro celé území státu?
- Datové vrstvy: DigitalniModelReliefu, StatPolygon (AC)
- Výsledek: 4,82° -1,52°
4.
- Jaká je plocha území v ha, kde se sklon limitně blíží k nule?
- Datové vrstvy: DigitalniModelReliefu
- Výsledek: 1 197 567ha
Aspect("DigitalniModelReliefu") < 0
5.
- Vytvořte pomocí Raster Calculatoru rastr, který obsahuje hodnotu 1 pro území, kde je nadmořská výška nad 700m a sklon menší než 5°; hodnotu 2, kde je platí, že je nadmořská výška nad 700m a sklon je větší než 5°. Jaká je výměra takto určeného území v ha?
- Datové vrstvy: DigitalniModelReliefu
- Výsledek: 1 - 249 974ha; 2 - 361 191ha
Con( ("DigitalniModelReliefu" > 700) & (Slope("DigitalniModelReliefu") < 5), 1, Con(("DigitalniModelReliefu" > 700) & (Slope("DigitalniModelReliefu") > 5), 2) )
6.
- Pro území Ústeckého kraje vytvořte rastr s prostorovým rozlišením 100m, jehož buňky mají hodnoty s normálním rozdělením.
- Datové vrstvy: KrajePolygony
CreateNormalRaster(100, Extent(-842071.81,-1033413.16,-714001.93,-935236,59))
7.
- Vypočítejte DMT s využitím výškových kót, vrstevnic, vodních toků, vodních ploch a státní hranice ČR. Dále vypočítejte DMT pouze s využitím výškových kót, vrstevnic a státní hranice ČR. Oba rastry vytvořte s prostorovým rozlišením 1km. Minimální Z hodnotu nastavte na 0. Jaká je výměra území v ha, pro kterou je absolutní hodnota rozdílu výšek DMT větší než 100 metrů?
- Datové vrstvy: VyskoveKoty, Vrstevnice, VodniToky, VodniPlochy, StatPolygon (AC)
- Výsledek: 11 500ha
Další úlohy najdete v cvičení k předmětu 153GIS2 (2.cvičení) a 153GIS2 (3.cvičení).