[GRASS-user] Setting region to add map to location

Rich Shepard rshepard at appl-ecosys.com
Mon Jul 9 10:44:37 PDT 2018


On Mon, 9 Jul 2018, Markus Metz wrote:

> first set the region to match the raster:
> g.region raster=elwood_dem2013
> now proceed as before:

Markus,

   Okay. I thought the vector box matched that of the raster. I'll make that
change.

> this:
>> nsres=0.99999517
>> ewres=1.00000368
>
> is probably not what you want, it is recommended to align the region to
> exactly 1 map unit:
> g.region -p -a res=1
> before reprojection

   I think this is what I did:

>> g.region vect=elwood_reproj res=1 -a
>> r.proj loc=elwood map=PERMANENT in=elwood_dem2013 res=1 method=lanczos

> You should notice a difference when using the suggested g.region calls.

   I'll set the source region to the raster map prior to creating the vector
box. Will post results.

Thanks,

Rich


More information about the grass-user mailing list