[gdal-dev] Gdal_calc.py with more than 26 input files

[Mike Toews]

On 17 July 2015 at 11:58, Even Rouault <even.rouault at spatialys.com> wrote:
> In your above example, only the first band of var_stacked.vrt will be used. So
> basically you're computing the average of all lines of a block (possibly a
> single line depending on the dataset), and this average line will be written
> in the first line of each block.
>
> Perhaps there should be a check on the shape of the result to verify it
> matches the shape of the input blocks and error out if not, if that's what you
> meant ?
> Unless there would be some new option to mean that a multiband raster should
> be loaded as a 3d numpy array. The current behaviour of using a single band is
> a documented one, so we cannot change that without breaking existing valid
> uses.

Thanks, that's clear to me now. I see that calculations are currently
done per band, and that 3D array support is not implemented. I've
started an enhancement ticket for a -3D option:
https://trac.osgeo.org/gdal/ticket/6040

-Mike