<div dir="ltr"><br><br>On Sun, Jul 8, 2018 at 1:41 AM, Rich Shepard <<a href="mailto:rshepard@appl-ecosys.com">rshepard@appl-ecosys.com</a>> wrote:<br>><br>> On Sun, 8 Jul 2018, Markus Neteler wrote:<br>><br>>> r.patch does exactly that, no flags needed.<br>><br>><br>> Markus,<br>><br>>   A related issue: when I want to reproject raster DEMs from their source<br>> location to a common target location r.proj fails because the source's<br><div>> region is outside the target's region.</div><div><br></div><div>The source's region does not matter, what matters are the extents of the source raster map.<br></div><div><br></div><div>> How this can be makes no sense as the</div>> target's region is an entire state and the DEMs are 7.5' topographic quads.<br>> But, that's what r.proj tells me.<br>><br><div>>   I've not found how to resolve this issue.</div><div><br></div><div>Apparently the input raster, when reprojected to the current region in the location/mapset, does not overlap with the current region.</div><div><br></div><div>You can check the reprojected extents of the input raster with r.proj -g.<br></div><div><br></div><div>Markus M</div><div><br></div>><br>>   Example: the target location's (OR) default region is:<br>><br>> g.region -l<br>> north-west corner: long: 125:01:06.208515W lat: 46:54:29.991261N<br>> north-east corner: long: 115:30:24.620534W lat: 46:53:20.948674N<br>> south-east corner: long: 116:00:36.987691W lat: 40:41:04.56057N<br>> south-west corner: long: 124:33:45.521802W lat: 40:42:06.559349N<br>> center longitude:  120:16:28.077102W<br>> center latitude:   43:53:09.60569N<br>> rows:              691017<br>> cols:              724617<br>><br>> g.region -p<br>> projection: 99 (unnamed)<br>> zone:       0<br>> datum:      ** unknown (default: WGS84) **<br>> ellipsoid:  grs80<br>> north:      369907.96287773<br>> south:      -321108.56796525<br>> west:       2155817.84893823<br>> east:       2880434.93778086<br>> nsres:      0.99999932<br>> ewres:      1.00000012<br>> rows:       691017<br>> cols:       724617<br>> cells:      500722665489<br>><br>> and a source's smaller region (Elwood) is:<br>><br>> g.region -l<br>> north-west corner: long: 122:22:30.177355W lat: 45:14:59.998015N<br>> north-east corner: long: 122:15:00.21557W lat: 45:15:07.129823N<br>> south-east corner: long: 122:14:59.838614W lat: 45:14:54.839228N<br>> south-west corner: long: 122:22:29.773496W lat: 45:14:47.707842N<br>> center longitude:  122:18:45.005157W<br>> center latitude:   45:14:57.48019N<br>> rows:              415<br>> cols:              10732<br>><br>> g.region -p<br>> projection: 99 (unnamed)<br>> zone:       0<br>> datum:      ** unknown (default: WGS84) **<br>> ellipsoid:  grs80<br>> north:      582841.5<br>> south:      581596.5<br>> west:       7719262.5<br>> east:       7751458.5<br>> nsres:      3<br>> ewres:      3<br>> rows:       415<br>> cols:       10732<br>> cells:      4453780<br>><br>>   I need to understand why r.proj fails and how I can fix these situations.<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>