发明名称 基于大气反射-散射原理构建模型复原夜间交通视频车灯的方法
摘要 本发明涉及的是基于大气反射‑散射原理构建复原模型,提取夜间交通视频中车灯的方法。首先根据漫反射原理抑制路面的漫反射光,应用暗原色先验理论估计环境光和透射率,得到只含有车灯与路面强反射光的复原图像,对复原图像再利用暗原色先验理论重新估计环境光,进一步抑制其中的强光光晕,最后根据车灯的几何特征筛选出复原图像中的车灯,有效检测车灯。验证实验结果表明,本方法能够适应多种复杂光照下的车灯检测,能够用于挡风玻璃较暗且前灯照明功率过大这类特殊情况下的车灯提取,能够有效降低路面强反射光区域以及车身镜面反射对车灯检测造成的影响。
申请公布号 CN105844595A 申请公布日期 2016.08.10
申请号 CN201610146282.5 申请日期 2016.03.14
申请人 天津工业大学 发明人 汤春明;曹志升;涂丽平;林祥清;于翔
分类号 G06T5/00(2006.01)I;G06K9/00(2006.01)I 主分类号 G06T5/00(2006.01)I
代理机构 代理人
主权项 一种根据大气反射‑散射原理构建复原模型,提取夜间交通视频车灯的方法,所述方法包括以下步骤:A.抑制漫反射输入彩色视频图像I(x)的最大漫反射色度图D(x)可表示为:<maths num="0001"><math><![CDATA[<mrow><mi>D</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><mi>I</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>-</mo><mfrac><mrow><mover><mi>&lambda;</mi><mo>~</mo></mover><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>I</mi><mi>r</mi></msub><mo>(</mo><mi>x</mi><mo>)</mo><mo>+</mo><msub><mi>I</mi><mi>g</mi></msub><mo>(</mo><mi>x</mi><mo>)</mo><mo>+</mo><msub><mi>I</mi><mi>b</mi></msub><mo>(</mo><mi>x</mi><mo>)</mo><mo>)</mo></mrow><mo>-</mo><mover><mi>I</mi><mo>~</mo></mover><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow></mrow><mrow><mn>3</mn><mover><mi>&lambda;</mi><mo>~</mo></mover><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>-</mo><mn>1</mn></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FSA0000127683270000011.GIF" wi="1400" he="124" /></maths>其中,<img file="FSA0000127683270000012.GIF" wi="1018" he="65" />是最大漫反射色度,d<sub>r</sub>(x)、d<sub>g</sub>(x)、d<sub>b</sub>(x)分别是R、G、B三个颜色通道中的漫反射分量,I<sub>r</sub>(x),I<sub>g</sub>(x),I<sub>b</sub>(x)是I(x)中的三个颜色通道,<img file="FSA0000127683270000013.GIF" wi="270" he="54" /><img file="FSA0000127683270000014.GIF" wi="204" he="48" />令E(x)=I(x)‑D(x);B.构建模型计算机视觉及计算机图形学领域广泛采用的大气散射模型表示为:E(x)=R(x)t(x)+E<sub>∞</sub>(x)(1‑t(x))           (2)其中,t(x)=e<sup>‑βd(x)</sup>为透射率(β是散射系数,d(x)是场景点x的深度值),E<sub>∞</sub>(x)是参与成像的大气光,R(x)是清晰图像,E(x)是经过漫反射抑制后图像中点x的亮度,等号右端第一项为衰减项,第二项为天空光;由于参与成像的不仅含有大气光,而且还有环境光,将公式(2)扩展成式(3):E(x)=R(x)t(x)+(E<sub>∞</sub>(x)+A(x))(1‑t(x))          (3)其中,A(x)是参与成像的环境光,E<sub>∞</sub>(x)是参与成像的大气光;C.估计参量公式(3)复原R(x)时,需要对未知参量t(x)、E<sub>∞</sub>(x)和A(x)进行估计:首先根据暗原色先验理论估计环境光A(x),然后根据大气反射‑散射模型重新定义透射率t(x)并估计大气光E<sub>∞</sub>(x);①图像E的暗通道E<sup>dark</sup>定义如公式(4):<maths num="0002"><math><![CDATA[<mrow><msup><mi>E</mi><mrow><mi>d</mi><mi>a</mi><mi>r</mi><mi>k</mi></mrow></msup><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><munder><mi>min</mi><mrow><mi>c</mi><mo>&Element;</mo><mo>{</mo><mi>r</mi><mo>,</mo><mi>g</mi><mo>,</mo><mi>b</mi><mo>}</mo></mrow></munder><mrow><mo>(</mo><munder><mi>min</mi><mrow><mi>y</mi><mo>&Element;</mo><mi>&Omega;</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow></mrow></munder><mo>(</mo><mrow><msup><mi>E</mi><mi>c</mi></msup><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow></mrow><mo>)</mo><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FSA0000127683270000015.GIF" wi="1379" he="84" /></maths>其中,y是邻域变量,c是颜色通道变量,r、g、b分别代表红、绿、监三个颜色通道,Ω(x)是以x为中心的邻域,这里取15*15邻域,E<sup>c</sup>代表E的某个颜色通道;环境光A(x)可以借助暗通道图估计,求解过程如下:Step1原始图像I(x)应用漫反射抑制,得到漫反射抑制图E(x);Step2对漫反射抑制图E(x)取暗通道E<sup>dark</sup>;Step3将暗通道处理后的图像E<sup>dark</sup>中像素的亮度值按从大到小顺序排列,取前0.5%最亮的像素点,这些点对应于原图像I(x)中最大亮度值区域,用它们在图像中的灰度估计环境光A(x);②由于夜间大气光强度已经降的很低,认为大气光E<sub>∞</sub>(x)接近于0;③透射率为环境光占所有成像光源的比例,可表示为:t(x)=ω*A(x)/E(x)           (5)其中,ω是调节系数,防止得到的复原图像远景区域变模糊;D.对t(x)应用导向滤波平滑处理;E.复原图像由公式(3)可以得到复原后的图像R(x)如公式(6)所示:<maths num="0003"><math><![CDATA[<mrow><mi>R</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>E</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>-</mo><mrow><mo>(</mo><msub><mi>E</mi><mi>&infin;</mi></msub><mo>(</mo><mi>x</mi><mo>)</mo><mo>+</mo><mi>A</mi><mo>(</mo><mi>x</mi><mo>)</mo><mo>)</mo></mrow></mrow><mrow><mi>m</mi><mi>a</mi><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>(</mo><mi>x</mi><mo>)</mo><mo>,</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow></mrow></mfrac><mo>+</mo><mrow><mo>(</mo><msub><mi>E</mi><mi>&infin;</mi></msub><mo>(</mo><mi>x</mi><mo>)</mo><mo>+</mo><mi>A</mi><mo>(</mo><mi>x</mi><mo>)</mo><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FSA0000127683270000021.GIF" wi="1441" he="141" /></maths>其中,t<sub>0</sub>为防止分母为零而设;F.为了抑制R(x)图中产生的强光光晕,利用公式(4),对R(x)遍历取15*15的暗通道,重新估计环境光A(x),并假定同一场景中的透射率t(x)不变;G.用公式(6)去掉光晕,得到最终的复原图R’(x);H.车灯筛选复原图R’(x)中存在非车灯噪声,需要进一步滤除,具体方法如下:根据所拍摄的视频,车头灯的形状大致分为以下四类:(1)圆形车灯(2)椭圆形车灯(3)近似圆形车灯(4)方形车灯首先根据面积大小和长宽比例,去除面积较大且形状不规则的反射光区域:将面积小于阈值Th的亮斑保留,这里Th取300个像素,一般车灯近似圆形,这里取长宽比的范围在(0.8,1.2);然后根据亮斑周长去除明显小亮斑:去除大面积反射区域以后,取每一帧视频中亮斑周长P的平均值P<sub>mean</sub>,如果P<0.5P<sub>mean</sub>,这类亮斑将被抑制掉;最后根据亮斑圆度M作进一步细化处理:由于车头灯的M大多近似为1,通过公式(7)定义圆度偏心率temp,如果|1‑M|<temp,说明亮斑近似圆形,将被保留;<maths num="0004"><math><![CDATA[<mrow><mi>t</mi><mi>e</mi><mi>m</mi><mi>p</mi><mo>=</mo><msqrt><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mrow><mo>(</mo><mo>|</mo><mrow><mn>1</mn><mo>-</mo><msub><mi>M</mi><mi>i</mi></msub></mrow><mo>|</mo><mo>-</mo><mi>m</mi><mi>e</mi><mi>a</mi><mi>n</mi><mo>(</mo><mrow><mo>|</mo><mrow><mn>1</mn><mo>-</mo><msub><mi>M</mi><mi>i</mi></msub></mrow><mo>|</mo></mrow><mo>)</mo><mo>)</mo></mrow><mo>/</mo><mi>N</mi></mrow></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FSA0000127683270000022.GIF" wi="1369" he="159" /></maths>其中,temp是圆度偏心率,N是当前图像中总车灯数,M是第i个亮斑的圆度,mean(|1‑M<sub>i</sub>|)表示取|1‑M<sub>i</sub>|的均值。
地址 300387 天津市西青区宾水西道399号