Matplotlib Basemap圆形无法跨球面另一半及投影边缘异常问题求助

嘿，这个问题我之前踩过坑！Basemap在处理跨180°经线（国际日期变更线）的圆形时，确实会出现这种“被投影边界截断”的情况——本质是它默认把经度范围限制在-180°到180°之间，超过这个范围的坐标点会被直接丢弃，导致圆形无法延伸到球面的另一半。下面给你几个实用的解决办法：

1. 手动拆分圆形为两段圆弧 #

最直接的思路是把跨180°的圆拆成两部分：一部分在0°_{180°区间，另一部分转换为-180°}0°区间（也就是把大于180°的经度减去360°），然后分别绘制。示例代码如下：

from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt
import numpy as np

# 初始化圆柱投影的Basemap
m = Basemap(projection='cyl', llcrnrlon=-180, urcrnrlon=180, 
            llcrnrlat=-90, urcrnrlat=90)
m.drawcoastlines()

# 目标圆心和半径（单位：度）
center_lon, center_lat = 180, 0
radius_deg = 10

# 生成圆形的原始坐标点
theta = np.linspace(0, 2 * np.pi, 100)
lons = center_lon + radius_deg * np.cos(theta)
lats = center_lat + radius_deg * np.sin(theta)

# 拆分坐标：把>180°的经度转成-180°~0°范围
mask = lons > 180
# 第一部分：经度在-180~180之间的点
lon_part1, lat_part1 = lons[~mask], lats[~mask]
# 第二部分：经度>180的点，减去360转成负经度
lon_part2, lat_part2 = lons[mask] - 360, lats[mask]

# 分别绘制两段圆弧
x1, y1 = m(lon_part1, lat_part1)
m.plot(x1, y1, 'blue', linewidth=2)
x2, y2 = m(lon_part2, lat_part2)
m.plot(x2, y2, 'blue', linewidth=2)

plt.show()

这样处理后，圆形就能在180°经线两侧完整显示了。

2. 用大圆绘制替代（仅适用于球面大圆） #

如果你的“圆形”是球面上的大圆（比如等角距离的圆弧），可以直接用Basemap的drawgreatcircle()方法，同样需要拆分两段绘制：

# 接上面的Basemap初始化代码
# 绘制跨180°的大圆圆弧
m.drawgreatcircle(170, 0, -170, 0, delon=10, color='red', linewidth=2)
m.drawgreatcircle(-170, 0, 170, 0, delon=10, color='red', linewidth=2)

不过这个方法只适合大圆场景，如果你需要的是固定半径的小圆圈，还是第一种方法更通用。

3. 切换到Cartopy（长期推荐方案） #

Basemap已经被官方标记为弃用了，现在更推荐用Cartopy来处理地理可视化——它对跨180°经线的图形支持非常友好，不需要手动拆分就能自动渲染完整的圆形：

import cartopy.crs as ccrs
import matplotlib.pyplot as plt

# 初始化PlateCarree投影
ax = plt.axes(projection=ccrs.PlateCarree())
ax.coastlines()
ax.set_global()

# 直接绘制跨180°的圆形，Cartopy会自动处理边界
circle = plt.Circle((180, 0), 10, color='green', alpha=0.5, 
                    transform=ccrs.PlateCarree())
ax.add_artist(circle)

plt.show()

这段代码里的圆形会自动在180°经线两侧显示完整的形状，完全不需要手动调整坐标，省心很多。

为什么会出现这个问题？ #

Basemap的投影坐标范围默认被限制在-180°~180°经度之间，当你生成的圆形坐标点超过180°时，这些点会被判定为超出投影范围而被截断，所以看起来像是“碰到了边界”，但实际上是球面的另一侧。通过坐标转换拆分图形，或者换用更现代的Cartopy库，就能完美解决这个问题。

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