153ZODH / 13. cvičení: Porovnání verzí
| Řádek 44: | Řádek 44: | ||
<source lang="bash"> | <source lang="bash"> | ||
g.region rast=lsat7_2002_1 | g.region rast=lsat7_2002_1 | ||
r.out.gdal input=lsat7_2002 format=GTiff type=Byte output=lsat7_2002 | r.out.gdal input=lsat7_2002 format=GTiff type=Byte output=lsat7_2002.tif | ||
</source> | </source> | ||
| Řádek 50: | Řádek 50: | ||
<source lang="bash"> | <source lang="bash"> | ||
gdalinfo lsat7_2002 | grep Band | gdalinfo lsat7_2002.tif | grep Band | ||
</source> | </source> | ||
Verze z 30. 11. 2009, 22:14
< Stránky předmětu • Předchozí cvičení • Další cvičení
Úvod do GRASS/R
Osnova
Toto cvičení je zaměřeno propojení GRASS GIS a prostředí pro statistickou analýzu dat R.
Seznam příkazů
Rozhraní GRASS/R
Základní informace o rozhraní R pro GRASS GIS jsou dostupné zde.
Výpočet NDVI
Spustíme GRASS s lokací nc_spm_08. K dispozici máme snímky LandSat 7 z roku 2002.
g.mlist type=rast pattern='lsat7*'
lsat7_2002_10 lsat7_2002_20 lsat7_2002_30 lsat7_2002_40 lsat7_2002_50 lsat7_2002_61 lsat7_2002_62 lsat7_2002_70 lsat7_2002_80
Vytvoříme obrazovou skupinu, která obsahuje 1-3, 4-5 a 7 pásmo.
i.group group=lsat7_2002 input=lsat7_2002_10,lsat7_2002_20,lsat7_2002_30,lsat7_2002_40,lsat7_2002_50,lsat7_2002_70
Nastavíme region a data vyexportujeme do formátu GeoTIFF.
g.region rast=lsat7_2002_1
r.out.gdal input=lsat7_2002 format=GTiff type=Byte output=lsat7_2002.tif
Soubor by měl obsahovat 6 kanálů.
gdalinfo lsat7_2002.tif | grep Band
Band 1 Block=527x2 Type=Byte, ColorInterp=Gray Band 2 Block=527x2 Type=Byte, ColorInterp=Undefined Band 3 Block=527x2 Type=Byte, ColorInterp=Undefined Band 4 Block=527x2 Type=Byte, ColorInterp=Undefined Band 5 Block=527x2 Type=Byte, ColorInterp=Undefined Band 6 Block=527x2 Type=Byte, ColorInterp=Undefined
Z příkazové řádky GRASSu spustíme prostředí R a nahrajeme balíček 'spgrass6'.
GRASS (nc_spm_08):~ > R > library(spgrass6)