[发明专利]一种利用结构保持滤波的遥感图像融合方法

权利要求书

1.一种利用结构保持滤波的遥感图像融合方法，其特征在于包括下述步骤：

第一步，对于包含N个波段MS₁，MS₂，…，MS_N的原始多光谱图像，通过双三次插值将N个波段图像上采样到和全色图像PAN尺寸相同的图像LRM₁，LRM₂，…，LRM_N；

第二步，将图像LRM₁，LRM₂，…，LRM_N进行超球彩色变换，得到变换后的一个强度分量I和N-1个角度分量θ₁，θ₂，…，θ_N-1，

第三步，分别计算强度分量I与全色图像PAN的均值与方差，将强度分量I与全色图像PAN做直方图匹配得到匹配图像式中μ_I与μ_P分别为强度分量I和全色图像PAN的均值，σ_I与σ_P分别为强度分量I和全色图像PAN的标准差；

第四步，使用结构保持滤波器对强度分量I进行滤波，得到滤波输出图像SI，任滤波参数r、k、σ中，r是计算用的滑动窗口大小，该窗口自左上至右下遍历全图，k决定了计算区域协方差的窗口大小，σ控制着滤波过程中的平滑程度；使用结构保持滤波器对匹配图像MP进行两级滤波操作，其中第一级滤波输出图像记为S₁，再将S₁作为第二级滤波的输入，第二级滤波输出图像记为S₂，第一级的滤波参数记为r₁、k₁、σ₁；第二级的滤波参数记为r₂、k₂、σ₂；记SPF_r,k,σ(.)为结构保持滤波函数，SI＝SPF_r,k,σ(I)，

第五步，分别计算I与MP对应的细节图像DI＝I-SI与DP＝MP-S₂；将细节图像DI与DP注入到匹配图像MP中，得到新强度分量I_new＝MP+DI+DP；

第六步，对新强度分量I_new和角度分量θ₁，θ₂，…，θ_N-1作超球彩色反变换，得到N个波段的融合结果HRM₁，HRM₂，…，HRM_N，