<div dir="ltr"><div>Dear all,</div><div><br></div><div>I am trying to project a DSM raster layer from Belgium (<a href="https://overheid.vlaanderen.be/en/producten-diensten/digital-elevation-model" target="_blank">https://overheid.vlaanderen.be/en/producten-diensten/digital-elevation-model</a>; EPSG 31370 to EPSG 28992 (Amersfoort / RD New).</div><div><br></div><div>The result is a DEM raster layer with the features shifted roughly 100 m as compared to the same features in the AHN (Dutch DEM, <a href="https://www.ahn.nl/ahn-4" target="_blank">https://www.ahn.nl/ahn-4</a>). See screenshot and code below.<br></div><div><br></div><div>Reprojecting the same in QGIS using the warp function (gdalwarp) results in the Belgium DSM layer to be fairly closely aligned to the Dutch one.</div><div><br></div><div>I tried to set the bounds before, align to the target raster, and other suggestions in the help file of r.proj, to no avail. Any idea what I might be doing wrong?</div><div><br></div><div>Cheers,</div><div><br></div><div>Paulo<br></div><div><br></div><div><img src="cid:ii_liy51mbd0" alt="image795.png" width="558" height="253"><br></div><div><br></div><div>AHN (Netherlands)<br><br>EPSG 28992<br>projection: 99 (Amersfoort / RD New)<br>zone:       0<br>datum:      rd18<br>ellipsoid:  bessel<br>north:      397562.15979713<br>south:      373497.15979713<br>west:       107708.49748491<br>east:       133913.49748491<br>nsres:      5<br>ewres:      5<br>rows:       4813<br>cols:       5241<br>cells:      25224933<br><br>Belgium DSM<br><br>EPSG 31370<br>projection: 99 (BD72 / Belgian Lambert 72)<br>zone:       0<br>datum:      bel72<br>ellipsoid:  international<br>north:      250000<br>south:      238000<br>west:       162000<br>east:       194000<br>nsres:      5<br>ewres:      5<br>rows:       2400<br>cols:       6400<br>cells:      15360000<br><br><br>Reproject<br><br>(Thu Jun 15 17:11:08 2023)                                                      <br>r.proj location=Merkske_Belgium mapset=PERMANENT input=DHMVIIDSMRAS5m_k02 output=test memory=3000 resolution=5<br>Selected PROJ pipeline:<br>+proj=pipeline +step +inv +proj=sterea +lat_0=52.1561605555556 +lon_0=5.38763888888889 +k=0.9999079 +x_0=155000 +y_0=463000 +ellps=bessel +step +proj=lcc +lat_0=90 +lon_0=4.36748666666667 +lat_1=51.1666672333333 +lat_2=49.8333339 +x_0=150000.013 +y_0=5400088.438 +ellps=intl<br>************************<br>Selected PROJ pipeline:<br>+proj=pipeline +step +inv +proj=lcc +lat_0=90 +lon_0=4.36748666666667 +lat_1=51.1666672333333 +lat_2=49.8333339 +x_0=150000.013 +y_0=5400088.438 +ellps=intl +step +proj=sterea +lat_0=52.1561605555556 +lon_0=5.38763888888889 +k=0.9999079 +x_0=155000 +y_0=463000 +ellps=bessel<br>************************<br><br>Input:<br>Cols: 4114 (original: 6400)<br>Rows: 2400 (original: 2400)<br>North: 250000.000000 (original: 250000.000000)<br>South: 238000.000000 (original: 238000.000000)<br>West: 173430.000000 (original: 162000.000000)<br>East: 194000.000000 (original: 194000.000000)<br>EW-res: 5.000000<br>NS-res: 5.000000<br><br>Output:<br>Cols: 4108 (original: 5241)<br>Rows: 2456 (original: 4813)<br>North: 396852.159797 (original: 397562.159797)<br>South: 384572.159797 (original: 373497.159797)<br>West: 107708.497485 (original: 107708.497485)<br>East: 128248.497485 (original: 133913.497485)<br>EW-res: 5.000000<br>NS-res: 5.000000<br><br>Allocating memory and reading input raster map...<br>Selected PROJ pipeline:<br>+proj=pipeline +step +inv +proj=sterea +lat_0=52.1561605555556 +lon_0=5.38763888888889 +k=0.9999079 +x_0=155000 +y_0=463000 +ellps=bessel +step +proj=lcc +lat_0=90 +lon_0=4.36748666666667 +lat_1=51.1666672333333 +lat_2=49.8333339 +x_0=150000.013 +y_0=5400088.438 +ellps=intl<br>************************</div></div>