<div><div dir="auto">That’s a problem for us. With gdal 3 we started to force values for the towgs84 so that proj would not decide.</div></div><div dir="auto">This is wanted because our front end will place markers and it does coordinates transform. Without access to proj 6. Meaning it needs to have the transform Value.</div><div dir="auto"><br></div><div dir="auto">I’ll check what you mentioned (on Monday though I think)</div><div dir="auto"><br></div><div dir="auto">Thanks for having a look</div><div dir="auto"><br></div><div dir="auto">Greg</div><div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, 8 Nov 2019 at 12:30, Even Rouault <<a href="mailto:even.rouault@spatialys.com">even.rouault@spatialys.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">On vendredi 8 novembre 2019 08:18:11 CET Grégory Bataille wrote:<br>
> Hey all,<br>
> <br>
> So I have a strange issue that I somewhat isolated, but I don't know how to<br>
> move further.<br>
> I have dataset in EPSG 31468-15949<br>
<br>
I would check if that wouldn't be related to availability of the BETA2007.gsb <br>
grid. osgeo/gdal:alpine-normal-3.0.2 has it. Does your MAC install have it ?<br>
<br>
That said, your TIFF file has a hardcoded TOWGS84=598.1,73.7,418.2,0.202           <br>
,0.045,-2.455,6.7 . When using<br>
<br>
gdalwarp Investigation_tile_alignment/ClippedProject.tif out.tif -overwrite -<br>
t_srs EPSG:4326 --debug on<br>
<br>
I see locally that BETA2007.gsb is used. This is not so bad, but a bit <br>
unexpected since I'd thought that forced TOWGS84 would have take the <br>
precedence.<br>
<br>
Checking further in ogrct.cpp, I see that as the GeoTIFF has also the EPSG:<br>
31468 code, and when instanciating the PROJ coordinate transformation, this is <br>
the preferred initializer, and thus the hardcoded TOWGS84 is ignored. So this <br>
explains what I observe. I'm not completely clear of what would be the desired <br>
behaviour here. In some use cases, one might prefer to honour the possibly <br>
suboptimal TOWGS84 from the geotiff info (the simple fact of doing <br>
gdal_translate -a_srs EPSG:31468 hardcodes it. Even with GDAL 3, which is <br>
probably no longer a good idea...). In other cases, just trust PROJ to use the <br>
best transformation method.<br>
<br>
Even<br>
<br>
-- <br>
Spatialys - Geospatial professional services<br>
<a href="http://www.spatialys.com" rel="noreferrer" target="_blank">http://www.spatialys.com</a><br>
</blockquote></div></div>-- <br><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature">Greg</div>