绘制基于纬度/经度值的等高线

Question

我正在尝试使用 mpl_toolkits.basemap 在俄克拉荷马州的地图上绘制天气变量，但在弄清楚如何插入数据以绘制在地图上时遇到问题。

这是我当前代码的总体思路：

lons = [-97.9547, -97.9747, -97.4256]
lats = [35.5322, 35.864, 35.4111]
data = [2,2,2]
map = Basemap(llcrnrlon = -103.068237, llcrnrlat = 33.610045, urcrnrlon = -94.359076, urcrnrlat = 37.040928, resolution = 'i')  


CS = map.contour(X, Y, data)
map.drawstates()
plt.show()

我试图完成的是根据 lons/lats 列表中的相关参考索引在地图上绘制数据值，然后勾勒出数据变量的值。

现在这显然行不通，因为我需要插入数据。有没有办法可以使用 griddata 函数来完成这个任务？鉴于纬度和经度值不是线性间隔的，我对如何建立网格边界感到非常困惑。

有没有我想念的更简单的方法来做到这一点？

任何帮助和/或提示将不胜感激，这阻碍了我继续研究项目的下一个主要部分！

Answer 1

我没有在这台机器上安装 python，所以无法测试。但是这样的事情应该可以为您提供计数图所需的输入......

import numpy as np

lons = [-97.9547, -97.9747, -97.4256]
lats = [35.5322, 35.864, 35.4111]
data = [2,2,2]

xs, ys = np.meshgrid(lons, lats)

dataMesh = np.empty_like(xs)
for i, j, d in zip(lons, lats, data):
    dataMesh[lons.index(i), lats.index(j)] = d

map = Basemap(llcrnrlon = -103.068237, llcrnrlat = 33.610045, urcrnrlon = -94.359076, urcrnrlat = 37.040928, resolution = 'i')  


CS = map.contour(xs, ys, dataMesh)
map.drawstates()
plt.show()

就像我说的那样，我还没有测试过这个。我不知道如果您尝试绘制单元化值会发生什么。您可能需要使用不同的 numpy 数组初始化。