一种基于似物性估计的道线检测方法与流程

本发明涉及图像处理领域，特别涉及一种基于似物性估计的道线检测方法。

背景技术：

随着传感器技术和电子技术的发展，先进驾驶辅助系统(adas)成为汽车工业发展的一个重要方向。在adas系统中，道线检测是车道偏离预警、车道跟随、前车防撞预警等应用的重要基础。

目前已有的道线检测方法可以粗略地分成两大类：基于低层次图像特征的方法和基于模型匹配的方法。在使用低层次特征的方法中，中国专利cn100424724c、cn103617412b、cn105426868a、cn105046198a、cn104573618a以边缘点检测为算法基础；在基于模型匹配的方法中，中国专利cn106228125a以学习的方法建立包含道线的图像块所对应的模型。

除了道线标记之外，adas系统的摄像机所观察的视域中经常出现如各式车辆、车道分隔护栏、路旁的绿化等物体。无论是基于低层次特征的方法还是基于模型匹配的方法或多或少都会受到这些对象所引入的虚假信息的困扰。例如，基于边缘点检测的方法通常不加区分地将车辆、路旁植物的边缘点与真实道线的边缘点一起输入到hough变换的投票模块；对基于模型匹配的方法而言，由于图像中存在与道线结构类似的区域，增加检测率与降低虚检率始终是一个两难的问题。另外，由车辆等非道线对象引入的虚假信息不可避免地会增加adas系统花费在道线检测中的计算代价。

技术实现要素：

本发明实施例提供一种基于似物性估计的道线检测方法，所提供的技术以车载摄像机采集的图像序列为输入，通过似物性估计自动地确定候选道线区域中每个像素的投票权值，以hough投票的方式检测道线在图像中的位置。

一种基于似物性估计的道线检测方法，如下：

以块扫描的形式扫描输入图像，在每个扫描位置计算块内图像数据的似物性值c；

以如下的分段函数对似物性值作映射处理，增强较大的似物性值，抑制较小的值

其中η是一个实数，且0＜η≤0.5，t1和t2是两个预设的阈值，且t1＞t2；

以处理后的似物性值cm和像素的梯度强度的乘积作为后续hough投票的权值；

在图像平面预设两条直线，第一条直线lver平行于图像平面的y轴且经过前一帧图像检测所得的消失点(x0，y0)，第二条直线lhor平行于图像平面的x轴且接近图像的底部，其y坐标为h-δh，其中h为图像的高度，δh是一个预设的常数；

假设所检测的直线与所述第一条直线lver存在交点(x0，ym)，且交点的y坐标ym的取值范围为y0-δy≤ym≤y0+δy，其中δy为一个预设的常数；若点(x，y)有投票权值，则由点(x0，ym)和(x，y)可确定一条直线，该直线与所述第二条直线lhor的交点可按下式计算：

对ym取值范围内的所有整数，由上述公式计算对应的xb，对每一个计算所得的(xb，ym)，在参数平面的相应位置累加投票权值；

对所有具有大于阈值的投票权值的点实施上述投票过程；

在xb-ym参数平面中搜索极大值点，由极大值点的ym和xb坐标对应的直线确定两侧的道线。

附图说明

图1示出了方框标示的四个以图像块表示的局部区域，其中的区域①和②包含了部分道线，③和④没有包含道线；

图2为对图1的4个图像块计算方向梯度分布直方图所得的结果；

图3示意了在前一帧图像计算所得的消失点设置一条垂直线，在接近图像底部设置一条水平线，由这两条直线上的交点可确定道线所在的直线；

图4为本发明实施例基于似物性估计道线检测方法的流程示意图。

具体实施方式

本发明提供一种基于似物性估计的道线检测方法。所提供的方法对输入图像以基于似物性估计的算法自动地为候选区域内的边缘点确定进行后续hough投票的权值，由hough投票检测得到的直线确定道线。

adas系统的摄像机所采集的道路图像经常包含车辆、路面破损和坑洼、路旁树木和建筑物的阴影等引入的边缘。如果不加区别地让这些区域的边缘信息参与hough投票，则势必会引入虚假的线段，对道线检测的结果产生不良影响。

观察包含道线的道路图像，不难发现在道线区域，不仅其边缘两侧具有明显的梯度变化，而且在一条道线同侧的边缘像素具有接近的梯度方向。图1示出了一个示意图，图中示出了4个图像块，其中的块①包含了道线的左侧部分，块②则在水平方向跨越了道线区域，块③位于路面具有明显纹理的区域，块④则包含车身的边缘部分；图2是对这4个图像块计算方向梯度分布直方图所得的结果，其中的梯度方向以公式(2)计算的弧度值转换成角度表示，且将[-90，90]范围内的整数均匀量化成18个量化等级，对具有相应区间范围梯度方向的像素的梯度强度值进行累加计算形成图2所示的直方图。由图中例子容易发现，与位于非道线区域的块相比，①块和②块的直方图的分布更为集中。

基于上述观察，本发明提供了一种以如下方法计算的针对道线检测的似物性度量：首先对输入图像计算x和y方向的梯度图像gx和gy，按公式(1)和公式(2)计算像素(x，y)的梯度强度和方向：

α(x，y)＝tan^-1(gy(x，y)/gx(x，y))(2)

其次，按公式(2)计算的值为[-π/2，π/2]范围内实数表示的弧度，将其换算成对应的[-90，90]范围内的角度值，且作均匀量化；再次，统计预设大小的块的方向梯度直方图h，其横坐标表示角度经均匀量化后的量化等级，纵坐标是块内具有该等级对应角度范围的所有像素的梯度强度之和；最后，按以下公式(3)计算以(x，y)为中心，大小为m×m的图像块的似物性度量：

其中hs代表对直方图h按由大到小排序后形成的数组，hs(1)和hs(2)是排序后数组中具有最大和次大值的元素。

注意到按上述公式(3)计算的似物性值更多地取决于区域是否包含道线的边缘部分，而对包含的区域面积不敏感，所以本发明的一个实施例对较大分辨率的图像首先进行缩小处理。具体地，对于分辨率不小于1280×720的图像，x和y方向各缩小到其原始尺寸的二分之一；对于不小于1920×1080的图像，缩小到其原始尺寸的四分之一。

以块扫描的形式处理输入图像，若某块具有较大的c值，则该块有较大的可能性位于道线边缘区域；反之，较小的c值则对应较小的可能性。本发明的一个实施例进一步地按公式(4)的分段函数增强较大的似物性值，抑制较小的似物性值，且将处理以后的似物性值赋值给块内的所有像素：

其中η是一个实数，且0＜η≤0.5，t1和t2是两个预设的阈值，本发明的一个实施例取t1＝0.7，t2＝0.4。

对图像的每个像素，以该像素的梯度值乘以cm值作为该像素的投票权值，从而得到权值图r，具体地，计算公式如公式(5)：

r(x，y)＝cm(x，y)×g(x，y)(5)

对于大量正常安装的驾驶辅助系统，地平线通常出现在图像的中心附近，平行的道线相交于消失点没入地平面而消失不见。参照图3，设前一帧图像检测所得的消失点为(x0，y0)，则当前帧图像中左右道线的交点应该在(x0，y0)附近的位置。本发明的实施例在图像平面中设定两条直线，第一条直线lver平行于图像平面的y轴且经过点(x0，y0)；第二条直线lhor平行于图像平面的x轴且接近图像的底部，其y坐标为h-δh。如果检测到了道线上的一点(x，y)，设道线与lver的交点e1的y坐标为ym，连接e1和p成一条直线，由公式(6)可以得到该直线与lhor的交点e2的x坐标值xb。

以xb和ym为参数，公式(6)表示了xb和ym之间的函数关系，图像平面中的任意一点(x，y)由公式(5)将映射为xb-ym平面中的一条曲线。进一步地，图像平面中位于一条直线上的点，由公式(5)将映射成xb-ym平面中的多条曲线，由于图像平面中的直线上的点有相同的xb和ym值，所以这些曲线将在xb-ym平面中交于一点。由此，若对于图像平面中具有投票权值的点，在xb-ym平面由该点确定的曲线经过的点上累加其权值，然后，通过检测xb-ym平面中的极大值点可以确定直线与lver的交点的y坐标ym，与lbot的交点的x坐标xb，从而确定图像平面中的一条直线。

根据上述分析，本发明提供了一种基于似物性估计的道线检测方法，为了更具体地描述本发明，以下结合附图以及具体实施方式对本发明的技术方案进行详细说明。

如图4所示，本发明基于似物性估计的道线检测方法的流程可以包括：

步骤401，以一维梯度算子[-1，0，+1]和[-1，0，+1]^t分别沿x和y两个方向与图像作卷积运算计算x和y方向的梯度图像gx和gy，按公式(1)计算每个像素(x，y)的梯度强度，按公式(2)计算梯度方向，且将计算所得以弧度表示的梯度方向转换成以角度表示。

步骤402，以δx和δy为x和y方向的步长扫描输入图像，在每个扫描位置计算以该位置为中心，m×m大小块的方向梯度分布直方图h，其横坐标表示按步骤401计算的角度经均匀量化后的量化等级，纵坐标是块内具有该等级对应角度范围的所有像素的梯度强度之和，其中，本发明的一个实施例预设量化等级数为18。

步骤403，根据块的方向梯度分布直方图计算块所在位置的似物性值，具体地，对直方图h，按其元素从大到小排序，按公式(3)计算块所在位置的似物性值。

步骤404，按公式(4)的分段函数增强较大的似物性值，抑制较小的似物性值，且将映射处理后的似物性值记作cm。

对图像进行块扫描的时候，相邻块可以存在重叠区域，对位于重叠区域的像素，取所有覆盖所述像素的块计算所得的cm值的最大值作为该像素的cm值。

步骤405，对图像的每个像素，以该像素的梯度强度值乘以cm值作为该像素的投票权值，从而得到权值图r，具体地，计算公式如公式(5)。

步骤406，以hough变换的方式将图像平面中每个投票权值大于预设阈值的点(x，y)以公式(6)映射成xb-ym参数平面的一条曲线，且在曲线经过的点累加权值r(x，y)。

具体地，设定参数平面的一个坐标轴对应参数xb，另一个坐标轴对应参数ym。本发明的一个实施例使用一个二维矩阵q记录参数平面中的累加权值，且初始化所有的矩阵元素为0。签于在辅助驾驶系统中，容易保证图像中消失点位置不会在较小的时间间隔内有很大的变化，可以通过减少ym的搜索范围来减小计算代价，本发明的一个实施例设定ym的取值范围为y0-δy≤ym≤y0+δy，其中(x0，y0)为前一帧图像计算所得的消失点，δy是根据图像大小预设的常数，对于不大于1920×1080的图像，本发明的一个实施例取δy＝32；对权值图r的每个权值大于预设阈值rth的位置，由图像平面坐标(x，y)，依次取所述ym取值范围内的所有整数，根据公式(6)计算相应的xb，在矩阵q的(xb，ym)位置按公式(7)累加权值：

q(xb，ym)＝q(xb，ym)+r(x，y)(7)

步骤407，按公式(8)以窗口扫描的方式在矩阵q中搜索满足条件的极大值位置，

q(m，b)＝max{q(m+u，b+v)，-bth≤u，v≤bth}且q(m，b)≥qth(8)

其中(m，b)表示矩阵中的一个搜索位置，u和v为整数，表示两个方向与搜索位置的偏离，bth是表示窗口大小的预设值，本发明的一个实施例取bth＝16，qth是一个预设的阈值，本发明的一个实施例取qth＝64。对于满足公式(8)的极大值位置，将其坐标与q值构成一个元组(m，b，q(m，b))，加入集合z。

步骤408，若z中的元素数目不多于一个预设的阈值nth，则保留所有的极大值位置；否则对z中的元素按其q值从大到小排序，只保留集合中较大的nth个元素进行后续运算。

步骤409，集合z中的每一个坐标可确定一条直线，具体地，若有坐标(m0，b0)，则图像平面中对应的直线由两点(x0，m0)和(b0，h-δh)确定，其中x0为所述前一帧检测所得消失点的x方向坐标，h为图像的高度，δh是一个预设的常数。

步骤410，将前一个步骤确定的直线按其斜率是否为负数划分成两个集合，在斜率为负的直线集合中，选择累积了最大权值的一条直线作为左道线；在斜率为非负的直线集合中，选择累积了最大权值的一条直线作为右道线。

步骤411，更新消失点，具体地，设当前帧计算所得的左右道线的交点为(xc，yc)，前一帧的消失点为(xp，yp)，则按公式(9)更新消失点为(x0，y0)：

其中q是一个预设的常数，本发明的一个实施例取0.7≤q≤0.9。

以上所述仅为本发明的较佳实施例，但本发明的保护范围并不局限于此，凡在本发明的精神和原则之内，所做的任何修改或替换等，都应涵盖在本发明的保护范围内。