[GRASS-user] regular grid

Hamish hamish_nospam at yahoo.com
Thu Nov 9 19:57:54 EST 2006


Tomas Lanczos wrote:
> 
> how can I create a regular grid of proints within a given area defined
> bz its border?

g.region to set up the border and point resolution, then

r.mapcalc surface=1
r.to.vect feature=point in=surface out=dots

if you will have more than a million or three grid points you might want
to use the "-b" flag.

note this will put grid points at the middle of each raster cell
 (x= bound + 1/2 res), adjust you bounds by half a cell if this isn't
what you want.

if you just want the coords,
  r.stats -1gn surface | cut -f1-2 -d' '


Hamish




More information about the grass-user mailing list