如何使用空间polygonsdataframe对象掩盖栅格？

Question

我想使用“空间polygonsdataframe”对象中包含的所有多边形掩盖栅格。 到目前为止，我尝试了此操作：

r=raster("myraster.asc")

然后我得到了带有6个多边形列表的对象（称为“ vor_spdf”）对象，所以我运行了：

typeof(vor_spdf) # It results in [1] "S4"
r1=crop(r,extent(vor_spdf))
r2=mask(r1,vor_spdf@polygons[[1]]) #Trying to mask the raster on the geometry of the first polygon

这是错误... 错误... ）：无法找到用于函数'掩码'签名'rasterlayer的函数'''rastlayer'的方法，“多边形”'

您是否有有关如何处理 这？ 谢谢

Answer 1

示例数据

library(raster)
library(terra)
p <- shapefile(system.file("external/lux.shp", package="raster"))
r <- rast(system.file("ex/elev.tif", package="terra")) |> raster()

解决方案

m <- mask(r, p[1,])

---

，但您正在使用过时的软件包。您可能会使用“ Terra”更好。

library(terra)
v <- vect(system.file("ex/lux.shp", package="terra"))
r <- rast(system.file("ex/elev.tif", package="terra"))

m <- mask(r, v[1,])