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

Question

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

raster::crop(rasterFile, SpatialPolygonsObject)

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