153GIS2 / 4. cvičení: Porovnání verzí
m →Úkoly |
mBez shrnutí editace |
||
(Není zobrazeno 48 mezilehlých verzí od 3 dalších uživatelů.) | |||
Řádek 1: | Řádek 1: | ||
{{Zastaralé|155GIS2}} | |||
{{Cvičení|153GIS2|4|Spatial analyst, Map Algebra, Raster calculator, extraction, raster creation}} | |||
== Zadání == | |||
[http://josef.fsv.cvut.cz/~cajthaml/vyuka/gis2/gis2_cviceni4.zip Data] | |||
; Náplň | |||
[http:// | * [http://web.pdx.edu/~jduh/courses/geog475f09/Students/W5_Raster%20presentation%2010-27-09.pdf Raster Operations] | ||
* [http://gis.fzp.ujep.cz/files/8.pdf Rastrový datový model a mapová algebra] | |||
* [http://web.natur.cuni.cz/IGP/main/staff/martinek/DPZdata/2-RastrovaData.pdf Rastrová data] | |||
* Velikostní stupnice (poměrová - ratio, intervalová - interval, pořadová, kategoriální) | |||
** [http://en.wikipedia.org/wiki/Level_of_measurement Level of measurement] | |||
* Mapová algebra | |||
** ''Lokální'' funkce (Cell Statistics) | |||
*** jedna rastrová vrstva (lokální sinus, kosinus, reklasifikace, ...) | |||
*** více rastrových vrstev (min, max, mean, sum, product, majority, minority, variety, ...) | |||
*** [http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#/How_Cell_Statistics_works/009z00000088000000/ How Cell Statistics works] | |||
** ''Fokální'' funkce (Neighborhood Statistics) | |||
*** Fokalní min, max, mean, sum,... | |||
*** Kernelové'' operace (density, slope, aspect) | |||
*** [http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//009z000000r7000000.htm How Focal Statistics works] | |||
** ''Inkrementální'' funkce (globalní funkce) | |||
*** Např. vzdálenostní funkce, obalová zóna, cost distance, cast path | |||
** ''Zonální'' funkce (Zonal Statistics) | |||
*** Zonalní sum, mean, min, max, ... | |||
*** [http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#/How_Zonal_Statistics_works/009z000000wt000000/ How Zonal Statistics works] | |||
* Spatial Analyst - Raster Calculator | |||
* Vyhledání funkcí v ArcTooboxu | |||
** ArcToolbox - Tool Help - Syntax | |||
* Environments, nastavení Extent pro danou funkci | |||
* Extrahování dat (Spatial Analyst - Extraction) | |||
** Exctract By Mask | |||
* Vytvoření nové rastrové vrstvy (Spatial Analyst - Raster Creation) | |||
** Konstantní rastr, rastr s normálním rozdělením a náhodný rastr | |||
==Úkoly== | ==Úkoly== | ||
Řádek 13: | Řádek 50: | ||
# Jaká je plocha území, kde je sklon větší nebo roven jak 20 gonů a menší nebo roven jak 25 gonů? | # Jaká je plocha území, kde je sklon větší nebo roven jak 20 gonů a menší nebo roven jak 25 gonů? | ||
# Jaký průměrný sklon mají svahy, které jsou vzdáleny do 100 metrů od silnic? | # Jaký průměrný sklon mají svahy, které jsou vzdáleny do 100 metrů od silnic? | ||
# Vytvořte rastr, který bude obsahovat pro každý pixel minimální hodnotu z interpolací IDW, Kriging, Spline (výchozí nastavení, pixel 10m). Jakou teplotu má bod | # Vytvořte rastr, který bude obsahovat pro každý pixel minimální hodnotu z interpolací IDW, Kriging, Spline (výchozí nastavení, pixel 10m). Jakou teplotu má bod o souřadnicích 564 300, 1 057 100? | ||
o souřadnicích 564 300, 1 057 100? | # Vytvořte rastr, který vznikne z interpolace IDW (výchozí nastavení, pixel 10m) a následně fokální funkcí jako průměrná hodnota (Focal Mean) z oblasti 15x15 pixelů. Jakou teplotu má bod o souřadnicích 564 300, 1 057 100? | ||
# Vytvořte rastr, který vznikne z interpolace IDW (výchozí nastavení, pixel 10m) a následně fokální funkcí jako průměrná hodnota (Focal Mean) z oblasti | |||
Jakou teplotu má bod o souřadnicích 564 300, 1 057 100? | |||
# Vyzkoušejte si vytvořit pomocí Raster Calculatoru vrstvu stínovaného reliéfu s azimutem osvětlení 315 stupňů a výškou 45 stupnů. | # Vyzkoušejte si vytvořit pomocí Raster Calculatoru vrstvu stínovaného reliéfu s azimutem osvětlení 315 stupňů a výškou 45 stupnů. | ||
# Vyzkoušejte si vytvořit shapefile s pomocí Raster Calculatoru, který obsahuje pouze plochu území s nadmořskou výškou do 600 metrů. | # Vyzkoušejte si vytvořit shapefile s pomocí Raster Calculatoru, který obsahuje pouze plochu území s nadmořskou výškou do 600 metrů. | ||
Řádek 23: | Řádek 58: | ||
==Výsledky== | ==Výsledky== | ||
# 2068 ha | # 2068 ha | ||
# 5,79 °C | #:<code>("tingrid" >= 650) & ("tingrid" <= 750)</code>, kde tingrid je rastr DEM | ||
# | # 5,79 °C | ||
# | #:<code>(Spline("stare_mesto_teploty", "Teplota", 10) + Kriging("stare_mesto_teploty", "Teplota", "Spherical", 10) + Idw("stare_mesto_teploty", "Teplota", 10)) / 3</code>, na vytvořenou vrstvu aplikovat fci "Add Surface Information" + vytvořit si bod o daných souřadnicích | ||
# | # 1039 ha | ||
# | #:<code>((Slope("tingrid") * 10 / 9 ) >= 20) & ((Slope("tingrid") * 10 / 9 ) <= 25)</code> | ||
# 11,11° | |||
#:buffer, Exctract by Mask, na nově vytvořený rastr fce Slope, a mean v atributové tabulce | |||
# 5,510 °C | |||
#:<code>(CellStatistics([Idw("stare_mesto_teploty", "Teplota", 10), Kriging("stare_mesto_teploty", "Teplota", "Spherical", 10), Spline("stare_mesto_teploty", "Teplota", 10)], "MINIMUM"))</code>, na vytvořenou vrstvu aplikovat fci "Add Surface Information" + vytvořit si bod o daných souřadnicích | |||
# 5,511 °C | |||
#:<code>(FocalStatistics(Idw("stare_mesto_teploty", "Teplota", 10), NbrRectangle(15, 15, "CELL"), "MEAN"))</code>, na vytvořenou vrstvu aplikovat fci "Add Surface Information" + vytvořit si bod o daných souřadnicích | |||
# <code>Hillshade("tingrid", 315, 45)</code> | |||
# <code>RasterToPolygon_conversion("tingrid" <= 600, "I:/Pracovni/gis2/vysledek.shp", "NO_SIMPLIFY")</code> | |||
# <code>CreateNormalRaster(10, Extent(-568952,-1061606,-558277,-1049760)) *15 + 120</code> | |||
{{GIS}} | {{GIS}} |
Aktuální verze z 16. 2. 2015, 15:41
Spatial analyst, Map Algebra, Raster calculator, extraction, raster creation
Zadání
- Náplň
- Velikostní stupnice (poměrová - ratio, intervalová - interval, pořadová, kategoriální)
- Mapová algebra
- Lokální funkce (Cell Statistics)
- jedna rastrová vrstva (lokální sinus, kosinus, reklasifikace, ...)
- více rastrových vrstev (min, max, mean, sum, product, majority, minority, variety, ...)
- How Cell Statistics works
- Fokální funkce (Neighborhood Statistics)
- Fokalní min, max, mean, sum,...
- Kernelové operace (density, slope, aspect)
- How Focal Statistics works
- Inkrementální funkce (globalní funkce)
- Např. vzdálenostní funkce, obalová zóna, cost distance, cast path
- Zonální funkce (Zonal Statistics)
- Zonalní sum, mean, min, max, ...
- How Zonal Statistics works
- Lokální funkce (Cell Statistics)
- Spatial Analyst - Raster Calculator
- Vyhledání funkcí v ArcTooboxu
- ArcToolbox - Tool Help - Syntax
- Environments, nastavení Extent pro danou funkci
- Extrahování dat (Spatial Analyst - Extraction)
- Exctract By Mask
- Vytvoření nové rastrové vrstvy (Spatial Analyst - Raster Creation)
- Konstantní rastr, rastr s normálním rozdělením a náhodný rastr
Úkoly
(pro rastry používejte pixel 10 m, pokud není uvedeno jinak)
- Jaká je plocha území o nadmořské výšce 650 až 750 metrů? Plochu určete pomocí Raster Calculatoru.
- Jakou teplotu má bod o souřadnicích 564 300, 1 057 100? Rastr, ze kterého budete teplotu určovat vypočítejte jako průměr z metod IDW, Kriging a Spline. U všech metod použijte výchozí nastavení a velikost pixelu 10 metrů. (K interpolaci teplot využijte shapefile stare_mesto_teploty).
- Jaká je plocha území, kde je sklon větší nebo roven jak 20 gonů a menší nebo roven jak 25 gonů?
- Jaký průměrný sklon mají svahy, které jsou vzdáleny do 100 metrů od silnic?
- Vytvořte rastr, který bude obsahovat pro každý pixel minimální hodnotu z interpolací IDW, Kriging, Spline (výchozí nastavení, pixel 10m). Jakou teplotu má bod o souřadnicích 564 300, 1 057 100?
- Vytvořte rastr, který vznikne z interpolace IDW (výchozí nastavení, pixel 10m) a následně fokální funkcí jako průměrná hodnota (Focal Mean) z oblasti 15x15 pixelů. Jakou teplotu má bod o souřadnicích 564 300, 1 057 100?
- Vyzkoušejte si vytvořit pomocí Raster Calculatoru vrstvu stínovaného reliéfu s azimutem osvětlení 315 stupňů a výškou 45 stupnů.
- Vyzkoušejte si vytvořit shapefile s pomocí Raster Calculatoru, který obsahuje pouze plochu území s nadmořskou výškou do 600 metrů.
- Vyzkoušejte si vytvořit rastr, jehož buňky mají hodnoty s normálním rozdělením o střední hodnotě 120 a směrodatné odchylce 15. Jako Output Extent používejte vrstvu obrys, velikost pixelu 10 metrů. (funkce Create Normal Raster vytváří rastr o N(0,1)
Výsledky
- 2068 ha
("tingrid" >= 650) & ("tingrid" <= 750)
, kde tingrid je rastr DEM
- 5,79 °C
(Spline("stare_mesto_teploty", "Teplota", 10) + Kriging("stare_mesto_teploty", "Teplota", "Spherical", 10) + Idw("stare_mesto_teploty", "Teplota", 10)) / 3
, na vytvořenou vrstvu aplikovat fci "Add Surface Information" + vytvořit si bod o daných souřadnicích
- 1039 ha
((Slope("tingrid") * 10 / 9 ) >= 20) & ((Slope("tingrid") * 10 / 9 ) <= 25)
- 11,11°
- buffer, Exctract by Mask, na nově vytvořený rastr fce Slope, a mean v atributové tabulce
- 5,510 °C
(CellStatistics([Idw("stare_mesto_teploty", "Teplota", 10), Kriging("stare_mesto_teploty", "Teplota", "Spherical", 10), Spline("stare_mesto_teploty", "Teplota", 10)], "MINIMUM"))
, na vytvořenou vrstvu aplikovat fci "Add Surface Information" + vytvořit si bod o daných souřadnicích
- 5,511 °C
(FocalStatistics(Idw("stare_mesto_teploty", "Teplota", 10), NbrRectangle(15, 15, "CELL"), "MEAN"))
, na vytvořenou vrstvu aplikovat fci "Add Surface Information" + vytvořit si bod o daných souřadnicích
Hillshade("tingrid", 315, 45)
RasterToPolygon_conversion("tingrid" <= 600, "I:/Pracovni/gis2/vysledek.shp", "NO_SIMPLIFY")
CreateNormalRaster(10, Extent(-568952,-1061606,-558277,-1049760)) *15 + 120