[发明专利]基于非降采样Laplacian金字塔和BEMD的遥感图像融合处理方法

权利要求书

1.一种基于非降采样Laplacian金字塔和BEMD的遥感图像融合处理方法，其特征在于包括以下步骤：

1)对低分辨率多光谱影像进行HIS变换，得到色调H、强度I、饱和度S三个分量；以及采用非降采样Laplacian金字塔对高分辨率全色影像进行J层分解，得到J个细节分量d_j和一个低频分量A_J，其中1≤j≤J；

2)对所述强度分量I进行K层BEMD分解，得到二维内蕴模态函数和余项函数低频分量A_J进行K层BEMD分解，得到二维内蕴模式函数和余项函数其中1≤k≤K；

3)利用二维希尔伯特变换提取和图像的瞬时特征；

4)利用所述瞬时特征，采用加权方法对和进行融合，得到

5)将和强度分量I的余项进行BEMD逆变换，得到高频融合分量bimfH=RKI+Σk=1KbimfkF;]]>

6)将所述高频融合分量结合d_j进行非降采样Laplacian金字塔逆变换，得到融合后图像的强度分量I′，即

7)结合色调H、融合后图像的强度分量I′和饱和度S，进行HIS逆变换，得到融合后的高分辨率多光谱图像。

2.根据权利要求1所述基于非降采样Laplacian金字塔和BEMD的遥感图像融合处理方法，其特征在于步骤1)中采用非降采样Laplacian金字塔对高分辨率全色影像进行J层分解包括：首先将高分辨率全色影像和非降采样滤波器组做卷积，得到原始图像的一层分解高频图像PAN^H和低频图像PAN^L，然后对非降采样滤波器做插值，插值后的非降采样滤波器再和低频图像PAN^L卷积，重复上述步骤，即完成对高分辨率全色影像的J层非降采样Laplacian金字塔分解。

3.根据权利要求1或2所述基于非降采样Laplacian金字塔和BEMD的遥感图像融合处理方法，其特征在于：所述分解层数J由确定，其中，R_pan为高分辨率全色影像的分辨率，R_ms为低分辨率多光谱影像的分辨率。

4.根据权利要求1所述基于非降采样Laplacian金字塔和BEMD的遥感图像融合处理方法，其特征在于，步骤3)中的瞬时特征包括：内蕴模态函数水平、垂直、45度和135度四个方向的瞬时频率ω_i和频率能量E_ω，Eω(x,y)=14Σi=14ωi2,]]>其中1≤i≤4。

5.根据权利要求4所述基于非降采样Laplacian金字塔和BEMD的遥感图像融合处理方法，其特征在于：以ω_i和E_ω作为高分辨率全色影像的近似分量A_J与低分辨率多光谱影像的强度分量I的融合标准，采用加权方法对和图像的高频部分进行融合：

D_F(x，y)表示融合后的像素值，D_A(x，y)表示第一幅待融合影像的像素值，D_B(x，y)表示第二幅待融合影像的像素值，表示第一幅影像的频率能量，表示第二幅影像的频率能量；

当频率能量和像素点D_A(x，y)、D_B(x，y)的绝对值的大小关系不一致时，权系数α和β选取的规则如下：

(1)当|DA|+EωA+|DB|+EωB≠0]]>时，令α′=|DA|+EωA|DA|+EωA+|DB|+EωB,]]>β′＝1-α′，当|DA|+EωA+|DB|+EωB=0,]]>令α′=12,]]>β′＝1-α′；

(2)当D_A(x，y)和D_B(x，y)同号时，α＝α′，β＝β′；当D_A(x，y)和D_B(x，y)异号时：

如果|D_A(x，y)|≥|D_B(x，y)|，则α＝α′，β＝-β′；

如果|D_A(x，y)|＜|D_B(x，y)|，则α＝-α′，β＝β′。