# 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$