Hi<br>If I run g.region rast=rainfall -p, I get:<br><br>projection: 3 (Latitude-Longitude)<br>zone: 0<br>datum: wgs84<br>ellipsoid: wgs84<br>north: 33:40:46.49916S<br>south: 34:20:55.49928S<br>west: 18:17:55.50036E<br>
east: 19:10:16.50036E<br>nsres: 0:00:05.01875<br>ewres: 0:00:02.18125<br>rows: 480<br>cols: 1440<br>cells: 691200<br><br>and If I run r.univar rainfall, I get the following output:<br><br>
100%<br>total null and non-null cells: 691200<br>total null cells: 0<br><br>Of the non-null cells:<br>----------------------<br>n: 691200<br>minimum: 0<br>maximum: 3094<br>range: 3094<br>mean: 22.0228<br>mean of absolute values: 22.0228<br>
standard deviation: 76.1639<br>variance: 5800.94<br>variation coefficient: 345.841 %<br>sum: 15222164<br><br><br><br><div class="gmail_quote">On Tue, Jun 15, 2010 at 12:22 PM, Hamish <span dir="ltr"><<a href="mailto:hamish_b@yahoo.com">hamish_b@yahoo.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div class="im">Micha wrote:<br>
> The only unusual thing I notice above is that the resolution settings<br>
> for the raster are different N-S and E-W. This came from the original<br>
> tiff (see below) which also has rectangular pixels,<br>
<br>
</div>that is perfectly normal for a lat/lon map away from the equator.<br>
longitude scales a cos(lat).<br>
<div class="im"><br>
<br>
> (the v.rast.stats module creates a temp raster at the *current region's<br>
> resolution* settings, which might be different from this rainfall<br>
> raster's rectangular resolution...)<br>
<br>
</div>the results of:<br>
<br>
g.region -p rast=mapname<br>
r.univar mapname<br>
<br>
<br>
could help.<br>
<font color="#888888"><br>
<br>
Hamish<br>
<br>
<br>
<br>
</font></blockquote></div><br><br clear="all"><br>-- <br>Kind Regards<br>TS Gumede<br>CSIR, Meraka Institute<br>072 258 1650<br><br>