Gdalwarp no utiliza el re muestreo de lanczos en la imagen con canal alfa

Question

Quiero envolver la imagen Geotiff de relieve sombreado (RGBA) en la proyección EPSG:8353 a la proyección EPSG:3857 para usarla en un mapa web en mosaico de Mapnik. Quiero usar un re muestreo lanczos, pero la imagen resultante se ve "dentada", como si se hubiera utilizado un re muestreo de vecinos más cercanos. Estoy usando este comando:

gdalwarp \
  -of GTiff -overwrite -r lanczos  \
  -tr 1.194328566968441 1.194328566968441 -tap \
  -t_srs EPSG:3857 shading-sjtsk.tif shading-web.tif

La resolución objetivo en unidades georreferenciadas se establece para ser la misma que el nivel de zoom 17 en la proyección mercator web del mapa en mosaico.

Por favor, ¿podrías aconsejarme por qué no se está utilizando el re muestreo lanczos? También probé otros métodos (como cúbico), pero la imagen siempre parece ser indistinguible de la más cercana.

Imagen de entrada:

Imagen de salida:

Imágenes tif originales se comparten aquí.

Answer 1

El problema estaba con la banda alfa, la cual no se deformó correctamente. Necesitamos evitar que la banda alfa de la imagen fuente se considere como tal con el parámetro -nosrcalpha. Luego se deformará como una banda regular. El siguiente comando produce una imagen suave y bonita:

gdalwarp \
  -of GTiff -overwrite -r lanczos -nosrcalpha \
  -tr 1.194328566968441 1.194328566968441 -tap \
  -t_srs EPSG:3857 shading-sjtsk.tif shading-web.tif

Imagen de resultado: