python - 在matplotlib中将经/纬度转换为x/y以进行imshow

真是愚蠢的问题。我有一个全局网格数据集，我想将其子集到某个纬度/经度范围（50N-80N，185E-230E）并将其放入墨卡托投影中。由于我正在使用imshow，因此我需要正确设置范围。我找不到我需要的示例。我显然不明白该如何设置。任何帮助表示赞赏。

mapproj = Basemap(projection = 'merc', resolution = 'l', lat_0=65,lon_0=207.5)

x0, y0 = mapproj(185.0,50.0)
x1, y1 = mapproj(230.0,80.0)

myplot = plt.imshow(datamasked.T, interpolation = 'nearest', cmap = my_cmap, \
norm = matplotlib.colors.Normalize(vmin = 0.0, vmax = 10.0, clip = False), \
extent = (x0, x1, y0, y1))

最佳答案

这看起来不错，唯一的事情是您需要使用底图imshow方法。也就是说，应使用plt.imshow()而不是mapproj.imshow()

这是一个完整的工作示例：

from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt
import numpy as np

m = Basemap(llcrnrlon=-10,llcrnrlat=35,urcrnrlon=35,urcrnrlat=60.,
             resolution='i', projection='tmerc', lat_0 = 48.9, lon_0 = 15.3)

m.drawcoastlines()

#Madrid
x,y = m([-3.703889],[40.4125])
m.plot(x,y, marker="o", color="blue", label="Madrid", ls="")

# image background
x0, y0 = m(-10, 35)
x1, y1 = m(35, 60)
extent= (x0, y0, x1, y1) # lonmin, lonmax, latmin, latmax
m.imshow(np.random.rand(23,56), interpolation = 'nearest',
          extent = extent, alpha=0.4, cmap="terrain")


plt.legend()
plt.show()

关于python - 在matplotlib中将经/纬度转换为x/y以进行imshow，我们在Stack Overflow上找到一个类似的问题：https://stackoverflow.com/questions/42933002/