Change projection from Unknown to given EPSG with GDAL

Question

I've got a TIFF in projection (proj4: +proj=stere +lat_0=90 +lat_ts=60 +lon_0=-80 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs). Would like to convert it to EPSG:3857. ArcMap does it correctly, but GDAL does not. I tried gdalwarp:

gdalwarp -s_srs "+proj=stere +lat_0=90 +lat_ts=60 +lon_0=-80 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs" -t_srs EPSG:3857 -co "COMPRESS=LZW"

Gdalwarp generates more or less empty output.

So I tried gdal_translate:

gdal_translate -projwin_srs "+proj=stere +lat_0=90 +lat_ts=60 +lon_0=-80 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs" -a_srs EPSG:3857 -co "COMPRESS=LZW"

It gives some result, but it's wrong. Both results are the same when doing them in Python.

ArcMap doesn't have a problem with that. Any ideas how to make it automatically with Python (not using ArcPy!)?

The results and input raster are here (around 1-3MB each): https://drive.google.com/drive/folders/1NyXjYqSJI8KUSgA4AwTQbfLHCQhxILPG?usp=sharing

Answer 1

Seems to work fine for me with a simple:

gdalwarp -t_srs EPSG:3857 -r average -of GTiff input.tif out.tif

Your file already has a CRS set so no need to override it.

GDAL 3.8.1, released 2023/11/28