使用Cartopy绘制跨越国际日期线的线

Question

我试图绘制全球风暴轨迹，但当风暴穿过日期线（经度从~360到~0）时，该线会在绘图空间中一路循环。 Here's what the plot looks like. See the weird straight lines near the top. 这是我的密码：

ax = plt.axes(projection=ccrs.Robinson())
ax.set_global()

ax.coastlines()

for i in range(nstorms-1): #loop through each TC
        bidx = start_idx[i]
        eidx = start_idx[i+1]
        plt.plot(clons[bidx:eidx],clats[bidx:eidx],transform=ccrs.PlateCarree())

如果我尝试将变换更改为大地坐标，它看起来如下：

Answer 1

要绘制跨越日期线的多段线，需要正确清理经度。例如，值359到2应调整为359到362。在下面的演示代码中，sanitize_lonlist()用于在绘制红色之字形线之前对经度值列表进行清理

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

def sanitize_lonlist(lons):
    new_list = []
    oldval = 0
    treshold = 10  # used to compare adjacent longitudes
    for ix,ea in enumerate(lons):
        diff = oldval - ea
        if (ix>0):
            if (diff>treshold):
                ea = ea+360
        oldval = ea
        new_list.append(ea)
    return new_list

ax = plt.axes(projection=ccrs.Robinson()) 
ax.set_global() 
ax.coastlines(alpha=0.3)

# sample long/lat data for demo purposes
# xdateline: list of longitudes that cross dateline several times
xdateline = [347,349,352,358,4,7,8,3,359,358,360,3,5,359,1,357,0,8,12,6,357,349]
# ydateline: list of accompanying latitudes
ydateline = range(len(xdateline))

# plot the line crossing dateline using `sanitized` values of longitudes
plt.plot(sanitize_lonlist(xdateline), ydateline, transform=ccrs.PlateCarree(), color='red') 

plt.show()

使用xdateline的原始值绘制代码行：-

plt.plot(xdateline, ydateline, transform=ccrs.PlateCarree(), color='red')

该地块将为：-