[GRASS-user] some questions on r.tile
Blumentrath, Stefan
Stefan.Blumentrath at nina.no
Thu Aug 11 02:21:57 PDT 2016
In the meantime you could try the following as a workaround:
Create a grid (v.mkgrid) for your tiles,
then loop over the tiles and set the region to your tile plus / minus the overlap you like to have at the boundaries after checking that current boundary does not exceed north or south limits.
It seems that g.region does not wrap latitudinal boundaries in the north, while it does for longitudinal limits and rounds latitudinal boundaries in the south.
Compare output from:
g.region -p n=90 s=-90 w=-180 e=180
g.region -p n=n+1 s=s-1 w=w-1 e=e+1
g.region -p n=n s=s-1 w=w-1 e=e+1
in a lat/lon location.
Cheers
Stefan
More information about the grass-user
mailing list