使用 spdplyr 对空间数据进行分组？

Question

嗨，感谢阅读我的文章，我正在尝试对包含城市多边形的数据集进行分组，以避免出现重复值，但是当我尝试使用 group_by（用 spdplyr 包屏蔽）时，数据未分组。 有谁知道我怎么能做到这一点？ 谢谢您的帮助

我使用的代码如下：

# DATA:

 nombre_alcaldia                                                                     n
   <chr>                                                                       <POLYGON>
 1 Álvaro Obregón  ((-99.18871 19.39475, -99.18868 19.39468, -99.18857 19.39442, -99.18…
 2 Álvaro Obregón  ((-99.18871 19.39475, -99.18868 19.39468, -99.18857 19.39442, -99.18…
 3 Álvaro Obregón  ((-99.18871 19.39475, -99.18868 19.39468, -99.18857 19.39442, -99.18…
 4 Álvaro Obregón  ((-99.18871 19.39475, -99.18868 19.39468, -99.18857 19.39442, -99.18…
 5 Álvaro Obregón  ((-99.18871 19.39475, -99.18868 19.39468, -99.18857 19.39442, -99.18…
 6 Álvaro Obregón  ((-99.18871 19.39475, -99.18868 19.39468, -99.18857 19.39442, -99.18…
 7 Álvaro Obregón  ((-99.18871 19.39475, -99.18868 19.39468, -99.18857 19.39442, -99.18…
 8 Álvaro Obregón  ((-99.18871 19.39475, -99.18868 19.39468, -99.18857 19.39442, -99.18…
 9 Álvaro Obregón  ((-99.18871 19.39475, -99.18868 19.39468, -99.18857 19.39442, -99.18…
10 Álvaro Obregón  ((-99.18871 19.39475, -99.18868 19.39468, -99.18857 19.39442, -99.18…

# grouping

datos_esp1 |> 
  group_by(nombre_alcaldia) |> 
  summarise(n = geometry )

Answer 1

请检查此链接

https://github.com/r-spatial/sf/issues/915

spdplyr 不能 group_by 几何。 有两种选择，

1. 您将获得独特的几何图形，而不是根据实际范围进行分组。

2. 您 group_by X 列，稍后您将加入一组独特的几何图形。