[Qgis-user] Qgis : Calculat the median of multiple bands of raster (90 bands)

[Karlo no Dabas maajas]

Hi.

Here comes a bit more polished script version. The results are joined 
into a multiband GeoTiff file.

________________________________
INPUT=MultidimensionGrid
INPUT_EXT=tif
TMP_SAGA=MultiDim
FOCAL_CMD=MEDIAN

for NR in {0..5} ; do
	echo -ne "\033]0;Working with Band $NR\007"
	saga_cmd io_gdal 0 \
		-GRIDS $TMP_SAGA$NR -MULTIPLE 0 \
		-SELECTION  $NR \
		-FILES $INPUT.$INPUT_EXT
	saga_cmd statistics_grid 1 \
		-GRID $TMP_SAGA$NR.sgrd  \
		-$FOCAL_CMD $FOCAL_CMD$NR
	rm $TMP_SAGA$NR.*
	JOIN=$JOIN$FOCAL_CMD$NR'.sgrd'\;
done

JOIN=${JOIN%;}
saga_cmd io_gdal 1 \
	-GRIDS $JOIN \
	-MULTIPLE 1 -FILE $INPUT'_'$FOCAL_CMD.tiff \
	-FORMAT GeoTIFF \
	-OPTIONS COMPRESS=DEFLATE PREDICTOR=2

rm $FOCAL_CMD*

________________________________

B.r.
K.


01.01.26 15:50, Karlo no Dabas maajas via QGIS-User rakstīja:
> Hi again!
> 
> „Focal Statistics” works only with one layer. But it's easy to overcome 
> this with small script.
> 
> Look at this example. As the source, I have a GeoTIFF with 6 bands:
> 
> [..]
> 
> 
> So I create a script to calculate median for all bands:
> 
> [..]