Question

我试图绘制两点之间的大圆距离。我在cartopy docs中找到了一个（introductionctory_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

此问题是由投影中的硬编码阈值引起的。目前，这不是用户可控制的参数，但您可以通过定义自己的子类来解决这个问题：

class LowerThresholdRobinson(ccrs.Robinson):

    @property
    def threshold(self):
        return 1e3

如果在定义轴时使用LowerThresholdRobinson()而不是ccrs.Robinson()，则可以消除此问题。

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

1 个答案: