cartopy：大圆距离线的更高分辨率

Question

我正在尝试绘制两点之间的大圆距离。我在 cartopy 文档中找到了一个 (introductory_examples/01.great_circle.html):

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

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

ax.set_global()

ax.coastlines()

plt.plot([-0.08, 132], [51.53, 43.17], color='red',      transform=ccrs.Geodetic())
plt.plot([-0.08, 132], [51.53, 43.17], color='blue', transform=ccrs.PlateCarree())

plt.show()

生成下图：

great circle example

问题是，在我自己的工作中，这两个点靠得更近，并且在不同的投影中（尽管我认为这在这里并不重要）。如果我将此代码更改为较小区域中的一条线，如下所示：

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

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

ax.set_extent([-5, 55, 40, 55])

ax.coastlines()

plt.plot([-0.08, 50], [51.53, 43.17], color='red',      transform=ccrs.Geodetic())
plt.plot([-0.08, 50], [51.53, 43.17], color='blue', transform=ccrs.PlateCarree())

plt.show()

这样就形成了下图：shorter line

本例中的红色大圆线看起来很蹩脚，看起来是因为分辨率太低。如何增加构成大圆线的点数？

Answer 1

此问题是由于投影中的硬编码阈值造成的。目前这不是用户可控的参数，但您可以通过定义自己的 sub-class:

来解决这个问题

class LowerThresholdRobinson(ccrs.Robinson):

    @property
    def threshold(self):
        return 1e3

如果您在定义轴时使用 LowerThresholdRobinson() 而不是 ccrs.Robinson()，这应该可以解决问题。

cartopy：大圆距离线的更高分辨率

cartopy: higher resolution for great circle distance line

python

matplotlib

cartopy