<div dir="ltr"><div dir="ltr"><div>Stefan,<br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Le ven. 22 mars 2019 à 10:44, Stefan Moebius <<a href="mailto:Stefan.Moebius@amdocs.com">Stefan.Moebius@amdocs.com</a>> a écrit :<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi again,<br>
<br>
We recently encountered an issue with failing coordinate transformations.<br>
The scenario is that a current (Java-based) software using GDAL 2.4.0/PROJ<br>
5.2.0 reads HFA files generated a long time ago by another software using<br>
GDAL 1.11.1/PROJ 4.8.0. The HFA files happen to come in a DHDN projection<br>
system. The current software then transforms the files into a UTM projection<br>
system and fails.<br>
<br>
Investigating this further, we've come up with the following minimal test<br>
showing the issue:<br>
<br>
The current version of gdalinfo returns this coordinate system:<br>
PROJCS["DHDN_VfD2_Germany_zone_3",<br>
    GEOGCS["GCS_DHDN",<br>
        DATUM["Deutsche_Hauptdreiecksnetz",<br>
            SPHEROID["Bessel_1841",6377397.16,299.15281],<br>
            TOWGS84[583,68,399.5,-0,-0,578614.3160276701,11300000]],<br></blockquote><div><br></div></div><div>--> OK, the TOWGS84 is wrong here. The rotation (4th,5th,6th) and scale factor difference (7th) terms of the TOWGS84 are obviously nonsense.</div><div>They should be respectively in arc-seconds and in PPM.</div><div>From what you mention, I believe this is due to<br></div><div><a href="https://github.com/OSGeo/gdal/commit/e3d18d2070d8d73eef37ff3b844b7e441ffb5149">https://github.com/OSGeo/gdal/commit/e3d18d2070d8d73eef37ff3b844b7e441ffb5149</a></div><div>which is likely correct by itself.</div><div>The issue here is that as you mentionned the file was produced by GDAL 1.11.1 which incorrectly stored the values</div><div>of the rotation and sale factor difference as arc-seconds and PPM whereas Erdas expected radians and unity-based scale factor difference.</div><div>Thus when reading this with GDAL >=2.2, the file is wrongly interpretated.</div><div>You should override the TOWGS84 values of this WKT with the corrected ones.</div><div><br></div><div>Even<br></div><div><br></div></div></div>