[发明专利]基于自动调节引导滤波器参数的遥感图像融合方法

权利要求书

1.一种基于自动调节引导滤波器参数的遥感图像融合方法，其特征在于，包括如下步骤：

(1)获取待融合遥感图像：

将从数据库中选取对同一目标拍摄的一幅多光谱MS图像和一幅全色图PAN图像作为待融合遥感图像；

(2)对待融合遥感图像进行多尺度变换：

(2a)对多光谱MS图像进行RGB-IHS变换，得到多光谱MS图像的I分量、H分量和S分量；

(2b)对多光谱MS图像的I分量和全色图PAN图像分别进行多尺度变换，得到多光谱MS图像I分量的高频系数和低频系数，以及全色图PAN图像的高频系数和低频系数；

(3)获取待融合遥感图像的低频权重图：

将多光谱MS图像I分量的低频系数和全色图PAN图像的低频系数基于每个像素进行比较，得到待融合遥感图像的低频权重图W_1L和W_2L；

(4)获取待融合遥感图像的高频权重图：

将多光谱MS图像I分量的高频系数和全色图PAN图像的高频系数基于每个像素进行比较，得到待融合遥感图像的高频权重图W_1H和W_2H；

(5)设置引导滤波器调整参数ε_l1和ε_l2，以及计算引导滤波器窗口半径r的迭代次数t：

(5a)设置用于对低频权重图W_1L和W_2L滤波的引导滤波器的调整参数为ε_l1；

(5b)设置用于对高频权重图W_1H和W_2H滤波的引导滤波器的调整参数为ε_l2；

(5c)设引导滤波器引导图像行方向上的像素个数为width，引导滤波器引导图像列方向上的像素个数为height，则计算引导滤波器窗口半径r的迭代次数其中min为取小值运算；

(6)获取待融合遥感图像的重定义低频权重图：

令引导滤波器的窗口半径r＝1，将待融合遥感图像的低频权重图W_1L作为引导滤波器的输入图像，将多光谱MS图像I分量的低频系数作为引导滤波器的引导图像，对输入图像和引导图像进行滤波，得到待融合遥感图像的重定义低频权重图M_1L，同时将待融合遥感图像的低频权重图W_2L作为引导滤波器的输入图像，将全色图PAN图像的低频系数作为引导滤波器的引导图像，对输入图像和引导图像进行滤波，得到待融合遥感图像的重定义低频权重图M_2L；

(7)利用待融合遥感图像的重定义低频权重图M_1L、M_2L对与其对应的低频系数分别进行加权融合：

利用待融合遥感图像的重定义低频权重图M_1L、M_2L对与其对应的低频系数分别进行加权融合，得到融合后遥感图像的低频系数L_N；

(8)获取待融合遥感图像的重定义高频权重图：

令引导滤波器的窗口半径r＝1，将待融合遥感图像的高频权重图W_1H作为引导滤波器的输入图像，将多光谱MS图像I分量的高频系数作为引导滤波器的引导图像，对输入图像和引导图像进行滤波，得到待融合遥感图像的重定义高频权重图M_1H，同时将待融合遥感图像的高频权重图W_2H作为引导滤波器的输入图像，将全色图PAN图像的高频系数作为引导滤波器的引导图像，对输入图像和引导图像进行滤波，得到待融合遥感图像的重定义高频权重图M_2H；

(9)利用待融合遥感图像的重定义高频权重图M_1H、M_2H对与其对应的高频系数分别进行加权融合：

利用待融合遥感图像的重定义高频权重图M_1H、M_2H对与其对应的高频系数分别进行加权融合，得到融合后遥感图像的高频系数H_N；

(10)获取遥感融合图像：

(10a)对融合后遥感图像的低频系数L_N和高频系数H_N进行多尺度逆变换，得到多光谱MS图像I分量和全色图PAN图像的过渡遥感融合图像；

(10b)对过渡遥感融合图像，以及其对应的多光谱MS图像的H分量和S分量进行IHS-RGB变换，得到多光谱MS图像和全色图PAN图像的遥感融合图像；

(11)获取最佳遥感融合图像：

令引导滤波器的窗口半径r＝2,3,...,t+1，并执行步骤(6)，将得到的t幅多光谱MS图像和全色图PAN图像的遥感融合图像和步骤(10)中得到遥感融合图像中每幅图像的综合评价参数C_r最大值对应的遥感融合图像，作为最佳遥感融合图像。