[gdal-dev] Zonal statistics with GDAL/StarSpan

[Brent Fraser]

Greg,

  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.

Brent Fraser


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)?
> 
> Thanks!
> 
> Greg
>