反距离加权水平内插，附matlab代码(ERA5和GNSS站点不并址的处理方法之水平补偿)

1.内插方法

我在学习过程，内插方法为反距离加权水平内插，分享我的方法和公式，以及matlab代码。

2.使用该内插法的原因

GNSS与ERA5格网位置不并址，需要进行水平方向和垂直方向的补偿的补偿获得。水平方向不并址如第3节图所示，使用的是反距离加权水平内插得到GNSS站点P的气象数据。

3.具体方法

假定GNSS站经纬度为P(φ,λ)，附近四个格网点上的气象参数依次为 P1、P2、P3 和 P4，如下图所示：

4.测站处的气象参数P可表示为：

w1，w2，w3和w4为权重。

为P到P1、P2、P3 和 P4的距离（j=1、2、3、4），为地球半径。

其中用下面公式求得：

5.详细免费的代码请关注公众号WZZHHH，下面是部分代码展示。

functionPdata=interpo_Inverse(lat_lon,lon,lat,P)%反距离加权水平内插获得GNSS站的气象数据%假定GNSS站经纬度为P(lon,lat)%四个格网点上的气象参数依次为P1、P2、P3和P4P1=P(1,1);P2=P(1,2);P3=P(1,3);P4=P(1,4);%lat_lon为P1、P2、P3和P4的经纬度，按照行排列%例如P1经纬度为lat_lon(1,1)、lat_lon(1,2)