# register GeoTIFF file to be used in current mapset: r.external input=terra_lst1km20030314.LST_Day.tif output=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) # this stores the output map directly as GeoTIFF: r.mapcalc "warm.tif = if(modis_celsius > 20.0, modis_celsius, null() )" # cease GDAL output connection and turn back to write GRASS raster files: r.external.out -r # now use the resulting file elsewhere gdalinfo $HOME/gisoutput/warm.tif
Last changed: $Date$