[gdal-dev] How to overlay a raster with shapes to calc statistics?

[Andreas H.]

Hi all,

I have a GeoTiff file covering the whole globe at approx. 1km
resolution. And I have a shapefile containing say 20 of polygons.

For each of the 20 polygons, I want to know the mean value of the raster
data falling inside the polygon, and possibly more statistics (max, min,
medeian, std, ...).

How can I achieve this using GDAL?

Thanks for your insight,
Andreas.