¿Cómo importar correctamente un archivo raster con estos metadatos en Tilemill?

Question

Tengo un archivo GeoTIFF para el que gdalsrsinfo informa de esto:

PROJ.4 : '+proj=utm +zone=32 +datum=WGS84 +units=m +no_defs '

OGC WKT :
PROJCS["WGS 84 / UTM zone 32N",
GEOGCS["WGS 84",
    DATUM["WGS_1984",
        SPHEROID["WGS 84",6378137,298.257223563,
            AUTHORITY["EPSG","7030"]],
        AUTHORITY["EPSG","6326"]],
    PRIMEM["Greenwich",0,
        AUTHORITY["EPSG","8901"]],
    UNIT["degree",0.0174532925199433,
        AUTHORITY["EPSG","9122"]],
    AUTHORITY["EPSG","4326"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",9],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",0],
UNIT["metre",1,
    AUTHORITY["EPSG","9001"]],
AXIS["Easting",EAST],
AXIS["Northing",NORTH],
AUTHORITY["EPSG","32632"]]

Si intento importar ese archivo a Tilemill, no puedo verlo correctamente. Si elijo dejar que Tilemill autodetecte el SRS, falla:

Unable to determine SRS for layer "antani2" at /home/strappanai/antani2.tiff

Si en Tilemill especifico en su lugar una cadena PROJ.4 personalizada, es decir, la que obtengo de gdalsrsinfo (y también de QGIS), esto es

+proj=utm +zone=32 +datum=WGS84 +units=m +no_defs

entonces la capa realmente consigue importado, pero se muestra significativamente inclinado hacia la derecha - no como se esperaba.

No sé si esto puede importar en alguna medida, pero QGIS importa y visualiza el archivo perfectamente.

También intenté reproyectar el archivo fuente como se indica aquí: https://www.mapbox.com/tilemill/docs/guides/reprojecting-geotiff/

Pero gdalwarp devuelve este error mil veces y luego falla:

ERROR 1: Too many points (441 out of 441) failed to transform,
unable to compute output bounds.
Warning 1: Unable to compute source region for output window     3593619,2225230,6394,6394, skipping.
More than 1000 errors or warnings have been reported. No more will be reported from now.

EDITAR: Esta es la extensión de la capa dada por QGIS.

Extensión de la capa (proyección de la fuente original de la capa) 638153.2011202700668946,5115644.8502587005496025 : 662908.4373886700486764,5132649.3878203006461263

Aquí hay un archivo de comportamiento idéntico: archivo de muestra

Y así es como se ve en QGIS y en Tilemill (QGIS lee el SRS del archivo, yo especifico la cadena PROJ.4 en Tilemill como se indica arriba ya que no puede leerla del TIFF).

Answer 1

He probado un poco con tu archivo de ejemplo. La solución fue la siguiente. Reproyecté el archivo a WGS84 usando QGIS y luego lo añadí a mi proyecto Tilemill definiendo la proyección como WGS84.

El comando gdal es el siguiente:

gdalwarp -overwrite -s_srs EPSG:32632 -t_srs EPSG:4326 -of GTiff antanired.tiff antanired_warp_WGS84.tiff

Como puede ver en el ejemplo, funciona con los datos proporcionados. Pero falló dejando la configuración en auto detectar srs.

La inclinación de la imagen hacia la derecha sólo se debe a las diferencias de srs. Aunque la imagen es perfectamente rectangular en UTM, estará un poco apretada cuando se muestre en WGS84. Lo mismo ocurrirá si lo haces al revés.

ACTUALIZACIÓN

Puedes deshacerte del fondo negro pasando

nodata="0"

en el campo Avanzado durante la carga de la trama. Como se muestra para los datos Landsat aquí

Para darle la respuesta final a su pregunta: ¡La importación es absolutamente correcta! Parte de su error fue la falta de comprensión de los archivos raster y las proyecciones. Espero que el truco de la definición de los nódulos funcione.