R：如何通过SpatialPolygonsDataFrame裁剪和屏蔽栅格列表？

Question

我想通过SpatialPolygonsDataFrame裁剪和屏蔽栅格列表。 spdf是GADM的一个国家，这里有一个伊朗的例子，我将如何仅为一个栅格做到这一点：

library(raster)

iran <- getData("GADM", country = "IRN", level = 1)
r <- raster(iran, nrow=100, ncol=100)
set.seed(-99)
values(r) <- runif(ncell(r))
r2 <- crop(r, iran)
res <- mask(r2, iran)

但是，我有很多栅格（4个不同国家的66个栅格），我不想复制粘贴我的代码。任何想法如何解决这个问题？谢谢