¿Intercalar geotiff con otro utilizando el canal alfa (o el valor de los nódulos) en GDAL?

Question

Tengo dos geotiffs relacionados con una ubicación pero con diferentes áreas cubiertas:

1.

2.

¿Cómo puedo intersecar sus partes opacas en un solo geotiff así?

Answer 1

Nota: Esta solución sólo funciona para la imagen de una sola banda

Si tiene varias bandas, puede repetir esto para cada banda añadiendo -A_Band y -B_Band opciones y luego gdal_merge.py -separate las bandas.

---

Deberías ser capaz de hacerlo:

(Asumimos que el valor de NoData de sus rasters es 0 )

gdal_calc.py -A file1.tif -B file2.tif --NoDataValue=0 \
--calc="A*(A>0)+B*logical_and(A==0,B>0)" \
--outfile=outfile.tif

Explicación:

• A*(A>0) Si A no es Nodata, será A*1 , si no A*0
• B*logical_and(A==0,B>0) Si A es NoData y B no es Nodata, será B*1 , si no B*0 .

(Esto siempre preferirá A sobre B).