obtener imágenes rgb de sentinel-2 mediante gdal

Question

Tengo algunos problemas para obtener una imagen en colores reales RGB de las bandas del satélite Sentinel

que producen esto

a mí me parece que no es nada natural

también he probado gdal para hacer la fusión con :

gdal_merge -o merged.tif -of GTIFF B04.jp2 B03.jp2 B02.jp2

pero me dio un archivo inutilizable, al menos en microDEM .

Así que me pregunto qué es lo que está mal aquí?

Answer 1

Es necesario añadir el -separate con el fin de colocar cada archivo de entrada en una banda separada y (opcionalmente) el -co PHOTOMETRIC=RGB opción de creación para forzar la interpretación fotométrica (para evitar, por ejemplo, la ColorInterp=undefined y establecer la interpretación del color adecuado para cada banda):

gdal_merge -separate -co PHOTOMETRIC=RGB -o merged.tif B04.jp2 B03.jp2 B02.jp2

Nota: -of GTIFF no es necesario porque GeoTIFF es el formato de salida por defecto.

Otro enfoque consiste en escalar las bandas de entrada antes de gdal_merging (pero no es estrictamente necesario en los clientes, donde las bandas se estiran):

for %%i in (4,3,2) do gdal_translate -scale B0%%i.jp2 B0%%i.tif
gdal_merge -separate -co PHOTOMETRIC=RGB -o merged.tif B04.tif B03.tif B02.tif

Por último, la salida necesita algunas operaciones de ajuste fino, como el corte del histograma (por ejemplo, un corte de recuento acumulativo) y la calibración del brillo, el contraste y la saturación, con el fin de ser una representación RGB aceptable para utilizar en clientes GIS/RS o en cualquier otro lugar. A continuación se presenta una muestra rápida y sucia de Sentinel-2 RGB creada con gdal_merge y luego se visualiza en QGIS: