发明名称 基于生物视觉机理的图像去雾方法
摘要 本发明公开了一种基于生物视觉机理的图像去雾方法,本发明的方法通过非经典感受野单拮抗模型,利用不同敏感系数下外周的抑制作用和亚区的去抑制作用的相互作用,将有雾图像的弥漫光照去除,增强原始图像色彩对比度,实现有雾图像的自动增强,本发明提供了一种符合生理特性的方法,在保持细节特征的同时,准确、高效地实现了单幅有雾彩色图像的颜色去雾增强。
申请公布号 CN103955900A 申请公布日期 2014.07.30
申请号 CN201410191124.2 申请日期 2014.05.07
申请人 电子科技大学 发明人 张显石;李永杰;李朝义
分类号 G06T5/00(2006.01)I 主分类号 G06T5/00(2006.01)I
代理机构 成都宏顺专利代理事务所(普通合伙) 51227 代理人 周永宏
主权项 一种基于生物视觉机理的图像去雾方法,包括以下步骤:步骤1:设定感受野大小及确定相应的模型参数;步骤2:对色偏图像每一个像素点分别提取红色分量I<sub>R</sub>、绿色分量I<sub>G</sub>、蓝色分量I<sub>B</sub>,计算出平均亮度分量I:I=(I<sub>R</sub>+I<sub>G</sub>+I<sub>B</sub>)/3;步骤3:设定去雾深度th和兴奋抑制比K;步骤4:将红、绿、蓝三通道像素值均值最小者与去雾深度比较,若该值小于去雾深度,即min(mean(I<sub>R</sub>),mean(I<sub>G</sub>),mean(I<sub>B</sub>)))&lt;th,图像无需进行去雾,处理完成;大于等于去雾深度则执行步骤5;步骤5:根据兴奋抑制比计算出红通道非经典感受野外周敏感系数A<sub>R2</sub>,亚区敏感系数A<sub>R3</sub>,绿通道非经典感受野外周敏感系数A<sub>G2</sub>,亚区敏感系数A<sub>G3</sub>,蓝通道非经典感受野外周敏感系数A<sub>B2</sub>,亚区敏感系数A<sub>B3</sub>,A<sub>R2</sub>=K3   A<sub>R3</sub>=A<sub>R2</sub>3A<sub>G2</sub>=K3   A<sub>G3</sub>=A<sub>G2</sub>3A<sub>B2</sub>=K3   A<sub>B3</sub>=A<sub>B2</sub>3步骤7:根据公式<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>R</mi><mrow><mi>R</mi><mn>3</mn></mrow></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>;</mo><msub><mi>&sigma;</mi><mi>u</mi></msub><mo>)</mo></mrow><mo>=</mo><mi>MAX</mi><mo>[</mo><mn>0</mn><mo>,</mo><mi>I</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>A</mi><mrow><mi>R</mi><mn>3</mn></mrow></msub><mo>&times;</mo><mi>I</mi><mrow><mo>(</mo><mi>p</mi><mo>,</mo><mi>q</mi><mo>)</mo></mrow><mo>&CircleTimes;</mo><mi>g</mi><mrow><mo>(</mo><mo>|</mo><mi>p</mi><mo>-</mo><mi>x</mi><mo>|</mo><mo>,</mo><mo>|</mo><mi>q</mi><mo>-</mo><mi>y</mi><mo>|</mo><mo>;</mo><msub><mi>&sigma;</mi><mi>u</mi></msub><mo>)</mo></mrow><mo>]</mo></mrow>]]></math><img file="FDA0000501911380000011.GIF" wi="1428" he="84" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>R</mi><mrow><mi>G</mi><mn>3</mn></mrow></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>;</mo><msub><mi>&sigma;</mi><mi>u</mi></msub><mo>)</mo></mrow><mo>=</mo><mi>MAX</mi><mo>[</mo><mn>0</mn><mo>,</mo><mi>I</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>A</mi><mrow><mi>G</mi><mn>3</mn></mrow></msub><mo>&times;</mo><mi>I</mi><mrow><mo>(</mo><mi>p</mi><mo>,</mo><mi>q</mi><mo>)</mo></mrow><mo>&CircleTimes;</mo><mi>g</mi><mrow><mo>(</mo><mo>|</mo><mi>p</mi><mo>-</mo><mi>x</mi><mo>|</mo><mo>,</mo><mo>|</mo><mi>q</mi><mo>-</mo><mi>y</mi><mo>|</mo><mo>;</mo><msub><mi>&sigma;</mi><mi>u</mi></msub><mo>)</mo></mrow><mo>]</mo></mrow>]]></math><img file="FDA0000501911380000012.GIF" wi="1432" he="80" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>R</mi><mrow><mi>B</mi><mn>3</mn></mrow></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>;</mo><msub><mi>&sigma;</mi><mi>u</mi></msub><mo>)</mo></mrow><mo>=</mo><mi>MAX</mi><mo>[</mo><mn>0</mn><mo>,</mo><mi>I</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>A</mi><mrow><mi>B</mi><mn>3</mn></mrow></msub><mo>&times;</mo><mi>I</mi><mrow><mo>(</mo><mi>p</mi><mo>,</mo><mi>q</mi><mo>)</mo></mrow><mo>&CircleTimes;</mo><mi>g</mi><mrow><mo>(</mo><mo>|</mo><mi>p</mi><mo>-</mo><mi>x</mi><mo>|</mo><mo>,</mo><mo>|</mo><mi>q</mi><mo>-</mo><mi>y</mi><mo>|</mo><mo>;</mo><msub><mi>&sigma;</mi><mi>u</mi></msub><mo>)</mo></mrow><mo>]</mo></mrow>]]></math><img file="FDA0000501911380000013.GIF" wi="1428" he="90" /></maths>计算出红通道亚区去抑制后响应R<sub>R3</sub>(x,y;σ<sub>u</sub>),绿通道亚区去抑制后响应R<sub>G3</sub>(x,y;σ<sub>u</sub>),蓝通道亚区去抑制后响应R<sub>B3</sub>(x,y;σ<sub>u</sub>),其中,<img file="FDA0000501911380000017.GIF" wi="44" he="44" />表示卷积,I(x,y)表示中心点亮度,I(p,q)表示亚区内除中心外的点的亮度,(p,q)为落在亚区内除中心外的点,MAX表示取两者中较大值;步骤8:根据公式<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>R</mi><mrow><mi>R</mi><mn>2</mn></mrow></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>;</mo><msub><mi>&sigma;</mi><mi>s</mi></msub><mo>)</mo></mrow><mo>=</mo><msub><mi>A</mi><mrow><mi>R</mi><mn>2</mn></mrow></msub><mo>&times;</mo><msub><mi>R</mi><mrow><mi>R</mi><mn>3</mn></mrow></msub><mrow><mo>(</mo><mi>p</mi><mo>,</mo><mi>q</mi><mo>;</mo><msub><mi>&sigma;</mi><mi>u</mi></msub><mo>)</mo></mrow><mo>&CircleTimes;</mo><mi>g</mi><mrow><mo>(</mo><mo>|</mo><mi>p</mi><mo>-</mo><mi>x</mi><mo>|</mo><mo>,</mo><mo>|</mo><mi>q</mi><mo>-</mo><mi>y</mi><mo>|</mo><mo>;</mo><msub><mi>&sigma;</mi><mi>s</mi></msub><mo>)</mo></mrow><mo>]</mo></mrow>]]></math><img file="FDA0000501911380000014.GIF" wi="1060" he="71" /></maths><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>R</mi><mrow><mi>G</mi><mn>2</mn></mrow></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>;</mo><msub><mi>&sigma;</mi><mi>s</mi></msub><mo>)</mo></mrow><mo>=</mo><msub><mi>A</mi><mrow><mi>G</mi><mn>2</mn></mrow></msub><mo>&times;</mo><msub><mi>R</mi><mrow><mi>G</mi><mn>3</mn></mrow></msub><mrow><mo>(</mo><mi>p</mi><mo>,</mo><mi>q</mi><mo>;</mo><msub><mi>&sigma;</mi><mi>u</mi></msub><mo>)</mo></mrow><mo>&CircleTimes;</mo><mi>g</mi><mrow><mo>(</mo><mo>|</mo><mi>p</mi><mo>-</mo><mi>x</mi><mo>|</mo><mo>,</mo><mo>|</mo><mi>q</mi><mo>-</mo><mi>y</mi><mo>|</mo><mo>;</mo><msub><mi>&sigma;</mi><mi>s</mi></msub><mo>)</mo></mrow><mo>]</mo></mrow>]]></math><img file="FDA0000501911380000015.GIF" wi="1066" he="72" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>R</mi><mrow><mi>B</mi><mn>2</mn></mrow></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>;</mo><msub><mi>&sigma;</mi><mi>s</mi></msub><mo>)</mo></mrow><mo>=</mo><msub><mi>A</mi><mrow><mi>B</mi><mn>2</mn></mrow></msub><mo>&times;</mo><msub><mi>R</mi><mrow><mi>B</mi><mn>3</mn></mrow></msub><mrow><mo>(</mo><mi>p</mi><mo>,</mo><mi>q</mi><mo>;</mo><msub><mi>&sigma;</mi><mi>u</mi></msub><mo>)</mo></mrow><mo>&CircleTimes;</mo><mi>g</mi><mrow><mo>(</mo><mo>|</mo><mi>p</mi><mo>-</mo><mi>x</mi><mo>|</mo><mo>,</mo><mo>|</mo><mi>q</mi><mo>-</mo><mi>y</mi><mo>|</mo><mo>;</mo><msub><mi>&sigma;</mi><mi>s</mi></msub><mo>)</mo></mrow><mo>]</mo></mrow>]]></math><img file="FDA0000501911380000016.GIF" wi="1060" he="73" /></maths>计算出红通道外周抑制R<sub>R2</sub>(x,y;σ<sub>s</sub>),绿通道外周抑制R<sub>G2</sub>(x,y;σ<sub>s</sub>),蓝通道外周抑制R<sub>B2</sub>(x,y;σ<sub>s</sub>),其中,(p,q)为落在外周内的点;步骤9:根据公式<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>R</mi><mrow><mi>R</mi><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>;</mo><msub><mi>&sigma;</mi><mi>c</mi></msub><mo>)</mo></mrow><mo>=</mo><mi>MAX</mi><mo>[</mo><mn>0</mn><mo>,</mo><msub><mi>I</mi><mi>R</mi></msub><mrow><mo>(</mo><mi>p</mi><mo>,</mo><mi>q</mi><mo>)</mo></mrow><mo>&CircleTimes;</mo><mi>g</mi><mrow><mo>(</mo><mo>|</mo><mi>p</mi><mo>-</mo><mi>x</mi><mo>|</mo><mo>,</mo><mo>|</mo><mi>q</mi><mo>-</mo><mi>y</mi><mo>|</mo><mo>;</mo><msub><mi>&sigma;</mi><mi>c</mi></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>R</mi><mrow><mi>R</mi><mn>2</mn></mrow></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>;</mo><msub><mi>&sigma;</mi><mi>s</mi></msub><mo>)</mo></mrow><mo>]</mo></mrow>]]></math><img file="FDA0000501911380000021.GIF" wi="1419" he="90" /></maths><maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msub><mi>R</mi><mrow><mi>G</mi><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>;</mo><msub><mi>&sigma;</mi><mi>c</mi></msub><mo>)</mo></mrow><mo>=</mo><mi>MAX</mi><mo>[</mo><mn>0</mn><mo>,</mo><msub><mi>I</mi><mi>G</mi></msub><mrow><mo>(</mo><mi>p</mi><mo>,</mo><mi>q</mi><mo>)</mo></mrow><mo>&CircleTimes;</mo><mi>g</mi><mrow><mo>(</mo><mo>|</mo><mi>p</mi><mo>-</mo><mi>x</mi><mo>|</mo><mo>,</mo><mo>|</mo><mi>q</mi><mo>-</mo><mi>y</mi><mo>|</mo><mo>;</mo><msub><mi>&sigma;</mi><mi>c</mi></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>R</mi><mrow><mi>G</mi><mn>2</mn></mrow></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>;</mo><msub><mi>&sigma;</mi><mi>s</mi></msub><mo>)</mo></mrow><mo>]</mo></mrow>]]></math><img file="FDA0000501911380000022.GIF" wi="1425" he="82" /></maths><maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><msub><mi>R</mi><mrow><mi>B</mi><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>;</mo><msub><mi>&sigma;</mi><mi>c</mi></msub><mo>)</mo></mrow><mo>=</mo><mi>MAX</mi><mo>[</mo><mn>0</mn><mo>,</mo><msub><mi>I</mi><mi>B</mi></msub><mrow><mo>(</mo><mi>p</mi><mo>,</mo><mi>q</mi><mo>)</mo></mrow><mo>&CircleTimes;</mo><mi>g</mi><mrow><mo>(</mo><mo>|</mo><mi>p</mi><mo>-</mo><mi>x</mi><mo>|</mo><mo>,</mo><mo>|</mo><mi>q</mi><mo>-</mo><mi>y</mi><mo>|</mo><mo>;</mo><msub><mi>&sigma;</mi><mi>c</mi></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>R</mi><mrow><mi>B</mi><mn>2</mn></mrow></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>;</mo><msub><mi>&sigma;</mi><mi>s</mi></msub><mo>)</mo></mrow><mo>]</mo></mrow>]]></math><img file="FDA0000501911380000023.GIF" wi="1419" he="90" /></maths>计算出红通道中心区抑制后响应R<sub>R1</sub>(x,y;σ<sub>c</sub>),绿通道中心区抑制后响应R<sub>G1</sub>(x,y;σ<sub>c</sub>),蓝通道中心区抑制后响应R<sub>B1</sub>(x,y;σ<sub>c</sub>),其中,(p,q)为落在中心区Center内的点,取R<sub>R1</sub>(x,y;σ<sub>c</sub>)、R<sub>G1</sub>(x,y;σ<sub>c</sub>)、R<sub>B1</sub>(x,y;σ<sub>c</sub>)作为像素点(x,y)新的红、绿、蓝分量;步骤10:在红、绿、蓝通道上分别计算,以红通道输出为红色分量I<sub>R</sub>,绿通道输出为绿色分量I<sub>G</sub>,蓝通道输出为蓝色分量I<sub>B</sub>合成彩色图像,若满足步骤4的条件,则去雾停止,否则,兴奋抑制比K加1,返回步骤5。
地址 611731 四川省成都市高新区(西区)西源大道2006号