osmextract提取欧洲绿地多边形时遇vectortranslate错误及特征缺失

问题 #

使用R语言osmextract包的oe_get()函数从geofabrik_europe-latest.osm.pbf提取绿地相关multipolygons，生成.gpkg文件时遇到以下问题：

• 触发报错：Error in sf::gdal_utils(util = "vectortranslate", ...) : gdal_utils vectortranslate: an error occured
• 弹出50+条重复警告：GDAL Message 1: Non closed ring detected. To avoid accepting it, set the OGR_GEOMETRY_ACCEPT_UNCLOSED_RING configuration option to NO
• 提取的特征数仅约80万，但用相同参数提取奥地利单区域.pbf时已得到约60万，欧洲全域结果存在大量多边形缺失
• 尝试添加-nlt GEOMETRY参数后，特征数量未改善，反而新增了大量空几何无效特征；同时不清楚如何在osmextract/vectortranslate中设置OGR_GEOMETRY_ACCEPT_UNCLOSED_RING=NO

使用的代码如下：

library(osmextract)
poly_amenities_green_low_vectortranslate = c(
  "-t_srs", "EPSG:3035",
  "-select", "osm_id, name, leisure, landuse, natural, access, garden_type", 
  "-where", "(leisure IN ('park', 'nature_reserve') AND access NOT IN ('private', 'no', 'customers', 'permit', 'license', 'restricted', 'agricultural', 'forestry')) OR (leisure = 'garden' AND access NOT IN ('private', 'no', 'customers', 'permit', 'license', 'restricted', 'agricultural', 'forestry') AND landuse NOT IN ('allotments', 'residential') AND garden_type IN ('community', 'botanical', 'public', 'flowerbed', 'municipal', 'street_side')) OR (landuse IN ('grass', 'forest', 'meadow', 'flowerbed', 'village_green') AND access NOT IN ('private', 'no', 'customers', 'permit', 'license', 'restricted', 'agricultural', 'forestry')) OR (natural IN ('beach', 'heath', 'wood', 'fell', 'grassland', 'scrub', 'tundra') AND access NOT IN ('private', 'no', 'customers', 'permit', 'license', 'restricted', 'agricultural', 'forestry'))"
)
oe_get("Europe", layer = "multipolygons", provider = "geofabrik", force_download = FALSE, vectortranslate_options = poly_amenities_green_low_vectortranslate, extra_tags = c("leisure", "landuse", "access", "natural", "garden:type"), download_only = TRUE, skip_vectortranslate = FALSE, never_skip_vectortranslate = TRUE)    

解答 #

1. 如何设置OGR_GEOMETRY_ACCEPT_UNCLOSED_RING参数 #

GDAL的配置参数可以通过--config选项传入vectortranslate_options向量，直接添加以下内容即可启用该限制：

"--config", "OGR_GEOMETRY_ACCEPT_UNCLOSED_RING", "NO"

同时注意修正标签名错误：OSM标准标签是garden:type而非garden_type，之前的写法会导致该字段无法被正确识别，直接过滤掉大量符合条件的花园类绿地。修改后的完整参数向量示例：

poly_amenities_green_low_vectortranslate = c(
  "--config", "OGR_GEOMETRY_ACCEPT_UNCLOSED_RING", "NO",
  "-t_srs", "EPSG:3035",
  "-select", "osm_id, name, leisure, landuse, natural, access, \"garden:type\"", 
  "-where", "(leisure IN ('park', 'nature_reserve') AND access NOT IN ('private', 'no', 'customers', 'permit', 'license', 'restricted', 'agricultural', 'forestry')) OR (leisure = 'garden' AND access NOT IN ('private', 'no', 'customers', 'permit', 'license', 'restricted', 'agricultural', 'forestry') AND landuse NOT IN ('allotments', 'residential') AND \"garden:type\" IN ('community', 'botanical', 'public', 'flowerbed', 'municipal', 'street_side')) OR (landuse IN ('grass', 'forest', 'meadow', 'flowerbed', 'village_green') AND access NOT IN ('private', 'no', 'customers', 'permit', 'license', 'restricted', 'agricultural', 'forestry')) OR (natural IN ('beach', 'heath', 'wood', 'fell', 'grassland', 'scrub', 'tundra') AND access NOT IN ('private', 'no', 'customers', 'permit', 'license', 'restricted', 'agricultural', 'forestry'))"
)

2. 特征缺失的核心原因 #

• GDAL转换提前中断：大量未闭合环的存在触发了GDAL的错误，导致vectortranslate进程提前终止，未完成全部数据的提取，这是特征缺失的主要原因。
• 标签名不匹配：代码中错误使用garden_type替代OSM标准标签garden:type，导致符合条件的社区花园、市政花园等数据被遗漏。
• 全域数据处理压力：欧洲全域PBF文件体积远大于奥地利单区域，处理时更容易触发内存不足或GDAL的默认资源限制，进一步加剧数据提取不完整的问题。
• -nlt GEOMETRY的副作用：该参数强制将所有几何类型转为通用GEOMETRY，与OSM multipolygons的原生类型冲突，生成大量无效空几何，占用处理资源，干扰正常提取流程。

3. 后续优化建议 #

• 先应用上述修正后的参数（包含OGR配置和标签名修正）重新运行提取，观察错误是否消失、特征数量是否提升。
• 调整内存分配：处理欧洲全域数据时，建议给R分配更大的内存（Windows系统可使用memory.limit(size=16384)设置16G内存；Linux/macOS可在启动R时添加--max-mem-size=16G参数）。
• 分区域提取：如果全域处理仍有问题，可将欧洲拆分为多个子区域（如按国家）分别提取，再用sf::rbind()合并结果，降低单次处理的数据规模。
• 验证过滤条件：检查access字段的排除列表是否过于严格，是否误排除了部分开放绿地；同时可逐步简化-where条件，定位是否存在其他过滤逻辑导致的数据遗漏。

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