基于轨迹特征挖掘的追尾风险路段动态识别与管控方法

1.本发明属于道路驾驶安全监管技术领域，具体涉及一种基于轨迹特征挖掘的追尾风险路段动态识别与管控方法。


背景技术：

2.根据美国公路交通安全局(nhtsa)公布的数据显示，高速公路上车辆之间产生的碰撞事故中，追尾事故占比达到42.61％，且死亡率高达42.83％。事故中，车辆在静止停放时如果被追尾，造成人员伤亡的后果更为严重。
3.现有技术应对车辆追尾风险主要以加固车身(如采用防撞车身)、车载雷达主动预警等思路展开，缺少全局视野，导致头痛医头脚痛医脚，存在盲目性及时空局限性。本方案从全局视野出发，借助轨迹大数据研判追尾风险行为及特征路段，为追尾风险防范提供了新的管控方法和技术框架。


技术实现要素：

4.针对现有技术中的上述不足，本发明提供的基于轨迹特征挖掘的追尾风险路段动态识别与管控方法解决了传统管控方法仅着眼于车辆本身追尾行为导致管控效率低、车辆传感设备依赖性高的问题。
5.为了达到上述发明目的，本发明采用的技术方案为：基于轨迹特征挖掘的追尾风险路段动态识别与管控方法，包括以下步骤：
6.s1、采集目标区域的路网数据，并进行追尾风险路段语义划分；
7.s2、采集目标区域的动态/历史车辆轨迹数据，进行轨迹出行路径匹配，得到追尾风险对应的分向路段；
8.s3、根据目标区域的动态/历史车辆轨迹数据，分别计算匹配分向路段的疲劳指数、流量指数以及急减速指数；
9.s4、根据计算的疲劳指数、流量指数以及急减速指数，计算追尾风险路段的加权追尾风险指数；
10.s5、根据加权追尾风险指数，分级动态识别追尾风险路段的追尾风险等级，并根据风险等级进行分类管控。
11.进一步地，所述步骤s1具体为：
12.s11、采集目标区域的路网数据，并修复异常路网数据，确定路网信息图；
13.s12、在路网信息图中，将100米以内的所有交汇口端点的形心识别为逻辑道路交叉口；
14.s13、对相邻两个逻辑道路交叉口之间的路段划分逻辑分向路段，实现追尾风险路段的语义划分。
15.进一步地，所述步骤s2具体为：
16.s21、采集目标区域的动态/历史车辆轨迹数据，对其处理并分别构建各车辆的历
史出行轨迹集合和动态出行轨迹集合；
17.s22、遍历各车辆的历史/动态出行轨迹集合中的轨迹点数据；
18.s23、对历史/动态出行轨迹集合中的轨迹点，获取其临近范围内的临近连接线集合；
19.s24、对临近连接线集合不为空的轨迹点，遍历其对应的临近连接线集合，并确定该轨迹点的方位角射线与各临近连接线切线的夹角；
20.s25、根据夹角大小，匹配分向路段：
21.当存在小于45度的夹角时，将最小夹角的临近连接线作为匹配连接线，其对应的虚拟分向路段作为匹配的分向路段；
22.当所有夹角均大于45度时，则若临近连接线最大夹角大于135度，则以该连接线对应分向路段的反向路段作为匹配的分向路段。
23.进一步地，所述步骤s3中，计算疲劳指数的方法具体为：
24.sa3-1、根据目标区域的路网数据，进行疲劳驾驶路段语义划分；
25.sa3-2、根据目标区域的动态/历史车辆轨迹数据，进行轨迹出行路径匹配，得到疲劳驾驶对应的分向路段；
26.sa3-3、根据疲劳驾驶的分向路段和追尾风险的分向路段的对应关系，对于每个追尾风险的分向路段，在疲劳驾驶的分向路段中划分与追尾风险对应的路段；
27.sa3-4、计算在疲劳驾驶的分向路段中与追尾风险对应的路段的疲劳驾驶指数，作为该追尾风险分向路段的疲劳指数；其中，疲劳指数包括相对疲劳指数和绝对疲劳指数。
28.进一步地，所述步骤s3中，计算流量指数的方法具体为：
29.sb3-1、获取历史轨迹出行集合中各轨迹点对应的匹配连接线集合；
30.sb3-2、依次确定各轨迹点处的车速，对于车速大于0的轨迹点，参考其匹配连接线，统计对应分向路段各时段各车道流量；
31.sb3-3、将分向路段内动态时段的平均车道流量作为该分向路段的流量指数。
32.进一步地，所述步骤s3中，计算急减速指数的方法为：
33.sc3-1、遍历各车辆的历史轨迹点出行集合中的轨迹点；
34.sc3-2、根据前向连续三个轨迹点之间的速度关系，识别减速事件；
35.sc3-3、计算减速事件中中间轨迹点的减速速度；
36.sc3-4、将各分向路段中所属连接线中轨迹点减速速度大于预设值的次数与减速事件总次数的比值，作为该分向路段的急减速指数。
37.进一步地，所述步骤sc3-2中，当轨迹点p
i+1
的速度v
i+1
＝0，且0《vi《v
i-1
时，识别为出现减速事件；其中，下标i为轨迹点序号；
38.所述步骤sc3-3中，中间轨迹点pi的减速速度di计算公式为：
[0039][0040]
式中，difv为轨迹点pi和p
i-1
的速度差，refdist1为第一参考距离，refdist2为第二参考距离；
[0041]
其中，refdist1＝v
i-1
*diftime，diftime为轨迹点pi和p
i-1
间的时间间隔，
difdist为轨迹点pi和p
i-1
的欧氏距离，diffh为(refdist1-refdist)》1米时，轨迹点p
i-1
和pi的高差。
[0042]
进一步地，所述步骤s4中，追尾风险路段的加权追尾风险指数reindi的计算公式为：
[0043]
reindi＝w1*w2*w3*dindi
[0044]
式中，w1为相对疲劳指数对应的相对疲劳风险等级权重，w2为绝对疲劳指数对应的绝对疲劳风险等级权重，w3为流量指数权重，dindi为急减速指数；
[0045]
其中，相对疲劳指数对应的风险等级根据分向路段的相对疲劳指数进行1～10级的风险等级划分，绝对疲劳风险等级根据分向路段的绝对疲劳风险指数，按自然间断法分10级进行间断排序得到，流量指数权重根据动态时段平均车道流量所属范围确定。
[0046]
进一步地，所述步骤s5中，动态识别追尾风险等级方法为：
[0047]
根据当前时刻，筛选上一时段车辆采样数大于50次的分向路段，根据其对应的加权追尾风险指数，按自然间断法分10级进行间断排序，获得各分向路段的追尾风险等级。
[0048]
进一步地，所述步骤s5中，进行分类管控的方法为：
[0049]
对追尾风险等级小于4的分向路段，不进行管控；
[0050]
对追尾风险等级为5～7的分向路段，联动电子导航地图进行提醒；
[0051]
对追尾风险等级为8～10的分向路段，对分向路段间隔1公里设置提醒标识，同时联动电子导航地图进行提醒；
[0052]
对加权追尾风险指数为目标区域内前100的分向路段，对分向路段间隔1公里设置提醒标识、联动电子导航地图进行提醒并设置临时休息区。
[0053]
本发明的有益效果为：
[0054]
(1)本发明实现了基于车辆实时轨迹数据分析的追尾风险路段动态识别管控，包括基于道路结构特征的车辆追尾风险分析路段语义分割、基于实时轨迹特征分析的加权追尾风险指数统计(疲劳指数+流量指数+急减速指数)、追尾风险路段分类动态识别，风险路段分类管控；
[0055]
(2)本发明方法突破了传统管控方法仅着眼于车辆本身追尾行为导致管控效率低、车辆传感设备依赖性高的局限；
[0056]
(3)本发明方法可显著提高道路网络整体追尾风险路段动态识别管控水平，降低相关事故风险。
附图说明
[0057]
图1为本发明提供的追尾风险路段识别与管控方法流程图。
[0058]
图2为本发明提供的路段语义划分前示意图。
[0059]
图3为本发明提供的路段语义划分后示意图。
[0060]
图4为本发明提供的异常轨迹点数据剔除示意图。
[0061]
图5为本发明提供的轨迹点匹配临近连接线示意图。
[0062]
图6为本发明提供的轨迹点减速速度估计示意图。
具体实施方式
[0063]
下面对本发明的具体实施方式进行描述，以便于本技术领域的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。
[0064]
本发明实施例提供了一种基于轨迹特征挖掘的追尾风险路段动态识别与管控方法，如图1所示，包括以下步骤：
[0065]
s1、采集目标区域的路网数据，并进行追尾风险路段语义划分；
[0066]
s2、采集目标区域的动态/历史车辆轨迹数据，进行轨迹出行路径匹配，得到追尾风险对应的分向路段；
[0067]
s3、根据目标区域的动态/历史车辆轨迹数据，分别计算匹配分向路段的疲劳指数、流量指数以及急减速指数；
[0068]
s4、根据计算的疲劳指数、流量指数以及急减速指数，计算追尾风险路段的加权追尾风险指数；
[0069]
s5、根据加权追尾风险指数，分级动态识别追尾风险路段的追尾风险等级，并根据风险等级进行分类管控。
[0070]
在本发明实施例的步骤s1中，原始地理信息gis路网数据由于具有不规则的底层数据结构，其中道路连接线数据无法直接作为研究统计分析对象，需首先根据连接线数据特征进行语义划分，划分后的分向路段逻辑单元作为研究统计对象单元。基于此，本发明实施例的步骤s1具体为：
[0071]
s11、采集目标区域的路网数据，并修复异常路网数据，确定路网信息图；
[0072]
s12、在路网信息图中，将100米以内的所有交汇口端点的形心识别为逻辑道路交叉口；
[0073]
s13、对相邻两个逻辑道路交叉口之间的路段划分逻辑分向路段，实现追尾风险路段的语义划分。
[0074]
本实施例的步骤s11中，获取目标区域内的路网gis底层数据，并人工修复异常数据，如悬挂线和数据逻辑错误的。
[0075]
本实施例的步骤s12中，如图2所示，定义路网信息图中两条连接线断点距离小于1米为相交状态，对于三条连接线及以上具有相交状态的端点，识别为交汇点，即图2中的gis连线端点a～h；对于100米以内的所有交汇口端点的形心，识别为逻辑道路交叉口，即图2中的虚拟点ab。
[0076]
在本实施例的步骤s13中，以图2为例，对于识别出的临近两个逻辑道路交叉口a,b之间的连线进行语义路段划分，划分方法具体为：
[0077]
s13-1、连接相邻两个逻辑道路交叉口，获得射线ab；
[0078]
s13-2、根据射线ab的长度，划分逻辑分向路段，并进行路段编号；
[0079]
当射线ab小于400米时，将射线ab作为逻辑分向路段，并进行路段编号；
[0080]
当射线ab大于400米时，以逻辑道路交叉口a为起点，200米为间隔对射线ab依次进行逻辑分向路段划分，直到最后一个划分点c至逻辑道路交叉口b的距离小于400米时，以射线cb中点进行切割划分得到对应两段逻辑分向路段，对每个划分出的逻辑分向路段，依次
进行路段编号；
[0081]
其中，进行路段编号的方法为：
[0082]
将逻辑道路交叉口a到b方向的逻辑分向路段依次编号为1,2,3,
…
,s，将逻辑道路交叉口b到a方向的逻辑分向路段依次编号为s+1,s+2,s+3,
…
,s+s；s为划分出的逻辑分向路段数量，s为正整数。
[0083]
划分出的逻辑分向路段如图3所示，具体地：
[0084]
(1)当射线ab小于400米时，则对应与ab同向的连接线如cd、dg、gh，赋同一分向路段编号s，与ab反向的连接线如fe、eb、ba，赋另一分向路段编号s+1；
[0085]
(2)当射线ab大于400米时，从a出发，按1公里为单位打断ab，直至剩余段落距离不足400米时，在中点打断，如图3所示，划分对应逻辑分向路段并编号。
[0086]
本发明实施例的步骤s2具体为：
[0087]
s21、采集目标区域的动态/历史车辆轨迹数据，对其处理并分别构建各车辆的历史出行轨迹集合和动态出行轨迹集合；
[0088]
s22、遍历各车辆的历史/动态出行轨迹集合中的轨迹点数据；
[0089]
s23、对历史/动态出行轨迹集合中的轨迹点，获取其临近范围内的临近连接线集合；
[0090]
s24、对临近连接线集合不为空的轨迹点，遍历其对应的临近连接线集合，并确定该轨迹点的方位角射线与各临近连接线切线的夹角；
[0091]
s25、根据夹角大小，匹配分向路段：
[0092]
当存在小于45度的夹角时，将最小夹角的临近连接线作为匹配连接线，其对应的虚拟分向路段作为匹配的分向路段；
[0093]
当所有夹角均大于45度时，则若临近连接线最大夹角大于135度，则以该连接线对应分向路段的反向路段作为匹配的分向路段。
[0094]
本实施例的步骤s21中，构建历史出行轨迹集合的方法为：
[0095]
s21-a1、采集目标区域的历史车辆轨迹数据，并按车辆对轨迹点数据分组；
[0096]
s21-a2、对各车辆轨迹点数据按采样时间顺序排列，并剔除异常轨迹点数据，获得轨迹序列；
[0097]
s21-a3、识别轨迹序列中的出行中断点，并构建中断点集合p；
[0098]
s21-a4、根据构建的中断点集合p，识别车辆出行段落，进而获得各车辆的出行轨迹集合；
[0099]
在步骤s21-a1中，通过采样车辆卫星定位、国际货运监管平台、网络货运平台等数据源获取货车轨迹数据，其中获取的轨迹点数据格式如表1所示；
[0100]
表1轨迹点字段属性
[0101][0102]
在步骤s21-a2中，剔除异常轨迹点数据的方法为：
[0103]
如图4所示，按轨迹点数据的时间顺序遍历起始轨迹点到到达轨迹点之间的所有轨迹点p(i)，若p(i)到p(i-1)与p(i+1)距离均超过阈值1公里，同时p(i-1)到p(i+1)距离d(i-1,i+1)小于1公里，则将p(i)视为离群的异常轨迹点数据并剔除。
[0104]
在步骤s21-a3中，将对前向车速为0的轨迹点，若累积停止时间超过30分钟的轨迹点识别为出行中断点，进而获得中断点集合。
[0105]
在步骤s21-a4中，将每个出行中断点后第一个速度不为0的轨迹点视为起点，到下一个出行中断点间连续过个轨迹点，识别为一个出行段落，进而获得各车辆的出行轨迹集合。
[0106]
本实施例的步骤s21中，构建动态出行轨迹集合的方法为：
[0107]
s21-b1、设置时段间隔，按当前时刻汇总上一时段间隔内的轨迹点数据；
[0108]
s21-b2、按车辆对汇总时段的轨迹点数据分组；
[0109]
s21-b3、对各车辆的轨迹点数据集合进行遍历，并剔除异常轨迹点数据，构建出动态出行轨迹集合。
[0110]
在步骤s21-b1中，通过采样车辆卫星定位、国际货运监管平台、网络货运平台等数据源获取车辆实时轨迹数据，其中获取的轨迹点数据格式同表1所示；以10分钟时段间隔，按当前时刻上一时段汇总上传轨迹点，如当前时刻为早8点11分，则所处时段为早8点10分-8点20分时段，上一时段即为早8点-8点10分时段。
[0111]
在步骤s21-b3中，剔除异常轨迹点数据的方法同步骤s21-a2。
[0112]
在步骤s23中，临近连接线是指地图gis数据中的最底层连接线单元，每个连接线由一系列道路节点连接组成，有对应道路属性、等级、反向(单向或双向)、起点编号、终点编号以及前面确定的分向路段编号；对于每一个轨迹点，获取临近范围内临近连接线集合，首选搜索30米范围内临近连接线集合，若为空集则进一步扩大搜索范围为40米，依次类推，直至60米，若60米内仍无临近连接线，则输出该记录并进行路网检查。
[0113]
在步骤s25中，虚拟分向路段是指通过逻辑道路交叉口a指向b的射线划分的逻辑路段；在匹配分向路段的过程中，当所有角均大于45度但小于135度时，表示连接线匹配失
败，即不认为轨迹点从该连接线经过。
[0114]
在步骤s25中，对轨迹点匹配连接线如图5所示，其中，β1和β2为均为45度以下的夹角。
[0115]
在本发明实施例的步骤s3中，计算疲劳指数的方法具体为：
[0116]
sa3-1、根据目标区域的路网数据，进行疲劳驾驶路段语义划分；
[0117]
sa3-2、根据目标区域的动态/历史车辆轨迹数据，进行轨迹出行路径匹配，得到疲劳驾驶对应的分向路段；
[0118]
sa3-3、根据疲劳驾驶的分向路段和追尾风险的分向路段的对应关系，对于每个追尾风险的分向路段，在疲劳驾驶的分向路段中划分与追尾风险对应的路段；
[0119]
sa3-4、计算在疲劳驾驶的分向路段中与追尾风险对应的路段的疲劳驾驶指数，作为该追尾风险分向路段的疲劳指数；其中，疲劳指数包括相对疲劳指数和绝对疲劳指数。
[0120]
在上述步骤sa3-1中，进行疲劳路段的语义划分方法与进行追尾风险路段语义划分的方法类似，在具体划分逻辑道分向路段时，以图2为例，其划分方法为：
[0121]
对于相邻两个逻辑道路交叉口确定的射线ab，根据射线ab的长度，划分逻辑分向路段，并进行路段编号；
[0122]
当射线ab小于1公里时，将射线ab作为逻辑分向路段，并进行路段编号；
[0123]
当射线ab大于1公里但小于2公里时，以射线ab中点进行切割划分，获得两段逻辑分向路段，并进行路段编号；
[0124]
当射线ab大于2公里时，以逻辑道路交叉口a为起点，1公里为间隔对射线ab依次进行逻辑分向路段划分，直到最后一个划分点c至逻辑道路交叉口b的距离小于2公里时，以射线cb中点进行切割划分得到对应两段逻辑分向路段，对每个划分出的逻辑分向路段，依次进行路段编号。其中，路段编号方式与上述划分追尾风险路段时的编号方法相同，在此不再赘述。
[0125]
在上述步骤sa3-2中，进行轨迹出行路径匹配得到疲劳驾驶对应分向路段的方法相同，在此不再赘述。
[0126]
在上述步骤sa3-3中，对各追尾风险分向路段所属连接线，分析其在疲劳驾驶分向路段中的对应划分情况，若全部划分至同一分向路段，则该分向路段为对应的路段，若划分至多于1个疲劳驾驶分向路段，则选择划分长度比例最大的疲劳驾驶分向路段作为对应的路段。
[0127]
在上述步骤sa3-4中，计算疲劳指数的方法为：
[0128]
sa3-41、对于确定的对应路段，在设定时段内，统计设定时间间隔内的出行次数，以及平均连续驾驶时长；
[0129]
其中，在统计出行次数时，属于同一辆车的一个或多个轨迹点在该时段内匹配连接线对应分向路段出行次数+1；所有车辆轨迹遍历后，即可统计出个分向路段对应各时段经过车辆数及平均连续驾驶时长；
[0130]
sa3-42、根据统计的出行次数及平均连续驾驶时长，计算各分向路段的疲劳指数；
[0131]
其中，疲劳指数包括相对疲劳指数和绝对疲劳指数；
[0132]
其中，对于任一时段内，路段连续驾驶t小时的相对疲劳指数ri
t
为：
[0133][0134]
式中，λ为预设的时段权重，n为连续驾驶时间大于t小时的出行次数，n为出行总次数；
[0135]
对于任一时段内，路段连续驾驶t小时的绝对疲劳指数ai
t
为：
[0136]
ai
t
＝λ∑icnti×
wi[0137]
式中，cnti为连续驾驶时段i的出行采样次数，wi为cnti对应的权重。
[0138]
具体地，对于每一路段，按20分钟时间间隔统计该间隔内的出行次数，若同一出行轨迹在同一路段上有连续多个轨迹点，则以最后轨迹点对应驾驶时间为准，不重复统计；并统计预设时段内，每一时间间隔下的平均连续驾驶时长，统计结果如表2所示；
[0139]
表2连续驾驶时间统计
[0140][0141]
在计算相对疲劳指数时，时段1权重为0.5，时段2权重为1.5，时段3权重为1.0，时段4权重为3.0。
[0142]
在计算绝对疲劳指数时，wi设置如表3所示；
[0143]
表3连续驾驶时间权重
[0144]
连续驾驶时段权重2小时20分1
2小时40分1.23小时1.4......9小时40分5.410小时5.6大于10小时10
[0145]
在本发明实施例的步骤s3中，计算流量指数的方法具体为：
[0146]
sb3-1、获取历史轨迹出行集合中各轨迹点对应的匹配连接线集合；
[0147]
sb3-2、依次确定各轨迹点处的车速，对于车速大于0的轨迹点，参考其匹配连接线，统计对应分向路段各时段各车道流量；
[0148]
具体地，对于同一车辆同一追尾风险分向路段仅统计一次，不重复统计；
[0149]
sb3-3、将分向路段内动态时段的平均车道流量作为该分向路段的流量指数。
[0150]
在本发明实施例中步骤s3中，计算急减速指数的方法为：
[0151]
sc3-1、遍历各车辆的历史轨迹点出行集合中的轨迹点；
[0152]
sc3-2、根据前向连续三个轨迹点之间的速度关系，识别减速事件；
[0153]
sc3-3、计算减速事件中中间轨迹点的减速速度；
[0154]
sc3-4、将各分向路段中所属连接线中轨迹点减速速度大于预设值的次数与减速事件总次数的比值，作为该分向路段的急减速指数。
[0155]
在上述步骤sc3-2中，当轨迹点p
i+1
的速度v
i+1
＝0，且0《vi《v
i-1
时，识别为出现减速事件；其中，下标i为轨迹点序号；
[0156]
在上述步骤sc3-3中，如图6所示，中间轨迹点pi的减速速度di计算公式为：
[0157][0158]
式中，difv为轨迹点pi和p
i-1
的速度差，refdist1为第一参考距离，refdist2为第二参考距离；
[0159]
其中，refdist1＝v
i-1
*diftime，diftime为轨迹点pi和p
i-1
间的时间间隔，difdist为轨迹点pi和p
i-1
的欧氏距离，diffh为(refdist1-refdist)》1米时，轨迹点p
i-1
和pi的高差。
[0160]
在上述步骤sc3-4中，减速速度预设值为3.0m/s2。
[0161]
在本发明实施例的步骤s4中，追尾风险路段的加权追尾风险指数reindi的计算公式为：
[0162]
reindi＝w1*w2*w3*dindi
[0163]
式中，w1为相对疲劳指数对应的相对疲劳风险等级权重，w2为绝对疲劳指数对应的绝对疲劳风险等级权重，w3为流量指数权重，dindi为急减速指数；
[0164]
其中，相对疲劳指数对应的风险等级根据分向路段的相对疲劳指数进行1～10级的风险等级划分，绝对疲劳风险等级根据分向路段的绝对疲劳风险指数，按自然间断法分10级进行间断排序得到，流量指数权重根据动态时段平均车道流量所属范围确定。
[0165]
具体地，疲劳风险等级权重参考如表4所示；
[0166]
表4疲劳风险等级权重参考
[0167]
相对疲劳风险等级w1绝对疲劳风险等级w20级0.50级0.51级0.61级0.72级0.72级0.93级0.83级1.14级0.94级1.35级15级1.56级1.16级1.77级1.27级1.98级1.38级2.19级1.49级2.310级1.510级2.5
[0168]
流量指数权重参考如表5所示；
[0169]
表5动态流量权重参考
[0170]
动态时段平均车道流量w3《1000.1100-1500.5150-2001200-2501.5》＝2503
[0171]
在本发明实施例的步骤s5中，动态识别追尾风险等级方法为：
[0172]
根据当前时刻，筛选上一时段车辆采样数大于50次的分向路段，根据其对应的加权追尾风险指数，按自然间断法分10级进行间断排序，获得各分向路段的追尾风险等级。
[0173]
基于上述追尾风险等级的划分方式，对各分向路段进行分类管控的方法为：
[0174]
对追尾风险等级小于4的分向路段，不进行管控；
[0175]
对追尾风险等级为5～7的分向路段，联动电子导航地图进行提醒；
[0176]
对追尾风险等级为8～10的分向路段，对分向路段间隔1公里设置提醒标识，同时联动电子导航地图进行提醒；
[0177]
对加权追尾风险指数为目标区域内前100的分向路段，对分向路段间隔1公里设置提醒标识、联动电子导航地图进行提醒并设置临时休息区。
[0178]
本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。
[0179]
本领域的普通技术人员将会意识到，这里所述的实施例是为了帮助读者理解本发明的原理，应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合，这些变形和组合仍然在本发明的保护范围内。

技术特征：
1.基于轨迹特征挖掘的追尾风险路段动态识别与管控方法，其特征在于，包括以下步骤：s1、采集目标区域的路网数据，并进行追尾风险路段语义划分；s2、采集目标区域的动态/历史车辆轨迹数据，进行轨迹出行路径匹配，得到追尾风险对应的分向路段；s3、根据目标区域的动态/历史车辆轨迹数据，分别计算匹配分向路段的疲劳指数、流量指数以及急减速指数；s4、根据计算的疲劳指数、流量指数以及急减速指数，计算追尾风险路段的加权追尾风险指数；s5、根据加权追尾风险指数，分级动态识别追尾风险路段的追尾风险等级，并根据风险等级进行分类管控。2.根据权利要求1所述的基于轨迹特征挖掘的追尾风险路段动态识别与管控方法，其特征在于，所述步骤s1具体为：s11、采集目标区域的路网数据，并修复异常路网数据，确定路网信息图；s12、在路网信息图中，将100米以内的所有交汇口端点的形心识别为逻辑道路交叉口；s13、对相邻两个逻辑道路交叉口之间的路段划分逻辑分向路段，实现追尾风险路段的语义划分。3.根据权利要求1所述的基于轨迹特征挖掘的追尾风险路段动态识别方法，其特征在于，所述步骤s2具体为：s21、采集目标区域的动态/历史车辆轨迹数据，对其处理并分别构建各车辆的历史出行轨迹集合和动态出行轨迹集合；s22、遍历各车辆的历史/动态出行轨迹集合中的轨迹点数据；s23、对历史/动态出行轨迹集合中的轨迹点，获取其临近范围内的临近连接线集合；s24、对临近连接线集合不为空的轨迹点，遍历其对应的临近连接线集合，并确定该轨迹点的方位角射线与各临近连接线切线的夹角；s25、根据夹角大小，匹配分向路段：当存在小于45度的夹角时，将最小夹角的临近连接线作为匹配连接线，其对应的虚拟分向路段作为匹配的分向路段；当所有夹角均大于45度时，则若临近连接线最大夹角大于135度，则以该连接线对应分向路段的反向路段作为匹配的分向路段。4.根据权利要求3所述的基于轨迹特征挖掘的追尾风险路段动态识别与管控方法，其特征在于，所述步骤s3中，计算疲劳指数的方法具体为：sa3-1、根据目标区域的路网数据，进行疲劳驾驶路段语义划分；sa3-2、根据目标区域的动态/历史车辆轨迹数据，进行轨迹出行路径匹配，得到疲劳驾驶对应的分向路段；sa3-3、根据疲劳驾驶的分向路段和追尾风险的分向路段的对应关系，对于每个追尾风险的分向路段，在疲劳驾驶的分向路段中划分与追尾风险对应的路段；sa3-4、计算在疲劳驾驶的分向路段中与追尾风险对应的路段的疲劳驾驶指数，作为该追尾风险分向路段的疲劳指数；其中，疲劳指数包括相对疲劳指数和绝对疲劳指数。
5.根据权利要求4所述的基于轨迹特征挖掘的追尾风险路段动态识别与管控方法，其特征在于，所述步骤s3中，计算流量指数的方法具体为：sb3-1、获取历史轨迹出行集合中各轨迹点对应的匹配连接线集合；sb3-2、依次确定各轨迹点处的车速，对于车速大于0的轨迹点，参考其匹配连接线，统计对应分向路段各时段各车道流量；sb3-3、将分向路段内动态时段的平均车道流量作为该分向路段的流量指数。6.根据权利要求5所述的基于轨迹特征挖掘的追尾风险路段动态识别与管控方法，其特征在于，所述步骤s3中，计算急减速指数的方法为：sc3-1、遍历各车辆的历史轨迹点出行集合中的轨迹点；sc3-2、根据前向连续三个轨迹点之间的速度关系，识别减速事件；sc3-3、计算减速事件中中间轨迹点的减速速度；sc3-4、将各分向路段中所属连接线中轨迹点减速速度大于预设值的次数与减速事件总次数的比值，作为该分向路段的急减速指数。7.根据权利要求6所述的基于轨迹特征挖掘的追尾风险路段动态识别与管控方法，其特征在于，所述步骤sc3-2中，当轨迹点p
i+1
的速度v
i+1
＝0，且0<v
i
<v
i-1
时，识别为出现减速事件；其中，下标i为轨迹点序号；所述步骤sc3-3中，中间轨迹点p
i
的减速速度d
i
计算公式为：式中，difv为轨迹点p
i
和p
i-1
的速度差，refdist1为第一参考距离，refdist2为第二参考距离；其中，refdist1＝v
i-1
*diftime，diftime为轨迹点p
i
和p
i-1
间的时间间隔，difdist为轨迹点p
i
和p
i-1
的欧氏距离，diffh为(refdist1-refdist)>1米时，轨迹点p
i-1
和p
i
的高差。8.根据权利要求6所述的基于轨迹特征挖掘的追尾风险路段动态识别与管控方法，其特征在于，所述步骤s4中，追尾风险路段的加权追尾风险指数reindi的计算公式为：reindi＝w1*w2*w3*dindi式中，w1为相对疲劳指数对应的相对疲劳风险等级权重，w2为绝对疲劳指数对应的绝对疲劳风险等级权重，w3为流量指数权重，dindi为急减速指数；其中，相对疲劳指数对应的风险等级根据分向路段的相对疲劳指数进行1～10级的风险等级划分，绝对疲劳风险等级根据分向路段的绝对疲劳风险指数，按自然间断法分10级进行间断排序得到，流量指数权重根据动态时段平均车道流量所属范围确定。9.根据权利要求8所述的基于轨迹特征挖掘的追尾风险路段动态识别与管控方法，其特征在于，所述步骤s5中，动态识别追尾风险等级方法为：根据当前时刻，筛选上一时段车辆采样数大于50次的分向路段，根据其对应的加权追尾风险指数，按自然间断法分10级进行间断排序，获得各分向路段的追尾风险等级。10.根据权利要求9所述的基于轨迹特征挖掘的追尾风险路段动态识别与管控方法，其特征在于，所述步骤s5中，进行分类管控的方法为：
对追尾风险等级小于4的分向路段，不进行管控；对追尾风险等级为5～7的分向路段，联动电子导航地图进行提醒；对追尾风险等级为8～10的分向路段，对分向路段间隔1公里设置提醒标识，同时联动电子导航地图进行提醒；对加权追尾风险指数为目标区域内前100的分向路段，对分向路段间隔1公里设置提醒标识、联动电子导航地图进行提醒并设置临时休息区。

技术总结
本发明公开了一种基于轨迹特征挖掘的追尾风险路段动态识别与管控方法，实现了基于车辆实时轨迹数据分析的追尾风险路段动态识别管控，包括基于道路结构特征的车辆追尾风险分析路段语义分割、基于实时轨迹特征分析的加权追尾风险指数统计(疲劳指数+流量指数+急减速指数)、追尾风险路段分类动态识别，风险路段分类管控，该方法突破了传统管控方法仅着眼于车辆本身追尾行为导致管控效率低、车辆传感设备依赖性高的局限，可显著提高道路网络整体追尾风险路段动态识别管控水平，降低相关事故风险。险。险。


技术研发人员：吴洲豪 党晓旭 周晓东 王玉娇 汪泽 冯套柱 刘圆圆 王春祥 李玲 金书鑫 王秋玉 张越 康秋会 杨柳
受保护的技术使用者：西安科技大学
技术研发日：2023.03.14
技术公布日：2023/6/12