¿Establecer la banda 4 del raster en alfa con GDAL?

Question

Tengo una trama con 4 bandas, pero la banda 4 no se interpreta como un canal alfa. Cuando ejecuto gdalinfo obtengo la siguiente información sobre las bandas:

Band 1 Block=360000x1 Type=Byte, ColorInterp=Gray
Band 2 Block=360000x1 Type=Byte, ColorInterp=Undefined
Band 3 Block=360000x1 Type=Byte, ColorInterp=Undefined
Band 4 Block=360000x1 Type=Byte, ColorInterp=Undefined

Me gustaría que tuviera este aspecto:

Band 1 Block=360000x1 Type=Byte, ColorInterp=Red
  Mask Flags: PER_DATASET ALPHA
Band 2 Block=360000x1 Type=Byte, ColorInterp=Green
  Mask Flags: PER_DATASET ALPHA
Band 3 Block=360000x1 Type=Byte, ColorInterp=Blue
 Mask Flags: PER_DATASET ALPHA
Band 4 Block=360000x1 Type=Byte, ColorInterp=Alpha

Soy capaz de hacer esto en QGIS mediante el establecimiento de la banda 4 a transparente y, a continuación, la representación de la trama.

¿Hay alguna forma de cambiar el "ColorInterp" en GDAL ya que quiero automatizar el proceso?

Answer 1

Debería ser posible con gdal_edit.py http://www.gdal.org/gdal_edit.html

-colorinterp_X rojo|verde|azul|alfa|indefinido:

(GDAL >= 2.3) Cambiar la interpretación del color de la banda X (donde X es un número de banda válido, empezando por 1)

El comando completo sería el siguiente

gdal_edit -colorinterp_1 red -colorinterp_2 green -colorinterp_3 blue -colorinterp_4 alpha image.tif