发明名称 一种基于杜鹃搜索算法的自适应图像增强方法
摘要 本发明公开了一种基于杜鹃搜索算法的自适应图像增强方法,本发明利用定义的图像质量评价函数,综合考虑增强后的图像的空间统计特征,熵等信息,通过利用杜鹃搜索算法对归一化的非完全Beta函数图像增强应用中最优参数问题优化求解,从而可以快速的获得最优增强参数,可用于数字图像处理相关技术领域中;本发明能够快速的获得归一化的非完全Beta函数图像增强的最优参数,所定义的图像质量评价函数能够客观的评价图像的质量,能够用于偏亮或者偏暗图像增强以后的质量评价,能够满足图像自动增强的要求。
申请公布号 CN104376543A 申请公布日期 2015.02.25
申请号 CN201410711875.2 申请日期 2014.11.28
申请人 湖北工业大学 发明人 叶志伟;王明威;赵伟;尹宇洁;王春枝;刘伟;陈宏伟;徐慧;宗欣露
分类号 G06T5/40(2006.01)I 主分类号 G06T5/40(2006.01)I
代理机构 武汉科皓知识产权代理事务所(特殊普通合伙) 42222 代理人 张火春
主权项 一种基于杜鹃搜索算法的自适应图像增强方法,其特征在于:步骤1:输入待增强原始图像,用f(i,j)表示原始图像在像素点(i,j)处的灰度值,读取原始图像各个像素点的灰度值,统计各个灰度级k,k值范围为0‑255;k出现的次数为G(k),得到原始图像的灰度直方图;步骤2:利用已经获得图像的灰度直方图,对于k从0到255,扫描G(k),得到原始图像的最大灰度值G<sub>max</sub>和最小灰度值G<sub>min</sub>,并将原始图像进行归一化处理记为I’,将图像的灰度值变换到[0,1]区间;步骤3:采用自定义的图像质量评价函数作为原图像增强后图像质量的适应度评价函数,将归一化后的图像灰度值代入归一化的非完全Beta函数中,并利用图像质量评价函数求得归一化的非完全Beta函数最优α,β参数;自定义的图像质量评价函数为:Fitness=log(E(I<sub>1e</sub>)×N<sub>T</sub>/Δh)×(sum(h<sub>T</sub>)/(M×N))<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>E</mi><mrow><mo>(</mo><msub><mi>I</mi><mi>le</mi></msub><mo>)</mo></mrow><mo>=</mo><mo>-</mo><munderover><mi>&Sigma;</mi><mn>0</mn><mn>255</mn></munderover><msub><mi>e</mi><mi>i</mi></msub></mrow>]]></math><img file="FDA0000620316760000011.GIF" wi="314" he="141" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>e</mi><mi>i</mi></msub><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>h</mi><mi>i</mi></msub><msub><mi>log</mi><mn>2</mn></msub><mrow><mo>(</mo><msub><mi>h</mi><mi>i</mi></msub><mo>)</mo></mrow></mtd><mtd><mi>while</mi><msub><mi>h</mi><mi>i</mi></msub><mo>&NotEqual;</mo><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mi>other</mi></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000620316760000012.GIF" wi="636" he="158" /></maths>其中,M为原始图像的宽度,N为原始图像高度,单位是像素,M×N为原始图像的大小;I<sub>1e</sub>是对原始图像的灰度值进行归一化处理后,利用归一化的非完全Beta函数对该灰度值进行变换,增强变换后的图像,E(I<sub>1e</sub>)为图像I<sub>1e</sub>的熵值,e<sub>i</sub>是灰度级i的熵;h<sub>i</sub>表示第i级灰度级概率密度,N<sub>T</sub>为灰度级概率密度大于阈值T的像素个数;sum(h<sub>T</sub>)表示灰度级概率密度大于阈值T的所有灰度级概率密度之和;Δh表示灰度级概率密度h<sub>i</sub>的均方差;归一化的非完全Beta函数定义为下式:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>F</mi><mrow><mo>(</mo><mi>u</mi><mo>,</mo><mi>&alpha;</mi><mo>,</mo><mi>&beta;</mi><mo>)</mo></mrow><mo>=</mo><msup><mi>B</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><mi>&alpha;</mi><mo>,</mo><mi>&beta;</mi><mo>)</mo></mrow><msubsup><mo>&Integral;</mo><mn>0</mn><mi>u</mi></msubsup><msup><mi>t</mi><mrow><mi>&alpha;</mi><mo>-</mo><mn>1</mn></mrow></msup><msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>t</mi><mo>)</mo></mrow><mrow><mi>&beta;</mi><mo>-</mo><mn>1</mn></mrow></msup><mi>dt</mi></mrow>]]></math><img file="FDA0000620316760000013.GIF" wi="785" he="101" /></maths>其中,u在这里表示归一化后图像的灰度值变量,是步骤2中的f'(i,j),t是积分变量,B<sup>‑1</sup>(α,β)是B(α,β)的逆变换,β的函数B(α,β)由下式给出:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>B</mi><mrow><mo>(</mo><mi>&alpha;</mi><mo>,</mo><mi>&beta;</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mo>&Integral;</mo><mn>0</mn><mn>1</mn></msubsup><msup><mi>t</mi><mrow><mi>&alpha;</mi><mo>-</mo><mn>1</mn></mrow></msup><msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>t</mi><mo>)</mo></mrow><mrow><mi>&beta;</mi><mo>-</mo><mn>1</mn></mrow></msup><mi>dt</mi></mrow>]]></math><img file="FDA0000620316760000021.GIF" wi="597" he="133" /></maths>上式中t是积分变量,α,β是归一化的非完全Beta函数参数,步骤4,利用步骤3中求得最优α,β参数,采用归一化的非完全Beta变换函数F(u,α,β)(0≤u≤1),对归一化图像进行变换,所述的变换公式如下:g'(i,j)=F(f'(i,j),α,β)上式中,f'(i,j)表示归一化图像在像素点(i,j)处的灰度值,就是归一化的非完全Beta函数的输入变量u,g'(i,j)表示变换后图像在像素点(i,j)处的灰度值,其中f'(i,j)和g'(i,j)均满足0≤g'(i,j)≤1,0≤f'(i,j)≤1;步骤5,对归一化增强变换后的图像进行反变换处理,其采用的公式为f"(i,j)=(G'<sub>max</sub>‑G'<sub>min</sub>)g'(i,j)+G'<sub>min</sub>式中,G'<sub>max</sub>和G'<sub>min</sub>分别为变换后图像的最大和最小灰度值,对于8位灰度图像,G'<sub>max</sub>=255,G'<sub>min</sub>=0,g'(i,j)是指采用归一化的非完全Beta函数增强后像素点(i,j)的灰度值,f″(i,j)为进行反归一化后像素点(i,j)的灰度值;步骤6,输出增强后图像。
地址 430068 湖北省武汉市武昌区南湖李家墩1村1号