发明名称 一种单幅图像去雾方法
摘要 本发明一种单幅图像去雾方法,利用待去雾RGB图像I的区域颜色均值向量及其L2范数来获取传播图,再根据传播图和大气光值来计算去雾后的图像;该方法可以直接利用单幅RGB彩色图像,无需外界提供任何其他信息;同时本发明计算去雾图像过程中利用系数可以直接得到去雾图像,无需在去雾后进一步增加图像亮度以满足显示要求。
申请公布号 CN104008527B 申请公布日期 2017.02.01
申请号 CN201410153443.4 申请日期 2014.04.16
申请人 南京航空航天大学 发明人 丁萌;魏丽;王帮峰;刘中杰
分类号 G06T5/00(2006.01)I 主分类号 G06T5/00(2006.01)I
代理机构 南京经纬专利商标代理有限公司 32200 代理人 许方
主权项 一种单幅图像去雾方法,该方法通过有雾RGB彩色图像的大气光值及反映图像深度信息的传播图来计算去雾后的图像,其特征在于,根据下式计算去雾后图像J:<maths num="0001"><math><![CDATA[<mrow><mi>J</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>I</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><mi>A</mi></mrow><mrow><mi>m</mi><mi>a</mi><mi>x</mi><mrow><mo>(</mo><mi>T</mi><mo>(</mo><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mo>)</mo><mo>,</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow></mrow></mfrac><mo>+</mo><mi>&alpha;</mi><mi>A</mi></mrow>]]></math><img file="FDA0001112871200000011.GIF" wi="692" he="149" /></maths>式中,t<sub>0</sub>=0.3,α∈[0.6,0.9],T(i,j)为传播图,I(i,j)为待去雾原图;max(T(i,j),t<sub>0</sub>)表示当T(i,j)≥t<sub>0</sub>时,max(T(i,j),t<sub>0</sub>)=T(i,j);否则max(T(i,j),t<sub>0</sub>)=t<sub>0</sub>;该方法中的传播图通过以下方法获取:步骤A、获取待去雾RGB图像I的区域颜色均值向量及其L2范数;步骤A‑1、输入大小为m×n的待去雾RGB图像I,图像I中每个像素点对应一个1×3颜色向量[I<sup>R</sup>(i,j),I<sup>G</sup>(i,j),I<sup>B</sup>(i,j)],其中,i,j表示像素点坐标,i∈[1,m],j∈[1,n],其中,m、n均为正整数;步骤A‑2:获取图像I中满足以下条件的所有区域窗口Ω(i,j):以第i行第j列为中心像素点,构建大小为(2r+1)×(2r+1)的区域窗口;其中,i的取值范围是r+1到m‑r之间的整数;j的取值范围是r+1到n‑r之间的整数,r区域窗口半径;步骤A‑3:利用下式获取每个区域窗口Ω(i,j)对应的图像I区域颜色均值向量E(i,j)E(i,j)=[e<sup>R</sup>(i,j),e<sup>G</sup>(i,j),e<sup>B</sup>(i,j)]<sup>T</sup>,其中<img file="FDA0001112871200000012.GIF" wi="702" he="226" />C∈{R,G,B};步骤A‑4:根据每个E(i,j)得到对应的图像I区域颜色均值向量的L2范数u(i,j),其表达形式:u(i,j)=||E(i,j)||<sub>2</sub>;步骤B、获取反映图像深度信息的传播图T;步骤B‑1、利用下式计算传播图T的初始值<img file="FDA0001112871200000013.GIF" wi="62" he="62" /><maths num="0002"><math><![CDATA[<mrow><mover><mi>T</mi><mo>~</mo></mover><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mn>1</mn><mo>-</mo><mfrac><mrow><msup><mi>A</mi><mi>T</mi></msup><mi>E</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><msup><mi>u</mi><mn>2</mn></msup><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mrow><mrow><mo>|</mo><mo>|</mo><mi>A</mi><mo>|</mo><msubsup><mo>|</mo><mn>2</mn><mn>2</mn></msubsup><mo>-</mo><msup><mi>A</mi><mi>T</mi></msup><mi>E</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0001112871200000014.GIF" wi="581" he="150" /></maths>式中,A为大气光值,A<sup>T</sup>表示对A进行转置运算;步骤B‑2、利用引导滤波算法进一步优化传播图初始值<img file="FDA0001112871200000021.GIF" wi="67" he="61" />得到反映图像深度信息的传播图T。
地址 210016 江苏省南京市秦淮区御道街29号