如何将每个MultiLineString转换为单个LineString要素

问题说明 #

待处理的河段shapefile中，仅ID为8168547、8171738、8170616、8169920的4条河段为MultiLineString类型，其余要素均为LineString类型。直接调用sf包st_cast("LINESTRING")会将每个MultiLineString拆分为多条独立LineString，造成数据行数增加，无法满足每个MultiLineString对应输出单个LineString、保持原数据行数不变的需求。

实现逻辑 #

st_cast拆分MultiLineString的本质是默认把MultiLineString下的每个子线段作为独立要素输出，要得到单条LineString，需要先把MultiLineString下的所有子线段按节点顺序拼接为连续线，再转换类型即可，全程不会改变原数据的属性表结构和行数。

具体代码实现 #

1. 加载依赖、读取数据

library(sf)
# 替换为本地shapefile的实际存储路径
river_data <- st_read("river_shapefile.shp")

2. 编写几何转换函数

ml_to_ls <- function(geom) {
  # 本身就是LineString的要素直接返回，不做处理
  if (st_geometry_type(geom) == "LINESTRING") return(geom)
  # 提取MultiLineString所有节点坐标，按存储顺序拼接
  line_coords <- st_coordinates(geom)[, c("X", "Y")]
  # 生成新的单条LineString几何
  st_linestring(as.matrix(line_coords))
}

3. 批量替换几何列完成转换

# 逐行处理所有几何要素
processed_geom <- lapply(st_geometry(river_data), ml_to_ls)
# 替换原数据几何列，保留原有坐标系和属性表
st_geometry(river_data) <- st_sfc(processed_geom, crs = st_crs(river_data))

异常场景兼容处理 #

如果遇到MultiLineString的子线段是连续但未做节点连接的情况，可以先用st_line_merge()做线合并，再兜底处理，避免直接拼接节点产生错误：

library(dplyr)
river_data <- river_data |> 
  mutate(
    geom = st_sfc(
      lapply(st_geometry(geom), function(g) {
        if (st_geometry_type(g) != "MULTILINESTRING") return(g)
        merged_line <- st_line_merge(g)
        # 合并后已转为LineString直接返回，否则提取坐标拼接
        if (st_geometry_type(merged_line) == "LINESTRING") {
          return(merged_line)
        } else {
          coords <- st_coordinates(g)[, c("X", "Y")]
          return(st_linestring(as.matrix(coords)))
        }
      }),
      crs = st_crs(geom)
    )
  )

处理完成后可运行table(st_geometry_type(river_data))做校验，输出结果应全部为LineString，且数据行数和原文件完全一致。

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