# register GeoTIFF file in GRASS database: r.external terra_lst1km20030314.LST_Day.tif out=modis_celsius # define output directory for files resulting from GRASS calculation: r.external.out directory=$HOME/gisoutput/ format="GTiff" # perform GRASS calculation (here: extract pixels > 20 deg C) # write output directly as GeoTIFF: r.mapcalc "warm.tif = if(modis_celsius > 20, modis_celsius, null() )" # cease GDAL output connection and turn back to write GRASS raster files: r.external.out -r # use the result elsewhere qgis $HOME/gisoutput/warm.tif
Last changed: $Date: 2014-05-29 06:23:58 +0200 (Thu, 29 May 2014) $