GeoPandas转换坐标系(to_crs)导致Shapefile反转问题求助

将Natural Earth提供的10m分辨率海洋Shapefile转换为以尼日利亚为中心的投影（ESRI:102022）时，出现了“海洋干涸”的反转现象：

原始代码及正常效果 #

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

此时海洋显示为正常蓝色：

转换投影后代码及异常效果 #

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

原因分析 #

• 多边形缠绕方向与投影不兼容：Natural Earth的海洋图层是用一个覆盖全球的多边形表示“海洋区域”（本质是陆地的补集），它的顶点顺序符合WGS84（EPSG:4326）坐标系的逆时针规则。但ESRI:102022是非洲中心的非全局投影，转换后多边形的顶点顺序在该投影下被渲染引擎判定为顺时针，而大多数GIS工具默认逆时针多边形代表内部区域，顺时针代表外部区域，因此原本的海洋区域被识别为外部，反而填充了陆地部分。
• 简单来说：投影坐标系的范围限制导致多边形的“内外”判定逻辑反转了。

解决方法 #

方法1：手动修复多边形缠绕方向 #

通过调整多边形顶点顺序，确保转换后为逆时针方向：

import geopandas as gpd
from shapely.geometry import Polygon

# 读取海洋数据
oceans = gpd.read_file('./ne_10m_ocean')

# 定义修复函数：将顺时针多边形转为逆时针
def fix_winding_direction(geom):
    if isinstance(geom, Polygon):
        # 检查多边形外部环是否为逆时针
        if not geom.exterior.is_ccw:
            # 反转顶点顺序
            return Polygon(list(reversed(geom.exterior.coords)))
    return geom

# 应用修复
oceans['geometry'] = oceans['geometry'].apply(fix_winding_direction)

# 转换投影并绘图
oceans.to_crs("ESRI:102022").plot()

方法2：用buffer(0)自动修复拓扑与缠绕方向 #

buffer(0)是GeoPandas中修复多边形拓扑错误的常用技巧，同时会自动调整缠绕方向：

import geopandas as gpd

oceans = gpd.read_file('./ne_10m_ocean')
# 先转换投影，再修复
oceans_proj = oceans.to_crs("ESRI:102022")
oceans_proj['geometry'] = oceans_proj['geometry'].buffer(0)
# 绘图验证
oceans_proj.plot()

方法3：换用陆地图层做补集（更可靠） #

由于海洋图层本质是陆地的补集，直接使用陆地图层转投影后，通过背景色模拟海洋，反而更不易出问题：

import geopandas as gpd
import matplotlib.pyplot as plt

# 读取陆地图层
land = gpd.read_file('./ne_10m_land')
# 转换投影
land_proj = land.to_crs("ESRI:102022")

# 创建绘图
fig, ax = plt.subplots()
# 设置背景色为海洋蓝
ax.set_facecolor('#87CEEB')
# 绘制陆地
land_proj.plot(ax=ax, color='#8B7355')
plt.show()

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