如何在R中基于空间多边形裁剪栅格数据

Question

我想基于SpatialPolygons对象剪裁栅格图像。我知道我们可以使用raster软件包中的crop函数。

raster::crop(rasterFile, SpatialPolygonsObject)

但是这个函数基于SpatialPolygons对象的extent，因此裁剪结果是一个矩形。然而，在某些情况下，SpatialPolygons对象不是矩形，如何处理这些情况？

- just_rookie

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，或者之后使用trim（r2，values = NA）。

- dww

---

谢谢您的回复。但是如何删除那些不包含在掩码中的单元格呢？ - just_rookie

但是在进行mask操作之前，您可以先进行crop操作，以尽可能消除更多的细胞。 - Robert Hijmans

或者您可以使用 trim 函数（来自 raster 包）来删除具有 NA 值的外部行和列：trim(r2, values = NA)。 - Claudia