R语言sf包计算多边形质心报错：Loop无效，多边形闭合问题

错误原因分析 #

你的报错核心有两点：

1. Loop 7 is not valid: Edge 185 crosses edge 187：第7个几何图形存在自相交问题，即便执行st_make_valid，也可能因几何类型混合（比如生成了GEOMETRYCOLLECTION）未彻底修复。
2. not implemented for objects of class sfc_GEOMETRY：st_centroid默认使用的S2引擎不支持混合几何类型的sfc_GEOMETRY对象，这类对象可能同时包含多边形、点、线等多种几何类型。

分步解决方案 #

1. 拆分混合几何类型 #

先将所有几何强制转换为单一多边形类型，拆分可能存在的GEOMETRYCOLLECTION：

# 强制转换为多边形，拆分几何集合，忽略类型转换警告
shapefile_data <- shapefile_data %>%
  st_cast("POLYGON", warn = FALSE) %>%
  st_make_valid()

2. 修复自相交与未闭合多边形 #

如果st_make_valid无法修复自相交，用微小缓冲+闭合操作处理：

# 用极小缓冲抹平自相交边缘，再强制闭合多边形
shapefile_data$geometry <- shapefile_data$geometry %>%
  st_buffer(dist = 0.000001) %>%  # 数值可根据数据精度调整
  st_make_valid() %>%
  st_close()  # 确保多边形首尾点重合，解决未闭合问题

3. 绕过S2引擎计算质心 #

如果S2引擎对复杂几何兼容性差，暂时关闭S2用平面几何计算：

# 关闭S2引擎，使用平面坐标系计算质心
sf_use_s2(FALSE)

shapefile_data <- shapefile_data %>%
  mutate(centroid = st_centroid(geometry),
         lon = st_coordinates(centroid)[, 1],
         lat = st_coordinates(centroid)[, 2])

# 完成后可重新开启S2（按需）
sf_use_s2(TRUE)

4. 验证修复结果 #

修复后检查几何有效性，确保无残留错误：

# 检查所有几何是否有效
invalid_count <- sum(!st_is_valid(shapefile_data))
if (invalid_count > 0) {
  warning(paste("仍有", invalid_count, "个无效几何"))
} else {
  message("所有几何已修复为有效状态")
}

完整修复代码 #

library(sf)
library(dplyr)

# 加载Shapefile数据
shapefile_path <- "/path/to/shapefile.shp"
shapefile_data <- st_read(shapefile_path)

# 1. 转换为单一多边形类型并修复有效性
shapefile_data <- shapefile_data %>%
  st_cast("POLYGON", warn = FALSE) %>%
  st_make_valid()

# 2. 修复自相交与未闭合问题
shapefile_data$geometry <- shapefile_data$geometry %>%
  st_buffer(dist = 0.000001) %>%
  st_make_valid() %>%
  st_close()

# 3. 转换到WGS84坐标系
shapefile_data <- st_transform(shapefile_data, crs = "+proj=longlat +ellps=WGS84 +datum=WGS84")

# 4. 关闭S2计算质心
sf_use_s2(FALSE)
shapefile_data <- shapefile_data %>%
  mutate(centroid = st_centroid(geometry),
         lon = st_coordinates(centroid)[, 1],
         lat = st_coordinates(centroid)[, 2])
sf_use_s2(TRUE)

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