如何修复GML中的GE_P_NON_PLANAR_POLYGON_DISTANCE_PLANE错误？

这个错误核心是建筑物的多边形几何（如墙面、屋顶面）顶点不在同一平面，超出了CityDoctor设定的平面距离阈值。以下是几种落地的解决思路：

手动修复（小批量数据适用） #

• 优化坐标精度与投影：多数非平面问题源于坐标精度丢失或投影转换误差。将GML转至本地投影坐标系（如UTM分区），减少经纬度转投影时的浮点损耗，再重新验证。
• GIS工具平面化处理：用QGIS、ArcGIS打开GML，选中报错多边形，使用「使平面化」工具自动调整顶点到同一平面；或手动编辑偏离平面的顶点，将其移至其他顶点构成的平面上。

批量自动化修复（大规模数据适用） #

• GDAL/OGR命令行处理：通过ogr2ogr结合几何修复函数批量修正：

  ogr2ogr -f GML repaired_buildings.gml original_buildings.gml -sql "SELECT ST_MakeValid(ST_Force2D(geometry)) AS geometry, * FROM buildings"
  

  若OGR版本支持，可替换为ST_Planeize函数针对性处理平面问题。
• Python脚本批量处理：用geopandas和shapely实现自动化修复：

  import geopandas as gpd
  from shapely.geometry import Polygon
  from shapely.ops import make_valid
  
  # 读取原始数据
  gdf = gpd.read_file("original_buildings.gml")
  
  # 定义平面化逻辑（以拟合前三点平面为例，可按需优化）
  def planarize_poly(poly):
      if not poly.is_valid:
          poly = make_valid(poly)
      coords = list(poly.exterior.coords)
      if len(coords) >= 3:
          # 将所有顶点投影到前三点构成的平面（取第一点的Z值简化处理）
          planar_coords = [(x, y, coords[0][2]) for x, y, z in coords]
          return Polygon(planar_coords)
      return poly
  
  # 批量处理并保存
  gdf['geometry'] = gdf['geometry'].apply(planarize_poly)
  gdf.to_file("repaired_buildings.gml", driver="GML")
  

• 开源CityGML工具修复：使用citygml-tools（3D GeoInfo Group开发），其内置几何修复模块可批量处理非平面多边形问题，通过命令行执行即可。

临时阈值调整（工程误差可接受时） #

若非平面偏差在允许的工程误差范围内，可修改CityDoctor的规则配置文件，增大GE_P_NON_PLANAR_POLYGON_DISTANCE_PLANE对应的距离阈值，避免误报。

内容的提问来源于stack exchange，提问作者Gülce