<div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div>Hi Even,</div></div><div><br></div><div><div class="gmail_quote"><br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">

What version of PROJ does your GDAL 2.3.2 build use ? PROJ 5 ?<br>
And your GDAL 2.2.2 build: PROJ 4 ?<br></blockquote><div><br></div><div>GDAL 2.3.2 uses PROJ 5.2.0.</div><div>GDAL 2.2.3 uses PROJ.4 version 493</div><div>GDAL 2.2.2 uses 
PROJ.4 version 480</div><div><br></div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">

If so, that might explain things. The ESRI:102164 translated as a proj string <br>
has no +towgs84, whereas EPSG:20790 has.<br>
<br>
With PROJ 4 and how GDAL uses it, due to the absence of a +towgs84 in one of <br>
the source/target CRS, no datum shift is performed, which is what is expected.<br>
<br>
With PROJ 5 and how GDAL 2.3.2 uses it, a unwanted datum shift is done (it is <br>
a bit like if the ESRI 102164 had an implicit +towgs84=0,0,0).<br>
I've traced this as <a href="https://github.com/OSGeo/gdal/issues/1156" rel="noreferrer" target="_blank">https://github.com/OSGeo/gdal/issues/1156</a><br>
</blockquote><div><br></div><div>I think this is exactly the problem!</div><div><br></div><div>C:\>gdalinfo --version<br>GDAL 2.3.2, released 2018/09/21<br><br>C:\>echo "0 0 0" | gdaltransform -s_srs "+proj=longlat +ellps=GRS80 +towgs84=100,100,100" -t_srs "+proj=longlat +ellps=GRS80"<br>0.000898301199977168 0.000904355202316605 100.001573113725<br><br><br>C:\>gdalinfo --version<br>GDAL 2.2.3, released 2017/11/20<br><br>C:\>echo "0 0 0" | gdaltransform -s_srs "+proj=longlat +ellps=GRS80 +towgs84=100,100,100" -t_srs "+proj=longlat +ellps=GRS80"<br>0 0 0<br></div><div><br></div><div><br>pedro@pedro-VirtualBox:~$ gdalinfo --version<br>GDAL 2.2.2, released 2017/09/15<br><br>pedro@pedro-VirtualBox:~$ echo "0 0 0" | gdaltransform -s_srs "+proj=longlat +ellps=GRS80 +towgs84=100,100,100" -t_srs "+proj=longlat +ellps=GRS80"<br>0 0 0<br></div><div><br></div><div><br></div>Do you think this can be fixed soon?</div><div class="gmail_quote"><br></div><div class="gmail_quote">Thank you very much Even!</div><div class="gmail_quote"><br></div><div class="gmail_quote">Best regards,</div><div class="gmail_quote">Pedro</div><div class="gmail_quote"><br></div><div class="gmail_quote"><br></div></div></div></div></div>