[gdal-dev] converting from EPSG:3857 to EPSG:25832 , gdaltransform works, gdalwarp not ?

Javier Jimenez Shaw j1 at jimenezshaw.com
Sun Jul 27 02:09:25 PDT 2025


w_info.txt says epsg:4326, not epsg:25832. Coordinates are wrong anyhow.

On Sun, 27 Jul 2025, 10:57 Peter Bennewitz via gdal-dev, <
gdal-dev at lists.osgeo.org> wrote:

> Dear all,
>
> Overall, I'm trying to convert ESRI SHAPE data, given in "Web Mercator"
> EPSG:3857 , to a TIF in   ETRS89 / UTM Zone 32U , EPSG:25832 .
>
> The GDAL Version is the one included in current stable Debian,
> https://packages.debian.org/bookworm/gdal-bin ,  gdal-bin (3.6.2+dfsg-1
> and others)
>
> The first step, rasterization, works flawless:
>     gdal_rasterize -a_srs EPSG:3857 -tr 10 10 -burn 1 q.shp q.tif
>
> gdalinfo lists the file as: (full text attached)
>
> ....
> Coordinate System is:
> PROJCRS["WGS 84 / Pseudo-Mercator",
>     BASEGEOGCRS["WGS 84",
>         ENSEMBLE["World Geodetic System 1984 ensemble",
>             MEMBER["World Geodetic System 1984 (Transit)"],
> ..
> Corner Coordinates:
> Upper Left  (  836370.121, 6080896.490) (  7d30'47.67"E, 47d50'37.43"N)
> Lower Left  (  836370.121, 6045996.490) (  7d30'47.67"E, 47d37'58.40"N)
> Upper Right (  887980.121, 6080896.490) (  7d58'36.70"E, 47d50'37.43"N)
> Lower Right (  887980.121, 6045996.490) (  7d58'36.70"E, 47d37'58.40"N)
> Center      (  862175.121, 6063446.490) (  7d44'42.18"E, 47d44'18.30"N)
> Band 1 Block=5161x1 Type=Float64, ColorInterp=Gray
>
> Running the lower left corner through gdaltransform gives plausible ETRS89
> coordinates :
>
> echo  836375.120800 6045997.312000 | gdaltransform -s_srs EPSG:3857 -t_srs
> EPSG:25832
> 388314.422838145 5276569.91979615 0
>
> However, running the TIF through gdalwarp, and checking the result with
> gdalinfo, does not show a plausible bounding box:
>
> gdalwarp  -s_srs EPSG:3857 -t_srs EPSG:25832 q.tif w.tif
> gdalinfo w.tif
> ...
> Corner Coordinates:
> Upper Left  (  836370.121, 6080896.490) (Invalid angle,Invalid angle)
> Lower Left  (  836370.121, 6045996.490) (Invalid angle,Invalid angle)
> Upper Right (  887980.121, 6080896.490) (Invalid angle,Invalid angle)
> Lower Right (  887980.121, 6045996.490) (Invalid angle,Invalid angle)
> Center      (  862175.121, 6063446.490) (Invalid angle,Invalid angle)
>
> Looks like it passes gdalwarp unchanged.
>
> Apparently, I'm overlooking something obvious.
> Any pointers much appreciated,
> thanks, best regards, have a nice Sunday and thanks for GDAL !
> Peter
> _______________________________________________
> 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/20250727/ef550419/attachment.htm>


More information about the gdal-dev mailing list