一种基于杜鹃搜索算法的Contourlet变换图像融合方法技术

本发明专利技术公开了一种基于杜鹃搜索算法的Contourlet变换图像融合方法，对配准后的源图像A、B，分别进行Contourlet变换，得到1个低频子带和不同尺度、不同方向的多个高频子带；按照定义的融合规则在所有尺度和方向上对两幅图像的变换系数进行融合，得到融合后的系数；融合图像为R，对于融合后的系数，按照其低频子带和高频子带的顺序，依次进行Contourlet逆变换；最后输出融合后的图像。本发明专利技术与传统的图像融合方法相比，对各尺寸各方向上的融合系数权重进行了加权处理，计算出最优权值，并以此提出一种优化的Contourlet变换图像融合方法。

【技术实现步骤摘要】
一种基于杜鹃搜索算法的Contourlet变换图像融合方法
本专利技术属于智能计算和数字图像处理的交叉应用
，涉及一种Contourlet变换图像融合方法，具体涉及一种基于杜鹃搜索算法的Contourlet变换图像融合方法。
技术介绍
数字图像是人们获取的数字信息的一个重要种类，对于数字图像的处理一直是一个热点。最早的图像融合工作可以追溯到20世纪80年代中期，BurtPJ最早使用拉普拉斯金字塔方法对双筒望远镜图像进行了融合。随着现代影像传感器技术的迅猛发展，影像由单一的可见光模式逐渐发展成多种传感器模式，不同传感器获取的目标信息因其成像方式、使用的电磁波段以及分辨率不同，所反映的目标信息也不同。比如，红外影像、可见光影像，深度影像、多光谱、高光普、高分辨率、合成孔径雷达卫星影像、医学图像中的MRI图和CT图等不同光谱与空间分辨率的图像提供物理特性的互补性，从不同的侧面来揭示着目标本质。融合的概念和技术产生于20世纪80年代，是指在多种信息集成过程中的任一步骤，由不同信息源的信息复合成一种表达形式。影像信息融合是一门综合了传感器、影像处理、信号处理、计算机和人工智能等多种学科的现代高新技术，它是一个对多源传感器(如CCD摄像机、红外探测器、激光探测器、雷达、声呐等)获取的对应同一目标的不同影像数据和其它信息的处理过程，着重于把那些在空间或时间上冗余或互补的多源数据，按一定的规则(或算法)进行运算处理，获得比任何单一数据更精确、更丰富的信息，生成一副具有新的空间、波普、时间特征的合成影像。影像信息融合一般分为三个层次，即像素级融合、特征级融合和决策级融合。像素级影像融合算法有两个关键地方，即选择合适的影像表达工具以及制定合适的影像融合规则。在图像融合中，基于多尺度分析的系数融合规则仍不成熟。系数融合规则是基于多尺度分析的图像融合算法中另一个重要因素，源图像经多尺度变换工具分解后，得到各尺度下的分量，低频子带分量代表源图像的背景信息，高频子带分量则代表了图像空间细节信息和边缘特征，如何选取合理的系数融合规则，将其直接影响融合图像的质量。因此，如今大部分的图像融合算法的研究重点均集中于系数融合规则的制定。因为需要在所有尺度和方向上对两幅图像的变换系数进行融合，得到融合后的系数。融合规则体现在变换后源图像的低频分量和高频分量的处理中，其合理与否直接影响最终的融合图像效果。在传统的图像融合算法中，幅值较大的变换系数通常被认为是灰度值变换显著的边缘、纹理等图像特征，因此最常见的融合规则是采用模值取大，即对比相同子带相同位置的变换系数，取模值较大者作为融合后的系数。或者单一权值的融合，即各尺寸各方向上的融合系数权重为同一个值，且此权值是根据大量实验得到的。而事实上图像的重要特征并不能由单个系数得到有效反映，本专利技术提出一种基于杜鹃搜索算法来计算各尺度各方向的系数加权权重的系数融合规则，并以此提出一种优化的Contourlet变换图像融合方法。
技术实现思路
为了解决上述技术问题，本专利技术提出了一种基于杜鹃搜索算法的Contourlet变换图像融合方法，可以达到提高图像融合质量的目的。本专利技术所采用的技术方案是：一种基于杜鹃搜索算法的Contourlet变换图像融合方法，用于对配准后的源图像A、B进行融合，其中源图像A、B均为灰度图像，且A,B∈RM×N，RM×N是大小为M×N的空间；其特征在于，所述方法包括以下步骤：步骤1：输入两幅源图像A、B，并分别进行Contourlet变换，得到1个低频子带和不同尺度、不同方向的多个高频子带；步骤2：Contourlet系数融合；按照定义的融合规则在所有尺度和方向上对两幅图像的变换系数进行融合，得到融合后的系数；步骤3：融合图像为R；对于融合后的系数，按照其低频子带和高频子带的顺序，依次进行Contourlet逆变换；步骤4：输出融合后的图像。本专利技术的有益效果是：利用杜鹃搜索算法对各尺度方向上的融合系数加权权重的系数进行优化求解，并采用信息熵作为适应度函数对选择的权重系数的效果进行评价，得到更为合理的系数融合规则，进一步完善了Contourlet变换的图像融合方法。本专利技术采用的基于杜鹃搜索算法计算各尺度各方向的系数加权权重的系数融合规则，较之以往简单融合规则的以模值取打，即对比相同子带相同位置的变换系数，取模值较大者作为融合后的系数；或者根据大量实验得到的单一权值的融合，即各尺寸各方向上的融合系数权重为同一个值，可以更好的、更客观的融合两幅源图像的细节信息，同时采用杜鹃搜索算法求解最优的权重系数，可以达到快速图像融合的目的。附图说明图1：本专利技术实施例的流程图.具体实施方式为了便于本领域普通技术人员理解和实施本专利技术，下面结合附图及实施例对本专利技术作进一步的详细描述，需要说明的是杜鹃搜索算法和Contourlet变换是已有优化技术，应当理解，此处所描述的实施示例仅用于说明和解释本专利技术，并不用于限定本专利技术。请见图1，本专利技术提供的一种基于杜鹃搜索算法的Contourlet变换图像融合方法，用于对配准后的源图像A、B进行融合，其中源图像A、B均为灰度图像，且A,B∈RM×N，RM×N是大小为M×N的空间；包括以下步骤：步骤1：输入两幅源图像A、B，并分别进行Contourlet变换，得到1个低频子带和不同尺度、不同方向的多个高频子带；本实施例首先假设对输入图像A、B进行L层Contourlet分解，第l(0≤l≤L-1)层上的方向数为k，则Contourlet变换过程可以表示为：A→(a0,a1,a2,...aL-1,alow)B→(b0,b1,b2,...bL-1,blow)al＝{Al,0,Al,1,Al,2,...Al,k-1},(0≤l≤L-1)bl＝{Bl,0,Bl,1,Bl,2,...Bl,k-1},(0≤l≤L-1)上式中，alow，blow分别为A，B的低频子带，al和bl分别为图像A，B第l层上的高频方向子带集合，Al,k-1为图像A第l层上的高频子带的k-1方向的分量，Bl,k-1为图像B第l层上的高频子带的k-1方向的分量。步骤2：Contourlet系数融合；按照定义的融合规则在所有尺度和方向上对两幅图像的变换系数进行融合，得到融合后的系数；具体实现包括以下子步骤：步骤2.1：根据步骤1进行Contourlet变换，分解层数为L，每层方向数为k，得到1个低频子带和不同尺度、不同方向的L×k个高频子带，则需要计算的各尺度各方向上的权重总数为L×k+1个；初始化杜鹃搜索算法所需的参数，得到算法的初始种群；其中，参数包括迭代次数Iter、算法的种群规模N、鸟蛋被发现的概率pa、速度步长控制参数η、列维飞行步长控制参数λ、鸟巢的初始空间位置；鸟巢的初始空间位置值向量是相应的每一维对应权重集合中的每一个权重；步骤2.2：利用融合后图像的信息熵计算每个鸟巢的适应度函数值，记录并保留当前最优鸟巢位置A；适应度函数计算公式：其中，P(l)为某灰度值l在图像中出现的概率，L为图像的灰度等级；步骤2.3：通过杜鹃搜索算法中的位置更新公式更新鸟巢的位置，得到一组新解并计算每个鸟巢的适应度函数值，比较这组解的适应度函数值并记其具有最优适应度函数值的位置为B；其中位置更新公式为：其中本文档来自技高网...

【技术保护点】
1.一种基于杜鹃搜索算法的Contourlet变换图像融合方法，用于对配准后的源图像A、B进行融合，其中源图像A、B均为灰度图像，且A,B∈RM×N，RM×N是大小为M×N的空间；其特征在于，所述方法包括以下步骤：步骤1：输入两幅源图像A、B，并分别进行Contourlet变换，得到1个低频子带和不同尺度、不同方向的多个高频子带；步骤2：Contourlet系数融合；按照定义的融合规则在所有尺度和方向上对两幅图像的变换系数进行融合，得到融合后的系数；步骤3：融合图像为R；对于融合后的系数，按照其低频子带和高频子带的顺序，依次进行Contourlet逆变换；步骤4：输出融合后的图像。

【技术特征摘要】
1.一种基于杜鹃搜索算法的Contourlet变换图像融合方法，用于对配准后的源图像A、B进行融合，其中源图像A、B均为灰度图像，且A,B∈RM×N，RM×N是大小为M×N的空间；其特征在于，所述方法包括以下步骤：步骤1：输入两幅源图像A、B，并分别进行Contourlet变换，得到1个低频子带和不同尺度、不同方向的多个高频子带；步骤2：Contourlet系数融合；按照定义的融合规则在所有尺度和方向上对两幅图像的变换系数进行融合，得到融合后的系数；步骤3：融合图像为R；对于融合后的系数，按照其低频子带和高频子带的顺序，依次进行Contourlet逆变换；步骤4：输出融合后的图像。2.根据权利要求1所述的基于杜鹃搜索算法的Contourlet变换图像融合方法，其特征在于：步骤1中，将图像A、B进行Contourlet变换，首先假设对输入图像A、B进行L层Contourlet分解，第l层上的方向数为k，其中0≤l≤L；则Contourlet变换过程为：A→(a0,a1,a2,...aL-1,alow)；B→(b0,b1,b2,...bL-1,blow)；al＝{Al,0,Al,1,Al,2,...Al,k-1},(0≤l≤L-1)；bl＝{Bl,0,Bl,1,Bl,2,...Bl,k-1},(0≤l≤L-1)；其中，alow、blow分别为A、B的低频子带，al和bl分别为图像A、B第l层上的高频方向子带集合，Al,k-1为图像A第l层上的高频子带的k-1方向的分量，Bl,k-1为图像B第l层上的高频子带的k-1方向的分量。3.根据权利要求1所述的基于杜鹃搜索算法的Contourlet变换图像融合方法，其特征在于，步骤2的具体实现包括以下子步骤：步骤2.1：根据步骤1进行Contourlet变换，分解层数为L，每层方向数为k，得到1个低频子带和不同尺度、不同方向的L×k个高频子带，则需要计算的各尺度各方向上的权重总数为L×k+1个；初始化杜鹃搜索算法所需的参数，得到算法的初始种群；其中，参数包括迭代次数Iter、算法的种群规模N、鸟蛋被发现的概率pa、速度步长控制参数η、列维飞行步长控制参数λ、鸟巢的初始空间位置；鸟巢的初始空间位置值向量是相应的每一维对应权重集合中的每一个权重；步骤2.2：利用融合后图像的信息熵计算每个鸟巢的适应度函数值，记录并保留当前最优鸟巢位置A；适应度函数计算公式：其中，P(l)为某灰度值l在图像中出现的概率，L为图像的灰度等级；步骤2.3：通过杜鹃搜索算法中的位置更新公式更新鸟巢的位置，得到一组新解并计算每个鸟巢的适应度函数值，比较这组解的适应度函数值并记其具有最优适应度函数值的位置为B；步骤2.4：用随机数r∈[0,1]与鸟巢的主人发现外来的鸟蛋概率pa对比，并判断：若r＞pa，则通过...

【专利技术属性】
技术研发人员：叶志伟，郑逍，王春枝，苏军，张旭，杨娟，孙爽，陈凤，孙一恒，金灿，
申请(专利权)人：湖北工业大学，
类型：发明
国别省市：湖北,42