如何在NorthPolarStereo投影中添加更多经度/纬度线

Question

我正在尝试增加在Python的NorthPolarStereo投影上出现经度线的间隔。我尝试增加/定义ticks / lim（对于x和y坐标），并且返回错误“具有多个元素的数组的真值不明确。请使用a.any（）或a.all（）” 。我还尝试过更改范围的界限。我找到了一篇有关使用底图的文章，但我认为不需要整体使用。我用于投影的基本代码的链接为here。和代码片段在这里：

# Generate axes, using Cartopy, drawing coastlines, and adding features
fig=plt.figure(figsize=(10,10))

projection = ccrs.NorthPolarStereo()
ax = plt.axes(projection=projection)
ax.coastlines(linewidths=0.5)


ax.add_feature(cfeature.LAND, facecolor='lightgray')
ax.set_extent([-180,180,90,0],ccrs.PlateCarree())
gl = ax.gridlines(crs=ccrs.PlateCarree(),linestyle="--", linewidth=1, color='k', alpha=0.5)


theta = np.linspace(0, 2*np.pi, 100)
center, radius = [0.5, 0.5], 0.5
verts = np.vstack([np.sin(theta), np.cos(theta)]).T
circle = mpath.Path(verts * radius + center)

ax.set_xticks([180,150,120,90,60,30])
ax.set_boundary(circle, transform=ax.transAxes)
p=wrap_U.plot.contour(ax=ax,vmin=-8,vmax=16, transform=ccrs.PlateCarree(),
                        levels = 7, linewidths=0.5, cmap='k')

ax.clabel(p, np.arange(-8,16,8),fmt='%d', inline=1, fontsize=14)


gvutil.set_titles_and_labels(ax, 
                                 lefttitle="Zonal Wind", righttitle="m/s")

  #This line leads to truth value error  
gvutil.set_axes_limits_and_ticks(ax, xticks=np.linspace(-180, 180), xticklabels=['180S', '90S', '0', '90N', '180N'],
                                        ylim=ax.get_ylim()[::-1], yticks=U)



# Show the plot
plt.show()

我确实设置了自定义的anaconda环境，因此，如果需要更多有关此方面的信息，我也可以链接到查找环境的位置。

Answer 1

似乎您正在尝试设置一组不匹配的xticks和xticklabels

gvutil.set_axes_limits_and_ticks(ax, xticks=np.linspace(-180, 180), xticklabels=['180S', '90S', '0', '90N', '180N'],
                                        ylim=ax.get_ylim()[::-1], yticks=U)

您的代码生成50个标记np.linspace(-180,180)（默认为numpy linspace）

但是只有5个标签xticklabels=['180S', '90S', '0', '90N', '180N']

作为替代，我建议尝试xticks=np.linspace(-180,180,5)或暂时删除特定的xticklabel