放大 Cartopy 正交图

Question

我正在尝试放大正交投影中的一些 Cartopy 图，但使用简单的 ax.set_extent() 方法似乎是不可能的。 我正在使用的代码：

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

fig = plt.figure()
ax = plt.axes(projection = ccrs.Orthographic())
ax.set_global()
ax.set_extent((-120,120,50,90), crs = ccrs.PlateCarree())

我收到以下错误： ValueError: Failed to determine the required bounds in projection coordinates.

实际上，我只需要设置一个较低的边界纬度来绘制极地区域的图，这是我以前使用 Basemap 所做的。但是我不知道如何继续使用 Cartopy。

有没有办法做到这一点？

Answer 1

这里的问题是您的经度范围太大，因为 -120 或 120 都不在中心经度/纬度为 0 的正交磁盘上（它们在“后面”）。您可以先查询全局经度范围，然后使用它而不是 (-120, 120)。

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

fig = plt.figure()
ax = plt.axes(projection = ccrs.Orthographic())
ax.set_global()

global_extent = ax.get_extent(crs=ccrs.PlateCarree())
ax.set_extent(global_extent[:2] + (50, 90), crs=ccrs.PlateCarree())
ax.coastlines()

Answer 2

我在想象具有受限坐标的所需正交投影的外观方面确实存在一些小问题 - 并且据说 cartopy 有类似的理解问题 :-)。
但是当然你可以使用不同的（并且可能更合适的）投影； LambertConformal 可能有意义。

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

fig = plt.figure()
proj = ccrs.LambertConformal(central_longitude=-96.0, central_latitude=39.0, )
ax = plt.axes(projection = proj)

ax.set_extent((-120,120,50,90), crs = ccrs.PlateCarree())

ax.coastlines(resolution='110m')
ax.gridlines()
plt.show()