gdal_translate KMLSUPEROVERLAY crea una imagen negra en blanco

Question

Estoy tratando de convertir un archivo geotiff a kmz para verlo en GE. El geotiff tiene 878.720 KB. Utilizando un "kml normal", la resolución es muy baja en GE. Después de haber leído, parece que la conversión a una superposición es el mejor método.

He optado por utilizar gdal_translate por encima de otros métodos (por ejemplo, MapTiler o la herramienta de Arc que cuesta) porque ya he estado utilizando la biblioteca de GDAL python y estoy tratando de automatizar todo el proceso.

Para probar, escribo en UNIX:

gdal_translate -of KMLSUPEROVERLAY filepath/2010_full_mosaic.tif filepath/2010_full_mosaic.kmz

Se crea un archivo kml....pero al cargarlo en GE aparece en blanco (negro) - aunque la extensión/proyección parecen correctas.

¿Alguien podría aconsejarme qué he hecho mal? Al utilizar el comando gdalinfo, se me proporcionan los siguientes datos:

Driver: GTiff/GeoTIFF
Files: Filepath/2010_full_mosaic.tif
       Filepath/2010_full_mosaic.tif.ovr
       Filepath/2010_full_mosaic.tif.aux.xml
Size is 10572, 22630
Coordinate System is:
PROJCS["WGS 84 / UTM zone 36S",
    GEOGCS["WGS 84",
         DATUM["WGS_1984",
            SPHEROID["WGS 84",6378137,298.257223563,
                AUTHORITY["EPSG","7030"]],
            AUTHORITY["EPSG","6326"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4326"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",33],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",10000000],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AUTHORITY["EPSG","32736"]]
Origin = (390313.896406000014395,8164584.107470000162721)
Pixel Size = (25.000000000000000,-25.000624509551976)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (  390313.896, 8164584.107) ( 31d58'18.28"E, 16d35'54.43"S)
Lower Left  (  390313.896, 7598819.975) ( 31d56'22.37"E, 21d42'38.47"S)
Upper Right (  654613.896, 8164584.107) ( 34d26'57.65"E, 16d35'45.40"S)
Lower Right (  654613.896, 7598819.975) ( 34d29'41.00"E, 21d42'26.42"S)
Center      (  522463.896, 7881702.041) ( 33d12'49.09"E, 19d 9'27.96"S)
Band 1 Block=10572x1 Type=Float32, ColorInterp=Gray
 Overviews: 5286x11315, 2643x5658, 1322x2829, 661x1415, 331x708, 166x354

Becky

Answer 1

Dependiendo del tipo de archivo tif de entrada, puede ser necesario añadir

-expand rgb

a la línea gdal_translate. Encontré que esto era necesario al convertir gráficos BSB que comienzan con una imagen indexada. El KMLSUPEROVERLAY (al menos en mi caso) estaba escribiendo imágenes JPG para cada mosaico que no no soportan mapas de colores indexados.