[GRASSLIST:187] Re: Sum for an area

Hamish hamish_nospam at yahoo.com
Thu Mar 16 20:25:43 EST 2006


> > This is probably a very dumb question but I'm stuck... I need to
> > calculate the total precipitation over a basin. I have the annual
> > precip raster and the basin raster. So, I'm issuing the command
> > 
> > r.statistics base=sto_ica cover=ppt.81 at cru method=sum output=test
> > 
> > And I get:
> > ERROR: This module currently only works for integer (CELL) maps
> 
> yes, an annoying problem...
>  
> > Of course I could easily convert the precip map (ppt.81) to integer
> > but I'm afraid that, since it's a big area, this could result in
> > some errors (sub-estimation in the final sum).
> > 
> > I'm sure there is another way to accomplish what I want, I just
> > don't know how. Or, should I not worry about errors
> > (sub-estimations)?
> 
> if you set a MASK for the basin, will
> r.sum
> do the job?

r.univar too; any thoughts from the list on if r.univar should report
sum in full instead of shortening to 1.23456e+07?


Hamish




More information about the grass-user mailing list