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