153GIS2 - 4. cvičení - GRASS GIS: Porovnání verzí
m (→4.) |
(→6.) |
||
Řádek 72: | Řádek 72: | ||
;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 15x15 pixelů. Jakou teplotu má bod o souřadnicích 564 300, 1 057 100? | ||
<source lang="bash"> | |||
r.neighbors input=dem_idw output=fok15 method=average size=15 | |||
v.drape input=obod type=point rast=fok15 output=obod_fok method=bilinear | |||
v.out.ascii obod_fok | |||
</source> | |||
===== 7. ===== | ===== 7. ===== |
Verze z 11. 3. 2010, 13:23
Stránka obsahuje řešení úloh z 4. cvičení 153GIS2 v GRASS GIS.
Pro rastry používejte pixel 10m, pokud není uvedeno jinak.
g.region rast=dem
1.
- Jaká je plocha území o nadmořské výšce 650 až 750 metrů? Plochu určete pomocí Raster Calculatoru.
r.mapcalc 'dem_650_750 = if(dem >= 650 && dem <= 750, 1, null())'
r.report -hn map=dem_650_750 units=h
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).
@TODO kriging
echo "-564300|-1057100" | v.in.ascii input=- output=obod
v.surf.idw input=stare_mesto_teploty column=Teplota output=dem_idw
v.surf.rst input=stare_mesto_teploty zcolumn=Teplota elev=dem_rst
r.mapcalc 'dem_3 = (dem_idw + dem_rst) / 2'
v.drape input=obod type=point rast=dem_3 output=obod_3 method=bilinear
v.out.ascii obod_3
3.
- Jaká je plocha území, kde je sklon větší nebo roven jak 20 gonů a menší nebo roven jak 25 gonů?
r.slope.aspect elevation=dem slope=slope
r.mapcalc 'slope_g = if(slope * 10/9 >= 20 && slope * 10/9 <=25, 1, null())'
r.report -hn map=slope_g units=h
4.
- Jaký průměrný sklon mají svahy, které jsou vzdáleny do 100 metrů od silnic?
v.buffer input=stare_mesto_silnice output=silnice_100 distance=100
v.to.rast input=silnice_100 output=MASK use=val
r.univar -g map=slope | grep 'mean='
r.mask -r
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?
r.mapcalc 'dem_min = min(dem_idw, dem_rst)'
v.drape input=obod type=point rast=dem_min output=obod_min method=bilinear
v.out.ascii obod_min
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?
r.neighbors input=dem_idw output=fok15 method=average size=15
v.drape input=obod type=point rast=fok15 output=obod_fok method=bilinear
v.out.ascii obod_fok
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)