[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