¿Fusión de raster en QGIS?

Question

Tengo un gran problema que no puedo resolver ni siquiera después de probar todo lo que se me ocurre.

Estoy intentando hacer un análisis de la cuenca del lago Chad en QGIS (2.16.3). He descargado los datos de Landsat 8 desde Earth Explorer, pero necesito 4 imágenes de satélite diferentes para cubrir toda la cuenca. Hice la corrección atmosférica y creé una pila de bandas, el problema que aparece ahora es que tengo un "borde negro" alrededor de las imágenes individuales.

Traté de arreglar esto (como se lee en otro post aquí) yendo a las propiedades de la capa -> transparencia -> adicional sin valor de datos = 0, pero no funcionó (también traté de raster -> análisis -> cerca de rojo), así que en su lugar tomé la imagen individual, e hice raster -> miscelánea -> fusionar -> y marcó "sin valor de datos". Esto funcionó para eliminar los bordes negros, pero ahora cuando intento fusionar las cuatro imágenes individuales en una sola, los bordes negros aparecen de nuevo.

¿Alguien sabe cómo arreglar esto?

Llevo días intentando hacerlo, pero realmente no sé cómo ayudarme a estas alturas. Voy a añadir algunas imágenes para aclarar lo que estoy hablando.

Answer 1

Si usted está estableciendo el valor NoData en las propiedades de la capa / transparencia, creo que es sólo para fines de visualización. Además, significa que sus datos rasterizados no tienen el valor NoData establecido. Se me ocurren dos soluciones:

1. Dentro de QGIS:
   Utilice r.null de la caja de herramientas de procesamiento para establecer correctamente su valor NoData. A continuación, realice la fusión, respetando el NoData

2. Desde la línea de comandos:
   Utilice gdal_merge.py desde la línea de comandos para que puedas controlar mejor tus entradas. Puede utilizar la función -n para indicar cuál es el valor de NoData. Su comando se vería algo parecido a esto:
   gdal_merge.py -o merged_all.tif -n 0 RT_LC....your list of inputs here...

Answer 2

La forma en que funcionó ahora para mí fue: Crear un nuevo shapefile (polígono) - alternar edición - añadir característica (hacer un cuadrado en la imagen del satélite, dejando los bordes negros fuera) - seguro. Luego ir a raster - extracción - clipper, poner el raster original que quería recortar como entrada, el nuevo shapefile como salida - ejecutar esto. ¡¡¡En la imagen recortada, ponga ningún valor de datos adicionales a 0; luego al ejecutar el misc - merge de las cuatro imágenes (de nuevo, poniendo ningún valor de datos a 0 en la ventana de merge) y entonces funcionó!!! Además, traté de usar r.null, pero de alguna manera siempre obtuvo un mensaje de error. De todos modos, ahora me ha funcionado.