# 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 $HOME/gisoutput/ # do something (here: extract pixels > 20°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$