Hi, There is also a tool called Starspan than uses GDAL and GEOS to do arbitrary zonal stats.
Dylan On Wed, Aug 26, 2009 at 7:13 AM, Jose Gomez-Dans<jgomezd...@gmail.com> wrote: > Hi, > > 2009/8/20 questions anon <questions.a...@gmail.com> >> >> 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 > > i've done a similar example to what you want here: > < http://sites.google.com/site/spatialpython/aggregating-data-to-grid-cells> > Only difference is that I calculate the histogram, rather than just the mean > and std dev. It is trivial to modify this for your own purposes. In fact, > when I get some time, I'll throw in a full "zonal stats" class, but for now > this should get you started. > > J > > _______________________________________________ > gdal-dev mailing list > gdal-dev@lists.osgeo.org > http://lists.osgeo.org/mailman/listinfo/gdal-dev > _______________________________________________ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev