时间戳

首页 > 解决方案 > 如何使用 sf::st_intersects 从其邻居列表中过滤掉每个多边形

r - 如何使用 sf::st_intersects 从其邻居列表中过滤掉每个多边形

问题描述

我需要为sf数据集中的每个多边形提取相邻多边形。

这是一个简单的例子：

library(tidyverse)
library(sf)

demo(nc, ask = FALSE, verbose = FALSE)
nc <- nc %>% 
  mutate(polygon_id = row_number())

我设法用sf::st_intersects

neighbors <- st_intersects(nc, nc)
neighbors[[5]]
[1]  5  6  9 16 28

问题是每个多边形（此处为 5）都包含在邻居列表中。只使用一个nc数据集给我同样的结果

neighbors <- st_intersects(nc)
neighbors[[5]]
[1]  5  6  9 16 28

关于如何从相邻/相邻多边形列表中过滤掉实际多边形的任何提示？

标签： rpolygonsf

解决方案

好问题。这个问题可以有很多解决方案。但是这个问题的简单答案："Any tips on how to filter out the actual polygon from the list of adjacent/neighboring polygons?"，是使用带有 R 内核的 Jupyter Lab IDE 完成的。以下代码提供了一种回答问题的方法。

nc 数据集中有 100 个县。此代码以颜色显示所选县并显示所有相邻县。此代码适用于北卡罗来纳州 100 个县中的任何一个县。百强县选在这里。

代码：

nc1 <- nc %>% mutate(c_id = 1:nrow(nc))        
n = 100                  
grp <- st_intersects(nc1, nc1[n,1]   , sparse = F ) 

neighborhood <- nc1[grp,]
neighborhood 

plot(neighborhood$geom)
plot(nc1[n,1], col = 'blue', add = TRUE)   #

这段代码很容易扩展。我写了一个快速的小函数来显示邻国的名字（这里没有显示），但这个问题似乎很可能是一个与绘图相关的问题。

该图显示在 链接

推荐阅读