<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<meta content="text/html; charset=windows-1255"
http-equiv="Content-Type">
</head>
<body bgcolor="#ffffff" text="#000000">
On 15/06/2010 14:35, Sandile Gumede wrote:
<blockquote
cite="mid:AANLkTin3QPP9teMqnmvVGlDtqXIyVJXjDbaoVFGdR_HV@mail.gmail.com"
type="cite">Hi<br>
If I run g.region rast=rainfall -p, I get:<br>
</blockquote>
OK, what you've done here is change the current region to match the
raster "rainfall".<br>
Can you now try: <br>
v.rast.stats -c vect=catchments rast=rainfall pref=precip<br>
<br>
<br>
<blockquote
cite="mid:AANLkTin3QPP9teMqnmvVGlDtqXIyVJXjDbaoVFGdR_HV@mail.gmail.com"
type="cite"><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 moz-do-not-send="true"
href="mailto:hamish_b@yahoo.com">hamish_b@yahoo.com</a>></span>
wrote:<br>
<blockquote class="gmail_quote"
style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; 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>
<br>
This mail was received via Mail-SeCure System.<br>
</blockquote>
<br>
<br>
<pre class="moz-signature" cols="62">--
Micha Silver
<a class="moz-txt-link-freetext" href="http://www.surfaces.co.il/">http://www.surfaces.co.il/</a>
Arava Development Co. +972-52-3665918
</pre>
</body>
</html>