<div dir="ltr"><br><br>On Sun, Sep 15, 2019 at 8:39 PM Rich Shepard <<a href="mailto:rshepard@appl-ecosys.com">rshepard@appl-ecosys.com</a>> wrote:<br>><br>> On Sat, 14 Sep 2019, Markus Metz wrote:<br>><br>> > gdalwarp estimates the target extents unless you use the -te option. In<br>> > contrast, r.proj uses the current region. r.proj can estimate the target<br>> > extents with the -p or -g flag. I suggest to first check you current<br>> > region in the target location/mapset and compare with r.proj -p or -g<br>><br>> Markus M, et al.:<br>><br>> I've worked on this for 2+ days and cannot get the DEM from the source<br>> location to the target location when both have the same projection and<br>> units.<br>><br>> Following the above still fails and I've no idea what to do now.<br>><br>> Source location:<br>> > r.proj -p loc=Oregon map=topography in=n45w121<br>> Input map <n45w121@topography> in location <Oregon>:<br>> Source cols: 80297<br>> Source rows: 111356<br>> Local north: 148340.8192959<br>> Local south: 36984.78211177<br>> Local west: 2459852.09770149<br>> Local east: 2540149.13365527<br>><br>> Current (target) region:<br>> > g.region -p<br>> projection: 99 (NAD83(HARN) / Oregon North)<br>> zone:       0<br>> datum:      nad83harn<br>> ellipsoid:  grs80<br>> north:      307220.87234928<br>> south:      -199204.29675669<br>> west:       2161358.9640483<br>> east:       2822158.25992759<br>> nsres:      1.00000033<br>> ewres:      1.00000045<br>> rows:       506425<br>> cols:       660799<br>> cells:      334645133575<br>><br>> Reset current region to the source:<br>> > g.region n=148340.8192959 s=36984.78211177 w=2459852.09770149 e=2540149.13365527 cols=80297 rows=111356<br><div><br></div><div>You should now use</div><div>r.proj loc=Oregon map=topography in=n45w121</div><div><br></div><div>and not r.unpack, that is a completely different workflow that has nothing to do with reprojection.</div><div><br></div><div>Markus M</div><div><br></div>>> Unpack the DEM in the current location/mapset;<br>> > r.unpack n45w121.pack<br>> WARNING: Difference between PROJ_INFO file of packed map and of current<br>>           location:<br>>           name: NAD83(HARN) / Oregon North<br>>           datum: nad83harn<br>>           ellps: grs80<br>>           proj: lcc<br>>           lat_1: 46<br>>           lat_2: 44.33333333333334<br>>           lat_0: 43.66666666666666<br>>           lon_0: -120.5<br>>           x_0: 2500000<br>>           y_0: 0<br>>           no_defs: defined<br>>           - towgs84: 0.000,0.000,0.000<br>> ERROR: Projection of dataset does not appear to match current location. In<br>>         case of no significant differences in the projection definitions,<br>>         use the -o flag to ignore them and use current location definition.<br>><br>> If I correctly understand gdalwarp it functions like r.proj in reprojecting<br>> a raster map from one location to another location. Here, both the source<br>> and target locations have the same projection. Will gdalwarp work in this<br>> case? Isn't r.pack/r.unpack the correct module pair?<br>><br>> Please advise me of what I should do now or what other information I can<br>> provide.<br>><br>> Regards,<br>><br>> Rich<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></div>