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)
