[gdal-dev] Zonal statistics with GDAL/StarSpan
bfraser at geoanalytic.com
Thu Dec 18 13:59:48 EST 2008
I've used StarSpan to do similar things; in my case I wanted "mode" not "avg". Using v1.2.01 (I don't think I made major mods to the source):
starspan2.exe --progress --RID none --vector boundaries.shp --raster big.tif --stats results.txt mode
I edited the resulting results.txt text (csv) file and used GDAL's ogr2ogr to do a join on the shapefile's DBF file based on the FID.
Greg Ederer wrote:
> Hello all,
> I'm a GDAL newbie. I have a shape file containing subnational
> boundaries for Africa. I have a raster containing rainfall for Africa.
> I need to pull out the spatial average of rainfall for each of the
> subnational units.
> I think that I might be able to pull this off using StarSpan. But, the
> docs on the StarSpan wiki do not line up with the current version of the
> program (the changelog mentions that certain command line switches have
> been eliminated, and others added; but it doesn't say what the new
> switches do, or how to use them). I am digging through the StarSpan
> sources to try to figure out how to use it.
> Does anyone have any suggestions as to how I might get the needed data
> using GDAL (with, or without StarSpan)?
More information about the gdal-dev