153GIS2 / 4. cvičení: Porovnání verzí
mBez shrnutí editace |
|||
(Není zobrazeno 15 mezilehlých verzí od 2 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] | |||
[http:// | |||
; Náplň | ; Náplň | ||
Řádek 13: | Řádek 12: | ||
* Velikostní stupnice (poměrová - ratio, intervalová - interval, pořadová, kategoriální) | * Velikostní stupnice (poměrová - ratio, intervalová - interval, pořadová, kategoriální) | ||
** [http://en.wikipedia.org/wiki/Level_of_measurement Level of measurement] | |||
* Mapová algebra | * Mapová algebra | ||
** ''Lokální'' funkce (Cell Statistics) | ** ''Lokální'' funkce (Cell Statistics) | ||
*** jedna rastrová vrstva (lokální sinus, kosinus, reklasifikace, ...) | *** jedna rastrová vrstva (lokální sinus, kosinus, reklasifikace, ...) | ||
*** více rastrových vrstev (min, max, mean, sum, product, majority, minority, variety, ...) | *** více rastrových vrstev (min, max, mean, sum, product, majority, minority, variety, ...) | ||
*** [http:// | *** [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) | ** ''Fokální'' funkce (Neighborhood Statistics) | ||
*** Fokalní min, max, mean, sum,... | *** Fokalní min, max, mean, sum,... | ||
*** Kernelové'' operace (density, slope, aspect) | *** Kernelové'' operace (density, slope, aspect) | ||
*** [http:// | *** [http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//009z000000r7000000.htm How Focal Statistics works] | ||
** ''Inkrementální'' funkce (globalní funkce) | ** ''Inkrementální'' funkce (globalní funkce) | ||
*** Např. vzdálenostní funkce, obalová zóna, cost distance, cast path | *** Např. vzdálenostní funkce, obalová zóna, cost distance, cast path | ||
** ''Zonální'' funkce (Zonal Statistics) | ** ''Zonální'' funkce (Zonal Statistics) | ||
*** Zonalní sum, mean, min, max, ... | *** Zonalní sum, mean, min, max, ... | ||
*** [http:// | *** [http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#/How_Zonal_Statistics_works/009z000000wt000000/ How Zonal Statistics works] | ||
* Spatial Analyst - Raster Calculator | * Spatial Analyst - Raster Calculator | ||
* Vyhledání funkcí v ArcTooboxu | * Vyhledání funkcí v ArcTooboxu | ||
** | ** ArcToolbox - Tool Help - Syntax | ||
* Environments, nastavení Extent pro danou funkci | * Environments, nastavení Extent pro danou funkci | ||
Řádek 65: | Řádek 62: | ||
# 5,79 °C | # 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 | #:<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> | #:<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 | #:buffer, Exctract by Mask, na nově vytvořený rastr fce Slope, a mean v atributové tabulce | ||
# 5,510 °C | # 5,510 °C | ||
Řádek 74: | Řádek 71: | ||
#:<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>(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>Hillshade("tingrid", 315, 45)</code> | ||
# <code>RasterToPolygon_conversion("tingrid" <= 600, I:/Pracovni/ | # <code>RasterToPolygon_conversion("tingrid" <= 600, "I:/Pracovni/gis2/vysledek.shp", "NO_SIMPLIFY")</code> | ||
# <code>CreateNormalRaster (10, Extent(-568952,-1061606,-558277,-1049760)) *15 + 120</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