本发明属于道路交通安全分析技术，具体为一种基于车辆加速度的碰撞风险评价方法。
背景技术：
传统道路交通安全评价方法主要基于事故统计模型。交通事故是一种小概率随机事件，因此建立科学可靠的事故统计模型需要多年的数据积累，一旦数据质量出现问题，统计学模型的结果将会产生偏差。
与之相反的，交通冲突技术是一种典型的基于非事故数据的交通安全分析与评价方法，以大样本、高时效、低社会成本、定量分析与评价交通安全的特点而异于传统的以事故数据统计分析为基础的方法。目前，学者对于机动车之间冲突的研究较为成熟，已验证交通冲突技术可用于道路交通安全评价方法。
近年来，在交通安全评价方面，有学者使用极值理论对信号交叉口的事故频率以及高速公路的事故频率进行预测，并取得了良好结果。总体来看，现阶段基于交通冲突的安全评价主要利用冲突数(率)，且冲突本身的可靠性也缺乏验证，存在可信度低的缺点；而极值统计方法则更多的被应用在建立交通冲突与交通事故的关系方面，没有面向交通安全评价的相关指标。
技术实现要素：
本发明的目的在于提出了一种基于车辆加速度的碰撞风险评价方法。
实现本发明目的的技术解决方案为：一种基于车辆加速度的碰撞风险评价方法，具体步骤为：
基于视频数据，获取目标车辆的速度与相邻车辆的速度，确定目标车辆与相邻车辆的间距；
根据间距阈值d判定相邻车辆与目标车辆是否存在冲突可能性：当目标车辆与相邻车辆的间距小于间距阈值d，相邻车辆与目标车辆存在冲突可能性，否则不存在冲突可能性；
根据目标车辆与冲突车辆运动状态计算避免碰撞的减速度drac；
根据避免碰撞的减速度drac，构建超阈值极值分布模型，对超阈值极值分布模型进行估计，得到超阈值极值分布模型中广义pareto极值分布的尺度参数及形状参数；
根据避免碰撞的减速度drac和超阈值极值分布模型，构建冲突极值模型对碰撞风险进行预测。
优选地，所述相邻车辆包括目标车辆的前车、后车、相邻车道变道车辆、相交路口与之垂直方向车辆以及相交路口右侧左转弯车辆。
优选地，间距阈值d具体为：
式中，vg(t)为目标车辆t时刻的速度；vb(t)为相邻车辆t时刻的速度；ag为目标车辆的加速度。
优选地，根据目标车辆与冲突车辆运动状态计算避免碰撞的减速度drac分别具体为：
①相邻车道后车变道为跟驰行驶
②相邻车道后车变道超车行驶
③相交路口与之垂直方向车辆冲突行驶
④相交路口右侧车辆左转弯行驶
式中，v1(t)为前车t时刻的速度；v2(t)为后车t时刻的速度；s12(t)为t时刻两车之间的距离，α为车辆行驶方向与路边夹角。
优选地，根据避免碰撞的减速度drac，构建超阈值极值分布模型的具体过程为：
利用平均剩余寿命图确定超阈值极值分布模型的阈值，具体为：
式中，u为超阈值极值分布模型的阈值；yj为超出阈值u的样本，j＝1,2,3……n，ymax为最大样本；
统计超过阈值的样本组成新的样本集合，具体为：对获得的交通冲突数据进行处理，将超过阈值的数据放入集合x内，组成可用于超阈值极值分布模型的集合：
步骤5-3：构建超阈值极值分布模型，具体为：
式中，h(xi,σ,ξ)为广义pareto分布函数；xi为超过阈值的样本，i＝1,2,3…；x为由xi组成的集合；μ为广义pareto极值分布函数的位置参数；σ为广义pareto极值分布函数的尺度参数；ξ为广义pareto极值分布函数的形状参数。
优选地，对超阈值极值分布模型进行估计，得到超阈值极值分布模型中广义pareto极值分布的尺度参数及形状参数的具体方法为：
令位置参数μ＝0；
采用极大似然法对超阈值极值分布模型进行估计得到尺度参数，具体为：
p{x＝x}＝p(x；θ)
p{x＝x}＝p(x；θ)为集合x的概率密度函数；θ＝(θ1，θ2，…θn)为模型参数；l(θ)为似然函数；
优选地，构建的冲突极值模型具体为：
事故风险指标，计算公式为：
r为事故风险；h(u)为广义pareto分布在阈值u处的取值；
事故数预测指标，计算公式为：
nt为预测期tp时段内的预测事故次数；tp为预测期；
事故重现水平指标，计算公式为：
zt为对应于重现期t的重现水平；t为重现期；n为观测得到的数据总量；
本发明与现有技术相比，其显著优点为：
本发明将交通冲突和交通安全相关联，通过分析交通冲突来评价道路交通安全水平，通过视频采集路段或交叉口高峰时段的运行状况，利用基于多目标追踪算法的图像处理技术可同时提取多辆车辆轨迹，即可得到目标车辆与周围前车、后车、相邻车道变道车辆、相交路口与之垂直方向车辆以及相交路口右侧左转弯车辆的距离；结合极值模型对交通冲突数据进行建模，构建高时效的交通安全分析指标，实现以分钟级别的的道路交通安全分析；
本发明所选取的交通冲突指标为避免碰撞的减速度drac，冲突过程中的每一刻都对应一个drac值，故通过短时间的观测即可获取大量的交通冲突样本，本发明提出了四种不同车辆行驶情况下的避免碰撞的减速度drac的计算方法，可对不同的场景进行应用；
本发明所提出的指标能够反映以分钟为级别的交通安全水平变化情况，并且该安全水平指标能够直接与交通事故风险相联系。
下面结合附图对本发明做进一步详细的描述。
附图说明
图1是本发明的流程图。
图2为跟驰情况drac计算示意图。
图3位相邻车道后车变道为跟驰行驶情况drac计算示意图。
图4相邻车道后车变道为引领行驶情况drac计算示意图。
图5相交路口与之垂直方向车辆冲突行驶情况drac计算示意图。
图6相交路口右侧车辆左转弯行驶情况drac计算示意图。
具体实施方式
一种基于车辆加速度的碰撞风险评价方法，具体过程为：
步骤1：利用无人机航拍得到目标视频，基于多目标追踪算法对目标视频内的车辆进行轨迹提取，可同时获取目标车辆的速度与相邻车辆的速度，所述相邻车辆包括目标车辆的前车、后车、相邻车道变道车辆、相交路口与之垂直方向车辆以及相交路口右侧左转弯车辆；
步骤2：计算目标车辆与相邻车辆的间距；
步骤3：根据间距阈值d判定相邻车辆与目标车辆是否存在冲突可能性：当目标车辆与相邻车辆的间距小于间距阈值d，相邻车辆与目标车辆存在冲突可能性，否则不存在冲突可能性，其中，间距阈值d的计算公式为：
式中，vg(t)为目标车辆t时刻的速度；vb(t)为相邻车辆t时刻的速度；ag为目标车辆的加速度；在某些实施例中，ag＝3.35m/s2，原因为如果给定车辆的drac阈值超过3.35m/s2，则该车辆处于交通冲突中。
步骤4：根据目标车辆与冲突车辆运动状态计算避免碰撞的减速度(decelerationratetoavoidacrash，drac)，具体包括为：
①相邻车道后车变道为跟驰行驶
②相邻车道后车变道超车行驶
③相交路口与之垂直方向车辆冲突行驶
④相交路口右侧车辆左转弯行驶
式中，v1(t)为前车t时刻的速度；v2(t)为后车t时刻的速度；s12(t)为t时刻两车之间的距离，α为车辆行驶方向与路边夹角。
根据上述公式得到目标车辆与相邻冲突车辆的drac。
本发明是利用时空接近型交通冲突指标drac，drac可以考虑到速度差异和减速在潜在的碰撞发生中的作用，碰撞减速度的大小可有用来度量发生追尾事故的风险度，若计算出的碰撞减速度超过了车辆本身可利用的最大减速度，发生追尾事故是难免的，如图2—6所示。
步骤5：根据步骤4获得的drac，以目标车辆与相邻冲突车辆计算得到的drac为变量构建超阈值极值分布模型，对超阈值极值分布模型进行估计，得到超阈值极值分布模型中广义pareto极值分布的尺度参数及形状参数，具体过程为：
步骤5-1：利用平均剩余寿命图确定超阈值极值分布模型的阈值，具体公式为：
式中，u为超阈值极值分布模型的阈值；yj为超出阈值u的样本，j＝1,2,3……n，ymax为最大样本。
公式(7)所描述的一系列点所构成的图即为平均剩余寿命图，如果存在一个阈值范围，使得超出该范围内任一阈值的样本均服从广义pareto分布，那么在该范围内的平均剩余寿命图应该是线性或近似线性的；选取该范围内最大值作为超阈值极值分布模型的合理阈值。
步骤5-2：统计超过阈值的样本组成新的样本集合，具体为：对获得的交通冲突数据进行处理，将超过阈值的数据放入集合x内，组成可用于超阈值极值分布模型的集合：
步骤5-3：构建超阈值极值分布模型，具体为：
式中，h(xi,σ,ξ)为广义pareto分布函数；xi为超过阈值的样本，i＝1,2,3…，满足公式(5)；x为由xi组成的集合；μ为广义pareto极值分布函数的位置参数；σ为广义pareto极值分布函数的尺度参数；ξ为广义pareto极值分布函数的形状参数；
超阈值极值分布模型服从广义pareto极值分布函数，由于位置参数μ的大小不影响广义pareto极值分布的形状，因此令公式(9)中位置参数μ＝0。
步骤5-4：采用极大似然方法对超阈值极值分布模型进行估计，得到超阈值极值分布模型中广义pareto极值分布的尺度参数及形状参数，具体过程为：
采用极大似然方法对超阈值极值分布模型进行估计，极大似然法是参数估计的方法之一，在参数θ的可能取值范围内，选取使l(θ)达到最大的参数值θ，作为参数θ的估计值。
公式为：
p{x＝x}＝p(x；θ)(11)
公式(11)为集合x的概率密度函数；θ＝(θ1，θ2，…θn)为模型参数；l(θ)为似然函数；
步骤6：根据避免碰撞的减速度drac和超阈值极值分布模型，构建冲突极值模型对碰撞风险进行预测，具体为：
所述冲突极值模型包括事故风险指标，事故数预测指标和事故重现水平指标，分别具体为：
构建事故风险指标，计算公式为：
r为事故风险；h(u)为广义pareto分布在阈值u处的取值；
构建事故数预测指标，计算公式为：
nt为预测期tp时段内的预测事故次数；tp为预测期；
构建事故重现水平指标，计算公式为：
zt为对应于重现期t的重现水平；t为重现期；n为观测得到的数据总量；
根据所构建的三个安全评价指标对道路设施的交通安全水平进行评价，具体评价标准为：事故风险指标能够直接反映发生交通事故的风险，与交通事故数(率)相关联，当该值为0时意为该时段不会发生交通事故，当该值大于0时即为有发生交通事故的风险；事故数预测指标同样与交通事故数(率)相关联，当其大于0时，则该时段会发生交通事故、安全水平低，并且事故数预测值越大，安全水平越低，否则不会发生交通事故、安全水平高；事故重现水平指标意为该时段发生交通事故的概率，得到的结果越大则意味着发生交通事故的概率越大、安全水平越低，且该指标与交通事故数(率)的直接关联较弱。
1.一种基于车辆加速度的碰撞风险评价方法，其特征在于，具体步骤为：
基于视频数据，获取目标车辆的速度与相邻车辆的速度，确定目标车辆与相邻车辆的间距；
根据间距阈值d判定相邻车辆与目标车辆是否存在冲突可能性：当目标车辆与相邻车辆的间距小于间距阈值d，相邻车辆与目标车辆存在冲突可能性，否则不存在冲突可能性；
根据目标车辆与冲突车辆运动状态计算避免碰撞的减速度drac；
根据避免碰撞的减速度drac，构建超阈值极值分布模型，对超阈值极值分布模型进行估计，得到超阈值极值分布模型中广义pareto极值分布的尺度参数及形状参数；
根据避免碰撞的减速度drac和超阈值极值分布模型，构建冲突极值模型对碰撞风险进行预测。
2.根据权利要求1所述的基于车辆加速度的碰撞风险评价方法，其特征在于，所述相邻车辆包括目标车辆的前车、后车、相邻车道变道车辆、相交路口与之垂直方向车辆以及相交路口右侧左转弯车辆。
3.根据权利要求1所述的基于车辆加速度的碰撞风险评价方法，其特征在于，间距阈值d具体为：
式中，vg(t)为目标车辆t时刻的速度；vb(t)为相邻车辆t时刻的速度；ag为目标车辆的加速度。
4.根据权利要求1所述的基于车辆加速度的碰撞风险评价方法，其特征在于，根据目标车辆与冲突车辆运动状态计算避免碰撞的减速度drac分别具体为：
①相邻车道后车变道为跟驰行驶
②相邻车道后车变道超车行驶
③相交路口与之垂直方向车辆冲突行驶
④相交路口右侧车辆左转弯行驶
式中，v1(t)为前车t时刻的速度；v2(t)为后车t时刻的速度；s12(t)为t时刻两车之间的距离，α为车辆行驶方向与路边夹角。
5.根据权利要求1所述的基于车辆加速度的碰撞风险评价方法，其特征在于，根据避免碰撞的减速度drac，构建超阈值极值分布模型的具体过程为：
利用平均剩余寿命图确定超阈值极值分布模型的阈值，具体为：
式中，u为超阈值极值分布模型的阈值；yj为超出阈值u的样本，j＝1,2,3……n，ymax为最大样本；
统计超过阈值的样本组成新的样本集合，具体为：对获得的交通冲突数据进行处理，将超过阈值的数据放入集合x内，组成可用于超阈值极值分布模型的集合：
步骤5-3：构建超阈值极值分布模型，具体为：
式中，h(xi,σ,ξ)为广义pareto分布函数；xi为超过阈值的样本，i＝1,2,3…；x为由xi组成的集合；μ为广义pareto极值分布函数的位置参数；σ为广义pareto极值分布函数的尺度参数；ξ为广义pareto极值分布函数的形状参数。
6.根据权利要求1所述的基于车辆加速度的碰撞风险评价方法，其特征在于，对超阈值极值分布模型进行估计，得到超阈值极值分布模型中广义pareto极值分布的尺度参数及形状参数的具体方法为：
令位置参数μ＝0；
采用极大似然法对超阈值极值分布模型进行估计得到尺度参数，具体为：
p{x＝x}＝p(x；θ)
p{x＝x}＝p(x；θ)为集合x的概率密度函数；θ＝(θ1，θ2，…θn)为模型参数；l(θ)为似然函数；
7.根据权利要求1所述的基于车辆加速度的碰撞风险评价方法，其特征在于，构建的冲突极值模型具体为：
事故风险指标，计算公式为：
r为事故风险；h(u)为广义pareto分布在阈值u处的取值；
事故数预测指标，计算公式为：
nt为预测期tp时段内的预测事故次数；tp为预测期；
事故重现水平指标，计算公式为：
zt为对应于重现期t的重现水平；t为重现期；n为观测得到的数据总量；
