<div dir="ltr">Thank to Jukka,<div><br></div><div>so the problem was with my gdal 1.11.4 on Centos 7.</div><div>When I translated it to EPSG:4647, it doesn't contain proper WKT for EPSG, e.g:</div><div><br></div><div>Coordinate System is:<br>LOCAL_CS["ETRS89 / UTM zone 32N (zE-N)",<br>    GEOGCS["ETRS89",<br>        DATUM["unknown",<br>            SPHEROID["unretrievable - using WGS84",6378137,298.257223563]],<br>        PRIMEM["Greenwich",0],<br>        UNIT["degree",0.0174532925199433]],<br>    AUTHORITY["EPSG","4647"],<br>    UNIT["metre",1]]<br></div><div><br></div><div>compared to output from GDAL 2.2.4 on Ubuntu 18.04:</div><div><br></div><div>Coordinate System is:<br>PROJCS["ETRS89 / UTM zone 32N (zE-N)",<br>    GEOGCS["ETRS89",<br>        DATUM["European_Terrestrial_Reference_System_1989",<br>            SPHEROID["GRS 1980",6378137,298.257222101,<br>                AUTHORITY["EPSG","7019"]],<br>            TOWGS84[0,0,0,0,0,0,0],<br>            AUTHORITY["EPSG","6258"]],<br>        PRIMEM["Greenwich",0,<br>            AUTHORITY["EPSG","8901"]],<br>        UNIT["degree",0.0174532925199433,<br>            AUTHORITY["EPSG","9122"]],<br>        AUTHORITY["EPSG","4258"]],<br>    PROJECTION["Transverse_Mercator"],<br>    PARAMETER["latitude_of_origin",0],<br>    PARAMETER["central_meridian",9],<br>    PARAMETER["scale_factor",0.9996],<br>    PARAMETER["false_easting",32500000],<br>    PARAMETER["false_northing",0],<br>    UNIT["metre",1,<br>        AUTHORITY["EPSG","9001"]],<br>    AXIS["Easting",EAST],<br>    AXIS["Northing",NORTH],<br>    AUTHORITY["EPSG","4647"]]<br></div><div><br></div><div>With gdal 1.11.4, I cannot warp it to EPSG:4326 with error:</div><div><br></div><div>ERROR 1: No PROJ.4 translation for source SRS, coordinate<br>transformation initialization has failed.<br></div><div><br></div><div>and it works with gdal 2.2.4 for gdalwarp.</div><div><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, 25 Nov 2020 at 12:05, jratike80 <<a href="mailto:jukka.rahkonen@maanmittauslaitos.fi">jukka.rahkonen@maanmittauslaitos.fi</a>> wrote:<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>
Sorry, I did not read properly what you did. That should work and it does<br>
work for me with GDAL 3.2. I use a copy of your test.tif as "test_00.tif".<br>
<br>
<br>
gdal_translate -a_srs epsg:4647 test_00.tif test_000.tif<br>
Input file size is 100, 100<br>
0...10...20...30...40...50...60...70...80...90...100 - done.<br>
<br>
gdalwarp -t_srs epsg:4326 test_000.tif test_4326_2.tif<br>
Creating output file that is 125P x 73L.<br>
Processing test_000.tif [1/1] : 0Using internal nodata values (e.g. -9999)<br>
for image test_000.tif.<br>
Copying nodata values from source test_000.tif to destination<br>
test_4326_2.tif.<br>
...10...20...30...40...50...60...70...80...90...100 - done.<br>
<br>
<br>
-Jukka-<br>
<br>
<br>
jratike80 wrote<br>
> Hi,<br>
> <br>
> What you have wrong is the input coordinate system, not the output as you<br>
> tested.<br>
<br>
<br>
<br>
<br>
<br>
--<br>
Sent from: <a href="http://osgeo-org.1560.x6.nabble.com/GDAL-Dev-f3742093.html" rel="noreferrer" target="_blank">http://osgeo-org.1560.x6.nabble.com/GDAL-Dev-f3742093.html</a><br>
_______________________________________________<br>
gdal-dev mailing list<br>
<a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/gdal-dev" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br>
</blockquote></div>