On 17 July 2015 at 11:58, Even Rouault <[email protected]> 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 _______________________________________________ gdal-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/gdal-dev
