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@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@lists.osgeo.org > https://lists.osgeo.org/mailman/listinfo/gdal-dev >
_______________________________________________ gdal-dev mailing list gdal-dev@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/gdal-dev
