在R terra包处理矢量图层遇错误：LinearRing点未形成闭合线串

问题背景 #

• 涉及数据集：GSHHS shp 2.3.7版本（以GSHHS_f_L1.shp为例，压缩包内多个图层均存在此问题）
• 当前环境：R 4.4.0、terra 1.7.78、Linux Ubuntu 22.04.4 LTS，近期仅从GitHub安装过gower包

问题演变 #

1. 一月前操作完全正常：

library(terra)
xx <- vect('~/GSHHS_f_L1.shp')
plot(xx)
is.valid(xx)
# 输出：TRUE

绘图及几何有效性检查均无异常。

2. 一周前有效性检查失效，但可通过makeValid修复：

is.valid(xx)
# 输出：FALSE

yy <- makeValid(xx)
is.valid(yy)
# 输出：TRUE

仍可正常执行后续操作。

3. 近日出现无法修复的致命错误：

is.valid(xx)
# 报错：Error: IllegalArgumentException: Points of LinearRing do not form a closed linestring

yy <- makeValid(xx)
# 触发相同报错

仅绘图功能可用，尝试创建0缓冲区的常规修复方案也触发相同错误。已更新terra包、重新下载数据集，均无改善；偏好使用terra包而非sf，寻求专属解决方案。

Terra包下的解决方案 #

方法1：手动闭合线性环 #

遍历每个要素的线性环，检查首尾点是否一致，不一致则手动添加首点作为尾点闭合环：

library(terra)
xx <- vect('~/GSHHS_f_L1.shp')

fix_closed_ring <- function(vec_obj) {
  # 提取几何信息
  geom_df <- geom(vec_obj)
  # 按要素分组处理
  for (geom_id in unique(geom_df$geom)) {
    geom_subset <- geom_df[geom_df$geom == geom_id, ]
    for (part_id in unique(geom_subset$part)) {
      ring_subset <- geom_subset[geom_subset$part == part_id, c("x", "y")]
      # 检查首尾坐标是否完全一致
      if (!all(ring_subset[1, ] == ring_subset[nrow(ring_subset), ])) {
        # 复制首行作为新的尾行，闭合环
        new_row <- geom_subset[geom_subset$part == part_id, ][1, ]
        new_row$x <- ring_subset$x[1]
        new_row$y <- ring_subset$y[1]
        geom_df <- rbind(geom_df, new_row)
      }
    }
  }
  # 重新构建矢量对象
  vect(geom_df, type = "polygons", crs = crs(vec_obj))
}

# 修复几何
yy <- fix_closed_ring(xx)
# 验证有效性
is.valid(yy)

方法2：预处理几何后修复 #

先通过simplify轻微简化几何，再尝试makeValid修复，可能解决闭合问题：

xx_simplified <- simplify(xx, tolerance = 0.001) # tolerance可根据数据精度调整
yy <- makeValid(xx_simplified)
is.valid(yy)

方法3：降级terra版本 #

问题可能源于terra新版本的兼容性问题，尝试降级到之前正常工作的版本（例如1.7-55）：

# 从CRAN安装指定旧版本
install.packages("terra", version = "1.7-55", repos = "https://cran.r-project.org")

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