if("fasterize" %in% rownames(installed.packages()) == FALSE) {install.packages("fasterize")} if("raster" %in% rownames(installed.packages()) == FALSE) {install.packages("raster")} clip_raster_with_sf <- function(rast, sf){ r_mask <- fasterize::fasterize(sf, rast) r_masked <- raster::mask(rast, r_mask) return(r_masked) }