sf包st_read导入GDB文件报错：OGR error: cannot export to WKT

可行修复方案 #

1. 定向导入图层并跳过无效几何 #

更新后的GDAL/GEOS对几何格式校验更严格，GDB中可能存在个别异常图层或对象导致导入失败。先排查具体图层：

# 列出GDB内所有图层
st_layers("ProtectedConservedArea_2023.gdb")

# 单独导入目标图层，同时跳过无法解析的几何
st_read("ProtectedConservedArea_2023.gdb", layer = "你的目标图层名", options = "SKIPFAILURES=YES")

SKIPFAILURES=YES参数会让GDAL自动跳过损坏的几何对象，避免整个导入流程中断。

2. 调整GDAL的WKT输出配置 #

新版本GDAL的WKT格式默认设置可能与旧数据不兼容，强制切换为WKT1格式并关闭S2几何计算：

# 关闭S2球面几何支持
sf::sf_use_s2(FALSE)
# 设置GDAL使用WKT1格式
Sys.setenv(GDAL_WKT_FORMAT = "WKT1")

# 重新尝试导入
st_read("ProtectedConservedArea_2023.gdb")

3. 用GDAL命令行预处理数据 #

如果R内的方法无效，先通过GDAL命令行将GDB转换为兼容性更好的格式（如GeoPackage）：

# 转换整个GDB到GeoPackage
ogr2ogr -f GPKG protected_areas.gpkg ProtectedConservedArea_2023.gdb

# 或仅转换指定图层
ogr2ogr -f GPKG protected_areas.gpkg ProtectedConservedArea_2023.gdb 目标图层名

之后在R中导入转换后的文件：

st_read("protected_areas.gpkg")

4. 检查并修复异常几何 #

用GDAL工具定位并修复损坏的几何对象：

# 查找图层内的无效几何
ogrinfo -sql "SELECT * FROM 目标图层名 WHERE ST_IsValid(geometry) = 0" ProtectedConservedArea_2023.gdb

# 修复几何并生成新的GDB文件
ogr2ogr -f GDB fixed_protected_areas.gdb ProtectedConservedArea_2023.gdb -sql "SELECT ST_MakeValid(geometry) AS geometry, * FROM 目标图层名"

再导入修复后的GDB文件即可。

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