CN104732235B - 一种消除夜间道路反光干扰的车辆检测方法 - Google Patents
一种消除夜间道路反光干扰的车辆检测方法 Download PDFInfo
- Publication number
- CN104732235B CN104732235B CN201510121596.5A CN201510121596A CN104732235B CN 104732235 B CN104732235 B CN 104732235B CN 201510121596 A CN201510121596 A CN 201510121596A CN 104732235 B CN104732235 B CN 104732235B
- Authority
- CN
- China
- Prior art keywords
- mrow
- value
- image
- vehicle
- light
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000001514 detection method Methods 0.000 claims abstract description 42
- 238000012544 monitoring process Methods 0.000 claims abstract description 25
- 238000001914 filtration Methods 0.000 claims abstract description 9
- 230000011218 segmentation Effects 0.000 claims abstract description 8
- 230000008030 elimination Effects 0.000 claims abstract description 5
- 238000003379 elimination reaction Methods 0.000 claims abstract description 5
- 238000012545 processing Methods 0.000 claims description 13
- 238000006243 chemical reaction Methods 0.000 claims description 9
- 238000007781 pre-processing Methods 0.000 claims description 8
- 238000012216 screening Methods 0.000 claims description 7
- 238000003708 edge detection Methods 0.000 claims description 6
- 238000003709 image segmentation Methods 0.000 claims description 6
- 238000004891 communication Methods 0.000 claims description 5
- 238000000605 extraction Methods 0.000 claims description 5
- 238000009499 grossing Methods 0.000 claims description 3
- 230000005855 radiation Effects 0.000 claims description 2
- 238000002405 diagnostic procedure Methods 0.000 abstract 1
- 238000005286 illumination Methods 0.000 description 4
- 238000012850 discrimination method Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
Landscapes
- Image Processing (AREA)
- Image Analysis (AREA)
- Traffic Control Systems (AREA)
Abstract
本发明公开了一种消除夜间道路反光干扰的车辆检测方法。本发明利用大气散射原理和LOG滤波原理对道路监控图像进行预处理,得到反射强光图和反射弱光图,通过对车灯和反射光的样本统计获得区分车灯和反射光的Fisher统计特征,利用Fisher判别法区分图像中的车灯与反射光,对图像中的反射光进行分割并消除;然后,利用车灯形状特征从高亮连通区域中提取出车灯区域,建立车灯匹配准则,对车灯进行匹配,完成车辆的检测和统计。本发明采用的Fisher分类方法区分反射光并进行分割消除,能有效排除路面反射光对夜间车辆检测的干扰,从而提高了车辆检测率;本方法的车辆检测不需要依赖其他参照物,具有很好的夜间道路环境的适应性。
Description
技术领域
本发明属于图像处理技术领域,尤其涉及一种消除夜间道路反光干扰的车辆检测方法。
背景技术
近些年来,针对基于监控视频的白天道路车辆检测技术已经较为成熟;而针对夜间道路监控的车辆检测还存在检测准确率低,易受到道路环境光照影响。
夜间道路监控的车辆测算法主要分为:(1)基于车辆特征分类的方法,这种方法虽然受噪声影响较小,但是在外部光照条件较暗的情况下,车辆颜色、轮廓、形状等特征不明显,导致检测率低、鲁棒性差;(2)基于车尾灯颜色的车辆检测方法,这种方法易受周围环境光照和摄像头过度曝光影响,车尾灯的颜色有时候会变成单一的白色,因此,环境适应性和鲁棒性较差;(3)基于车前灯的车辆检测方法,由于车前灯是车辆最显著的特征,而且车灯检测算法复杂度上比特征分类方法低,比车尾灯检测方法的环境光照适应性更好。
但是,基于车前灯的车辆检测算法易受到车灯的路面反射光干扰,尤其是雨天道路场景中的干扰更加明显,影响到夜间车辆检测的准确率和鲁棒性。对于这种情况,已有的常用处理方法是根据车灯与反射光的形状、大小以及车辆之间的最小距离来区分车灯和车灯反射光。此类方法虽然能在一定程度上提高车辆检测率,但仍然不能从根本上解决反射光对夜间车辆检测的干扰问题。
因此,我们提出利用光线大气散射原理和LOG滤波对道路监控图像进行处理,根据图像中的车灯和反射光的亮度梯度的统计特征差别,采用基于Fisher判别法区分和分割消除图像中的反射光,提高夜间车辆检测准确率。
发明内容
针对上述问题,本发明公开了针对夜间道路监控图像的一种基于Fisher判别法的反射光区分、分割消除处理,以及车灯检测、匹配和跟踪的方法。
本发明解决其技术问题采用的技术方案步骤如下:
一种消除夜间道路反光干扰的车辆检测方法,包括如下步骤:
步骤1、对每帧道路监控图像进行预处理,得到预处理后的二值图像;
步骤1所述的预处理包括感兴趣区域提取、图像分割以及反射光消除,具体如下:
1-1.感兴趣区域提取和灰度转换;
从每帧道路监控图像中截取一段区域图像,对截取的区域图像进行车辆检测;然后对区域图像进行灰度转换得到灰度图I,具体转换公式如下:
Gray=R×0.299+G×0.587+B×0.114
其中,彩色位图由R/G/B组成;Gray为灰度值,其范围在0到255之间;
1-2.利用基于大气散射的原理对灰度图I进行处理,得到反射强光图RI,大气散射E(d,λ)如下定义:
E(d,λ)=I0(λ)×γ(λ)×exp(-d);
其中,I0为光源辐射强度,γ(λ)为波长λ的总散射系数,d为光源到摄像头屏幕的距离;
1-2-1.将灰度图I的像素点(x,y)划分成内部区域和外部区域其中划分半径r=2;
1-2-2.在内部区域中寻找像素最小值和最大值散射系数γ(x,y)如下定义:
其中,为最小值和最大值之间的距离;
1-2-3.在外部区域中找到像素最小值和最大值同时计算最小值和最大值之间的距离
1-2-4.最后得到反射强光图RI:
1-3.利用高斯拉普拉斯滤波对灰度图I进行处理,得到反射弱光图RS;
1-3-1.对灰度图I取反色,然后利用高斯滤波对取反色后的灰度图进行平滑处理,其中高斯函数为:
其中σ为标准偏差,取值为1.0;
1-3-2.利用拉普拉斯算子对平滑后的灰度图进行边缘检测;
1-3-3.最后通过扫面线填充算法对边缘检测后的灰度图进行填充,得到最终的反射弱光图RS;
1-4.选取并建立车灯和反射光的样本;
选取多张具有代表性的有反射光和车灯的夜间道路监控图像作为样本;每张图像都进行步骤1-2和1-3操作,得到灰度图I、反射强光图RI、反射弱光图RS三个位置对应的样本图,以灰度图I中车灯所在位置的像素值作为正样本,以灰度图I反射光所在位置的像素值作为负样本;
1-5.对所有样本进行统计,建立Fisher分类方程;
利用Fisher准则和步骤1-4的样本,统计计算出Fisher分类的直线方程y=wT·x中的参数:w和判别临界值;
1-6.根据亮度直方图双峰方法对每帧道路监控图像进行阈值分割处理,得到图像中的高亮区域;
根据直方图双峰方法,对需要检测的每帧道路监控图像提取其灰度直方图,首先找到灰度直方图中第一个峰值和第二个峰值,再找到第一个峰值和第二个峰值之间的谷值T,谷值T就是图像分割的阈值;再用如下公式对每帧道路监控图像进行分割;
f(x,y)=0,if(f(x,y)<T)
f(x,y)=255,if(f(x,y)≥T)
其中,f(x,y)为像素值;
1-7.利用Fisher判别方对分割后得到的高亮区域进行反射光分类和消除;
利用步骤1-5得到的Fisher分类参数,对大于谷值T的像素进行反射光分类区分,将判别为负样本的反射光区域像素值设为零,消除反射光,得到消除反射光干扰后的二值图像;对判别为正样本的车灯区域像素值不作处理;
步骤2、基于车灯形态的车灯检测处理,对车辆车灯进行检测具体采用如下步骤:
2-1.对二值图像采用八连通区域提取,在二值图像中的感兴趣区域内查找大于阈值T的像素点,找到后以这个像素点为目标像素,将目标像素的x值赋给这个区域的最左、最右值,即left和right,将目标像素的y值赋给这个区域的最上和最下值,即top和bottom;对目标像素进行上、下、左、右、左上角、左下角、右下角、右上角8个方向的像素点搜索,如果找到比left更小的值,则将这个更小的值赋值给left,同理,找到比top更小的值,则将这个更小的值赋值给top;而如果找到比right更大的值,则将这个更大的值赋值给right,同理,如果找到比bottom更大的值,则将这个更大的值赋值给bottom;最后,会检测出一个最左值、最右值分别为left、right,最上值和最下值分别为top和bottom的连通矩形区域;一副图像检测完之后,会检测出N个连通矩形区域C{C1,C2,…,Ci,…,Cn},每个连通矩形区域均具有宽度W和高度H;
2-2.对连通矩形区域进行筛选,假设有两个连通矩形区域Ci和Cj,具体筛选方式如下:
2-2-1.根据车灯区域的饱满度得到了如下的判断条件:
Si′≥Si/2Si>t;
其中Si为连通区域形成的矩形面积,Si′为连通区域的实际面积,t=12;
2-2-2.同一车辆的两个车灯处于相似水平线上时:
Ci,t≤Cj,b,Ci,b≤Cj,t,其中,Ci,t、Ci,b为Ci的顶部、底部;Cj,b和Cj,t为Cj的顶部、底部;
2-2-3.同一车辆的两个车灯其宽度与高度相似如下:T1<h<T2,T1<w<T2,其中,h、w分别为两区域的高度比例和宽度比例;T1和T2为判断阈值,当T1和T2分别为0.7和1.3最佳;
2-2-4.同一车辆两车灯的面积相近:
其中Si'、Sj'为两区域的面积;
2-2-5.满足以上四个条件并检测出来的为车灯;
步骤3、车辆检测
步骤2中已经加入了两车灯之间的高度、大小对比,根据车辆上两个车灯之间具有相似的高度、相似的大小以及一定的距离,在步骤2的条件上加入车灯距离的判断就能完成对车灯的匹配;当两个车灯满足步骤2-2的条件并且两个车灯距离在一个车灯区域宽度的2倍到8倍之间,则这两个车灯属于同一车辆,从而完成对车辆的检测;
步骤4、对检测到的车灯进行跟踪,具体步骤如下:
4-1.计算属于同一辆车中两个车灯在当前帧与上一帧的运动偏移量
其中,i为帧数;(x,y)作为车灯的中心点;
4-2.建立以()为中心坐标,宽和高分别为1.5×Wi-1和3×Hi-1的窗口进行搜索;对窗口中搜索到的连通矩形区域以步骤2-2车灯提取的判断条件对车灯进行判断,若是车灯,则转入步骤4-3;若不是车灯则结束;
4-3.判断两帧中的车灯是否为同一个车灯,若满足A1/2≤A2≤2×A1,则两帧中的车灯为同一个车灯;
其中,A2为在搜索区域中检测到的连通块面积,A1为上一帧的车灯连通块面积。
本发明有益效果如下:
本发明以固定的道路监控摄像机为基础,通过利用光线大气散射原理和LOG滤波对道路监控图像进行预处理,根据图像中的车灯和反射光的亮度梯度的统计特征差别,采用基于Fisher判别法区分图像中的反射光和车灯,并消除图像中的反射光;然后,进行连通区域提取并利用形态学对车灯进行提取和匹配,实现夜间车辆的检测和跟踪。
本发明提出了一种有效的解决夜间车灯的路面反射光干扰问题,能够有效排除路面反射光对夜间车辆检测的干扰,有效提高车辆检测率,与其他方法相比具有更好的夜间道路环境适应性。
本发明是在摄像头拍摄到的视频大小为640×480,而且从拍摄的角度来看,车辆的经过的一般都是从上往下运动,因此本发明提出一个判断公式:Topi-1<T。其中,Topi-1为在搜索区域中检测到的连通块的上边,当Bottomi-1>355时,T=Bottomi-1+10。当Bottomi-1<355时,T=Bottomi-1+5。这种方法只需要计算一个小范围内的连通区域,减少了运算量以及排除了窗口之外的很多干扰,增加了车灯检测的准确率。而且,由于视频拍摄的角度的问题,车辆离摄像头越近,偏移量越大,因此本发明采用的这个判断公式是分段判断的,这样能提高跟踪的检测准确率。
在车灯跟踪时可能会出现短暂的目标丢失。因此,本发明采用了在图像中选择一小段区域的计数方法,为检测到的车辆标记为0,当车辆经过这段区域时如标记为0则计数器加一,并将其标记为1,如标记为1,则不予处理。这样可以将经过的车辆统计出来,进行车流量统计。
附图说明
图1为本发明进行车辆检测和跟踪系统的流程图。
图2为夜间图像预处理的流程图。
图3为本发明对车辆车灯进行筛选的流程图。
图4为车灯跟踪的流程图。
具体实施方式
下面结合附图,对本发明的具体实施方案作进一步详细描述。
如图1所示,一种消除夜间道路反光干扰的车辆检测方法,具体步骤如下:
步骤101:对每帧道路监控图像进行预处理,得到预处理后的二值图像。
如图2所示,步骤101所述的预处理包括感兴趣区域提取、图像分割以及反射光消除,具体步骤如下:
步骤201:感兴趣区域提取和灰度转换。从每帧道路监控图像中截取一段区域图像,对截取的区域图像进行车辆检测;然后对区域图像进行灰度转换得到灰度图I,具体转换公式如下:
Gray=R×0.299+G×0.587+B×0.114
其中,彩色位图由R/G/B组成;Gray为灰度值,其范围在0到255之间。
步骤202:利用基于大气散射的原理对灰度图I进行处理,得到反射强光图RI,大气散射E(d,λ)如下定义:
E(d,λ)=I0(λ)×γ(λ)×exp(-d);
其中,I0为光源辐射强度,γ(λ)为波长λ的总散射系数,d为光源到摄像头屏幕的距离。
将灰度图I的像素点(x,y)划分成内部区域和外部区域其中划分半径r=2。
在内部区域中寻找像素最小值和最大值散射系数γ(x,y)如下定义:
其中,为最小值和最大值之间的距离。
在外部区域中找到像素最小值和最大值同时计算最小值和最大值之间的距离
最后得到反射强光图RI:
步骤203:利用高斯拉普拉斯(Laplacian of Gaussian,LOG)滤波对灰度图I进行处理,得到反射弱光图RS。
首先对灰度图I取反色,然后利用高斯滤波对取反色后的灰度图进行平滑处理,其中高斯函数为:其中σ为标准偏差,取值为1.0。
接着利用拉普拉斯算子对平滑后的灰度图进行边缘检测;
最后通过扫面线填充算法对边缘检测后的灰度图进行填充,得到最终的反射弱光图RS。
步骤204:选取并建立车灯和反射光的样本。
选取多张具有代表性的有反射光和车灯的夜间道路监控图像作为样本。每张图像都进行步骤202和203操作,得到灰度图I、反射强光图RI、反射弱光图RS三个位置对应的样本图,以灰度图I中车灯所在位置的像素值作为正样本,以灰度图I反射光所在位置的像素值作为负样本。
步骤205:对所有样本进行统计,建立Fisher分类方程。
利用Fisher准则和步骤204的样本,统计计算出Fisher分类的直线方程y=wT·x中的参数:w和判别临界值。
步骤206:根据亮度直方图双峰方法对每帧道路监控图像进行阈值分割处理,得到图像中的高亮区域。
根据直方图双峰方法,对需要检测的每帧道路监控图像提取其灰度直方图,首先找到灰度直方图中第一个峰值和第二个峰值,再找到第一个峰值和第二个峰值之间的谷值T,谷值T就是图像分割的阈值。再用如下公式对每帧道路监控图像进行分割。
f(x,y)=0,if(f(x,y)<T)
f(x,y)=255,if(f(x,y)≥T)
其中,f(x,y)为像素值;
步骤207:利用Fisher判别方对分割后得到的高亮区域进行反射光分类和消除。
利用步骤205得到的Fisher分类参数,对大于谷值T的像素进行反射光分类区分,将判别为负样本的反射光区域像素值设为零,消除反射光,得到消除反射光干扰后的二值图像;对判别为正样本的车灯区域像素值不作处理。
步骤102:基于车灯形态的车灯检测处理
如图3所示,对车辆车灯进行检测具体采用如下步骤:
步骤301:对二值图像采用八连通区域提取,在二值图像中的感兴趣区域内查找大于阈值T的像素点,找到后以这个像素点为目标像素,将目标像素的x值赋给这个区域的最左、最右值,即left和right,将目标像素的y值赋给这个区域的最上和最下值,即top和bottom。对目标像素进行上、下、左、右、左上角、左下角、右下角、右上角8个方向的像素点搜索,如果找到比left更小的值,则将这个更小的值赋值给left,同理,找到比top更小的值,则将这个更小的值赋值给top;而如果找到比right更大的值,则将这个更大的值赋值给right,同理,如果找到比bottom更大的值,则将这个更大的值赋值给bottom;最后,会检测出一个最左值、最右值分别为left、right,最上值和最下值分别为top和bottom的连通矩形区域。一副图像检测完之后,会检测出N个连通矩形区域C{C1,C2,…,Ci,…,Cn},每个连通矩形区域均具有宽度W和高度H。
步骤302:对连通矩形区域进行筛选,假设有两个连通矩形区域Ci和Cj,具体筛选方式如下:
第一,根据车灯区域的饱满度得到了如下的判断条件:
Si′≥Si/2Si>t;
其中Si为连通区域形成的矩形面积,Si′为连通区域的实际面积,t=12;
第二,同一车辆的两个车灯处于相似水平线上时:
Ci,t≤Cj,b,Ci,b≤Cj,t,其中,Ci,t、Ci,b为Ci的顶部、底部;Cj,b和Cj,t为Cj的顶部、底部;
第三,同一车辆的两个车灯其宽度与高度相似如下:T1<h<T2,T1<w<T2,其中,h、w分别为两区域的高度比例和宽度比例。T1和T2为判断阈值,当T1和T2分别为0.7和1.3最佳;
第四,同一车辆两车灯的面积相近:T1和T2 Si'、Sj'为两区域的面积。
满足以上四个条件并检测出来的为车灯。
步骤103:车辆检测
根据车辆上两个车灯之间具有相似的高度、相似的大小以及一定的距离,在步骤102中已经加入了两车灯之间的高度、大小对比。因此,本发明在步骤102的条件上加入车灯距离的判断就能完成对车灯的匹配。
当两个车灯满足步骤302的条件并且两个车灯距离在一个车灯区域宽度的2倍到8倍之间,则这两个车灯属于同一车辆,从而完成对车辆的检测。
步骤104:对检测到的车灯进行跟踪。
本发明采用如图4所示的跟踪过程进行跟踪,具体步骤如下:
步骤401:计算属于同一辆车中两个车灯在当前帧与上一帧的运动偏移量
其中,i为帧数。x,y均作为车灯的中心点。
步骤402:建立以()为中心坐标,宽和高分别为1.5×Wi-1和3×Hi-1的窗口进行搜索。对窗口中搜索到的连通区域以步骤302车灯提取的判断条件对车灯进行判断,若是车灯,则转入步骤403;若不是车灯则结束。
步骤403:判断两帧中的车灯是否为同一个车灯,若满足A1/2≤A2≤2×A1,则两帧中的车灯为同一个车灯。
其中,A2为在搜索区域中检测到的连通块面积,A1为上一帧的车灯连通块面积。
本发明提出了一种有效的解决夜间车灯的路面反射光干扰问题,可以有效排除路面反射光对夜间车辆检测的干扰,有效提高车辆检测率,与其他方法相比具有更好的夜间道路环境适应性。
本发明是在摄像头拍摄到的视频大小为640×480,而且从拍摄的角度来看,车辆的经过的一般都是从上往下运动,因此本发明提出一个判断公式:Topi-1<T。其中,Topi-1为在搜索区域中检测到的连通块的上边,当Bottomi-1>355时,T=Bottomi-1+10。当Bottomi-1<355时,T=Bottomi-1+5。这种方法只需要计算一个小范围内的连通区域,减少了运算量以及排除了窗口之外的很多干扰,增加了车灯检测的准确率。而且,由于视频拍摄的角度的问题,车辆离摄像头越近,偏移量越大,因此本发明采用的这个判断公式是分段判断的,这样能提高跟踪的检测准确率。
在车灯跟踪时可能会出现短暂的目标丢失。因此,本发明采用了在图像中选择一小段区域的计数方法,为检测到的车辆标记为0,当车辆经过这段区域时如标记为0则计数器加一,并将其标记为1,如标记为1,则不予处理。这样可以将经过的车辆统计出来,进行车流量统计。
Claims (1)
1.一种消除夜间道路反光干扰的车辆检测方法,其特征在于包括如下步骤:
步骤1、对每帧道路监控图像进行预处理,得到预处理后的二值图像;
步骤1所述的预处理包括感兴趣区域提取、图像分割以及反射光消除,具体如下:
1-1.感兴趣区域提取和灰度转换;
从每帧道路监控图像中截取一段区域图像,对截取的区域图像进行车辆检测;然后对区域图像进行灰度转换得到灰度图I,具体转换公式如下:
Gray=R×0.299+G×0.587+B×0.114
其中,彩色位图由R/G/B组成;Gray为灰度值,其范围在0到255之间;
1-2.利用基于大气散射的原理对灰度图I进行处理,得到反射强光图RI,大气散射E(d,λ)如下定义:
E(d,λ)=I0(λ)×γ(λ)×exp(-d);
其中,I0为光源辐射强度,γ(λ)为波长λ的总散射系数,d为光源到摄像头屏幕的距离;
1-2-1.将灰度图I的像素点(x,y)划分成内部区域和外部区域其中划分半径r=2;
<mrow>
<msubsup>
<mi>&Theta;</mi>
<mrow>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
</mrow>
<mi>i</mi>
</msubsup>
<mo>=</mo>
<mo>{</mo>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>+</mo>
<mi>u</mi>
<mo>,</mo>
<mi>y</mi>
<mo>+</mo>
<mi>v</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mo>-</mo>
<mi>r</mi>
<mo>&le;</mo>
<mi>u</mi>
<mo>&le;</mo>
<mi>r</mi>
<mo>,</mo>
<mo>-</mo>
<mi>r</mi>
<mo>&le;</mo>
<mi>v</mi>
<mo>&le;</mo>
<mi>r</mi>
<mo>}</mo>
</mrow>
<mrow>
<msubsup>
<mi>&Theta;</mi>
<mrow>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
</mrow>
<mi>e</mi>
</msubsup>
<mo>=</mo>
<mo>{</mo>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>+</mo>
<mi>u</mi>
<mo>,</mo>
<mi>y</mi>
<mo>+</mo>
<mi>v</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mo>-</mo>
<mn>2</mn>
<mi>r</mi>
<mo>&le;</mo>
<mi>u</mi>
<mo>&le;</mo>
<mn>2</mn>
<mi>r</mi>
<mo>,</mo>
<mo>-</mo>
<mn>2</mn>
<mi>r</mi>
<mo>&le;</mo>
<mi>v</mi>
<mo>&le;</mo>
<mn>2</mn>
<mi>r</mi>
<mo>}</mo>
</mrow>
1-2-2.在内部区域中寻找像素最小值和最大值散射系数γ(x,y)如下定义:
<mrow>
<mi>&gamma;</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msubsup>
<mi>MI</mi>
<mrow>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
</mrow>
<mi>i</mi>
</msubsup>
</mrow>
<mrow>
<msubsup>
<mi>MA</mi>
<mrow>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
</mrow>
<mi>i</mi>
</msubsup>
<mo>&times;</mo>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msubsup>
<mi>&epsiv;</mi>
<mrow>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
</mrow>
<mi>i</mi>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
其中,为最小值和最大值之间的距离;
1-2-3.在外部区域中找到像素最小值和最大值同时计算最小值和最大值之间的距离
1-2-4.最后得到反射强光图RI:
<mrow>
<mi>R</mi>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>|</mo>
<msubsup>
<mi>MI</mi>
<mrow>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
</mrow>
<mi>e</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>MA</mi>
<mrow>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
</mrow>
<mi>e</mi>
</msubsup>
<mo>&times;</mo>
<mi>&gamma;</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msubsup>
<mi>&epsiv;</mi>
<mrow>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
</mrow>
<mi>e</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mo>|</mo>
</mrow>
1-3.利用高斯拉普拉斯滤波对灰度图I进行处理,得到反射弱光图RS;
1-3-1.对灰度图I取反色,然后利用高斯滤波对取反色后的灰度图进行平滑处理,其中高斯函数为:
<mrow>
<mi>G</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>y</mi>
<mn>2</mn>
</msup>
<mo>-</mo>
<mn>2</mn>
<msup>
<mi>&sigma;</mi>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
<mi>&sigma;</mi>
</mrow>
</mfrac>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mfrac>
<mrow>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>y</mi>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mn>2</mn>
<msup>
<mi>&sigma;</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
其中σ为标准偏差,取值为1.0;
1-3-2.利用拉普拉斯算子对平滑后的灰度图进行边缘检测;
1-3-3.最后通过扫面线填充算法对边缘检测后的灰度图进行填充,得到最终的反射弱光图RS;
1-4.选取并建立车灯和反射光的样本;
选取多张具有代表性的有反射光和车灯的夜间道路监控图像作为样本;每张图像都进行步骤1-2和1-3操作,得到灰度图I、反射强光图RI、反射弱光图RS三个位置对应的样本图,以灰度图I中车灯所在位置的像素值作为正样本,以灰度图I反射光所在位置的像素值作为负样本;
1-5.对所有样本进行统计,建立Fisher分类方程;
利用Fisher准则和步骤1-4的样本,统计计算出Fisher分类的直线方程y=wT·x中的参数:w和判别临界值;
1-6.根据亮度直方图双峰方法对每帧道路监控图像进行阈值分割处理,得到图像中的高亮区域;
根据直方图双峰方法,对需要检测的每帧道路监控图像提取其灰度直方图,首先找到灰度直方图中第一个峰值和第二个峰值,再找到第一个峰值和第二个峰值之间的谷值T,谷值T就是图像分割的阈值;再用如下公式对每帧道路监控图像进行分割;
f(x,y)=0,if(f(x,y)<T)
f(x,y)=255,if(f(x,y)≥T)
其中,f(x,y)为像素值;
1-7.利用Fisher判别方对分割后得到的高亮区域进行反射光分类和消除;
利用步骤1-5得到的Fisher分类参数,对大于谷值T的像素进行反射光分类区分,将判别为负样本的反射光区域像素值设为零,消除反射光,得到消除反射光干扰后的二值图像;对判别为正样本的车灯区域像素值不作处理;
步骤2、基于车灯形态的车灯检测处理,对车辆车灯进行检测具体采用如下步骤:
2-1.对二值图像采用八连通区域提取,在二值图像中的感兴趣区域内查找大于阈值T的像素点,找到后以这个像素点为目标像素,将目标像素的x值赋给这个区域的最左、最右值,即left和right,将目标像素的y值赋给这个区域的最上和最下值,即top和bottom;对目标像素进行上、下、左、右、左上角、左下角、右下角、右上角8个方向的像素点搜索,如果找到比left更小的值,则将这个更小的值赋值给left,同理,找到比top更小的值,则将这个更小的值赋值给top;而如果找到比right更大的值,则将这个更大的值赋值给right,同理,如果找到比bottom更大的值,则将这个更大的值赋值给bottom;最后,会检测出一个最左值、最右值分别为left、right,最上值和最下值分别为top和bottom的连通矩形区域;一副图像检测完之后,会检测出N个连通矩形区域C{C1,C2,…,Ci,…,Cn},每个连通矩形区域均具有宽度W和高度H;
2-2.对连通矩形区域进行筛选,假设有两个连通矩形区域Ci和Cj,具体筛选方式如下:
2-2-1.根据车灯区域的饱满度得到了如下的判断条件:
Si′≥Si/2 Si>t;
其中Si为连通区域形成的矩形面积,Si′为连通区域的实际面积,t=12;
2-2-2.同一车辆的两个车灯处于相似水平线上时:
Ci,t≤Cj,b,Ci,b≤Cj,t,其中,Ci,t、Ci,b为Ci的顶部、底部;Cj,b和Cj,t为Cj的顶部、底部;
2-2-3.同一车辆的两个车灯其宽度与高度相似如下:T1<h<T2,T1<w<T2,其中,h、w分别为两区域的高度比例和宽度比例;T1和T2为判断阈值,且当T1和T2分别为0.7和1.3;
2-2-4.同一车辆两车灯的面积相近:
其中S′i、S′j为两区域的面积;
2-2-5.满足以上四个条件并检测出来的为车灯;
步骤3、车辆检测
步骤2中已经加入了两车灯之间的高度、大小对比,根据车辆上两个车灯之间具有相似的高度、相似的大小以及一定的距离,在步骤2的条件上加入车灯距离的判断就能完成对车灯的匹配;当两个车灯满足步骤2-2的条件并且两个车灯距离在一个车灯区域宽度的2倍到8倍之间,则这两个车灯属于同一车辆,从而完成对车辆的检测;
步骤4、对检测到的车灯进行跟踪,具体步骤如下:
4-1.计算属于同一辆车中两个车灯在当前帧与上一帧的运动偏移量
<mrow>
<msubsup>
<mi>&Delta;x</mi>
<mi>k</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>=</mo>
<msup>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>-</mo>
<msup>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msup>
</mrow>
<mrow>
<msubsup>
<mi>&Delta;y</mi>
<mi>k</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>=</mo>
<msup>
<mi>y</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>-</mo>
<msup>
<mi>y</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msup>
</mrow>
其中,i为帧数;(x,y)作为车灯的中心点;
4-2.建立以为中心坐标,宽和高分别为1.5×Wi-1和3×Hi-1的窗口进行搜索;对窗口中搜索到的连通矩形区域以步骤2-2车灯提取的判断条件对车灯进行判断,若是车灯,则转入步骤4-3;若不是车灯则结束;
4-3.判断两帧中的车灯是否为同一个车灯,若满足A1/2≤A2≤2×A1,则两帧中的车灯为同一个车灯;
其中,A2为在搜索区域中检测到的连通块面积,A1为上一帧的车灯连通块面积。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510121596.5A CN104732235B (zh) | 2015-03-19 | 2015-03-19 | 一种消除夜间道路反光干扰的车辆检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510121596.5A CN104732235B (zh) | 2015-03-19 | 2015-03-19 | 一种消除夜间道路反光干扰的车辆检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104732235A CN104732235A (zh) | 2015-06-24 |
CN104732235B true CN104732235B (zh) | 2017-10-31 |
Family
ID=53456109
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510121596.5A Active CN104732235B (zh) | 2015-03-19 | 2015-03-19 | 一种消除夜间道路反光干扰的车辆检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104732235B (zh) |
Families Citing this family (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105678240B (zh) * | 2015-12-30 | 2019-01-18 | 哈尔滨工业大学 | 一种用于道路去反光的图像处理方法 |
CN105893970A (zh) * | 2016-03-31 | 2016-08-24 | 杭州电子科技大学 | 基于亮度方差特征的夜间道路车辆检测方法 |
CN106651939B (zh) * | 2016-09-14 | 2019-07-09 | 天津工业大学 | 基于大气散射原理构建模型检测夜间交通视频多景深车灯的方法 |
CN106407951B (zh) * | 2016-09-30 | 2019-08-16 | 西安理工大学 | 一种基于单目视觉的夜间前方车辆检测方法 |
CN107067401A (zh) * | 2016-12-01 | 2017-08-18 | 浙江安谐智能科技有限公司 | 夜间行驶车辆的前照大灯精细分割方法 |
CN111368819B (zh) * | 2017-03-21 | 2023-05-30 | 北京中安未来科技有限公司 | 光斑检测方法和装置 |
CN108694699A (zh) * | 2017-04-10 | 2018-10-23 | 富士通株式会社 | 图像处理装置和方法 |
JP7114965B2 (ja) * | 2017-07-14 | 2022-08-09 | 富士通株式会社 | ターゲット検出方法、装置及び画像処理装置 |
CN107622494B (zh) * | 2017-08-28 | 2020-06-23 | 浙江工业大学 | 面向交通视频的夜间车辆检测与跟踪方法 |
CN108109138B (zh) * | 2017-12-15 | 2021-08-06 | 华南理工大学 | 一种针对类镜面物体的高光区域自适应匀光的方法 |
CN108052921B (zh) * | 2017-12-27 | 2020-12-04 | 海信集团有限公司 | 一种车道线检测方法、装置及终端 |
CN108538052A (zh) * | 2018-03-05 | 2018-09-14 | 华南理工大学 | 基于车头灯轨迹跟踪和动态配对的夜间交通流量检测方法 |
CN109063636B (zh) * | 2018-07-30 | 2021-05-07 | 安徽慧视金瞳科技有限公司 | 一种交互式课堂教学系统自动强光预警方法 |
CN109614859B (zh) * | 2018-11-01 | 2021-01-12 | 清华大学 | 视觉定位特征提取和匹配方法及装置 |
CN109753881B (zh) * | 2018-12-13 | 2023-05-12 | 小黄狗环保科技有限公司 | 一种防止强光干扰饮料瓶回收机识别投递物的方法 |
CN110197153B (zh) * | 2019-05-30 | 2023-05-02 | 南京维狸家智能科技有限公司 | 一种户型图中墙体自动识别方法 |
CN111723650A (zh) * | 2020-05-09 | 2020-09-29 | 华南师范大学 | 夜间车辆检测方法、装置、设备及存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102567965A (zh) * | 2011-12-20 | 2012-07-11 | 四川大学 | 基于亮度投影的车灯区域分割方法 |
CN102610104A (zh) * | 2012-03-16 | 2012-07-25 | 南京航空航天大学 | 一种车载的前方车辆检测方法 |
CN103455820A (zh) * | 2013-07-09 | 2013-12-18 | 河海大学 | 基于机器视觉技术的车辆检测和跟踪方法及系统 |
CN103914701A (zh) * | 2014-03-20 | 2014-07-09 | 燕山大学 | 一种基于图像的夜间车辆检测方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TW201002073A (en) * | 2008-06-23 | 2010-01-01 | Huper Lab Co Ltd | Method of vehicle segmentation and counting for nighttime video frames |
US8331623B2 (en) * | 2008-12-23 | 2012-12-11 | National Chiao Tung University | Method for tracking and processing image |
-
2015
- 2015-03-19 CN CN201510121596.5A patent/CN104732235B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102567965A (zh) * | 2011-12-20 | 2012-07-11 | 四川大学 | 基于亮度投影的车灯区域分割方法 |
CN102610104A (zh) * | 2012-03-16 | 2012-07-25 | 南京航空航天大学 | 一种车载的前方车辆检测方法 |
CN103455820A (zh) * | 2013-07-09 | 2013-12-18 | 河海大学 | 基于机器视觉技术的车辆检测和跟踪方法及系统 |
CN103914701A (zh) * | 2014-03-20 | 2014-07-09 | 燕山大学 | 一种基于图像的夜间车辆检测方法 |
Non-Patent Citations (5)
Title |
---|
On-Road Vehicle Detection during Dusk and at Night;Thomas Schamm 等;《2010 IEEE Intelligent Vehicles Symposium(IV)》;20100624;第418-423页 * |
Tracking and Pairing Vehicle Headlight in Night Scenes;Wei Zhang 等;《IEEE TRANSACTIONS ON INTELLIGENT TRANSPORTATION SYSTEMS》;20120331;第13卷(第1期);第140-153页 * |
Vision-Based Traffic Data Collection Sensor for Automotive Applications;David F. Llorca 等;《sensors》;20100122;第10卷(第1期);第860-875页 * |
基于单目视觉的夜间前方车辆检测方法;郭君斌 等;《汽车工程》;20140525;第36卷(第5期);第573-579,585页 * |
智能交通中视频车辆检测技术研究;陈琳;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20130715(第07期);C034-556 * |
Also Published As
Publication number | Publication date |
---|---|
CN104732235A (zh) | 2015-06-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104732235B (zh) | 一种消除夜间道路反光干扰的车辆检测方法 | |
US10127448B2 (en) | Method and system for dismount detection in low-resolution UAV imagery | |
CN102509098B (zh) | 一种鱼眼图像车辆识别方法 | |
CN106778551B (zh) | 一种高速路段及城市道路车道线识别方法 | |
CN109918971B (zh) | 监控视频中人数检测方法及装置 | |
CN107301405A (zh) | 自然场景下的交通标志检测方法 | |
CN106682586A (zh) | 一种复杂光照条件下基于视觉的实时车道线检测的方法 | |
CN109101924A (zh) | 一种基于机器学习的路面交通标志识别方法 | |
CN108256521B (zh) | 用于车身颜色识别的有效区域定位方法 | |
CN107092876A (zh) | 基于Retinex与S‑SIFT特征结合的低照度车型识别方法 | |
CN107133610B (zh) | 一种复杂路况下行车流量视觉检测与计数方法 | |
CN105893970A (zh) | 基于亮度方差特征的夜间道路车辆检测方法 | |
CN102842037A (zh) | 一种基于多特征融合的车辆阴影消除方法 | |
CN110060221B (zh) | 一种基于无人机航拍图像的桥梁车辆检测方法 | |
CN109886168B (zh) | 一种基于层阶的地面交通标志识别方法 | |
CN107563301A (zh) | 基于图像处理技术的红灯信号检测方法 | |
Cai et al. | Real-time arrow traffic light recognition system for intelligent vehicle | |
KR101026778B1 (ko) | 차량 영상 검지 장치 | |
CN106407951A (zh) | 一种基于单目视觉的夜间前方车辆检测方法 | |
CN111507391A (zh) | 一种基于机器视觉的有色金属破碎料智能识别方法 | |
Wali et al. | Shape matching and color segmentation based traffic sign detection system | |
CN111476804A (zh) | 托辊图像高效分割方法、装置、设备及存储介质 | |
CN109800693B (zh) | 一种基于颜色通道混合特征的夜间车辆检测方法 | |
CN113989771A (zh) | 一种基于数字图像处理的交通信号灯识别方法 | |
CN106529533A (zh) | 一种基于多尺度分析和匹配排序的复杂天气车牌定位方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
EE01 | Entry into force of recordation of patent licensing contract | ||
EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20150624 Assignee: HANGZHOU ZHISHU TECHNOLOGY CO.,LTD. Assignor: HANGZHOU DIANZI University Contract record no.: X2022330000062 Denomination of invention: A vehicle detection method for eliminating night road reflective interference Granted publication date: 20171031 License type: Common License Record date: 20220331 |