[gdal-dev] Temporal statistics

[Knut-Frode Dagestad]

Hi,

 From a time series of colocated images, I construct a VRT with one band
for each time step using gdalbuildvrt:

 > gdalbuildvrt -separate mtsat.vrt mtsat/*_IR*00.tif

Are there any tools that can be used to calculate some statistics (min,
max, mean, etc) versus time for each pixel of such a 3D Dataset?
The rasters are large and timeseries long, so reading everything into a
Python NumPy cube is not a good solution.

Thanks,
Knut-Frode