Python将地图绘制到球体上

Question

我有一张想要投影到球体上的地图，我按照另一个线程中的建议使用 cartopy。我在 Jupyter 中编程。

我的代码是：

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

os.chdir(os.getenv("HOME"))
os.chdir("Downloads/")
img = plt.imread("europamap1.jpg")

plt.figure(figsize=(3, 3))

ax = plt.axes(projection=ccrs.Orthographic(-10, 45))
ax.gridlines(color="black", linestyle="dotted")
ax.imshow(img, origin="upper", extent=(-180, 180, -90, 90),
          transform=ccrs.PlateCarree())  # Important

plt.show()

我得到的是一个白色实心圆圈enter image description here

我不知道为什么我看不到任何网格线或图像！我试图 plot.show(img) 并加载了图像！我基本上只是从this thread复制粘贴过来的

Answer 1

您使用plt.imread 加载的图像文件不是地理空间栅格格式，通常是以.tiff 扩展名结尾的GeoTIFF 文件。

问题在于您的图像缺少与cartopy 对齐的坐标参考系 (CRS)。例如print(img.crs) 很可能不会打印 CRS。

要投影您的数据，您需要找到一个 GeoTIFF 文件并将 CRS 与 cartopy 使用的投影对齐 - 例如，ccrs.PlateCarree()。

可以在此处找到一些测试卫星图像：https://gisgeography.com/usgs-earth-explorer-download-free-landsat-imagery/。

有关读取和写入栅格数据，请参阅rasterio:https://rasterio.readthedocs.io/en/latest/。

这是一个使用来自 cmets 中链接的源的 GeoTIFF 的示例。 https://photojournal.jpl.nasa.gov/tiff/PIA03526.tif:.

import matplotlib.pyplot as plt
import rasterio
from rasterio import plot

src = rasterio.open(r"../tests/data/PIA03526.tif")
plot.show(src)