Select 个多边形与另一个多边形相交(或完全在其中)[R sf]

Select polygons intersecting (or completely within) another polygon [R sf]

我想 select 来自一个简单要素集合的多边形与另一个 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)

现在,创建了多边形及其周围的鱼网后,我想要 select 那些与多边形相交或完全被多边形包含的鱼网单元。

感谢您的帮助!

马克

fishnet2 <- fishnet[poly,]

还有剧情

library(ggplot2)

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