选择与另一个多边形相交（或完全在另一个多边形内）的多边形[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)