使用Python Bokeh+Geopandas投影数据遇白令海峡显示问题求助

我之前处理过类似的俄罗斯远东地区投影错位问题，其实核心原因是墨卡托投影（EPSG:3857）对跨180°经线的几何体兼容性不好——俄罗斯的白令海峡区域刚好横跨了180°经线，直接转投影会让这部分几何体被错误地“拉”到地图的另一侧，导致填充颜色错位。不需要更换投影方式（因为Bokeh瓦片地图默认就是EPSG:3857），只需要先在WGS84（EPSG:4326）坐标系下处理跨经线的几何体，再转投影即可。

具体解决步骤： #

1. 拆分跨180°经线的几何体
   在EPSG:4326坐标系下，把横跨180°经线的多边形拆分成东、西两个部分，并将西半球的部分经度加360°，这样转墨卡托投影时就能正确映射到地图的右侧（白令海峡的正确位置）。

   代码示例：

   import geopandas as gpd
   from shapely.geometry import LineString
   from shapely.ops import split
   from shapely.affinity import translate
   
   # 创建180°经线的分割线
   split_line = LineString([(180, -90), (180, 90)])
   
   def process_cross_dateline_geom(geom):
       # 检查几何体是否跨180°经线
       if geom.intersects(split_line):
           # 拆分几何体
           split_geometries = list(split(geom, split_line).geoms)
           processed = []
           for g in split_geometries:
               # 根据质心判断几何体所在半球
               centroid = g.centroid
               if centroid.x < 0:
                   # 西半球部分，经度+360°，让它在转投影后处于地图右侧
                   g = translate(g, xoff=360)
               processed.append(g)
           return processed
       else:
           return [geom]
   
   # 应用处理函数到GeoDataFrame的geometry列
   geopdwd['geometry'] = geopdwd['geometry'].apply(process_cross_dateline_geom)
   # 拆分后的几何体展开为单独行
   geopdwd = geopdwd.explode('geometry', ignore_index=True)
   

2. 转换到EPSG:3857投影
   处理完跨经线的几何体后，再执行投影转换：

   geopdwd = geopdwd.to_crs(epsg=3857)
   

原理说明： #

EPSG:3857的坐标范围对应经度-180°到180°，直接转换跨180°的几何体时，投影算法会把跨线的部分“折叠”到地图的另一侧。通过拆分并调整西半球几何体的经度，我们让这部分在转投影后处于正确的位置，避免了错位问题。

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