发明名称 多聚焦图像融合方法
摘要 本发明公开了一种多聚焦图像融合方法,属于图像处理领域。在多聚焦图像融合过程中采用增强差分演化方法和扩块选择机制图像融合方法,其中相比于标准差分演化方法,该增强差分演化方法根据图像的大小,区间初始化种群,在每一次迭代后进行种群判断,在种群中各个体相同的前提下保留其中一个个体,再利用标准差分演化算法的初始化方法,在相应的区间上重新初始化其余个体,一起作为新一代种群,加强局部搜索,由此获得图像的最佳分块大小;另外,在图像块的适应度相等时,选择扩展块,比较获得融合的图像块。通过本发明,提高了图像融合中寻找最佳分块大小的效率,并且降低了计算的复杂度;另外,在图像融合中充分考虑了整体融合图像的清晰度。
申请公布号 CN102509280B 申请公布日期 2013.09.18
申请号 CN201110355867.5 申请日期 2011.11.10
申请人 重庆大学 发明人 冯永;李铁柱;钟将;周尚波;李季
分类号 G06T5/50(2006.01)I;G06N3/00(2006.01)I 主分类号 G06T5/50(2006.01)I
代理机构 重庆市前沿专利事务所(普通合伙) 50211 代理人 郭云
主权项 1.一种多聚焦图像融合方法,其特征在于包括以下步骤:S1、获取P张多聚焦图像,且每张多聚焦图像的图像大小为M像素×N像素,采用增强差分演化方法求取所述多聚焦图像的最佳分块大小:设定初始化种群由K个个体组成,个体为二维变量且表示为X<sub>i,G</sub>=(X<sub>1,i,G</sub>,X<sub>2,i,G</sub>),其中X<sub>i,G</sub>表示第G代种群中第i个个体,X<sub>1,i,G</sub>表示第G代种群中第i个个体的第一维变量,X<sub>2,i,G</sub>表示第G代种群中第i个个体的第二维变量,1≤i≤K,i、K均为整数,并且设定种群的变异因子F、交叉率CR以及迭代次数t的最大值N,其中所述变异因子F的取值范围为(0,2),所述交叉率CR的取值范围为[0,1],初始化迭代次数为0;S2、根据所述多聚焦图像将图像块大小的取值范围划分成K个区间,第j个区间表示为<maths num="0001"><![CDATA[<math><mrow><msub><mi>Z</mi><mi>j</mi></msub><mo>=</mo><mo>{</mo><mrow><mo>(</mo><mrow><mo>(</mo><mi>j</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>*</mo><mfrac><mi>M</mi><mi>K</mi></mfrac><mo>+</mo><mn>1</mn><mo>~</mo><mi>j</mi><mo>*</mo><mfrac><mi>M</mi><mi>K</mi></mfrac><mo>)</mo></mrow><mo>,</mo><mrow><mo>(</mo><mrow><mo>(</mo><mi>j</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>*</mo><mfrac><mi>N</mi><mi>K</mi></mfrac><mo>+</mo><mn>1</mn><mo>~</mo><mi>j</mi><mo>*</mo><mfrac><mi>N</mi><mi>K</mi></mfrac><mo>)</mo></mrow><mo>}</mo><mo>,</mo></mrow></math>]]></maths>其中j的取值范围为[1,K]且j为整数,第j个区间Z<sub>j</sub>中第一维变量的上界<img file="FDA00003362400000012.GIF" wi="384" he="127" />下界<img file="FDA00003362400000013.GIF" wi="527" he="126" />且第j个区间Z<sub>j</sub>中第二维变量的上界表示为<img file="FDA00003362400000014.GIF" wi="353" he="131" />下界表示为<maths num="0002"><![CDATA[<math><mrow><msub><mi>X</mi><mrow><mn>2</mn><mo>,</mo><mi>j</mi><mrow><mo>(</mo><mi>L</mi><mo>)</mo></mrow></mrow></msub><mo>=</mo><mrow><mo>(</mo><mi>j</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>*</mo><mfrac><mi>N</mi><mi>K</mi></mfrac><mo>+</mo><mn>1</mn><mo>;</mo></mrow></math>]]></maths>S3、根据公式X’<sub>j,G</sub>=(X<sub>1,j(L)</sub>,X<sub>2,j(L)</sub>)+rand[0,1]×{(X<sub>1,j(H)</sub>,X<sub>2,j(H)</sub>)-(X<sub>1,j(L)</sub>,X<sub>2,j(L)</sub>)},从而产生K个新的个体X’<sub>j,G</sub>,构成初始化种群X’<sub>G</sub>,其中X’<sub>j,G</sub>表示第j区间对应产生的新的个体,rand[0,1]表示取[0,1]范围内服从均匀分布的一个随机数,且对于不同的区间rand[0,1]的取值不同;S4、从所述初始化种群X’<sub>G</sub>中随机选取3个个体X<sub>r1,G</sub>、X<sub>r2,G</sub>、X<sub>r3,G</sub>,其中下标r1、r2、r3均为[1,K]范围内不相等的整数,根据公式V<sub>i,G</sub>=X<sub>r1,G</sub>+F×(X<sub>r2,G</sub>-X<sub>r3,G</sub>),获得变异后的个体V<sub>i,G</sub>,其中F为变异因子,重复执行K次此步骤,产生K个变异后的个体,重组获得变异种群V<sub>G</sub>;S5、对所述初始化种群X’<sub>G</sub>与所述变异种群V<sub>G</sub>执行交叉操作,生成交叉种群U<sub>G</sub>,使得所述交叉种群U<sub>G</sub>中每一个体至少有一个变量来自所述变异种群V<sub>G</sub>;S6、分别计算所述初始化种群X’<sub>G</sub>与交叉种群U<sub>G</sub>中每一个体的适应度f(X′<sub>i,G</sub>)、f(U<sub>i,G</sub>),其中f(X′<sub>i,G</sub>)表示所述初始化种群X’<sub>G</sub>中第i个个体的适应度,f(U<sub>i,G</sub>)表示所述交叉种群U<sub>G</sub>中第i个个体的适应度,并且根据公式<maths num="0003"><![CDATA[<math><mrow><msub><mi>X</mi><mrow><mi>i</mi><mo>,</mo><mi>G</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>U</mi><mrow><mi>i</mi><mo>,</mo><mi>G</mi></mrow></msub></mtd><mtd><mi>if</mi></mtd><mtd><mi>f</mi><mrow><mo>(</mo><msubsup><mi>X</mi><mrow><mi>i</mi><mo>,</mo><mi>G</mi></mrow><mo>&prime;</mo></msubsup><mo>)</mo></mrow><mo>&lt;</mo><mi>f</mi><mrow><mo>(</mo><msub><mi>U</mi><mrow><mi>i</mi><mo>,</mo><mi>G</mi></mrow></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><msup><mi>X</mi><mo>&prime;</mo></msup><mrow><mi>i</mi><mo>,</mo><mi>G</mi></mrow></msub></mtd><mtd></mtd><mtd><mi>otherwise</mi></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths>选择所述初始化种群X’<sub>G</sub>与交叉种群U<sub>G</sub>中对应的个体作为下一代种群的个体,其中X<sub>i,G+1</sub>表示第G+1代种群中的第i个个体:如果所述初始化种群X’<sub>G</sub>中第i个个体的适应度小于所述交叉种群U<sub>G</sub>中第i个个体的适应度,则选择所述交叉种群U<sub>G</sub>中第i个个体U<sub>i,G</sub>作为第G+1代种群的第i个个体X<sub>i,G+1</sub>,否则选择所述初始化种群X’<sub>G</sub>中第i个个体X'<sub>i,G</sub>作为第G+1代种群的第i个个体X<sub>i,G+1</sub>,从而形成第G+1代种群X<sub>G+1</sub>,并且所述迭代次数加1;S7、判断所述第G+1代种群X<sub>G+1</sub>中所有个体是否均相同:如果所有个体均相同则执行步骤S8,否则进一步判断所述迭代次数t是否大于迭代次数的最大值N:如果t≥N则表示获得最佳分块大小,将所述最佳分块大小输出,如果t&lt;N则重复执行步骤S4~S7;S8、保留所述第G+1代种群X<sub>G+1</sub>中的一个个体i,并且在i所在的区间上采用标准差分演化方法的初始化种群步骤重新初始化K-1个个体,与i一起形成新一代种群,重复执行步骤S4-S7;S9、所述多聚焦图像根据所述步骤S7中的所述最佳分块大小分割成多个图像块,并且对所述图像块进行图像融合,形成多聚焦融合图像。
地址 400045 重庆市沙坪坝区沙正街174号