153ZODH / 13. cvičení

Tento článek není dokončen nebo potřebuje upravit

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

< Stránky předmětu • Předchozí cvičení • Další cvičení