Python提取GeoJSON建筑足迹点X坐标并保留层级结构

1. 先确认坐标列表C的完整性 #

• 直接统计列表C的元素数量，若不等于9835，说明BIM软件导出环节丢了数据，检查导出设置：是否过滤了部分建筑、是否有坐标范围限制、导出时是否遗漏了某些楼层/类别。
• 随机抽选几个未被覆盖的建筑ID，对比原始GeoJSON文件和列表C，确认对应建筑的MultiPolygon数据是否存在。

2. 重构Extract函数的嵌套遍历逻辑 #

MultiPolygon的结构是**[ [ [ [x,y,z], [x,y,z], ... ], ... ], ... ]**，必须严格匹配三层嵌套（MultiPolygon → Polygon集合 → 每个Polygon的坐标环），不能跳过任何一层：
以下是通用逻辑示例（以Python为例，可对应转换为BIM可视化编程语言的节点逻辑）：

def extract_x_coords(multipolygon):
    building_x = []
    # 遍历MultiPolygon下的所有Polygon
    for polygon in multipolygon:
        polygon_x = []
        # 遍历单个Polygon下的所有坐标环（外环+可能的内环）
        for ring in polygon:
            # 提取当前坐标环里所有点的X值
            ring_x = [point[0] for point in ring]
            polygon_x.append(ring_x)
        building_x.append(polygon_x)
    return building_x

# 批量处理所有建筑的坐标列表
all_x_coords = [extract_x_coords(building_geom) for building_geom in C]

• 若用节点式编程（如Dynamo、Grasshopper），需搭建多层嵌套循环节点，确保每个层级的结构都被遍历到，避免仅做单层提取。

3. 处理异常数据 #

部分建筑可能存在空Polygon、空坐标环或不完整的点数据，需在函数中增加判断，避免中断遍历：

def extract_x_coords(multipolygon):
    building_x = []
    for polygon in multipolygon:
        if not polygon:
            continue
        polygon_x = []
        for ring in polygon:
            if not ring:
                continue
            # 确保点数据包含X坐标
            ring_x = [point[0] for point in ring if len(point) >= 1]
            polygon_x.append(ring_x)
        if polygon_x:
            building_x.append(polygon_x)
    return building_x

4. 验证结果 #

处理完成后，统计结果列表的长度，若等于9835则说明覆盖全部建筑；再随机选几个建筑，对比提取的X坐标层级结构与原始GeoJSON是否一致。

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