Below you’d find a simple R code that calculates the average value of a raster within polygons: