153ZODH / 13. cvičení: Porovnání verzí
| Řádek 39: | Řádek 39: | ||
i.group group=lsat7_2002 input=lsat7_2002_10,lsat7_2002_20,lsat7_2002_30,lsat7_2002_40,lsat7_2002_50,lsat7_2002_70 | i.group group=lsat7_2002 input=lsat7_2002_10,lsat7_2002_20,lsat7_2002_30,lsat7_2002_40,lsat7_2002_50,lsat7_2002_70 | ||
</source> | </source> | ||
Nastavíme region a data vyexportujeme do formátu [http://en.wikipedia.org/wiki/Geotiff GeoTIFF]. | |||
<source lang="bash"> | |||
g.region rast=lsat7_2002_1 | |||
r.out.gdal input=lsat7_2002 format=GTiff type=Byte output=lsat7_2002 | |||
</source> | |||
Soubor by měl obsahovat 6 kanálů. | |||
<source lang="bash"> | |||
gdalinfo lsat7_2002 | grep Band | |||
</source> | |||
<pre> | |||
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 | |||
</pre> | |||
Z příkazové řádky GRASSu spustíme prostředí R a nahrajeme balíček 'spgrass6'. | Z příkazové řádky GRASSu spustíme prostředí R a nahrajeme balíček 'spgrass6'. | ||
Verze z 30. 11. 2009, 22:11
< 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
Soubor by měl obsahovat 6 kanálů.
gdalinfo lsat7_2002 | 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)