crear una trama máxima con gdal_calc

Question

Tengo una serie NDVI-Time (24 tiffs de un año). Lo que me gustaría calcular es un único Raster con el valor máximo de NDVI para cada píxel de las 24 imágenes. He pensado en utilizar gdalcalc pero no encuentro la sintaxis adecuada para calcular el máximo de un mis archivos de entrada:

Ejemplo para 2 capas:

shell("E:/Programme/OSGeo4W.bat; gdal_calc -A E:/IVFL/Somalia/NDVI/201401_001_NDVI_NAfill.tif -B E:/IVFL/Somalia/NDVI/201409_121_NDVI_NAfill.tif --outfile E:/IVFL/Somalia/NDVI_maximum.tif --calc=(max(A,B))")

calc=max(A,B) siempre me da un error.

Y al final incluso tengo 24 archivos de entrada y no sólo 2. Alguien sabe la sintaxis (numpy) para calcular el máximo?

Answer 1

Debe utilizar máximo en lugar de max

--calc="maximum(A,B)"

Answer 2

Puede combinar varias funciones de máximo(A,B) por pares en una función larga.

--calc="maximum(maximum(maximum(A,B),maximum(C,D)),maximum(maximum(E,F),maximum(G,H)))"