在Python中将世界数据内插到局部网格

Question

1. 我有一个NC文件，它具有高分辨率(约25分钟)的全局数据，因此有很多坐标。
2. 我有一些区域的测量值坐标间隔不规则。
3. 为了在两者之间进行数学运算，我的目标是将全局数据插值到局部不规则网格中，而不是相反，以最小化不确定性。
4. 我该如何做到这一点?我想到了Scipy的griddata模块，但似乎无法启动这个过程。非常感谢。

编辑:

我尝试了下面的方法，虽然有点耗费资源，但我认为它是有效的。

从scipy.interpolate导入interp2d f = interp2d(长w,latw,Zw)

其中longw、latw、Zw分别是NETCDF文件和2D变量的世界坐标。

然后我申请了f(经度、纬度)对我所有的不规则观察y迭代插值。lon和lat是局部区域的观测坐标。

Answer 1

正如@gansub暗示的那样，通过定义输出网格也可以获得更好的结果。看看这个非常相似的问题:如何在Python中插入分散的数据到规则网格?

在评论之后，我意识到我错过了不规则的目标网格。也许这个答案会有所帮助:https://stackoverflow.com/questions/3242382/interpolation-over-an-irregular-grid