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

现在，我创建了多边形和周围的鱼网，我想选择那些与多边形相交的鱼网单元或完全被多边形包含的单元。

谢谢你的帮助

标记

Answer 1

fishnet2 <- fishnet[poly,]

还有一个情节

library(ggplot2)

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