如何筛选Polygon图层中与multiline要素邻接的多边形

按优先级依次排查以下问题，99%的场景都能解决：

• 坐标系不匹配（最高发）
  两个空间图层坐标系不一致时，所有空间关系计算都会返回错误结果，是最常见的踩坑点。
  先执行以下代码检查坐标系：

  st_crs(your_polygon_layer)
  st_crs(your_river_layer)
  

  如果两个返回的EPSG编码、WKT字符串不一致，先做坐标系统一，优先转成研究区域适用的投影坐标系（不建议直接用经纬度地理坐标系做邻接计算，精度误差大）：

  # 将河流图层坐标系转换为与地块图层完全一致
  your_river_layer <- st_transform(your_river_layer, st_crs(your_polygon_layer))
  

• 空间范围无实际重叠
  数据导入错误、坐标单位错配（比如把米单位的投影坐标当成经纬度）、坐标平移错误都会导致两个图层实际位置差了成百上千公里，根本没有交集。
  执行代码检查两个图层的边界范围：

  st_bbox(your_polygon_layer)
  st_bbox(your_river_layer)
  

  如果两个范围的坐标值差了2个数量级以上、没有重叠区间，先回溯数据导入流程修正坐标问题，再做后续计算。
• 几何间存在微小容差缝隙
  不同来源的矢量化数据、人工数字化误差，会导致看起来完全接壤的线和多边形，实际存在毫米到米级的缝隙，严格的相交判定会判定为不接触。
  这种情况给河流线生成一个极小的容差缓冲区（缓冲区大小取数据生产精度允许的最小值，比如投影坐标系下取0.1~1米，刚好覆盖缝隙又不会误选非邻接地块），再做相交筛选：

  # 注意dist参数的单位与坐标系单位一致
  river_tolerance_buffer <- st_buffer(your_river_layer, dist = 0.5)
  # 加sparse = F直接返回逻辑矩阵，可直接用于子集筛选
  adjacent_parcels <- your_polygon_layer[st_intersects(your_polygon_layer, river_tolerance_buffer, sparse = F), ]
  

• 几何存在拓扑错误
  多边形自交、多线断裂、几何闭合错误等无效拓扑问题，会导致空间关系计算返回异常空值。先修复几何再计算：

  your_polygon_layer <- st_make_valid(your_polygon_layer)
  your_river_layer <- st_make_valid(your_river_layer)
  # 修复后重新执行相交筛选
  adjacent_parcels <- your_polygon_layer[st_intersects(your_polygon_layer, your_river_layer, sparse = F), ]
  

筛选完成后可直接绘图验证结果正确性：

plot(st_geometry(adjacent_parcels), border = "red", lwd = 2)
plot(st_geometry(your_river_layer), col = "blue", add = T)

内容的提问来源于stack exchange，提问作者FreakyT