Geopandas叠加交集操作报错['geometry'] not found in axis 求助

问题分析 #

该报错源于gpd.overlay操作要求两个GeoDataFrame必须具备有效的活跃几何列，且内部处理依赖几何列的正确关联。尽管你已将wells1的几何列设为buffer，仍可能存在活跃几何列未正确识别、CRS不匹配或几何列命名冲突的问题。

具体解决步骤 #

1. 确认活跃几何列状态 #

先检查两个DataFrame的活跃几何列名称，确保它们被正确设置：

# 查看wells1的活跃几何列名称
print(wells1.geometry.name)
# 查看wells2的活跃几何列名称
print(wells2.geometry.name)

如果wells2的输出不是geometry，显式设置其活跃几何列：

wells2 = wells2.set_geometry("geometry")

2. 统一几何列名称（推荐） #

gpd.overlay对几何列命名的一致性更友好，可将wells1的活跃几何列重命名为geometry：

# 方法1：先删除原geometry列，再将buffer列重命名为geometry
wells1 = wells1.drop(columns="geometry").rename(columns={"buffer": "geometry"})

# 方法2：设置活跃几何列后，直接重命名几何列
wells1 = wells1.set_geometry("buffer")
wells1 = wells1.rename_geometry("geometry")

完成后执行交集操作：

res_intersection = gpd.overlay(wells2, wells1, how='intersection')

3. 检查并统一CRS #

两个GeoDataFrame的CRS必须完全一致，否则会引发异常：

# 打印当前CRS
print("wells1 CRS:", wells1.crs)
print("wells2 CRS:", wells2.crs)

# 统一CRS（将wells2转换为wells1的CRS）
if wells1.crs != wells2.crs:
    wells2 = wells2.to_crs(wells1.crs)

4. 修复无效几何对象 #

如果存在损坏的几何对象，也可能导致操作失败：

# 过滤无效几何
wells1 = wells1[wells1.is_valid]
wells2 = wells2[wells2.is_valid]

# 若仍有无效项，尝试修复（需shapely>=2.0）
from shapely.validation import make_valid
wells1['geometry'] = wells1['geometry'].apply(make_valid)
wells2['geometry'] = wells2['geometry'].apply(make_valid)

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