pandas合并geometry列触发KeyError索引不存在问题求解

• 直接触发KeyError的核心原因是数据集列索引异常：报错信息中geometry_x被拆分为单字符元组('g', 'e', 'o', 'm', 'e', 't', 'r', 'y', '_', 'x')，说明rpg15_19已经不是普通单级字符串列索引，大概率是之前做空间连接、表合并（sjoin/merge/join）或列重命名操作时误生成了多级列索引（MultiIndex），导致传入字符串列名时，pandas将字符串按字符迭代拆分为多级索引元组去匹配，自然找不到目标列。
• 原代码存在两层逻辑缺陷：一是用!= None判断地理列空值的写法无效，geopandas空几何、pd.NA/np.nan与None做普通等值比较时不会返回正确布尔值，会出现漏筛；二是逐列赋值的写法会出现覆盖问题，后执行的赋值会覆盖前一列已经写入的有效值，最终只会保留geometry列的非空值，达不到三列非空值整合的目的。

1. 先验证并修复列索引，运行以下代码确认列状态，拍平异常的多级列：

import pandas as pd
# 打印列名确认实际列名状态
print(rpg15_19.columns)
# 拍平多级列索引，将元组形式的列名转为正常字符串
if isinstance(rpg15_19.columns, pd.MultiIndex):
    rpg15_19.columns = [''.join(col).strip() if isinstance(col, tuple) else col for col in rpg15_19.columns]

2. 用combine_first方法按优先级合并三列非空值，该方法会自动保留前一列的非空值，仅用后一列的值补全前一列的空位置，不会出现覆盖问题：

import geopandas as gpd
# 按 geometry_x -> geometry_y -> geometry 的优先级合并非空几何值
rpg15_19['geom'] = (
    rpg15_19['geometry_x']
    .combine_first(rpg15_19['geometry_y'])
    .combine_first(rpg15_19['geometry'])
)

# 可选：过滤掉geom为空/空几何的无效行
rpg15_19 = rpg15_19[
    ~pd.isna(rpg15_19['geom']) 
    & ~gpd.GeoSeries(rpg15_19['geom']).is_empty
]

# 可选：将geom设为GeoDataFrame的活动几何列，保留原始坐标系
rpg15_19 = gpd.GeoDataFrame(
    rpg15_19, 
    geometry='geom', 
    crs=rpg15_19.crs if hasattr(rpg15_19, 'crs') else None
)

注意：如果拍平列索引后仍然找不到geometry_x/geometry_y/geometry列，对照第一步打印的实际列名，调整代码中对应的列名字符串即可。

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