GeoPandas overlay返回空GeoDataFrame，实际应有交集问题排查

你遇到的空交集问题，大概率是坐标系统设置错误导致的，以下是具体排查和修正步骤：

1. 核心错误：坐标系统初始化错误 #

如果你的CSV里的Coordinate1和Coordinate2是经纬度（WGS84，单位为度），直接设置coords_gdf.crs = "EPSG:26954"完全错误——EPSG:26954是科罗拉多州的平面投影坐标系（单位米），经纬度属于地理坐标系，两者不能直接绑定。这会导致你的点被错误定位，和Superfund图层完全不在同一空间范围，自然没有交集。

修正方法： #

先将点图层的CRS设置为经纬度对应的EPSG:4326，再转换到EPSG:26954：

# 转换坐标时先绑定地理坐标系（假设CSV里是WGS84经纬度）
coords_gdf = gpd.GeoDataFrame(coord, geometry=gpd.points_from_xy(hmsdf.Coordinate2, hmsdf.Coordinate1))
coords_gdf.crs = "EPSG:4326"  # 先设置为经纬度坐标系
coords_gdf = coords_gdf.to_crs("EPSG:26954")  # 再转换到目标投影

2. 简化缓冲区创建代码 #

你当前创建buffer的代码过于繁琐，可以简化为：

# 直接从coords_gdf生成buffer，自动继承CRS
buffer = coords_gdf.buffer(4828.03)
buffer_gdf = gpd.GeoDataFrame(geometry=buffer)

3. 验证图层是否真的重叠 #

分开绘图无法看出是否在同一空间，建议把两个图层画在同一张图里：

import matplotlib.pyplot as plt

fig, ax = plt.subplots()
buffer_gdf.plot(ax=ax, alpha=0.5)
shapefile.plot(ax=ax, color="red", alpha=0.5)
plt.show()

如果图里两个图层没有重叠，说明坐标设置还是有问题，需要确认：

• CSV里的坐标列是否搞反了（比如把纬度当x，经度当y）
• Shapefile原始CRS是否正确，转换后是否和点图层一致（打印shapefile.crs和coords_gdf.crs确认）

4. 其他排查点 #

• 检查几何有效性：运行print(buffer_gdf.is_valid.sum())和print(shapefile.is_valid.sum())，如果有无效几何，可能导致overlay失败
• 用空间关系提前验证：运行any(buffer_gdf.geometry.intersects(shapefile.geometry.unary_union))，如果返回False，说明确实没有空间交集，回到坐标设置环节排查

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