I'm trying to average 30 greyscale netCDF files pixel-by-pixel using
gdal_calc.py. However the script seems to require the input files to be enumerated with letters A-Z, so I'm four letters short.
The docs for
gdal_calc lead me to believe I can instead supply a single raster file with all the bands and run an expression on that, but I haven't been able to find the correct syntax.
The closest I've come has been with this command ...
gdal_calc.py -A stacked.vrt --outfile=mean.nc --calc="mean(A, axis=0)"
... which produced a weird blocky-looking result unfortunately:
Has anyone done this before? Ideally I'm looking for a solution that would also work if it were more than 30 input files/bands.