如何使用IMSHOW在地图上显示图像

Question

我有一个具有 9 种不同颜色值 （0-9） 的 geotif，并希望在地图上显示它。我正在尝试将其与 matplotlib 包中的底图一起使用。

from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt
im = plt.imread('a.tif')
extent = [-18, 52, -38, 38] # [left, right, bottom, top]
m = Basemap(projection='merc', llcrnrlon=extent[0], urcrnrlon=extent[1], llcrnrlat=extent[2], urcrnrlat=extent[3], resolution='c')
m.drawcountries() 
plt.imshow(im, extent=extent, alpha=0.6)
plt.show()

我只得到一张黑白图像。当我取消注释 drawcountries（） 时，我看到了来自 tif 的数据。

如何在地图上绘制tif的颜色并添加国家边界？

Answer 1

m.imshow(im, extent=extent, alpha=0.6)

Answer 2

使用范围并不能解决您的问题！它会产生一个更大的问题，因为它看起来没问题，但它是错误的。

该范围仅适用于圆柱投影和保持距离。它只是线性地将图像调整到某个框！墨卡托投影不是线性的。

尝试使用投影"cyl"。所以，答案是：非线性变换没有办法。可能的解决方法是将图像坐标转换为新投影中的距离，然后将 2D 插值（例如 matplotlib.tri）转换为某个新矩形（距离）。