[gdal-dev] calculate raster values within a vector

[questions anon]

Can anyone get me started on a script that can -

 -Open a shapefile

 -Open many rasters

 -Calculate the mean and standard deviation of those rasters within the
shapefile region, but is able to ignore areas of zeros or NoData

 -Output the mean and standard deviation to a table

 Otherwise any tutorials or books that I could follow to work this out would
be greatly appreciated! Thanks!
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/gdal-dev/attachments/20090820/a5275a81/attachment.html