Python中高斯过程/克里格插值法

Question

我试图在地图上插值风速和风向值，给出一些经纬坐标，然后将这些值与我的观测值进行比较。一些论文指出高斯过程/克里格法是解决这一问题的有效方法，但我对数学的理解不够好，无法直接实现它们的模型。在

我的数据框看起来像：

lons_   lats_   U2M_    V2M_
.
.
.

我可以参考SciPy上的示例，但不确定是否可以使用它们的方法，因为我正在尝试插值向量-风速和风向都有u&v分量：

^{pr2}$

其中U2M。。。。只是pd系列（…）的各行u、v分量。我可以分别对每个分量进行计算，然后从插值的u，v值重建速度和方向值吗？我想我的问题是，这是一个数学上合理的过程吗？在

Answer 1

风速可以直接从观测值中插值，但风向更复杂。标准的方法是转换为向量-计算了u和v组件，就完成了一半。在使用给定的公式转换回ws和dir变量之前，可以对这两个值进行平均和插值。这方面的背景比我在这本非常有用的资源中解释得更好：

田庄，斯图尔特。（2014年）。技术说明：平均风速和风向。10.13140/RG.2.1.3349.2006。 https://www.researchgate.net/publication/262766424_Technical_note_Averaging_wind_speeds_and_directions