使用 sf 包识别接触多边形边缘的线

Question

我正在使用 R 中的 sf 包进行空间连接。具体来说，我试图找到代表新泽西州每个人口普查区（多边形）边界的道路（线）。注意：对于新泽西州的许多人口普查区，人口普查区的边界都是道路。

我尝试使用 st_touches 在线条（道路）和多边形（人口普查区）之间进行连接，我相信它可以识别共享公共点但其内部不相交的几何图形。使用 st_touches 确实成功识别了构成人口普查区边界的道路，但它也识别了人口普查区内部接触人口普查区边缘的道路。

以下图为例，当我使用 st_touches 执行线-多边形连接时，被识别为接触的道路是红色的 4 条道路，它们代表人口普查区的实际边界。然而，它还以蓝色标识了 3 条道路，它们位于人口普查区多边形的内部。

这可能只是对多边形的误解，但我想知道多边形和线之间的空间连接是否将多边形本身视为一种“闭合”线。换句话说，我认为连接将构成多边形边缘的坐标视为线，这样接触多边形边缘的所有实际线（道路）中的任何一条都被标记为接触多边形的边缘。多边形，无论道路是在多边形内部还是在多边形顶部。当使用 st_touches 时，是否存在这种情况，多边形内部的线不应该被识别为接触多边形，因为它们位于内部多边形的？

如果这个问题有点（或非常）令人困惑，我很抱歉，所以如果需要一个可重现的示例，我可以制作一个。然而，我只是想提出这个问题，看看是否有人可以在没有样本数据的情况下回答它。谢谢！

Answer 1

看起来很可疑，所有三条不正确的道路都在多边形内部，没有一条在多边形外部。我期望相反 - 你应该得到边界道路，加上与边界共享一个点的外部道路（与边界正交）。

代码和数据示例可能会有所帮助。但我认为可能发生的事情：

• 您正在使用 sf 和球形 s2 模式（现在是默认模式）。
• 轨迹多边形在这里是倒置的 - 即它们描述了地球的所有表面除了人口普查轨迹。
• 具有倒置的多边形，结果在技术上是正确的，“内部”道路在多边形之外，并触摸它。