<div dir="ltr"><br><br>On Mon, Jul 9, 2018 at 6:11 PM, Rich Shepard <<a href="mailto:rshepard@appl-ecosys.com">rshepard@appl-ecosys.com</a>> wrote:<br>><br>> On Mon, 9 Jul 2018, Markus Metz wrote:<br>><br>>> The source's region does not matter, what matters are the extents of the<br>>> source raster map.<br>><br>><br>> Markus,<br>><br>>   I should have remembered that they're different.<br>><br>>> Apparently the input raster, when reprojected to the current region in the<br>>> location/mapset, does not overlap with the current region.<br>>><br>>> You can check the reprojected extents of the input raster with r.proj -g.<br>><br>><br>>   This I did:<br>><br>> r.proj loc=elwood map=PERMANENT in=elwood_dem2013 -g<br>> WARNING: Input and output locations are the same<br>> Input map <elwood_dem2013@PERMANENT> in location <elwood>:<br>> n=1281118.49999997 s=1235035.49999995 w=828457.49999998 e=861685.49999999 rows=15361 cols=11076<br>><br>>   Then I followed the 'box' map method of setting the target region as shown<br>> in the manual page's Notes, paragraph 4:<br>><br>> # In source location:<br><div>></div><div>first set the region to match the raster:</div><div><br></div><div>g.region raster=elwood_dem2013<br></div><div><br></div><div>now proceed as before:</div><div><br></div>> v.in.region -d out=elwood_region<br>><br>> # In target location:<br>><br>> v.proj in=elwood_region loc=elwood map=PERMANENT out=elwood_reproj<br>><br>> g.region -g<br>> projection=99<br>> zone=0<br>> n=201663.78224162<br>> s=137006.09426117<br>> w=2314182.92402788<br><div>> e=2388673.19838691</div><div><br></div><div>this:<br></div>> nsres=0.99999517<br><div>> ewres=1.00000368</div><div><br></div><div>is probably not what you want, it is recommended to align the region to exactly 1 map unit:</div><div>g.region -p -a res=1</div><div>before reprojection<br></div><div><br></div>> rows=64658<br>> cols=74490<br>> cells=4816374420<br>><br>> g.region vect=elwood_reproj res=1 -a<br>><br>> r.proj loc=elwood map=PERMANENT in=elwood_dem2013 res=1 method=lanczos<br>> ERROR: Input raster map is outside current region<br>><br><div>>   Please point out what I missed or did incorrectly.</div><div><br></div><div>You should notice a difference when using the suggested g.region calls.</div><div><br></div><div>Markus M<br></div>><br>><br>> Best regards,<br>><br>> Rich<br>><br>> _______________________________________________<br>> grass-user mailing list<br>> <a href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a><br>> <a href="https://lists.osgeo.org/mailman/listinfo/grass-user">https://lists.osgeo.org/mailman/listinfo/grass-user</a><br><br></div>