R语言sf包复杂LINESTRING转闭合POLYGON形状异常求解

问题背景 #

我和学生耗时近两周开展空间数据处理工作，目标是实现避开海岸线的两点GPS距离计算，基于R空间计算相关教程完成了基础逻辑搭建。
我们下载的目标区域shapefile仅包含LINESTRING类型几何数据，由于后续需要将地图转为网格计算，计划先做两项预处理：

• 添加边界"角点"避免网格计算出现边缘误差
• 将开放岸线闭合为多边形
  实际操作中闭合后的多边形形状完全不符合预期，测试所用的岸线数据为公开的行政区划岸线矢量文件。

已执行操作 #

首先提取线几何坐标、添加角点后绘制点线图，渲染形状正常，对应代码如下：

aoi_boundary_HARV <- st_read("coasts_subnational/coasts_subnational.shp")

# 闭合多边形预处理
coord <- as.numeric(as.character(aoi_boundary_HARV$geometry[[1]]))
coord <- data.frame(matrix(coord, ncol = 2))
colnames(coord) <- c("x","y")
corner1 <- c(x = coord$x[nrow(coord)], y= coord$y[nrow(coord)])
corner2 <- c(x = coord$x[1], y= coord$y[1])
corner <- c(x = corner2[1], y = corner1[2])
names(corner) <- names(cord_Table)
cord_Table_plusCorner <- rbind(cord_Table, corner)
cord_Table_plusCorner <- as.data.frame(cord_Table_plusCorner)
colnames(cord_Table_plusCorner) <- c("x","y")
polygone_land <- data.frame(id = 1, 
                            x=cord_Table_plusCorner$x, 
                            y=cord_Table_plusCorner$y)

plot(polygone_land$x, polygone_land$y, type = "b")

点线渲染效果：

确认点线形状无异常后，我们尝试将其转换为POLYGON类型的sf对象，对应代码如下：

CRS <- st_crs(aoi_boundary_HARV)
plot_aoi_boundary_HARV <- st_as_sf(polygone_land, 
                                   coords = c("x","y"), 
                                   crs = CRS)

polys = plot_aoi_boundary_HARV %>% 
  dplyr::group_by(id) %>% 
  dplyr::summarise() %>%
  st_cast("POLYGON")

plot(polys)

转换后多边形外轮廓大致符合预期，但内部拓扑结构完全错误，错误效果如下：

此前我们尝试用caveman方法修复拓扑无效，也参考了GIS社区内点转多边形、线转多边形的两类通用方案操作，均未解决问题。

问题根因 #

出现拓扑错误的核心原因是转换逻辑存在缺陷：

1. 直接将线几何拆分为独立点集后，用summarise()聚合转多边形时，程序会按点在数据框中的存储顺序依次连线，仅补单个角点会导致出现跨区域的对角连线，直接引发多边形自相交
2. 拆分点的过程丢失了原始线几何的连续拓扑关系，开放岸线的首尾没有沿规划的外边界连续衔接，强行转多边形会生成大量无效的内部面。

可行修复步骤 #

不要拆分为独立点集后重组，直接操作sf线几何对象完成闭合和转换，步骤如下：

1. 提取原始岸线的首尾端点，按外边界的连续走向拼接补充的边界线段，生成完整的闭合LINESTRING，不要跳点连线
   参考代码：

   library(sf)
   library(dplyr)
   # 提取原始岸线几何
   original_line <- aoi_boundary_HARV$geometry[[1]]
   # 提取线的首尾点坐标
   start_pt <- st_startpoint(original_line)
   end_pt <- st_endpoint(original_line)
   start_coord <- st_coordinates(start_pt)[,c("X","Y")]
   end_coord <- st_coordinates(end_pt)[,c("X","Y")]
   # 按连续顺序构造补充边界：从岸线终点到角点，再到岸线起点，保证连线不跨区域
   supplement_line <- st_linestring(
     rbind(
       end_coord,
       c(end_coord[,"X"], start_coord[,"Y"]), # 之前定义的角点
       start_coord
     )
   )
   # 拼接原始岸线和补充边界，合并为单条连续线
   closed_line <- st_union(original_line, supplement_line) %>% st_line_merge()
   # 校验是否为无断点的单条线
   stopifnot(st_geometry_type(closed_line) == "LINESTRING")
   

2. 对闭合的连续线直接做多边形化，不要走点聚合转多边形的逻辑：

   land_poly <- st_polygonize(closed_line) %>% 
     st_collection_extract("POLYGON")
   

3. 最后做拓扑校验和修复：用st_is_valid()检查多边形合法性，如果存在精度导致的微小自相交，直接调用st_make_valid()即可完成修复。

注意：如果原始岸线包含多个独立线段（如多个岛屿、多段行政边界），需要先按线段所属的地理单元分组，每组单独做闭合和转多边形操作，不要把所有点归到同一个分组里处理。

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