[gdal-dev] Temporal statistics

Etienne Tourigny etourigny.dev at gmail.com
Fri Jul 13 07:19:09 PDT 2012

Provided you can convert the format to netcdf, you could use a tool
like cdo which does exactly what you want

"CDO is a collection of command line Operators to manipulate and
analyse Climate and NWP model Data."

gdal supports netcdf output, you might have to do some fixing up of
the time axes though.

Basically I would do

- loop on each file:
 - convert to netcdf
 - set the time axis with cdo settime <def>
- merge all files: cdo mergetime <infiles> outfile.nc
- operate on the resulting file with cdo operators (cdo timmean, cdo
timmax, etc)


On Fri, Jul 13, 2012 at 11:01 AM, Knut-Frode Dagestad
<knutfrodesoppel at hotmail.com> wrote:
> 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
> _______________________________________________
> gdal-dev mailing list
> gdal-dev at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/gdal-dev

More information about the gdal-dev mailing list