Cartopy 在地图上绘制高/低海平面压力

Question

我正在从底图迁移到 cartopy。我想做的一件事是在地图上绘制高/低压，例如在底图中。这个页面上有一个很好的例子来说明如何做到这一点：https://matplotlib.org/basemap/users/examples.html（“绘制带有高点和低点标签的海平面压力天气图”）。我不打算从这个站点复制和粘贴代码，但想知道如何在 cartopy 中做同样的事情。我无法理解的主要事情是如何在cartopy中做m.xmax and x > m.xmin and y < m.ymax and y > m.ymin（我想像某种矢量变换。

我仔细看了看，并没有看到这个特定的例子被翻译成与 cartopy 兼容的东西。欢迎任何帮助！

Answer 1

要使用 cartopy 编写等效程序，您需要能够翻译两个概念。首先是找到投影的范围，这可以通过GeoAxes 的get_extent() 方法来完成：

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

my_proj = ccrs.Miller(central_longitude=180)
ax = plt.axes(projection=my_proj)
xmin, xmax, ymin, ymax = ax.get_extent()

还需要将坐标点从地理坐标转换为投影坐标，这是坐标参考系实例的transform_points()方法的功能：

import numpy as np

lons2d, lats2d = np.meshgrid(lons, lats) # lons lats are in degrees
transformed = my_proj.transform_points(ccrs.Geodetic(), lons2d, lats2d)
x = transformed[..., 0]  # lons in projection coordinates
y = transformed[..., 1]  # lats in projection coordinates

现在您可以使用与底图示例中相同的技术来过滤和绘制点，而不是 m.xmin 您使用 xmin 等。

当然还有其他替代方法，它们相对于底图示例各有利弊。如果你想出一些好的东西，你可以把它贡献给 Cartopy 画廊。