反距离加权水平内插，附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，下面是部分代码展示。

function?Pdata?=?interpo_Inverse(lat_lon,lon,lat,P)

%?反距离加权水平内插获得GNSS站的气象数据
%?假定GNSS站经纬度为P(lon,lat)

%?四个格网点上的气象参数依次为?P1、P2、P3?和?P4
P1?=?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)