问题描述
我有两个 shp 文件。
我只想知道那些完全 (100%) 落入 shpfile_2 的 shpfile_1。
library(sf)
shp_1 <- st_read(system.file("shape/nc.shp",package="sf"))
shp_1 <- st_transform(shp_1,2154);shp_1
shp_1 <- st_union(shp_1)
nc <- st_read(system.file("shape/nc.shp",package="sf"))
nc <- st_transform(nc,2154)
shp_2 <- nc[1,]
shp_2 <-st_make_grid(st_as_sfc(st_bbox(shp_2)),cellsize=5000)
within <- st_within(shp_2,shp_1)
问题是 st_within 为 600 万个多边形花费了太多时间。
有什么办法可以加快这个过程吗？
解决方法
暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！
如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@）