GeoPandas explode 怎么用：多部件几何展开、索引回写与面积质检
GeoPandas 读入的 MultiPolygon 在地图上没问题，可一做空间连接或逐地块统计，index 重复、属性回写错行的问题就出现了。explode 能把一个多部件几何展开为多行，但返回索引结构和后续主键策略需要先设计，否则导出后很难再追踪每个部件来自哪条记录。
GeoPandas explode：先确定可验证的结果边界
把 explode 看成数据表的规范化操作：空间被展开后，索引和业务 ID 必须同时重建。开始操作前，把输入批次、坐标参考、关键字段和成果用途写进处理记录；操作后至少比较数量、范围和一项业务统计。这样遇到异常时，才能定位是数据、参数还是规则出了问题。
索引不是业务主键
explode 默认可能保留多级索引，用于指出原行和部件序号；它适合临时定位，不适合作为跨文件关联键。导出前应把原始 ID 放入普通字段。
几何展开不会理解你的字段语义
每个子部件会继承同一组属性。若原字段是对象级总量，直接复制到每行后再求和会重复计算；应另建部件级指标或按面积权重分配。
投影决定面积核验是否可信
经纬度坐标下的 area 单位是平方度，不能拿来验收。应在适合研究区的投影 CRS 中比较 explode 前后的面积。
可执行实操流程
- 读取数据后检查 geom_type、is_valid、is_empty 和原始唯一 ID；先保存只读备份。
- 将原主键复制为 parent_id，在投影 CRS 下记录每个 parent_id 的原始面积。
- 执行 gdf.explode(index_parts=True)，reset_index 后把原行号和部件号明确改名。
- 生成新的 part_id，按 parent_id 汇总 exploded 面积与部件数，并比对原图层。
- 导出 GeoPackage 后重新读取，检查字段类型、索引是否已转普通列，并抽查空间位置。
parts = gdf.to_crs(4547).explode(index_parts=True).reset_index(names=['source_row','part_no'])
parts['parent_id'] = parts['parcel_id']
parts['part_id'] = parts['parent_id'].astype(str) + '_' + parts['part_no'].astype(str)项目避坑与质量检查
有些数据把人口或产值放在多部件对象上。explode 后如果把这些总量在部件层求和，会无声放大统计值。应该保留原对象总量列，并明确是否按面积、权重栅格或业务比例分配成 part_value。
| 检查节点 | 应保留的证据 | 异常时的第一动作 |
|---|---|---|
| 输入 | 来源、范围、CRS、字段统计 | 回看原始批次与空值/重复值 |
| 处理中 | 参数、日志与抽样结果 | 在最小样本复现而非直接重跑全量 |
| 输出 | 数量、范围与关键业务统计 | 与基线或人工判读样本比对 |
让流程可以被同事复跑
把这次选择的参数、拒绝的候选方案和抽样位置一并保存。GIS 项目最难交接的不是工具名称，而是某个阈值为何合理、某个异常为何可以接受。将这些判断写在处理日志中，数据更新后才能用同一把尺子复核。
复核记录应至少能回答三件事：输入数据从哪里来、为什么选这个参数、异常结果如何处置。把这三项与一份小样本结果绑定保存，比只留最终截图更能抵御人员交接、数据更新和审核追问。
FAQ
explode 后为什么出现 MultiIndex？
它用来保留原行与部件位置。导出或关联前可 reset_index，并把需要保留的层级命名为字段。
能否直接在 WGS84 下比面积？
不应使用平方度作为面积质量依据。先投影到合适的米制 CRS。
explode 会修复无效几何吗？
不会。应先识别并处理无效或空几何，再展开和统计。
总结
GeoPandas explode 的风险不在代码一行，而在展开后的索引、主键和总量语义是否仍然可追溯。可靠的 GIS 成果不只是一张看起来正确的图，而是输入边界、参数理由和复核证据都经得起追问的可复现过程。