在 Cartopy 中查看 Alaskan/East 北极的西伯利亚一侧

View Alaskan/East Siberian Side of Arctic in Cartopy

我使用以下代码在 Cartopy 的 NorthPolarStero 地图上绘制了一些点

fig, ax = plt.subplots(1,1,figsize=(8,8),  subplot_kw={'projection': ccrs.NorthPolarStereo()})

ax.scatter([x[0] for x in max_coords[:-1]],
           [x[1] for x in max_coords[:-1]],
           color='k',
           s=50,
           marker='x',transform=ccrs.PlateCarree(),zorder=5,
           label='BSH Center 1980-2020')

ax.set_extent([-180, 180, 90, 66], ccrs.PlateCarree())

ax.add_feature(cartopy.feature.LAND, edgecolor='black',zorder=1)

ax.legend(fontsize='x-large')

但实际上我只想绘制这张图片的上半部分,如下所示:

如果我将 ```ax.set_extent`` 行更改为:

[90, -90, 90, 66] 要么 [-90, 90, 90, 66]

它只给出了情节的下半部分,如下所示:

有谁知道如何获得上半部分?我也试过将 NorthPolarStereo 线的 central_longitude 关键字设置为 +/- 180,但它只是显示错误的一半。

要设置绘制所需区域的范围,您应该使用正在使用的投影坐标ccrs.NorthPolarStereo()。您需要的相关代码是:–

ax.set_extent([-2.633e+06, 2.696e+06, -6e+04, 2.9e+06], crs=ccrs.NorthPolarStereo()) 

你应该得到类似这样的情节: