如何用R语言SF库为Multipolygons添加相交Polylines列？

核心思路 #

借助sf包的空间连接能力，通过空间相交判断将多边形（Multipolygons）与对应的相交折线（Polylines）数据关联，再按需聚合匹配到的折线行数据，最终给多边形新增关联列。

具体实现步骤 #

1. 准备工作：加载库与确认数据格式 #

首先加载sf包，确保你的两个空间对象都是sf格式（如果不是，先转换）：

library(sf)

# 假设你的多边形对象名为poly_sf，折线对象名为line_sf
# 非sf格式转换示例：
# poly_sf <- st_as_sf(your_multipolygon_data)
# line_sf <- st_as_sf(your_polyline_data)

2. 空间连接：匹配相交的折线数据 #

使用st_join()函数，指定st_intersects作为空间谓词，只保留与多边形相交的折线数据：

# 左连接保留所有多边形，匹配到相交折线的行会关联对应数据
joined_data <- st_join(poly_sf, line_sf, join = st_intersects, left = TRUE)

此时如果一个多边形对应多条相交折线，会生成多行重复的多边形记录，每条记录对应一条折线数据。

3. 聚合多条折线数据（按需选择） #

如果需要将同一多边形匹配到的所有折线数据合并到单列中，结合dplyr进行聚合：

library(dplyr)

# 假设多边形有唯一标识列poly_id，折线的目标数据列是line_name、line_length
aggregated_data <- joined_data %>%
  group_by(poly_id) %>%
  summarise(
    # 合并折线名称，用逗号分隔
    matched_lines = paste(line_name, collapse = ", "),
    # 计算相交折线的总长度
    total_line_length = sum(line_length, na.rm = TRUE),
    # 保留多边形的几何信息，避免丢失空间属性
    geometry = first(geometry)
  ) %>%
  st_sf() # 确保结果仍为sf空间对象

4. 合并整行折线数据（进阶需求） #

如果需要把折线的整行非几何数据合并成一个列，可以这样处理：

aggregated_data <- joined_data %>%
  group_by(poly_id) %>%
  summarise(
    line_row_details = paste(
      do.call(paste, c(select(., -geometry, -poly_id), sep = "; ")),
      collapse = " | "
    ),
    geometry = first(geometry)
  ) %>%
  st_sf()

关键注意点 #

• 空间谓词可按需调整：如果需要匹配多边形包含折线的情况，用st_contains；匹配边界接触的情况，用st_touches
• 聚合时必须保留geometry列，否则结果会失去空间对象属性
• 如果不需要聚合，直接使用joined_data即可，每条多边形-相交折线的组合对应一行记录

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