As a follow up to my idea at the Spatial Ecology summer school in Matera I’d like to point you to my github repo: zonalbenchmark. The gist of it: “What tool can I use to calculate the zonal statistics for my data set?”

I hope to ignite some thoughts about different tools for this task. If this little python script is furthermore useful for someone to justify a decision, I am more than happy.

Sample usage:

1
2
3
python zonalStatBenchmark [tools] [input raster] [input mask / shape] [number of runs]
python zonalStatBenchmark.py 1-2-3-4 test_data/wc2.0_10m_tavg_07.tif test_data/mask.shp 1

There’s still a lot to to do; e.g.:

  • add more tools for benchmarking
  • improve error handling
  • add support for ‘per-feature-benchmarking’ (load one feature from shape, run benchmark, load the next one, run benchmark, etc.)
  • add multilanguage capabilities
  • add test data to repo
  • write unit tests

If you’re interested, please do not hesitate to contact me or to open an issue on github.