我正在使用cartopy来生成具有立体投影的北极地图,然后在顶部绘制一条线(以显示横截面的位置).如果我使用下面的代码,那么这条线不是直线穿过极点,而是沿着一条纬线.
import cartopy.crs as ccrs
import matplotlib.pyplot as plt
x=[180,0]
y=[50,50]
ax = plt.axes(projection=ccrs.NorthPolarStereo())
ax.set_extent([0, 360, 50, 90], crs=ccrs.PlateCarree())
ax.plot(x,y,transform=ccrs.PlateCarree())
plt.gca().stock_img()
plt.gca().coastlines()
plt.show()
Run Code Online (Sandbox Code Playgroud)
为了解决这个问题,我必须将x和y更改为:
x=[180,180,0,0]
y=[50,90,90,50]
Run Code Online (Sandbox Code Playgroud)
所以北极有两个数据点.有更好的解决方案吗?
编辑:附加图像
谢谢,
蒂姆