<div dir="ltr"><br><br>On Sat, Sep 14, 2019 at 1:15 AM Helmut Kudrnovsky <<a href="mailto:hellik@web.de">hellik@web.de</a>> wrote:<br>><br>> Rich Shepard wrote<br>> > Reprojecting a 10m DEM fails. Following are the command and input/output<br>> > projection and map description information. Screenshots of the input and<br>> > output maps are attached.<br>> ><br>> > I need to understand why the output is wrong.<br>> ><br>> > Command:<br>> > r.proj loc=Oregon map=topography in=n46w123 met=lanczos_f mem=16000 --o<br>> ><br>> > Input location:<br>> > name: unnamed<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>> > unit: Meter<br>> > units: Meters<br>> > meters: 1<br>> ><br>> > Cols: 502 (80961)<br>> > Rows: 39234 (113148)<br>> > North: 262349.035904 (262349.035904)<br>> > South: 223115.062541 (149201.112724)<br>> > West: 2302891.866970 (2302891.866970)<br>> > East: 2303393.867032 (2383852.876897)<br>> > EW-res: 1.000000<br>> > NS-res: 0.999999<br>> ><br>> > Output 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>> > unit: meter<br>> > units: meters<br>> > meters: 1<br>> ><br>> > Cols: 500 (71702)<br>> > Rows: 39232 (50757)<br>> > North: 262348.951018 (273873.984569)<br>> > South: 223116.836808 (223116.836808)<br>> > West: 2302892.017711 (2231689.553588)<br>> > East: 2303392.020970 (2303392.020970)<br>> > EW-res: 1.000007<br>> > NS-res: 1.000003Cols: 500 (71702)<br>> > Rows: 39232 (50757)<br>> > North: 262348.951018 (273873.984569)<br>> > South: 223116.836808 (223116.836808)<br>> > West: 2302892.017711 (2231689.553588)<br>> > East: 2303392.020970 (2303392.020970)<br>> > EW-res: 1.000007<br>> > NS-res: 1.000003<br>> ><br>> > Rich<br>> > _______________________________________________<br>> > grass-user mailing list<br>><br>> > grass-user@.osgeo<br>><br>> > <a href="https://lists.osgeo.org/mailman/listinfo/grass-user">https://lists.osgeo.org/mailman/listinfo/grass-user</a><br>> ><br>> > input-dem.png (296K)<br>> > &lt;<a href="http://osgeo-org.1560.x6.nabble.com/attachment/5416106/0/input-dem.png&gt">http://osgeo-org.1560.x6.nabble.com/attachment/5416106/0/input-dem.png&gt</a>;<br>> > output-dem.png (28K)<br>> > &lt;<a href="http://osgeo-org.1560.x6.nabble.com/attachment/5416106/1/output-dem.png&gt">http://osgeo-org.1560.x6.nabble.com/attachment/5416106/1/output-dem.png&gt</a>;<br>><br>> a good thing to see if reprojection of raster data works correctly, is<br>> comparing before/after r.proj with gdalwarp (1) results.<br><div><br></div><div>gdalwarp estimates the target extents unless you use the -te option. In contrast, r.proj uses the current region. r.proj can estimate the target extents with the -p or -g flag. I suggest to first check you current region in the target location/mapset and compare with r.proj -p or -g</div><div><br></div><div>Markus M<br></div><div>></div>> (1) <a href="https://gdal.org/programs/gdalwarp.html#gdalwarp">https://gdal.org/programs/gdalwarp.html#gdalwarp</a><br>><br>><br>><br>><br>> -----<br>> best regards<br>> Helmut<br>> --<br>> Sent from: <a href="http://osgeo-org.1560.x6.nabble.com/Grass-Users-f3884509.html">http://osgeo-org.1560.x6.nabble.com/Grass-Users-f3884509.html</a><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>