[发明专利]极轨高分辨率卫星全通道资料条纹噪声滤除方法

权利要求书

1.一种极轨高分辨率卫星全通道资料条纹噪声滤除方法，其特征在于，包括如下步骤：

1)由极轨高分辨率卫星资料形成观测亮温矩阵，再利用集合经验模态分解方法分解所述观测亮温矩阵，得到亮温大梯度矩阵C_M×N和剩余亮温矩阵T′_M×N；其中，M为极轨高分辨率卫星资料中单根扫描线的视场总数；N为极轨高分辨率卫星资料中扫描线总数；

2)利用主成分分析方法按下式对所述剩余亮温矩阵T′_M×N进行分解分析，得到各个主成分模态和对应的时间系数：

其中，E为模态矩阵；U为模态系数矩阵；为第i个主成分模态；为第i个主成分模态对应的时间系数；i为主成分模态序号，i＝1，2，...，M；

3)利用集合经验模态分解法对前N_pca个主成分模态的时间系数进行分解，每个分解为由高频到低频的多个本征模态分量；对于分解后的该N_pca个主成分模态的时间系数逐个去除前N_imf个高频本征模态分量，得到平滑后的该N_pca个主成分模态的时间系数；其中，N_pca和N_imf的取值范围均为[1，3]；

4)将步骤3中平滑后的主成分模态的时间系数和未进行分解的主成分模态的时间系数与步骤2)中的各个主成分模态进行重构，得到去除条纹噪声后的剩余亮温矩阵T′_reb，再按下式得到去掉条纹噪声后的观测亮温矩阵T_reb：

T_reb＝T_reb+C_M×N。

2.根据权利要求1所述的极轨高分辨率卫星全通道资料条纹噪声滤除方法，其特征在于：步骤1)中所述观测亮温矩阵为其中，第k个视场沿卫星行进方向的亮温序列记为经集合经验模态分解方法分解后得到的高频项和剩余项分别为和

所述亮温大梯度矩阵C_M×N为所述剩余亮温矩阵T′_M×N为

其中，k为极轨高分辨率卫星资料中视场序号；j为极轨高分辨率卫星资料中扫描线序号；index为大梯度阈值。