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

Question

我想根据 SpatialPolygons 对象裁剪光栅。我知道我们可以在 raster 包中使用 crop 函数，

raster::crop(rasterFile, SpatialPolygonsObject)

但是这个函数是基于SpatialPolygons对象的extent，所以裁剪后的结果是矩形。但是，在某些情况下，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)

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