[发明专利]一种卫星测高重力数据与重力卫星数据融合方法

权利要求书

1.一种卫星测高重力数据与重力卫星数据融合方法，其特征在于，包括以下步骤：

(1)将重力卫星数据进行插值处理，使其和卫星测高重力数据的空间分辨率保持一致；

(2)对卫星测高重力数据进行2层小波分解，得到一个低频系数分量aLL、两个水平高频系数分量aHL1和aHL2、两个垂直高频系数分量aLH1和aLH2、以及两个对角高频系数分量aHH1和aHH2；

(3)对重力卫星数据进行2层小波分解，得到一个低频系数分量sLL、两个水平高频系数分量sHL1和sHL2、两个垂直高频系数分量sLH1和sLH2、以及两个对角高频系数分量sHH1和sHH2；

(4)提取卫星测高重力数据分解结果中的低频系数分量aLL，作为融合数据的低频系数分量fLL，即fLL＝aLL；

(5)将分解结果中，卫星测高重力数据的第一层高频系数分量aHL1、aLH1、aHH1和重力卫星数据的第一层高频系数分量sHL1、sLH1、sHH1分别进行加权融合，得到融合数据的第一层高频系数分量fHL1、fLH1、fHH1，即

fHL1＝x₁·aHL1+y₁·sHL1，fLH1＝x₁·aLH1+y₁·sLH1，fHH1＝x₁·aHH1+y₁·sHH1，

其中，x₁和y₁分别为卫星测高重力数据和重力卫星数据第一层加权融合的权重因子；

(6)将分解结果中，卫星测高重力数据的第二层高频系数分量aHL2、aLH2、aHH2和重力卫星数据的第二层高频系数分量sHL2、sLH2、sHH2分别进行加权融合，得到融合数据的第二层高频系数分量fHL2、fLH2、fHH2，即

fHL2＝x₂·aHL2+y₂·sHL2，fLH2＝x₂·aLH2+y₂·sLH2，fHH2＝x₂·aHH2+y₂·sHH2

其中，x₂和y₂分别为卫星测高重力数据和重力卫星数据第二层加权融合的权重因子；

(7)利用步骤(4)中得到的融合数据的低频系数分量fLL，步骤(5)中得到的融合数据的第一层高频系数分量fHL1、fLH1、fHH1，以及步骤(6)中得到融合数据的第二层高频系数分量fHL2、fLH2、fHH2，使用这些系数分量进行小波重构，得到海洋重力融合数据；

其中，aLL表示卫星测高重力数据两层分解结果的低频系数分量；

aHL1表示卫星测高重力数据两层小波分解结果的第一层水平高频系数分量；

aLH1表示卫星测高重力数据两层小波分解结果的第一层垂直高频系数分量；

aHH1表示卫星测高重力数据两层小波分解结果的第一层对角高频系数分量；

aHL2表示卫星测高重力数据两层小波分解结果的第二层水平高频系数分量；

aLH2表示卫星测高重力数据两层小波分解结果的第二层垂直高频系数分量；

aHH2表示卫星测高重力数据两层小波分解结果的第二层对角高频系数分量；

sLL表示重力卫星数据两层分解结果的低频系数分量；

sHL1表示重力卫星数据两层小波分解结果的第一层水平高频系数分量；

sLH1表示重力卫星数据两层小波分解结果的第一层垂直高频系数分量；

sHH1表示重力卫星数据两层小波分解结果的第一层对角高频系数分量；

fLL表示融合数据的低频系数分量；

fHL1表示融合数据的第一层水平高频系数分量；

fLH1表示融合数据的第一层垂直高频系数分量；

fHH1表示融合数据的第一层对角高频系数分量；

fHL2表示融合数据的第二层水平高频系数分量；

fLH2表示融合数据的第二层垂直高频系数分量；

fHH2表示融合数据的第二层对角高频系数分量。

2.根据权利要求1所述的一种卫星测高重力数据与重力卫星数据融合方法，其特征在于：所述步骤(1)中插值处理的方法是三次多项式插值的方式。