CN104992557A - 一种城市交通警情等级预测方法 - Google Patents
一种城市交通警情等级预测方法 Download PDFInfo
- Publication number
- CN104992557A CN104992557A CN201510242669.6A CN201510242669A CN104992557A CN 104992557 A CN104992557 A CN 104992557A CN 201510242669 A CN201510242669 A CN 201510242669A CN 104992557 A CN104992557 A CN 104992557A
- Authority
- CN
- China
- Prior art keywords
- alert
- traffic
- attribute
- grade
- probability
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 40
- 239000003595 mist Substances 0.000 claims description 12
- 230000008569 process Effects 0.000 claims description 12
- 238000012549 training Methods 0.000 claims description 12
- 238000010219 correlation analysis Methods 0.000 claims description 7
- 238000012856 packing Methods 0.000 claims description 3
- 238000003672 processing method Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 description 12
- 238000012360 testing method Methods 0.000 description 7
- 230000008859 change Effects 0.000 description 4
- 206010039203 Road traffic accident Diseases 0.000 description 3
- 238000003066 decision tree Methods 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 206010019345 Heat stroke Diseases 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000004913 activation Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000012098 association analyses Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000000428 dust Substances 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000009440 infrastructure construction Methods 0.000 description 1
- 230000000977 initiatory effect Effects 0.000 description 1
- 230000031700 light absorption Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000001556 precipitation Methods 0.000 description 1
- 230000007115 recruitment Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000007790 scraping Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000007958 sleep Effects 0.000 description 1
- 238000009827 uniform distribution Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- 230000002618 waking effect Effects 0.000 description 1
Landscapes
- Traffic Control Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一种城市交通警情等级预测方法,包括如下步骤:步骤1:时间数据整理,时间数据包括5维特征属性,分别是第几单位时间、星期几、是否节假日、是否工作日和是否有重大活动;步骤2:天气数据与交通警情关联性分析及处理:天气数据包括温度、可见度、降雨量3个特征属性;步骤3:建立预测交通警情等级的属性集:根据数据库信息,将处理后的时间数据、天气数据、重大活动日、中小学放假特征属性整合成多维数据宽表,其中包括8个离散的特征属性维度;步骤4:辖区内交通警情等级划分;步骤5:构造朴素贝叶斯分类器;步骤6:预测各区域的交通警情等级。本发明提供了一种有效实现预测、准确性较好的城市交通警情等级预测方法。
Description
技术领域
本发明涉及智能交通领域,具体涉及一种城市辖区交通警情等级预测方法。
背景技术
随着我国经济的高速增长和现代交通的日益发达,与交通有关的警情数量不断增多,交通警情主要分为交通拥堵、交通事故、群众求助等,一般包含警情发生、警情上报、警情处理等过程,大量警情的发生给交通管理者造成很大的困扰。随着城市基础设施建设的不断完善和城市区域的扩大,给交警指挥部门合理调配警力带来很大挑战,将交警呈网格状均匀分配在区域内的部署方式需要耗费大量警力。
论文《基于时空数据应用模型的交通警情研判实现》通过阈值方法进行警情分级预警,其实质是当前警情超出预期以后的报警,并未建立预测模型,不能对未来多天各时段的警情情况进行预测和分级警示,其研究对实时指挥调度有支撑作用,但对勤务优化配置预案的制定无支撑作用;专利201410381532.4选取天气情况、道路情况、驾驶员反应时间、后车与前车距离、后车速度、后车相对于前车的速度差、后车加速度作为追尾事故的变量节点构建贝叶斯网络学习数据集,该方法仅侧重于针对追尾事故中的车辆建模,并采取预警措施。
发明内容
为了克服已有交通警情判别方式的无预测功能、准确性较差的不足,本发明提供了一种有效实现预测、准确性较好的城市交通警情等级预测方法。
本发明解决其技术问题所采用的技术方案是:
一种城市交通警情等级预测方法,所述预测方法包括如下步骤:
步骤1:时间数据整理
时间数据包括5维特征属性,分别是第几单位时间、星期几、是否节假日、是否工作日和是否有重大活动,所述第几单位时间属性处理方法为设定单位时间Δt=24h/N,按单位时间Δt将全天24小时分为N段;获取国务院颁布的法定节假日信息构成节假日属性;依据节假日和周末信息构成工作日属性;根据当天是 否有大型考试、当前辖区内学校是否放假、是否有大型场馆活动定义是否有重大活动属性;
步骤2:天气数据与交通警情关联性分析及处理
天气数据包括温度、可见度、降雨量3个特征属性;
步骤3:建立预测交通警情等级的属性集
根据数据库信息,将处理后的时间数据、天气数据、重大活动日、中小学放假等特征属性整合成多维数据宽表,其中包括8个离散的特征属性维度;
步骤4:辖区内交通警情等级划分
统计各辖区近期6个月内的警情发生次数,采用动态划分的方法将交通警情发生次数划分为三个等级,将各个交警支队在单位时间Δt内的警情数量采用等区间装箱方法求和,构成(x1,x2,…,xn)关于交通警情数量的一个样本空间,由n个样本点组成,其中可计算出当前交警支队的警情数量均值记为μ, 计算出警情数量的标准差记为σ,故警情数量属于[0,μ-σ)的标记警情等级为C1,警情数量属于[μ-σ,μ+σ)的标记警情等级为C2,警情数量属于[μ+σ,+∞)的标记警情等级为C3;
步骤5:构造朴素贝叶斯分类器
选取第几三小时X1、星期几X2、是否节假日X3、是否工作日X4、是否有重大活动X5、温度X6、可见度X7、降雨量X8构成预测交通警情等级的多维数据集T={X1,X2,X3,X4,X5,X6,X7,X8},将标记好的警情等级与多维数据集共同构成朴素贝叶斯分类的训练集R={X1,X2,X3,X4,X5,X6,X7,X8,Ci},其中Ci∈(C1,C2,C3),从而使每个训练样本都具有相应的警情等级类别标签。已知未来几天各个特征属性Ti={X1,X2,X3,X4,X5,X6,X7,X8}=d预测当前区域的交通警情等级,根据朴素贝叶斯分类器的原理可知需要比较在Ti=d的条件下属于C1、C2、C3三类警情等级的概率大小,即需要计算p(C1|Ti=d)、p(C2|Ti=d)、p(C3|Ti=d)三者的概率。由贝叶斯公式可知其 中Ci∈(C1,C2,C3),比较p(C1|Ti=d)、p(C2|Ti=d)、p(C3|Ti=d)三者的概率大小只需要比较p(C1)p(Ti=d|C1)、p(C2)p(Ti=d|C2)、p(C3)p(Ti=d|C3)三者的大小,分母p(Ti=d)对所有类别为常数;
通过统计训练集R中属于警情等级C1的样本个数可计算出p(C1)的概率,即 同理可计算出 p(Ti=d|Ci)即在类别Ci中属性Ti的取值为d的概率,具体p(Ti=d|Ci)=p(X1=d1|Ci)p(X2=d2|Ci)…p(X8=d8|Ci),其中d=(d1,d2,…,d8),8个离散特征属性的先验概率p(X1=d1|Ci)、…、p(X8=d8|Ci)通过计算在每个类别Ci中每个属性Xi取值为di的概率p(Xi=di|Ci)求得;
步骤6:预测各区域的交通警情等级
由步骤5形成各特征属性的概率表,将所有特征属性的概率相乘计算出p(C1|Ti)、p(C2|Ti)、p(C3|Ti)三者的概率大小,将所有特征属性的概率相乘时如果有某个特征属性的概率值为零,会消除乘积中涉及的所有其他属性概率的影响,导致最终p(Ci|Ti)的结果为零,每个特征属性从1开始计数,p(C1|Ti)、p(C2|Ti)、p(C3|Ti)三者中最大的概率值即为朴素贝叶斯分类器预测出的当前特征属性条件Ti下的交通警情等级。
进一步,所述步骤2中,将降雨量分为以下五个等级,降雨量0.1mm-9.9mm时为小雨,降雨量为10mm-24.9mm为中雨,降雨量25mm-49.9mm为大雨,降雨量50mm-99.9mm为暴雨,降雨量100mm以上为大暴雨;将可见度等级分为以下五类,可见度大于10km时为无雾,可见度为1km-10km时为轻雾,可见度为200m-1km时大雾,可见度为50m-200m时为浓雾,可见度小于50m时为强雾;将温度划分为小于5℃,5℃-28℃,29℃-35℃,大于35℃四个等级。
本发明的技术构思为:在数据来源和实现目标方面均与现有技术不同,本方法通过将天气因素、时间因素、重大活动、中小学放假作为特征属性,充分考虑了车、路、人、环境等因素对交通警情等级的影响,最终将整理后的多维数据通 过朴素贝叶斯分类器对不同辖区的交通警情等级进行预测。目前尚未发现以天气因素、时间因素、重大活动、中小学放假为特征属性的交通警情等级预测方法,针对全市各个交警支队所属辖区的交通警情等级预测结果可以优化警力部署,提高交通管理者的指挥效率。
本发明的有益效果主要表现在:(1)城市辖区交通警情等级预测是为了给所属辖区内的交警支队制定警力配置预案提供数据支撑,本方法动态划分各个辖区交通警情等级,对警情数量超出合理范围的辖区进行预警,交通警情等级展示结果简洁,便于交通管理者理解。
(2)特征属性的选择全面系统地揭示天气因素、时间因素、环境因素对交通警情等级的影响。通过特征属性与交通警情等级的关联性分析,并对连续属性值离散化,构造出符合朴素贝叶斯分类器的多维特征属性,利用概率统计知识进行分类,方法简单、分类准确率高、速度快,充分利用了贝叶斯分类器刻画外部环境变化如何影响交通警情等级变化,能够从原始样本数据中挖掘出天气数据、时间数据之间的耦合关系,为交警指挥部门调度提供数据支撑。
附图说明
图1是城市交通警情等级预测流程图。
图2是杭州市降雨等级与警情关系图。
图3是杭州市可见度与警情关系图。
图4是杭州市温度等级与警情关系图。
具体实施方式
下面结合附图对本发明作进一步描述。
参照图1~图4,一种城市交通警情等级预测方法,包括如下步骤:
步骤1:时间数据整理
时间数据包括5维特征属性,分别是第几单位时间、星期几、是否节假日、是否工作日和是否有重大活动,所述第几单位时间属性处理方法为设定单位时间Δt=24h/N,按单位时间Δt将全天24小时分为N段(例如N=8,单位时间为3h);获取国务院颁布的法定节假日信息构成节假日属性;依据节假日和周末信息构成工作日属性;根据当天是否有大型考试、当前辖区内学校是否放假、是否有大型 场馆活动定义是否有重大活动属性;
步骤2:天气数据与交通警情关联性分析及处理
天气数据包括温度、可见度、降雨量3个特征属性。降雨量往往会引起道路交通环境变差,导致路面摩擦系数下降,影响车辆制动,降雨量大时能见度变坏,影响视线,易发生交通事故,由降雨量与警情关联性分析结果可知,当天气为小到中雨或中雨时,交通事故量较多,当降雨量特别大时,事故量相对减少,因此降水和降雨量的大小与事故的关系较复杂,综合考虑以上因素结合天气预报标准可将降雨量分为以下五个等级,降雨量0.1mm-9.9mm时为小雨,降雨量为10mm-24.9mm为中雨,降雨量25mm-49.9mm为大雨,降雨量50mm-99.9mm为暴雨,降雨量100mm以上为大暴雨;雾对光的吸收和散射作用,使得道路能见度、视野下降,能见度有时降到30~40m甚至更低,驾驶员看不清前方和周围的情况,对前方车辆等识别困难,极易发生追尾等事故,雾中的水珠与路面上和空气中的灰尘混合,导致轮胎与路面附着系数减小,从而导致制动距离延长、行驶打滑等现象发生,进而引发道路交通事故,根据历史数据的相关性分析结果可知,可见度小于1km时与警情数呈负相关,可见度大于1km时与警情关系不明显,结合国家能见度等级和预报标准,将可见度等级分为以下五类,可见度大于10km时为无雾,可见度为1km-10km时为轻雾,可见度为200m-1km时大雾,可见度为50m-200m时为浓雾,可见度小于50m时为强雾;冬季温度突变不仅影响驾驶员、行人对环境反应的灵敏度,还会引起道路摩擦系数降低,车辆通过困难,春季气温回升快,此时驾驶员容易处于似睡似醒状态,同时路面湿滑极易引起车祸,相关研究表明在日平均相对湿度大于85%,平均气温大于17℃时,日均交通事故指数最大为4.69,有时春季也会出现倒春寒现象,雨雪交加的恶劣天气导致车祸频频发生,夏季气温大于35℃会出现汽车轮胎爆胎、驾驶员中暑等易引起交通事故的因素,气温与警情关联分析后可将温度划分为小于5℃,5℃-28℃,29℃-35℃,大于35℃四个等级。
步骤3:建立预测交通警情等级的属性集
由时间因素、天气因素、重大活动日、中小学放假等特征属性构成的多维数据间必然存在着直接或间接的规律和联系,例如:处于上下班的时间段肯定交通流数据比较大,可能引发的交通事故案件比较多,重大活动导致的人员出行需求剧增也有可能引发交通事故。各个特征属性间存在潜在的联系,根据数据库信息, 将处理后的时间数据、天气数据、重大活动日、中小学放假等特征属性整合成多维数据宽表,其中包括8个离散的特征属性维度,如表1所示,表1为预测交通警情等级的特征属性集。第几三小时、星期几、是否节假日、是否工作日、是否有重大活动、温度、可见度、降雨量等8个离散的特征属性维度。
表1
步骤4:辖区内交通警情等级划分
各个交警支队所属辖区内交通警情发生次数均有差异,核心辖区警情发生次数比偏远辖区多,采用固定阈值将警情发生次数等区间划分不符合实际情况,对各个交警支队指挥调度没有指导意义,统计各辖区近期6个月内的警情发生次数,采用动态划分的方法将交通警情发生次数划分为三个等级,将各个交警支队在单位时间Δt=3h内的警情数量采用等区间装箱方法求和,构成(x1,x2,…,xn)关于交通警情数量的一个样本空间,由n个样本点组成,其中可计算出当前交警支队的警情数量均值记为μ,计算出警情数量的标准差记为σ,故警情数量属于[0,μ-σ)的标记警情等级为C1,警情数量属于[μ-σ,μ+σ)的标记警情等级为C2,警情数量属于[μ+σ,+∞)的标记警情等级为C3。
步骤5:构造朴素贝叶斯分类器
选取第几三小时X1、星期几X2、是否节假日X3、是否工作日X4、是否有重大活动X5、温度X6、可见度X7、降雨量X8构成预测交通警情等级的多维数据集T={X1,X2,X3,X4,X5,X6,X7,X8}。将标记好的警情等级与多维数据集共同构成朴素贝叶斯分类的训练集R={X1,X2,X3,X4,X5,X6,X7,X8,Ci},其中Ci∈(C1,C2,C3),从而使每个训练样本都具有相应的警情等级类别标签。已知未来几天各个特征属性Ti={X1,X2,X3,X4,X5,X6,X7,X8}=d预测当前区域的交通警情等级,根据朴素贝叶斯分类器的原理可知需要比较在Ti=d的条件下属于C1、C2、C3三类警情 等级的概率大小,即需要计算p(C1|Ti=d)、p(C2|Ti=d)、p(C3|Ti=d)三者的概率。由贝叶斯公式可知其中Ci∈(C1,C2,C3),比较p(C1|Ti=d)、p(C2|Ti=d)、p(C3|Ti=d)三者的概率大小只需要比较p(C1)p(Ti=d|C1)、p(C2)p(Ti=d|C2)、p(C3)p(Ti=d|C3)三者的大小,分母p(Ti=d)对所有类别为常数,无需计算。
通过统计训练集R中属于警情等级C1的样本个数可计算出p(C1)的概率,即 同理可计算出 p(Ti=d|Ci)即在类别Ci中属性Ti的取值为d的概率,具体p(Ti=d|Ci)=p(X1=d1|Ci)p(X2=d2|Ci)…p(X8=d8|Ci),其中d=(d1,d2,…,d8),8个离散特征属性的先验概率p(X1=d1|Ci)、…、p(X8=d8|Ci)通过计算在每个类别Ci中每个属性Xi取值为di的概率p(Xi=di|Ci)即可求得。
步骤6:预测各区域的交通警情等级
由步骤5形成各特征属性的概率表,将所有特征属性的概率相乘计算出p(C1|Ti)、p(C2|Ti)、p(C3|Ti)三者的概率大小,将所有特征属性的概率相乘时如果有某个特征属性的概率值为零,会消除乘积中涉及的所有其他属性概率的影响,导致最终p(Ci|Ti)的结果为零,本方法假定训练数据R很大,以至于对每个特征属性计数加1造成的估计概率的变化可以忽略不计,故每个特征属性从1开始计数,此方法可以方便的避免概率值为零。p(C1|Ti)、p(C2|Ti)、p(C3|Ti)三者中最大的概率值即为朴素贝叶斯分类器预测出的当前特征属性条件Ti下的交通警情等级。
实例1:一种城市辖区交通警情等级预测方法,实施过程如下:
步骤1:时间数据整理
目前杭州几种比较典型的对交通影响较大的活动主要有:黄龙体育中心举办的大型演唱会、中超联赛,和平会展中心举办的大型招聘会和各类会展,石祥路 汽车城举办的大型车展,国际马拉松活动,钱江新城大型演出等。该表数据来源均为互联网。以下城区交警支队为例,下城区范围内主要的大型场馆有和平国际会展中心,从和平会展中心官方网站抓取历年的活动记录,在Oracle数据库中建立杭州市场馆活动情况表EVENT_HUODONG,各字段注释如表2所示,表2为杭州市场馆活动情况表EVENT_HUODONG:
PLACE | 场馆地点 | char(50) |
TYPE | 场馆活动类型 | char(50) |
TIME_BEGIN | 场馆活动起始时间 | DATE |
TIME_END | 场馆活动终止时间 | DATE |
AREA | 场馆所在的区域 | char(20) |
表2
杭州市场馆活动情况表内数据进行存储;
查询杭州市各大高校的寒暑假时间和中小学寒暑假时间,制成杭州市寒暑假时间记录表EVENT_SCHOOL。表内的各字段注释如表3所示,表3为杭州市寒暑假时间记录表EVENT_SCHOOL:
PLACE | 学校地点 | char(50) |
TYPE | 学校活动类型 | char(50) |
TIME_BEGIN | 学校活动起始时间 | DATE |
TIME_END | 学校活动终止时间 | DATE |
AREA | 学校所在的区域 | char(20) |
表3
杭州市寒暑假时间记录表内的数据进行存储。
步骤2:天气数据与交通警情关联性分析及处理
建立杭州市天气数据表HZTQ5,数据采集范围是从2012年1月份到2015年2月份,杭州市天气数据表格字段注释如表4所示,表4为杭州市天气数据表HZTQ5:
TIME | 时间,每隔三小时一次 | VARchar2(50) |
T | 温度,摄氏度 | NUMBER |
VV | 水平能见度(千米) | NUMBER |
RRR | 降雨量(毫米) | NUMBER |
表4
杭州市天气数据表数据进行存储。
如图2所示,分析降雨量与警情数量关系可知,天气为小到中雨或中雨时,交通事故量较多,当降雨量特别大时,事故量相对较少,进一步分析可知连续降雨5天或以上情况下,前两三天刮擦、堵路、事故警情呈现上升趋势,而后下降。
表5为可见度在1km以内与警情数之间的相关系数表,由表5和图3所示,分析可见度与警情数量关系可知,可见度小于1km时与警情数基本成负相关,其中在0.2-0.5km时负相关性最显著,相关系数大约为-0.76;而可见度大于1km时,两者没有明显的关系,浓雾天的能见度太低,发生事故量最多,轻雾天的交通事故量比雾天和大雾天的事故量都要多。
可见度 | 小于0.2km | 0.2-0.5km | 0.5-0.8km | 0.8-1km |
相关系数 | 0.009 | -0.76 | -0.38 | -0.05 |
表5
由图4所示,分析温度等级与警情数量关系可知,温度在28℃以下时与警情之间的关系不大,大于28℃时正相关性比较大,夏季最高气温大于35℃时日均交通事故量最多,车辆追尾类型在寒冷天易发生,可能与道路结冰有关。
步骤3:建立预测交通警情等级的属性集
预测交通警情等级的特征属性集中包括第几三小时、星期几、是否节假日、是否工作日、是否有重大活动、温度、可见度、降雨量等8个特征属性维度,将全天24小时划分为23:00--2:00为第1个三小时,2:00--5:00为第2个三小时,5:00--8:00为第3个三小时,8:00--11:00为第4个三小时,11:00--14:00为第5个三小时,14:00--17:00为第6个三小时,17:00--20:00为第7个三小时,20:00--23:00为第8个三小时,分别用[1,2,3,4,5,6,7,8]表示,星期一至星期日用[1,2,3,4,5,6,7]表示,工作日用1表示,非工作日用0表示,若当前区域内有大型活动、放假等情况用1表示,温度离散化为[1,2,3,4]四个等级,可见度离 散化为[1,2,3,4,5]五个等级,降雨量离散化为[1,2,3,4,5]五个等级。整理好的特征属性表如表6所示,表6为预测交通警情等级的特征属性集,数据库中存储形式,表6与杭州市温度等级与警情关系图中特征属性维度从左至右均一一对应。
表6
步骤4:辖区内交通警情等级划分
选取近6个月内各个交警支队的交通警情数据,将全天24小时分成8份,每份为3小时时间段,与采集到的天气数据时间段吻合。23:00--2:00为第1个三小时,2:00--5:00为第2个三小时,5:00--8:00为第3个三小时,8:00--11:00为第4个三小时,11:00--14:00为第5个三小时,14:00--17:00为第6个三小时,17:00--20:00为第7个三小时,20:00--23:00为第8个三小时,采用等区间装箱的方式将三小时时段内所有警情数量求和,构成(x1,x2,…,xn)关于当前辖区交通警情数量的一个样本空间,由n个样本点组成,其中可计算出当前交警支队的警情数量均值记为μ,计算出警情数量的方差记为σ,故警情数量属于[0,μ-σ)的标记警情等级为C1,警情数量属于[μ-σ,μ+σ)的标记警情等级为C2,警情数量属于[μ+σ,+∞)的标记警情等级为C3。
步骤5:构造朴素贝叶斯分类器并预测各个区域的交通警情等级
以事故类警情预警为例,将2013年1月至2015年2月数据分为两部分,2013年1月至2014年11月1日共4977条数据作为训练数据集,2014年11月1日至2015年2月3日共738条数据作为测试数据集。训练数据集中的4977条数据输 入分类器训练出一组与特征值有关的系数,去除738条测试集中的警情等级标签,然后将738条测试数据输入训练好的分类器中,记录下预测出的警情等级,将预测出的警情等级与测试集中原有的警情等级标签进行对比,测试预测结果的准确率。实验分别采用贝叶斯分类器和决策树分类器进行分类,并得到测试结果概况如表7所示,表7为贝叶斯分类器和决策树分类器预测准确率对比。
分类器 | 准确率 | 跨等级预测错误数据 |
贝叶斯分类器 | 77.6422764228%, | 15条 |
决策树分类器 | 73.8482384824% | 10条 |
表7
738条测试数据中预测准确573条数据,准确率达到77.6422764228%,预测错误165条数据,错误的165条数据中跨等级预测错误数据为15条,其余预测错误的150条警情等级均为相邻等级。
Claims (2)
1.一种城市交通警情等级预测方法,其特征在于:所述预测方法包括如下步骤:
步骤1:时间数据整理
时间数据包括5维特征属性,分别是第几单位时间、星期几、是否节假日、是否工作日和是否有重大活动,所述第几单位时间属性处理方法为设定单位时间Δt=24h/N,按单位时间Δt将全天24小时分为N段;获取国务院颁布的法定节假日信息构成节假日属性;依据节假日和周末信息构成工作日属性;根据当天是否有大型考试、当前辖区内学校是否放假、是否有大型场馆活动定义是否有重大活动属性;
步骤2:天气数据与交通警情关联性分析及处理
天气数据包括温度、可见度、降雨量3个特征属性;
步骤3:建立预测交通警情等级的属性集
根据数据库信息,将处理后的时间数据、天气数据、重大活动日、中小学放假特征属性整合成多维数据宽表,其中包括8个离散的特征属性维度;
步骤4:辖区内交通警情等级划分
统计各辖区近期6个月内的警情发生次数,采用动态划分的方法将交通警情发生次数划分为三个等级,将各个交警支队在单位时间Δt内的警情数量采用等区间装箱方法求和,构成(x1,x2,…,xn)关于交通警情数量的一个样本空间,由n个样本点组成,其中可计算出当前交警支队的警情数量均值记为μ,计算出警情数量的标准差记为σ,故警情数量属于[0,μ-σ)的标记警情等级为C1,警情数量属于[μ-σ,μ+σ)的标记警情等级为C2,警情数量属于[μ+σ,+∞)的标记警情等级为C3;
步骤5:构造朴素贝叶斯分类器
选取第几三小时X1、星期几X2、是否节假日X3、是否工作日X4、是否有重大活动X5、温度X6、可见度X7、降雨量X8构成预测交通警情等级的多维数据集T={X1,X2,X3,X4,X5,X6,X7,X8},将标记好的警情等级与多维数据集共同构成朴素贝叶斯分类的训练集R={X1,X2,X3,X4,X5,X6,X7,X8,Ci},其中Ci∈(C1,C2,C3),从而使每个训练样本都具有相应的警情等级类别标签。已知未来几天各个特征属性Ti={X1,X2,X3,X4,X5,X6,X7,X8}=d预测当前区域的交通警情等级,根据朴素贝叶斯分类器的原理可知需要比较在Ti=d的条件下属于C1、C2、C3三类警情等级的概率大小,即需要计算p(C1|Ti=d)、p(C2|Ti=d)、p(C3|Ti=d)三者的概率。由贝叶斯公式可知其中Ci∈(C1,C2,C3),比较p(C1|Ti=d)、p(C2|Ti=d)、p(C3|Ti=d)三者的概率大小只需要比较p(C1)p(Ti=d|C1)、p(C2)p(Ti=d|C2)、p(C3)p(Ti=d|C3)三者的大小,分母p(Ti=d)对所有类别为常数;
通过统计训练集R中属于警情等级C1的样本个数可计算出p(C1)的概率,即 同理可计算出 p(Ti=d|Ci)即在类别Ci中属性Ti的取值为d的概率,具体p(Ti=d|Ci)=p(X1=d1|Ci)p(X2=d2|Ci)…p(X8=d8|Ci),其中d=(d1,d2,…,d8),8个离散特征属性的先验概率p(X1=d1|Ci)、…、p(X8=d8|Ci)通过计算在每个类别Ci中每个属性Xi取值为di的概率p(Xi=di|Ci)求得;
步骤6:预测各区域的交通警情等级
由步骤5形成各特征属性的概率表,将所有特征属性的概率相乘计算出p(C1|Ti)、p(C2|Ti)、p(C3|Ti)三者的概率大小,将所有特征属性的概率相乘时如果有某个特征属性的概率值为零,会消除乘积中涉及的所有其他属性概率的影响,导致最终p(Ci|Ti)的结果为零,每个特征属性从1开始计数,p(C1|Ti)、p(C2|Ti)、p(C3|Ti)三者中最大的概率值即为朴素贝叶斯分类器预测出的当前特征属性条件Ti下的交通警情等级。
2.如权利要求1所述的一种城市交通警情等级预测方法,其特征在于:所述步骤2中,将降雨量分为以下五个等级,降雨量0.1mm-9.9mm时为小雨,降雨量为10mm-24.9mm为中雨,降雨量25mm-49.9mm为大雨,降雨量50mm-99.9mm为暴雨,降雨量100mm以上为大暴雨;将可见度等级分为以下五类,可见度大于10km时为无雾,可见度为1km-10km时为轻雾,可见度为200m-1km时大雾,可见度为50m-200m时为浓雾,可见度小于50m时为强雾;将温度划分为小于5℃,5℃-28℃,29℃-35℃,大于35℃四个等级。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510242669.6A CN104992557B (zh) | 2015-05-13 | 2015-05-13 | 一种城市交通警情等级预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510242669.6A CN104992557B (zh) | 2015-05-13 | 2015-05-13 | 一种城市交通警情等级预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104992557A true CN104992557A (zh) | 2015-10-21 |
CN104992557B CN104992557B (zh) | 2017-09-29 |
Family
ID=54304365
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510242669.6A Active CN104992557B (zh) | 2015-05-13 | 2015-05-13 | 一种城市交通警情等级预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104992557B (zh) |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106710209A (zh) * | 2016-11-09 | 2017-05-24 | 东南大学 | 一种大型活动期间交通意外事件的检测确认方法及系统 |
CN106845723A (zh) * | 2017-02-06 | 2017-06-13 | 安徽新华博信息技术股份有限公司 | 一种刑事案件的发生的预测方法 |
CN107103755A (zh) * | 2017-05-11 | 2017-08-29 | 厦门卫星定位应用股份有限公司 | 一种道路交通警情预测方法 |
CN107665577A (zh) * | 2017-09-14 | 2018-02-06 | 青岛海信网络科技股份有限公司 | 一种交通警情处理方法和装置 |
CN107783959A (zh) * | 2017-09-02 | 2018-03-09 | 南京中孚信息技术有限公司 | 一种基于贝叶斯预测的处警、接警信息评分方法 |
CN108417033A (zh) * | 2018-03-23 | 2018-08-17 | 四川高路交通信息工程有限公司 | 基于多维因素的高速路交通事故分析预测方法 |
CN108763242A (zh) * | 2018-03-26 | 2018-11-06 | 广州视源电子科技股份有限公司 | 标签生成方法及装置 |
CN109360421A (zh) * | 2018-11-28 | 2019-02-19 | 平安科技(深圳)有限公司 | 一种基于机器学习的交通信息预测方法及装置、电子终端 |
CN109767597A (zh) * | 2019-01-19 | 2019-05-17 | 跨越速运集团有限公司 | 一种车辆事故预警方法及系统 |
CN110197582A (zh) * | 2018-02-27 | 2019-09-03 | 沈阳美行科技有限公司 | 数据分析方法、装置以及交通预警方法、装置 |
CN110689131A (zh) * | 2019-09-26 | 2020-01-14 | 长安大学 | 一种基于朴素贝叶斯模型的车辆能耗影响分析方法 |
CN112019819A (zh) * | 2020-09-21 | 2020-12-01 | 富盛科技股份有限公司 | 一种适用于公安指挥中心监测并解决突发警情的方法 |
CN113570846A (zh) * | 2021-06-08 | 2021-10-29 | 北京交通大学 | 一种交通警情态势分析研判方法、设备及可读存储介质 |
CN113888850A (zh) * | 2021-12-03 | 2022-01-04 | 南京坤拓土木工程科技有限公司 | 沙尘气象市政道路拥堵态势多因素预警方法 |
CN115359657A (zh) * | 2022-08-16 | 2022-11-18 | 青岛海信网络科技股份有限公司 | 一种交通管理方法、装置、设备及介质 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012065196A (ja) * | 2010-09-16 | 2012-03-29 | Toshiba Corp | サービス情報配信システム及びサービス情報配信方法 |
KR20120115187A (ko) * | 2012-09-03 | 2012-10-17 | 여지홍 | 반(1/2) 입체 교차로를 이용한 논스톱 교통 시스템 및 이를 적용한 교통 시스템 |
CN103345208A (zh) * | 2013-06-04 | 2013-10-09 | 北京北排水务设计研究院有限责任公司 | 城市雨洪防治智能控制系统及控制方法 |
CN103578228A (zh) * | 2012-08-10 | 2014-02-12 | 陕西银河景天电子有限责任公司 | 一种暴雨洪涝灾害监测预警及应急响应系统与方法 |
CN104182618A (zh) * | 2014-08-06 | 2014-12-03 | 西安电子科技大学 | 一种基于贝叶斯网络的追尾预警方法 |
CN104268710A (zh) * | 2014-10-10 | 2015-01-07 | 北京交通大学 | 轨道交通路网动态应急处置方案生成方法 |
-
2015
- 2015-05-13 CN CN201510242669.6A patent/CN104992557B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012065196A (ja) * | 2010-09-16 | 2012-03-29 | Toshiba Corp | サービス情報配信システム及びサービス情報配信方法 |
CN103578228A (zh) * | 2012-08-10 | 2014-02-12 | 陕西银河景天电子有限责任公司 | 一种暴雨洪涝灾害监测预警及应急响应系统与方法 |
KR20120115187A (ko) * | 2012-09-03 | 2012-10-17 | 여지홍 | 반(1/2) 입체 교차로를 이용한 논스톱 교통 시스템 및 이를 적용한 교통 시스템 |
CN103345208A (zh) * | 2013-06-04 | 2013-10-09 | 北京北排水务设计研究院有限责任公司 | 城市雨洪防治智能控制系统及控制方法 |
CN104182618A (zh) * | 2014-08-06 | 2014-12-03 | 西安电子科技大学 | 一种基于贝叶斯网络的追尾预警方法 |
CN104268710A (zh) * | 2014-10-10 | 2015-01-07 | 北京交通大学 | 轨道交通路网动态应急处置方案生成方法 |
Non-Patent Citations (2)
Title |
---|
郑征 等: "基于时空数据应用模型的交通警情研判实现", 《警察技术》 * |
高爱霞: "基于可靠性分析的城市交通预警管理研究", 《北京人民警察学院学报》 * |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106710209A (zh) * | 2016-11-09 | 2017-05-24 | 东南大学 | 一种大型活动期间交通意外事件的检测确认方法及系统 |
CN106710209B (zh) * | 2016-11-09 | 2019-02-01 | 东南大学 | 一种大型活动期间交通意外事件的检测确认方法及系统 |
CN106845723A (zh) * | 2017-02-06 | 2017-06-13 | 安徽新华博信息技术股份有限公司 | 一种刑事案件的发生的预测方法 |
CN106845723B (zh) * | 2017-02-06 | 2020-12-29 | 安徽新华博信息技术股份有限公司 | 一种刑事案件的发生的预测方法 |
CN107103755A (zh) * | 2017-05-11 | 2017-08-29 | 厦门卫星定位应用股份有限公司 | 一种道路交通警情预测方法 |
CN107103755B (zh) * | 2017-05-11 | 2019-12-20 | 厦门卫星定位应用股份有限公司 | 一种道路交通警情预测方法 |
CN107783959A (zh) * | 2017-09-02 | 2018-03-09 | 南京中孚信息技术有限公司 | 一种基于贝叶斯预测的处警、接警信息评分方法 |
CN107665577A (zh) * | 2017-09-14 | 2018-02-06 | 青岛海信网络科技股份有限公司 | 一种交通警情处理方法和装置 |
CN107665577B (zh) * | 2017-09-14 | 2020-02-07 | 青岛海信网络科技股份有限公司 | 一种交通警情处理方法和装置 |
CN110197582A (zh) * | 2018-02-27 | 2019-09-03 | 沈阳美行科技有限公司 | 数据分析方法、装置以及交通预警方法、装置 |
CN108417033A (zh) * | 2018-03-23 | 2018-08-17 | 四川高路交通信息工程有限公司 | 基于多维因素的高速路交通事故分析预测方法 |
CN108763242A (zh) * | 2018-03-26 | 2018-11-06 | 广州视源电子科技股份有限公司 | 标签生成方法及装置 |
CN108763242B (zh) * | 2018-03-26 | 2022-03-08 | 广州视源电子科技股份有限公司 | 标签生成方法及装置 |
CN109360421A (zh) * | 2018-11-28 | 2019-02-19 | 平安科技(深圳)有限公司 | 一种基于机器学习的交通信息预测方法及装置、电子终端 |
CN109767597A (zh) * | 2019-01-19 | 2019-05-17 | 跨越速运集团有限公司 | 一种车辆事故预警方法及系统 |
CN110689131A (zh) * | 2019-09-26 | 2020-01-14 | 长安大学 | 一种基于朴素贝叶斯模型的车辆能耗影响分析方法 |
CN112019819A (zh) * | 2020-09-21 | 2020-12-01 | 富盛科技股份有限公司 | 一种适用于公安指挥中心监测并解决突发警情的方法 |
CN113570846A (zh) * | 2021-06-08 | 2021-10-29 | 北京交通大学 | 一种交通警情态势分析研判方法、设备及可读存储介质 |
CN113570846B (zh) * | 2021-06-08 | 2022-11-04 | 北京交通大学 | 一种交通警情态势分析研判方法、设备及可读存储介质 |
CN113888850A (zh) * | 2021-12-03 | 2022-01-04 | 南京坤拓土木工程科技有限公司 | 沙尘气象市政道路拥堵态势多因素预警方法 |
CN115359657A (zh) * | 2022-08-16 | 2022-11-18 | 青岛海信网络科技股份有限公司 | 一种交通管理方法、装置、设备及介质 |
CN115359657B (zh) * | 2022-08-16 | 2023-10-13 | 青岛海信网络科技股份有限公司 | 一种交通管理方法、装置、设备及介质 |
Also Published As
Publication number | Publication date |
---|---|
CN104992557B (zh) | 2017-09-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104992557A (zh) | 一种城市交通警情等级预测方法 | |
Liu et al. | Modelling the multilevel structure and mixed effects of the factors influencing the energy consumption of electric vehicles | |
CN102819955B (zh) | 基于车辆行程数据的道路网运行评价方法 | |
Ferreira et al. | The impact of driving styles on fuel consumption: A data-warehouse-and-data-mining-based discovery process | |
Jones et al. | Geographical variations in mortality and morbidity from road traffic accidents in England and Wales | |
CN109086246A (zh) | 一种道路交通源的排放量计算方法及装置 | |
CN108198425A (zh) | 一种电动汽车行驶工况的构建方法 | |
CN102592447B (zh) | 一种基于fcm的区域路网的道路交通状态判别方法 | |
CN102737504B (zh) | 一种基于驾驶特性的公交车到站时间实时估计方法 | |
CN103956050A (zh) | 基于车辆行程数据的道路网运行评价方法 | |
CN105825310A (zh) | 基于信息熵的出租车寻客路线推荐方法 | |
CN102044153B (zh) | 一种交通流通畅等级评价的方法与系统 | |
Ma et al. | Analysis of the logistic model for accident severity on urban road environment | |
CN105469602A (zh) | 一种基于ic卡数据的公交乘客候车时间范围的预测方法 | |
Ye et al. | Modeling level-of-safety for bus stops in China | |
CN113920732B (zh) | 一种针对特定驾驶人群的道路交通事故风险预警方法 | |
Haleem et al. | Identifying traditional and nontraditional predictors of crash injury severity on major urban roadways | |
CN104091440B (zh) | 基于道路等级和车辆类型的城市交通拥堵成本估算方法 | |
CN110867075B (zh) | 一种评估雨天条件下道路测速仪对驾驶人员反应行为影响的方法 | |
Ritchie | A statistical approach to statewide traffic counting | |
Kassu et al. | Identifying the principal factors influencing traffic safety on interstate highways | |
CN105608315A (zh) | 公交停靠站安全服务水平定量评价方法及系统 | |
CN114565155A (zh) | 一种公路安全测试方法、装置、电子设备及存储介质 | |
CN104118328A (zh) | 基于综合变量的汽车能耗率估计设备和方法 | |
CN113345252A (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 | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20231109 Address after: Room 301, 3rd Floor, Building 7, No. 2 Xiyuanba Road, Xihu District, Hangzhou City, Zhejiang Province, 310012 Patentee after: Chengbao (Hangzhou) Internet Technology Co.,Ltd. Address before: 310012 Room 101, floor 1, building 1, No. 223, Yile Road, Xihu District, Hangzhou City, Zhejiang Province Patentee before: ZHEJIANG ENJOYOR RESEARCH INSTITUTE CO.,LTD. |