GeoPandas 坐标系转换 (to_crs) 反转 Shapefile

Question

将 10m 分辨率海洋形状文件从自然地球转换为不同的投影会“反转”形状文件（“使海洋枯竭”）：

import geopandas as gpd
oceans = gpd.read_file('./ne_10m_ocean')

oceans.plot()

好吧，这里的大海又漂亮又蓝🐟

但是转换为另一个以尼日利亚为中心的投影（ESRI:102022），我们得到：

oceans.to_crs("ESRI:102022").plot()

这是怎么回事？我怎样才能重新填满海洋？

Answer 1

从跨越地图的文物中可以看出，你现在有一些跨越反子午线的特征，这导致你的一些多边形从里到外，因此错误的部分被着色。

您将需要吐出任何穿过“新”反子午线（180+/-25）的多边形 - 尽管事实上您选择的投影仅在非洲大陆上真正有效。请参阅此 GIS Stackexchange 问题，了解 QGis 中的流程示例。