发明名称 基于局部SoG的水下图像清晰化方法
摘要 本发明公开了一种基于局部SoG的水下图像清晰化方法,包括以下步骤:设水下图像为原图像,用该原图像的边界值进行边界扩充得到填充后的图像;然后,取一个局部域大小为N×N像素点的窗口,利用该窗口在边界扩充后的图像的各个通道中进行平移,估计每个局部域内的光源颜色,将局部域的光源颜色作为窗口中心像素点的光源颜色e<sub>λ</sub>;并通过光源颜色e<sub>λ</sub>与原始光源颜色<img file="DDA0000937960040000011.GIF" wi="50" he="54" />之间的衰减特性分别求出三个颜色通道的透射率t<sub>λ</sub>和其对应的方差σ<sub>λ</sub>;最后,对方差σ<sub>λ</sub>最大的透射率进行循环分块操作以获取背景光,从而得到最终的清晰的图像。通过大量实验证明,本发明提出的算法在提高图像清晰度、恢复图像颜色等方面有着非常好的性能。
申请公布号 CN105825481A 申请公布日期 2016.08.03
申请号 CN201610137128.1 申请日期 2016.03.09
申请人 天津大学 发明人 杨爱萍;曲畅;郑佳;王建
分类号 G06T5/00(2006.01)I 主分类号 G06T5/00(2006.01)I
代理机构 天津市北洋有限责任专利代理事务所 12201 代理人 李丽萍
主权项 一种基于局部SoG的水下图像清晰化方法,其特征在于,包括以下步骤:步骤一、输入水下图像为原图像I(x),用该原图像I(x)的边界值进行边界扩充得到填充后的图像I<sub>λ</sub>(x);I<sub>λ</sub>(x)=∫<sub>ω</sub>e(λ)S(x,λ)P(λ)dλ        (1)式(1)中,ω表示整个可见光波长范围,λ为入射光波长,x为像素点,e(λ)表示入射光的分布,S(x,λ)为图像I<sub>λ</sub>(x)中的某一点对入射光的反射率,P(λ)表示相机对入射光的感光特性;步骤二、在该图像中取一个大小为N×N像素点的窗口,利用该窗口在图像I<sub>λ</sub>(x)的各个通道上逐行的从左到右进行平移,每次移动一个像素形成一局部域Ω,估计每个局部域Ω的光源颜色,步骤是:设该图像I<sub>λ</sub>(x)的场景中所有物体表面的任意一个局部域Ω的平均反射是无色差的,即:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfrac><mrow><msub><mo>&Integral;</mo><mi>&Omega;</mi></msub><mi>S</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>&lambda;</mi><mo>)</mo></mrow><mi>d</mi><mi>x</mi></mrow><mrow><msub><mo>&Integral;</mo><mi>&Omega;</mi></msub><mrow><mi>d</mi><mi>x</mi></mrow></mrow></mfrac><mo>=</mo><mi>k</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000937960010000011.GIF" wi="1140" he="191" /></maths>式(2)中,Ω表示图像I<sub>λ</sub>(x)中某个局部域,k为[0,1]之间的一个常数;所以,局部域Ω的光源颜色如下:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msup><mrow><mo>(</mo><mfrac><mrow><msub><mo>&Integral;</mo><mi>&Omega;</mi></msub><mrow><msup><mrow><mo>(</mo><msub><mi>I</mi><mi>&lambda;</mi></msub><mo>(</mo><mi>x</mi><mo>)</mo><mo>)</mo></mrow><mi>p</mi></msup><mi>d</mi><mi>x</mi></mrow></mrow><mrow><msub><mo>&Integral;</mo><mi>&Omega;</mi></msub><mrow><mi>d</mi><mi>x</mi></mrow></mrow></mfrac><mo>)</mo></mrow><mfrac><mn>1</mn><mi>p</mi></mfrac></msup><mo>=</mo><msub><mi>ke</mi><mrow><mi>&Omega;</mi><mi>&lambda;</mi></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000937960010000012.GIF" wi="1118" he="230" /></maths>式(3)中,p为常数,取值范围为4~10,e<sub>Ωλ</sub>为局部域Ω的光源颜色,在一个局部域Ω中光源颜色的各通道分量为单值;将局部域Ω的光源颜色e<sub>Ωλ</sub>作为局部域中心像素点的光源颜色e<sub>λ</sub>:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msup><mrow><mo>(</mo><mfrac><mrow><msub><mo>&Integral;</mo><mi>&Omega;</mi></msub><mrow><msup><mrow><mo>(</mo><msub><mi>I</mi><mi>&lambda;</mi></msub><mo>(</mo><mi>x</mi><mo>)</mo><mo>)</mo></mrow><mi>p</mi></msup><mi>d</mi><mi>x</mi></mrow></mrow><mrow><msub><mo>&Integral;</mo><mi>&Omega;</mi></msub><mrow><mi>d</mi><mi>x</mi></mrow></mrow></mfrac><mo>)</mo></mrow><mfrac><mn>1</mn><mi>p</mi></mfrac></msup><mo>=</mo><msub><mi>ke</mi><mi>&lambda;</mi></msub><mo>=</mo><mi>k</mi><msub><mover><mi>e</mi><mo>&OverBar;</mo></mover><mi>&lambda;</mi></msub><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><msub><mi>c</mi><mi>&lambda;</mi></msub><mi>d</mi><mo>(</mo><mi>x</mi><mo>)</mo><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000937960010000013.GIF" wi="1133" he="222" /></maths>式(4)中,c<sub>λ</sub>是入射光在水中的衰减系数,d(x)为景深,e<sub>λ</sub>是局部域中心像素点的光源颜色,所述光源颜色e<sub>λ</sub>的元素值中包含了图像I<sub>λ</sub>(x)各通道的颜色衰减特征,<img file="FDA0000937960010000014.GIF" wi="61" he="54" />是图像I<sub>λ</sub>(x)的光源颜色,设<img file="FDA0000937960010000015.GIF" wi="51" he="55" />为常数;步骤三:利用局部域中心像素点的光源颜色e<sub>λ</sub>和图像I<sub>λ</sub>(x)的光源颜色<img file="FDA0000937960010000016.GIF" wi="48" he="54" />之间的颜色衰减特性求出图像I<sub>λ</sub>(x)的RGB三通道透射率t<sub>λ</sub>:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msup><mrow><mo>(</mo><mfrac><mrow><msub><mo>&Integral;</mo><mi>&Omega;</mi></msub><mrow><msup><mrow><mo>(</mo><msub><mi>I</mi><mi>&lambda;</mi></msub><mo>(</mo><mi>x</mi><mo>)</mo><mo>)</mo></mrow><mi>p</mi></msup><mi>d</mi><mi>x</mi></mrow></mrow><mrow><msub><mo>&Integral;</mo><mi>&Omega;</mi></msub><mrow><mi>d</mi><mi>x</mi></mrow></mrow></mfrac><mo>)</mo></mrow><mfrac><mn>1</mn><mi>p</mi></mfrac></msup><mo>=</mo><mi>k</mi><msub><mover><mi>e</mi><mo>&OverBar;</mo></mover><mi>&lambda;</mi></msub><msub><mi>t</mi><mi>&lambda;</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000937960010000021.GIF" wi="1126" he="222" /></maths>选择R通道的透射率t<sub>R</sub>作为待分块的透射率;步骤四:建立图像I<sub>λ</sub>(x)像素的索引index;与此同时,选择待分块的透射率t<sub>R</sub>上半部分中均值和标准差最小的部分进行循环分块操作,循环4次后得到透射率块,根据图像I<sub>λ</sub>(x)像素索引index,确定该透射率块在图像I<sub>λ</sub>(x)中对应的图像块,通过求取该对应的图像块RGB三通道的均值得到该图像I<sub>λ</sub>(x)的背景光B<sub>λ,∞</sub>;根据式(6)得到清晰图像J<sub>λ</sub>(x):<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>J</mi><mi>&lambda;</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mo>(</mo><msub><mi>I</mi><mi>&lambda;</mi></msub><mo>(</mo><mi>x</mi><mo>)</mo><mo>-</mo><msub><mi>B</mi><mrow><mi>&lambda;</mi><mo>,</mo><mi>&infin;</mi></mrow></msub><mo>)</mo></mrow><mrow><msub><mi>t</mi><mi>&lambda;</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow></mrow></mfrac><mo>+</mo><msub><mi>B</mi><mrow><mi>&lambda;</mi><mo>,</mo><mi>&infin;</mi></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mo>(</mo><mn>6</mn><mo>)</mo><mo>.</mo></mrow>]]></math><img file="FDA0000937960010000022.GIF" wi="1069" he="164" /></maths>
地址 300072 天津市南开区卫津路92号