Hi, I've tried to sort out the issue of the region resolution. Now when I'm running v.rast.stats it says:
ERROR: No categories found in raster map Here is the script I'm running thats giving me that error: #!/bin/sh #variable to customize: # path to GRASS software main directory GISBASE=/usr/lib/grass64 # path to GRASS database GISDBASE=$HOME/grassdata/Cape_ Town LOCATION_NAME=SRTMDEM MAPSET=PERMANENT # nothing to change below MAP=$1 LOCATION=$2 # generate temporal LOCATION: TEMPDIR=FLOODS mkdir -p $GISDBASE/$LOCATION_NAME/$MAPSET # save existing $HOME/.grassrc6 if test -e $HOME/.grassrc6 ; then mv $HOME/.grassrc6 /tmp/$TEMPDIR.grassrc6 fi echo "LOCATION_NAME: $LOCATION_NAME" > $HOME/.grassrc6 echo "MAPSET:$MAPSET" >> $HOME/.grassrc6 echo "DIGITIZER: none" >> $HOME/.grassrc6 echo "GISDBASE: $GISDBASE" >> $HOME/.grassrc6 export GISBASE=$GISBASE # Create a WIND file with minimal information and no projection: echo "proj: 3 zone: 0 north: 1 south: 0 east: 1 west: 0 cols: 1 rows: 1 e-w resol: 1 n-s resol: 1 top: 1 bottom: 0 cols3: 1 rows3: 1 depths: 1 e-w resol3: 1 n-s resol3: 1 t-b resol: 1 " > $GISDBASE/$LOCATION_NAME/$MAPSET/WIND # Copy WIND-file to DEFAULT_WIND: cp $GISDBASE/$LOCATION_NAME/$MAPSET/WIND \ $GISDBASE/$LOCATION_NAME/$MAPSET/DEFAULT_WIND echo "name: Latitude-Longitude datum: wgs84 towgs84: 0.000,0.000,0.000 proj: ll ellps: wgs84 "> $GISDBASE/$LOCATION_NAME/$MAPSET/PROJ_INFO echo "unit: degree ubits: degrees meters: 1.0 "> $GISDBASE/$LOCATION_NAME/$MAPSET/PROJ_UNITS export PATH=$GISBASE/bin:$GISBASE/scripts:$PATH export LD_LIBRARY_PATH=$GISBASE/lib:$LD_LIBRARY_PATH export GIS_LOCK=$$ export GISRC=$HOME/.grassrc6 # this should print GRASS version used: g.version # other calculations go here.... # import rainfall data set. cd /home/tgumede1/grassdata/Cape_Town # rainfall data set. r.in.gdal input=$HOME/grassdata/Cape_Town/TRMMLast1day.tif output=rainfall # DEM data set. r.in.gdal input=$HOME/grassdata/Cape_Town/Dem_CF.tif target=SRTMDEM output=dem g.region rast=dem -p # creating set of maps indicating flow acc, drainage dir, streams r.watershed --o elevation=...@permanent drainage=flow_direction basin=catch accumulation=acc threshold=1 memory=300 stream=str # convert catch raster to polygon vector r.to.vect in=ca...@permanent out=catchments feature=area g.region rast=rainfall -p # Calculate univariate statistics v.rast.stats -c vector=catchme...@permanent raster=rainf...@permanentcolpre=precp -- Kind Regards TS Gumede CSIR, Meraka Institute 072 258 1650
_______________________________________________ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user