主权项 |
一种基于杜鹃搜索算法的自适应图像增强方法,其特征在于:步骤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>Σ</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>≠</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>α</mi><mo>,</mo><mi>β</mi><mo>)</mo></mrow><mo>=</mo><msup><mi>B</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><mi>α</mi><mo>,</mo><mi>β</mi><mo>)</mo></mrow><msubsup><mo>∫</mo><mn>0</mn><mi>u</mi></msubsup><msup><mi>t</mi><mrow><mi>α</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>β</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>α</mi><mo>,</mo><mi>β</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mo>∫</mo><mn>0</mn><mn>1</mn></msubsup><msup><mi>t</mi><mrow><mi>α</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>β</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,输出增强后图像。 |