同CRS的两个GeoDataFrame无法拼接或合并的问题排查

排查与修复方案 #

• 强制统一CRS定义格式
  即便都是WGS84，若一个是EPSG代码格式（EPSG:4326）、一个是文字描述（WGS 84），geopandas可能无法识别为同一CRS。显式将两个数据集的CRS设置为EPSG代码：

  edges_with_maxspeed = edges_with_maxspeed.set_crs(epsg=4326, allow_override=True)
  edges_without_maxspeed = edges_without_maxspeed.set_crs(epsg=4326, allow_override=True)
  

• 规范使用pd.concat封装GeoDataFrame
  拼接后重新封装时，明确指定几何列和CRS：

  combined_df = pd.concat([edges_with_maxspeed, edges_without_maxspeed], ignore_index=True)
  edges = gpd.GeoDataFrame(combined_df, geometry='geometry', crs=4326)
  

• 清理无效/空几何数据
  空值或无效几何会干扰CRS检测，先过滤异常行：

  # 保留有效几何行
  edges_with_maxspeed = edges_with_maxspeed[edges_with_maxspeed.geometry.is_valid]
  edges_without_maxspeed = edges_without_maxspeed[edges_without_maxspeed.geometry.is_valid]
  # 删除几何为空的行
  edges_with_maxspeed = edges_with_maxspeed.dropna(subset=['geometry'])
  edges_without_maxspeed = edges_without_maxspeed.dropna(subset=['geometry'])
  

• 调整geopandas版本
  部分版本存在CRS字符串匹配bug，尝试升级到最新稳定版：

  pip install --upgrade geopandas
  

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