R中转换MultiPolygon Z列为sfc列时遇OGR数据损坏错误

根源分析 #

你遇到的OGR: Corrupt data错误，核心原因是CSV里的MultiPolygon Z类型WKT字符串被截断，导致OGR无法解析完整的空间数据结构。

解决步骤 #

1. 从源头避免WKT截断（QGIS导出环节） #

• 导出CSV时，在QGIS的导出矢量图层对话框中，切换到字段选项卡，将字段最大宽度设置为一个足够大的值（例如999999），覆盖默认的宽度限制，确保长WKT字符串能完整导出。
• 更稳妥的方式是放弃CSV，直接导出为GeoPackage格式：QGIS导出时选择GeoPackage，导入R时用st_read("your_file.gpkg")直接读取为sf对象，完全避免WKT截断问题。

2. 修复已截断的CSV数据 #

如果已经拿到截断的CSV，按以下步骤处理：

• 定位截断行：用字符长度排查异常行

  # 查看每行geometry字符串的长度
  all2$geom_length <- nchar(all2$geometry)
  # 找出长度明显短于其他行的记录
  subset(all2, geom_length < median(geom_length) * 0.5)
  

• 重新导入CSV，确保读取完整字符：
  使用readr包导入，明确指定geometry列为字符型，避免读取时截断

  library(readr)
  all2 <- read_csv("your_data.csv", col_types = cols(geometry = col_character()))
  

• 定位具体错误行并修复：
  单独尝试转换每个geometry字符串，找出失败的行，对比QGIS原始数据补全WKT

  library(sf)
  # 逐个转换，捕获错误
  convert_result <- lapply(all2$geometry, function(wkt) {
    tryCatch(st_as_sfc(wkt), error = function(e) e)
  })
  # 找出转换失败的行号
  error_indices <- which(sapply(convert_result, inherits, "error"))
  # 查看错误行的WKT内容
  all2$geometry[error_indices]
  

  补全这些行的WKT字符串（确保以MULTIPOLYGON Z (开头，以)闭合，坐标结构完整）后，再执行st_as_sf转换。

3. 最终转换验证 #

修复完成后，重新执行转换：

all_sf <- all2 %>% st_as_sf(wkt = "geometry", crs = 4326)
# 验证转换结果
st_crs(all_sf)
plot(st_geometry(all_sf))

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