使用Geopandas overlay(how="union")合并矢量数据遇多几何列报错

目标 #

将Natural Earth的国家多边形与争议区域合并，整合为同一个GeoDataFrame。

问题重现 #

执行代码 #

import geopandas as gpd

gpd.read_file("ne_50m_admin_0_countries.shp").overlay(gpd.read_file("ne_50m_admin_0_breakaway_disputed_areas.shp"), how="union")

报错信息 #

ValueError: GeoDataFrame does not support multiple columns using the geometry column name 'geometry'.

问题原因 #

这不是Geopandas的Bug，是操作上的疏漏。你使用的两个Natural Earth Shapefile都自带geometry列，overlay操作默认会保留两个数据集的所有属性列，导致结果里出现两个同名的geometry列——而GeoDataFrame的规则是不允许存在多个同名的几何列，因此触发报错。

解决方案 #

核心是避免几何列重名，有两种常用处理方式：

方式一：重命名其中一个数据集的几何列 #

import geopandas as gpd

# 分别读取两个Shapefile
countries = gpd.read_file("ne_50m_admin_0_countries.shp")
disputed_areas = gpd.read_file("ne_50m_admin_0_breakaway_disputed_areas.shp")

# 给争议区域的几何列重命名，避免冲突
disputed_areas = disputed_areas.rename_geometry("disputed_geom")

# 执行union合并操作
merged_gdf = countries.overlay(disputed_areas, how="union")

方式二：读取时只保留几何列（适合不需要其他属性的场景） #

import geopandas as gpd

# 读取时仅保留几何列，减少不必要的列冲突
countries = gpd.read_file("ne_50m_admin_0_countries.shp")[["geometry"]]
disputed_areas = gpd.read_file("ne_50m_admin_0_breakaway_disputed_areas.shp")[["geometry"]]

# 执行union合并
merged_gdf = countries.overlay(disputed_areas, how="union")

说明 #

overlay操作会将两个输入GeoDataFrame的所有属性列带入结果，只要存在同名列就会引发冲突。通过重命名几何列或筛选必要列，就能绕开这个限制，顺利完成无争议区域与争议区域的整合。

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