如何使用matplotlib在Cartopy极地立体图中获得平滑的网格线？

Question

使用matplotlib和cartopy的NorthPolarStereo投影进行绘图时，恒定纬度网格线（平行线）不是平滑圆。

这是一个最小的工作示例：

from matplotlib import pyplot
import cartopy.crs

ax = pyplot.axes(projection=cartopy.crs.NorthPolarStereo())
ax.set_extent((-180, 180, 60, 90), crs=cartopy.crs.PlateCarree())
ax.gridlines()

pyplot.show()

如您所见，平行线是点之间的一系列直线段。 每条线上没有足够的点使它们看起来光滑。 如何使平行线更圆？

Answer 1

问：如何使平行线更圆？

答：您需要更大的gridlines()创建的gridliner对象的n_steps属性值。 它的默认值是30.这是将值设置为90的相关代码，并生成更好的绘图。

gl = ax.gridlines(draw_labels=False, xlocs=None, ylocs=None)
gl.n_steps = 90

见这里的参考。

样本输出：

Answer 2

这是由于此用例的默认阈值太大。 这不是Cartopy投影界面的用户可控部分。 你可以使用private _threshold属性来破解它，将它除以100对你的用例来说是合理的：

from matplotlib import pyplot
import cartopy.crs

hacked_proj = cartopy.crs.NorthPolarStereo()
hacked_proj._threshold /= 100.
ax = pyplot.axes(projection=hacked_proj)
ax.set_extent((-180, 180, 60, 90), crs=cartopy.crs.PlateCarree())
ax.gridlines()

pyplot.show()

这不是一个很好的解决方案，因为它使用了Stereographic类的实现细节，最好在Cartopy中更好地设置它。 使用此用例在Cartopy Github仓库上打开票证将是一个很好的开始https://github.com/SciTools/cartopy/issues 。