r - 如何在 R 中基于 SpatialPolygons 裁剪栅格

Question

我想根据 SpatialPolygons 裁剪栅格目的。我知道我们可以使用 crop函数在 raster包裹，

raster::crop(rasterFile, SpatialPolygonsObject)

但此功能基于 extent的 SpatialPolygons对象，所以裁剪的结果是矩形。但是，在某些情况下， SpatialPolygons对象不是矩形，如何处理这些情况？

Answer 1

您可以使用 raster::mask .这是一个可重现的示例:

library(raster)
r = raster(vals = rnorm(400), nrows=20, ncols=20, ext= extent(c(0, 20, 0, 20)))
p = Polygon(matrix(5, 5, 15, 12, 7, 16, 3, 10), ncol=2, byrow = T))
p = SpatialPolygons(list(Polygons(list(p), "p")))

plot(r)
lines(p)

r2 = mask(r,p)
plot(r2)

如果您还需要裁剪栅格的范围以移除掩码周围的空行和列，那么您可以使用 crop申请前 mask ，或者您可以使用 trim(r2, values = NA)然后。