155GIS2 / 3. cvičení: Porovnání verzí

Z GeoWikiCZ
(Založena nová stránka s textem „{{GIS}} {{Cvičení|155GIS2|4|Interpolace rastrových dat}} == Základní pojmy == * [http://en.wikipedia.org/wiki/Inverse_distance_weighting IDW (metoda…“)
 
mBez shrnutí editace
Řádek 1: Řádek 1:
{{GIS}}
{{GIS}}
{{Cvičení|155GIS2|4|Interpolace rastrových dat}}
{{Cvičení|155GIS2|3|Mapová algebra}}


== Základní pojmy ==
== Základní pojmy ==
Řádek 45: Řádek 45:
* [http://training.gismentors.eu/grass-gis-zacatecnik/rastrova_data/rastrova-algebra.html Rastrová algebra]
* [http://training.gismentors.eu/grass-gis-zacatecnik/rastrova_data/rastrova-algebra.html Rastrová algebra]
-->
-->
== Úlohy ==
== Úlohy ==


{{arccr500}}
{{arccr500}}


Data ''meteorologických stanic'' byla převzata z
{{GISUloha|1|Jaká je plocha území v ha s nadmořskou výškou mezi 500 a 700m|DigitalniModelReliefu|197 106ha}}
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í '''[http://geo102.fsv.cvut.cz/~landa/vyuka/155GIS2/meteodata.zip 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.
 
{{GISUloha|1|Na základě naměřené teploty odvoďte rastr metodou IDW (výchozí hodnoty). Jaká je průměrná teplota na území ČR|meteodata, StatPolygon (AC)|3,5°C}}
 
{{GISUloha|2|Jaká je průměrná teplota v nadmořské výšce větší než 700 m při použití rastru vypočteného metodou Kriging (výchozí hodnoty)|meteodata, DigitalniModelReliefu|3,4°C}}


{{GISUloha|3|Jaká je průměrná teplota v nadmořské výšce větší než 700 m při použití rastru vypočteného metodou Spline (výchozí hodnoty)|meteodata, DigitalniModelReliefu|2,9°C}}
("DigitalniModelReliefu" > 500) & ("DigitalniModelReliefu" < 700)


{{GISUloha|4|Jaká je plocha území v ha, kde je teplota nižší než 3&deg;C (využijte interpolační metodu Natural Neighbor, prostorové rozlišení 100m) a je současně orientováno na jih. Kolik procent tohoto uzemí leží v nadmořské výšce větší než 1000m|meteodata, DigitalniModelReliefu|238 905ha; 40%}}
{{GISUloha|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ů|DigitalniModelReliefu|12 401ha}}


{{GISUloha|5|Jaká je interpolovaná hodnota teploty v reprezentačním bodě obce Peruc? (Použijte maximální hodnotu z interpolací Spline, IWD a Kriging, prostorové rozlišení 1000m, na 2 des. místa)|bod, meteodata|5,63}}
("DigitalniModelReliefu" > 700) & (Slope("DigitalniModelReliefu") * 10 / 9 > 25)


{{GISUloha|6|Vypočítejte pro reprezentační body obcí jejich teplotu. Rastr, ze kterého budete teplotu určovat vypočítejte jako průměr z metod IDW, Kriging a Spline (prostorové rozlišení 1km). Dále určete z takto vypočítaných hodnot průměrné teploty pro kraje. Který z krajů má nejnižší průměrnou teplotu a kolik to je|DigitalniModelReliefu, meteodata, ObceBody (AC)|Vysočina; 6,11&deg;}}
{{GISUloha|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|DigitalniModelReliefu, StatPolygon (AC)|4,82&deg; -1,52&deg;}}


Con("DigitalniModelReliefu", (Spline("meteodata", "teplota", 1000) +
{{GISUloha|4|Jaká je plocha území v ha, kde se sklon limitně blíží k nule|DigitalniModelReliefu|197 106ha}}
Kriging("meteodata", "teplota", "Spherical", 1000) + Idw("meteodata", "teplota", 1000)) / 3)


{{GISUloha|7|Vytvořte dva rastry teplot, které budou obsahovat pro každý pixel minimální, resp. maximální hodnotu z interpolací IDW, Kriging, Spline (výchozí nastavení, prostorové rozlišení 1km). Jaký je rozdíl takto odvozených teplot pro reprezentační bod obce Peruc|meteodata, ObceBody (AC)|0,30&deg;}}
Aspect("DigitalniModelReliefu") < 0


Abs(
{{GISUloha|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&deg;; hodnotu 2, kde je platí, že je nadmořská výška nad 700m a sklon je větší než 5&deg;. Jaká je výměra takto určeného území v ha|DigitalniModelReliefu|1 - 249 974ha; 2 - 361 191ha}}
(CellStatistics([Idw("meteodata", "teplota", 1000),
  Kriging("meteodata", "teplota", "Spherical", 1000),
  Spline("meteodata", "teplota", 1000)], "MAXIMUM"))
-
(CellStatistics([Idw("meteodata", "teplota", 1000),
  Kriging("meteodata", "teplota", "Spherical", 1000),
  Spline("meteodata", "teplota", 1000)], "MINIMUM")))
{{GISUloha|8|Vytvořte rastr teplot, který vznikne z interpolace IDW (výchozí nastavení, prostorové rozlišní 1km) a následně fokální funkcí jako průměrná hodnota (Focal Mean) z oblasti 5x5 pixelů. Jakou teplotu má oblast odpovídající reprezentačnímu bodu obce Peruc|meteodata, ObceBody (AC)|3,62}}


  (FocalStatistics(Idw("meteodata", "teplota", 1000), NbrRectangle(5, 5, "CELL"), "MEAN"))
  Con( ("DigitalniModelReliefu" > 700) & (Slope("DigitalniModelReliefu") < 5),
1,
Con(("DigitalniModelReliefu" > 700) & (Slope("DigitalniModelReliefu") > 5), 2) )

Verze z 12. 3. 2020, 10:18

Mapová algebra

Základní pojmy

Poznámky k systému ArcGIS

  • 3D Analyst Tools • Raster Interpolation • IDW
  • 3D Analyst Tools • Raster Interpolation • Kriging
  • 3D Analyst Tools • Raster Interpolation • Splines
  • 3D Analyst Tools • Raster Interpolation • Natural Neighbor
  • 3D Analyst Tools • Raster Interpolation • Topo To Raster
  • Data Management Tools • Raster • Raster Processing • Clip
  • Spatial Analyst Tools • Zona • Zonal Geometry
Rastrový kalkulátor v ArcGIS
Metoda IDW v ArcGIS

Odkazy


Ú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).

1.

Jaká je plocha území v ha s nadmořskou výškou mezi 500 a 700m?
Datové vrstvy: DigitalniModelReliefu
Výsledek: 197 106ha
("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: 197 106ha
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) )