On 15/06/2010 14:35, Sandile Gumede wrote:
Hi
If I run g.region rast=rainfall -p, I get:
OK, what you've done here is change the current region to match the
raster "rainfall".
Can you now try:
v.rast.stats -c vect=catchments rast=rainfall pref=precip
projection: 3 (Latitude-Longitude)
zone: 0
datum: wgs84
ellipsoid: wgs84
north: 33:40:46.49916S
south: 34:20:55.49928S
west: 18:17:55.50036E
east: 19:10:16.50036E
nsres: 0:00:05.01875
ewres: 0:00:02.18125
rows: 480
cols: 1440
cells: 691200
and If I run r.univar rainfall, I get the following output:
100%
total null and non-null cells: 691200
total null cells: 0
Of the non-null cells:
----------------------
n: 691200
minimum: 0
maximum: 3094
range: 3094
mean: 22.0228
mean of absolute values: 22.0228
standard deviation: 76.1639
variance: 5800.94
variation coefficient: 345.841 %
sum: 15222164
On Tue, Jun 15, 2010 at 12:22 PM, Hamish <[email protected]>
wrote:
Micha wrote:
> The only unusual thing I notice above is that the resolution
settings
> for the raster are different N-S and E-W. This came from the
original
> tiff (see below) which also has rectangular pixels,
that is perfectly normal for a lat/lon map away from the equator.
longitude scales a cos(lat).
> (the v.rast.stats module creates a temp raster at the *current
region's
> resolution* settings, which might be different from this rainfall
> raster's rectangular resolution...)
the results of:
g.region -p rast=mapname
r.univar mapname
could help.
Hamish
--
Kind Regards
TS Gumede
CSIR, Meraka Institute
072 258 1650
This mail was received via Mail-SeCure System.
--
Micha Silver
http://www.surfaces.co.il/
Arava Development Co. +972-52-3665918
|
_______________________________________________
grass-user mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/grass-user
