用正投影法再现cartopy地块的范围答案

【问题标题】：Reproducing extent of cartopy plots with orthographic projection用正投影法再现cartopy地块的范围
【发布时间】：2021-10-19 22:53:33
【问题描述】：

我正在生成一组不同的地图，我希望它们都具有相同的轴限制。但是，我很难弄清楚如何设置范围。

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

lon = -90
lat = 90
proj =  ccrs.Orthographic

projection = proj(lon, lat)
fig = plt.figure()
ax = fig.add_subplot(1, 1, 1, projection = projection)

transform = ccrs.Geodetic()

lat0 = 50
lon0 = 50

for x in range(-5,6):
    plt.plot(lon0 + x, lat0, 'bo', transform = transform)
    plt.plot(lon0, lat0 + x, 'bo', transform = transform)

我希望能够提取“地图”的范围，然后将其应用于以不同方式自动缩放的其他地图。但是，即使我尝试提取范围并将其重新应用于同一个图形，它也会失败。

a = ax.get_extent(crs = proj())
#a = (2769836.95350539, 3487499.8040009444, 4467034.2547478145, 5254224.255689873)
ax.set_extent(a, crs = proj())

【问题讨论】：

标签： python cartopy

【解决方案1】：

一旦我意识到这一点就很明显了，但是需要在范围描述中指定投影的lon 和lat。

a = ax.get_extent(crs = proj(lon, lat))
ax.set_extent(a, crs = proj(lon, lat))

【讨论】：