我想将CO(多边形)分成不被道路(线串)分割的部分(也包括多边形)。也就是说,我希望较小的多边形的部分以道路或州界为边界,而不要包含任何进入和退出多边形的道路。
我能够使用lwgeom::st_split
来生成几何图形集合,但是我不确定这是否对我有帮助;我坚持使用此解决方案,因为我不确定如何提取集合中的几何图形,例如,为其分配唯一的ID。
我的最终目标是确保我的点(单独的数据)不被道路分开。因此,如果您有一个可能更直接的解决方案,那么我也非常高兴。
library(tidyverse)
library(tigris)
library(sf)
library(lwgeom)
co <- states(cb = T) %>%
filter(NAME == "Colorado")
roads <- primary_secondary_roads(state = 'Colorado')
cosplit <- st_split(co,roads)
是否有人找到或看到了解决方案?
答案 0 :(得分:1)
我想我已经知道了...但是我绝对很想听听其他人的想法!!!
cosplitpoly <- cosplit %>%
st_collection_extract(c("POLYGON"))