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

[Rich Shepard]

On Mon, 9 Jul 2018, Markus Metz wrote:

> The source's region does not matter, what matters are the extents of the
> source raster map.

Markus,

   I should have remembered that they're different.

> Apparently the input raster, when reprojected to the current region in the
> location/mapset, does not overlap with the current region.
>
> You can check the reprojected extents of the input raster with r.proj -g.

   This I did:

r.proj loc=elwood map=PERMANENT in=elwood_dem2013 -g
WARNING: Input and output locations are the same
Input map <elwood_dem2013 at PERMANENT> in location <elwood>:
n=1281118.49999997 s=1235035.49999995 w=828457.49999998 e=861685.49999999 rows=15361 cols=11076

   Then I followed the 'box' map method of setting the target region as shown
in the manual page's Notes, paragraph 4:

# In source location:

v.in.region -d out=elwood_region

# In target location:

v.proj in=elwood_region loc=elwood map=PERMANENT out=elwood_reproj

g.region -g
projection=99
zone=0
n=201663.78224162
s=137006.09426117
w=2314182.92402788
e=2388673.19838691
nsres=0.99999517
ewres=1.00000368
rows=64658
cols=74490
cells=4816374420

g.region vect=elwood_reproj res=1 -a

r.proj loc=elwood map=PERMANENT in=elwood_dem2013 res=1 method=lanczos
ERROR: Input raster map is outside current region

   Please point out what I missed or did incorrectly.

Best regards,

Rich