Sampling based on distribution provided by Raster* object (does not have to be summed to one) over an sf object used as mask
brk_sampling(rasterStack, sf, sizeSite = 1)
Raster* object (RasterLayer or RasterStack)
An sf object
scalar or vector of the number of individual to sample per site. Default is 1. If vector, as to be the length of Raster* object