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

WZZHHH

已于 2024-07-30 14:50:35 修改

阅读量712

点赞数 7

分类专栏： ERA5数据处理文章标签： matlab 开发语言 ubuntu 数学建模 windows 数据库

于 2024-01-04 20:00:10 首次发布

本文链接：https://blog.csdn.net/WZZHHH/article/details/135394845

版权

ERA5数据处理专栏收录该内容

8 篇文章 4 订阅

订阅专栏

1.内插方法

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

2.使用该内插法的原因

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

3.具体方法

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

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

w1，w2，w3和w4为权重。

$w_{j}^{*}$ 为P到P1、P2、P3 和 P4的距离（j=1、2、3、4）， $R_{e}$ 为地球半径。

$\varnothing _{j}$ 其中用下面公式求得：

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)