153GIS2 / 4. cvičení: Porovnání verzí

Z GeoWikiCZ
← Porovnání se starší verzí
VizuálníWikitext
mBez shrnutí editace
 
(Není zobrazeno 48 mezilehlých verzí od 3 dalších uživatelů.)
Řádek 1: Řádek 1:
[[153GIS2 GIS 2| Hlavní stránka předmětu]] {{bullet}} [[153GIS2 - 3. cvičení|Předcházející cvičení]] {{bullet}} [[153GIS2 - 5. cvičení|Následující cvičení]]
{{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]


==Spatial analyst, Map Algebra, Raster calculator, extraction, raster creation== 
; Náplň


[http://maps.fsv.cvut.cz/~cajthaml/vyuka/gis2/gis2_cviceni4.zip Data]
* [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 5x5 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 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; [tingrid] >= 650 & [tingrid] <= 750; kde tingrid je rastr DEM
# 2068 ha
# 5,79 °C; ([Spline_stare1] + [Kriging_star1] + [Idw_stare_me1]) / 3; na vytvořenou vrstvu aplikovat fci Surface Spot + vytvořit si bod o daných souřadnicích
#:<code>("tingrid" >= 650) & ("tingrid" <= 750)</code>, kde tingrid je rastr DEM
# 1038 ha; (slope([tingrid]) * 10 / 9 ) >= 20 & ( slope([tingrid]) * 10 / 9 ) <= 25;
# 5,79 °C
# hillshade([tingrid], 315, 45)
#:<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
# gridshape(setnull([tingrid]  >= 600, 1), NOWEED)  
# 1039 ha
# [normal_rastr] * 15 + 120;
#:<code>((Slope("tingrid") * 10 / 9 ) >= 20) & ((Slope("tingrid") * 10 / 9 ) <= 25)</code>
# 10,5°; buffer, Exctract by Mask, na nově vytvořený rastr fce Slope, a mean v atributové tabulce
# 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


Tento text se vztahuje k předmětu, který již není součástí aktuálních studijních plánů.
Stránky navazujícího předmětu 155GIS2 jsou dostupné zde.
 


Spatial analyst, Map Algebra, Raster calculator, extraction, raster creation

Zadání

Data

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)
    • Inkrementální funkce (globalní funkce)
      • Např. vzdálenostní funkce, obalová zóna, cost distance, cast path
    • Zonální funkce (Zonal 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)

  1. Jaká je plocha území o nadmořské výšce 650 až 750 metrů? Plochu určete pomocí Raster Calculatoru.
  2. 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).
  3. Jaká je plocha území, kde je sklon větší nebo roven jak 20 gonů a menší nebo roven jak 25 gonů?
  4. Jaký průměrný sklon mají svahy, které jsou vzdáleny do 100 metrů od silnic?
  5. 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?
  6. 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?
  7. 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ů.
  8. Vyzkoušejte si vytvořit shapefile s pomocí Raster Calculatoru, který obsahuje pouze plochu území s nadmořskou výškou do 600 metrů.
  9. 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

  1. 2068 ha
    ("tingrid" >= 650) & ("tingrid" <= 750), kde tingrid je rastr DEM
  2. 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
  3. 1039 ha
    ((Slope("tingrid") * 10 / 9 ) >= 20) & ((Slope("tingrid") * 10 / 9 ) <= 25)
  4. 11,11°
    buffer, Exctract by Mask, na nově vytvořený rastr fce Slope, a mean v atributové tabulce
  5. 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
  6. 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
  7. Hillshade("tingrid", 315, 45)
  8. RasterToPolygon_conversion("tingrid" <= 600, "I:/Pracovni/gis2/vysledek.shp", "NO_SIMPLIFY")
  9. CreateNormalRaster(10, Extent(-568952,-1061606,-558277,-1049760)) *15 + 120
Citováno z „http://geo102.fsv.cvut.cz:8091/index.php?title=153GIS2_/_4._cvičení&oldid=30699“