[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