[GRASSLIST:5671] Re: sampling methods

Hamish hamish_nospam at yahoo.com
Mon Feb 7 21:11:54 EST 2005


> I wondering what is the best method to sample areas around survey
> sites. I'm using the version 6 beta.
..
> Method 2
> a. create a MASK from the raster buffer using r.mapcalc (in the hope
> to speed processing) b. run r.neighbors to calculate
> min,max,mean,variance, diversity, and mode using a window the size of
> my buffers, so that the center value becomes a measure for the entire
> buffer. c. extract to postgres using v.what.rast using the original
> points file
> 
> Most of my stuff is floating point so r.statistics doesn't work, and
> r.stats returns ranges (as far as I can get it to function) instead of
> actual values.
> 
> If any of the experience GRASS users has an opinion as to which is
> best and why or if there is an easier way to do this, I would
> appreciate the input.


what about something like

for site 1 2 3 4 5 ... ; do
  v.to.rast where='site=7'
  g.region (zoom to a little bigger than buffer dist)
  r.buffer out=MASK
  r.univar base_map
done


?
Hamish




More information about the grass-user mailing list