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

Z GeoWikiCZ
← Porovnání se starší verzí
VizuálníWikitext
mBez shrnutí editace
 
(Není zobrazeno 15 mezilehlých verzí od 2 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}}
==Spatial analyst, Map Algebra, Raster calculator, extraction, raster creation==
== Zadání ==
 
[http://josef.fsv.cvut.cz/~cajthaml/vyuka/gis2/gis2_cviceni4.zip Data]
[http://maps.fsv.cvut.cz/~cajthaml/vyuka/gis2/gis2_cviceni4.zip Data]


; 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://webhelp.esri.com/arcgisdesktop/9.2/index.cfm?TopicName=cell_statistics Cell statistics]
*** [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://webhelp.esri.com/arcgisdesktop/9.2/index.cfm?topicname=overlapping_neighborhood_statistics:focal_functions Overlapping neighborhood statistics: focal functions]
*** [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://webhelp.esri.com/arcgisdesktop/9.2/index.cfm?id=5240&pid=5234&topicname=Zonal_Statistics Zonal Statistics]
*** [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
** Usage (popis parametrů funkce)
 
** Recent Expressions


* Vyhledání funkcí v ArcTooboxu
* Vyhledání funkcí v ArcTooboxu
** Command Line syntax
** ArcToolbox - Tool Help - Syntax
** Scripting syntax
 
** Map Algebra 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
# 1038 ha
# 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>
# 10,5°
# 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/gis1/vysledek.shp, "NO_SIMPLIFY")</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>
# <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“