如何使用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

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

解决方法

好问题。这个问题可能有很多解决方案。但是，此问题的简单答案：iex -S mix是使用带有R内核的Jupyter Lab IDE完成的。以下代码提供了一种回答问题的方法。

nc数据集中有100个县。此代码以彩色显示所选县，并显示所有邻近县。此代码适用于nc的100个县。这里选择了县100。

代码：

"Any tips on how to filter out the actual polygon from the list of adjacent/neighboring polygons?"

此代码易于扩展。我编写了一个快速的小功能，用于显示邻国的名称（此处未显示），但是这个问题似乎很可能是与绘图有关的问题。

该图显示在Link