[gdal-dev] gdalwarp from EPSG:25832 to EPSG:4326 with error latitude or longitude exceeded limits?

[Pham Huu Bang]

Thank to Jukka,

so the problem was with my gdal 1.11.4 on Centos 7.
When I translated it to EPSG:4647, it doesn't contain proper WKT for EPSG,
e.g:

Coordinate System is:
LOCAL_CS["ETRS89 / UTM zone 32N (zE-N)",
    GEOGCS["ETRS89",
        DATUM["unknown",
            SPHEROID["unretrievable - using WGS84",6378137,298.257223563]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433]],
    AUTHORITY["EPSG","4647"],
    UNIT["metre",1]]

compared to output from GDAL 2.2.4 on Ubuntu 18.04:

Coordinate System is:
PROJCS["ETRS89 / UTM zone 32N (zE-N)",
    GEOGCS["ETRS89",
        DATUM["European_Terrestrial_Reference_System_1989",
            SPHEROID["GRS 1980",6378137,298.257222101,
                AUTHORITY["EPSG","7019"]],
            TOWGS84[0,0,0,0,0,0,0],
            AUTHORITY["EPSG","6258"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.0174532925199433,
            AUTHORITY["EPSG","9122"]],
        AUTHORITY["EPSG","4258"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",9],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",32500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AXIS["Easting",EAST],
    AXIS["Northing",NORTH],
    AUTHORITY["EPSG","4647"]]

With gdal 1.11.4, I cannot warp it to EPSG:4326 with error:

ERROR 1: No PROJ.4 translation for source SRS, coordinate
transformation initialization has failed.

and it works with gdal 2.2.4 for gdalwarp.



On Wed, 25 Nov 2020 at 12:05, jratike80 <jukka.rahkonen at maanmittauslaitos.fi>
wrote:

> Hi again,
>
> Sorry, I did not read properly what you did. That should work and it does
> work for me with GDAL 3.2. I use a copy of your test.tif as "test_00.tif".
>
>
> gdal_translate -a_srs epsg:4647 test_00.tif test_000.tif
> Input file size is 100, 100
> 0...10...20...30...40...50...60...70...80...90...100 - done.
>
> gdalwarp -t_srs epsg:4326 test_000.tif test_4326_2.tif
> Creating output file that is 125P x 73L.
> Processing test_000.tif [1/1] : 0Using internal nodata values (e.g. -9999)
> for image test_000.tif.
> Copying nodata values from source test_000.tif to destination
> test_4326_2.tif.
> ...10...20...30...40...50...60...70...80...90...100 - done.
>
>
> -Jukka-
>
>
> jratike80 wrote
> > Hi,
> >
> > What you have wrong is the input coordinate system, not the output as you
> > tested.
>
>
>
>
>
> --
> Sent from: http://osgeo-org.1560.x6.nabble.com/GDAL-Dev-f3742093.html
> _______________________________________________
> gdal-dev mailing list
> gdal-dev at lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/gdal-dev
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20201125/3b368dc8/attachment-0001.html>