R script used to calculate average of habitat suitability pixels under each buffer zone. This document contains the annotated code used to calculate the average pixel value under each buffer zone (overlapping areas act as a separate zone with their own averaged value) using R.