在R语言中将赞比西河多边形转为单一折线的方法咨询

针对你已筛选出的赞比西河两个多边形数据，以下提供两种常见的转换方案：提取岸线折线，或提取河道中线折线。

方案1：提取岸线折线 #

利用多边形的边界生成岸线，再合并为单一折线：

1. 加载依赖包并提取多边形边界
   多边形的轮廓即为岸线，用st_boundary()提取每个多边形的边界线：

library(sf)

# 提取两个多边形的边界
river_boundaries <- st_boundary(River1)

2. 合并边界并转为单一折线
   先合并所有边界为一个几何对象，再尝试转为单条折线（若边界连通则成功，否则为多段线）：

# 合并边界
united_boundaries <- st_union(river_boundaries)
# 尝试合并为单条折线
river_shoreline <- st_line_merge(united_boundaries)

3. 可选：简化折线（减少冗余点）
   若折线点数量过多，可通过st_simplify()简化，dTolerance为容差（单位与数据坐标系一致，建议先转投影坐标系）：

# 先转投影坐标系（以莫桑比克常用的UTM 36S为例）
River1_proj <- st_transform(River1, crs = 32736)
river_boundaries_proj <- st_boundary(River1_proj)
united_boundaries_proj <- st_union(river_boundaries_proj)
river_shoreline_simplified <- st_simplify(st_line_merge(united_boundaries_proj), dTolerance = 100)

4. 查看结果

plot(st_geometry(river_shoreline_simplified), col = "darkblue", lwd = 2)

方案2：提取河道中线折线 #

若需要河道中线，可使用centerline包生成：

1. 安装并加载依赖包

install.packages("centerline")
library(centerline)

2. 合并多边形并生成中线
   先将两个赞比西河多边形合并为一个，再生成中线：

# 合并多边形（转投影坐标系提升精度）
united_river_proj <- st_union(st_transform(River1, crs = 32736))
# 生成中线
river_centerline <- centerline(united_river_proj)

3. 查看结果

plot(st_geometry(united_river_proj), col = "lightblue")
plot(st_geometry(river_centerline), col = "red", lwd = 2, add = TRUE)

注意事项 #

• 若数据为地理坐标系（如WGS84），建议先转换为投影坐标系（如UTM），避免距离计算误差。
• 若合并后仍为多段线（MultiLineString），可根据需求保留或进一步拆分/合并。

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