r - 选择与另一个多边形相交(或完全在其中)的多边形 [R sf]

Question

我想从一个与另一个 sf 集合中的多边形相交的简单特征集合中选择多边形。这是一个可重现的示例:

library(sf)
nc <- st_read(system.file("shape/nc.shp", package="sf"))
nc <- st_transform(nc,32717);nc
poly <- nc[1,]
fishnet<-st_make_grid(st_as_sfc(st_bbox(poly)),cellsize=5000)

现在，创建了 poly 和围绕它的渔网，我想选择那些与 poly 相交或完全被 poly 包含的渔网单元。

谢谢你的帮助!

标记

Answer 1

fishnet2 <- fishnet[poly,]

还有一个情节

library(ggplot2)

ggplot() + 
  geom_sf(data = poly, fill = "blue") +
  geom_sf(data = fishnet2, color = "red", fill = NA)