rast.by.polys function

Summarise raster values by polygons