如何基于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，或者可以使用{{1}之后。