python - 如何使用 cartopy 绘制 SRTM 源栅格？

Question

我一直在研究 cartopy 提供的选项来制作快速自动 map 。我发现 SRTMSource 非常方便。但是，我还没有找到如何使用不同的投影(例如墨卡托投影)来绘制数据...

当我这样做时:

import cartopy.crs as ccrs
import matplotlib.pyplot as plt
from cartopy.io.srtm import SRTM3Source

source = SRTM3Source

plt.figure()
ax = plt.axes(projection=ccrs.PlateCarree())
ax.add_raster(source())
ax.set_extent([12, 13, 47, 48])
plt.show()()

它工作得很好，可以使用 PlateCarree 投影绘图

但是当我尝试将轴的投影更改为墨卡托时:

import cartopy.crs as ccrs
import matplotlib.pyplot as plt
from cartopy.io.srtm import SRTM3Source

source = SRTM3Source

plt.figure()
ax = plt.axes(projection=ccrs.Mercator())
ax.add_raster(source())
ax.set_extent([12, 13, 47, 48])
plt.show()()

我收到一条错误消息:

Unsupported projection for the SRTM3 source.

那么有没有办法保持 SRTMSource 的功能(自动下载)并使用任何投影绘制栅格？

我想简单地做一个 imshow 并更改投影，但是由于在绘图时检索了 RasterSource 的数据，我不知道该怎么做...

plt.figure()
ax = plt.axes(projection=ccrs.Mercator())
ax.imshow(source(), transform=ccrs.Mercator())
ax.set_extent([12, 13, 47, 48])
plt.show()

返回类型错误。

Answer 1

我知道这是一篇旧帖子，但我想了解更多关于 cartopy 的信息，我想我现在可以回答这个问题，也许这可以帮助其他有类似问题的人。

好的，我们开始吧。首先cartopy主页上对transform和projection的区别有很好的解释: https://scitools.org.uk/cartopy/docs/latest/tutorials/understanding_transform.html?highlight=understanding

projection关键字决定了plot的样子，所以这里针对这个问题，我们把projection改成Mercator:

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

你可以选择任何你喜欢的投影，但你必须添加 transform 关键字(在这种情况下这仍然是 PlateCarree，因为 DATA 在 PlateCarree 中提供):

ax.set_extent([12, 13, 47, 48],transform=ccrs.PlateCarree())

所以正确的代码应该是:

import cartopy.crs as ccrs
import matplotlib.pyplot as plt
from cartopy.io.srtm import SRTM3Source

source = SRTM3Source

plt.figure()
ax = plt.axes(projection=ccrs.Mercator())
ax.add_raster(source())
ax.set_extent([12, 13, 47, 48], transform=ccrs.PlateCarree())
plt.show()()

我希望这是正确的并且我的解释对您有所帮助。