CN103646533B - Traffic accident modeling and control method based on sparse multi-output regression - Google Patents

Traffic accident modeling and control method based on sparse multi-output regression Download PDF

Info

Publication number
CN103646533B
CN103646533B CN201310596073.7A CN201310596073A CN103646533B CN 103646533 B CN103646533 B CN 103646533B CN 201310596073 A CN201310596073 A CN 201310596073A CN 103646533 B CN103646533 B CN 103646533B
Authority
CN
China
Prior art keywords
traffic accident
accident
traffic
output regression
modeling
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.)
Expired - Fee Related
Application number
CN201310596073.7A
Other languages
Chinese (zh)
Other versions
CN103646533A (en
Inventor
陈小波
梁军
江浩斌
陈龙
张飞云
肖艳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jiangsu University
Original Assignee
Jiangsu University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Jiangsu University filed Critical Jiangsu University
Priority to CN201310596073.7A priority Critical patent/CN103646533B/en
Publication of CN103646533A publication Critical patent/CN103646533A/en
Application granted granted Critical
Publication of CN103646533B publication Critical patent/CN103646533B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Traffic Control Systems (AREA)

Abstract

The present invention discloses a kind of traffic accident modeling and control method based on sparse multi-output regression, can be used to the situation of certain selected area generation traffic accident to set up model. Multi-output regression hazard model is set up in certain area, can automatically detect the characteristic parameter irrelevant with accident and distinguish the impact of parameters, bring the characteristic parameter relevant to traffic into model, can also predict the following situation that traffic accident occurs in this area. The present invention utilizes the history and current situation data of traffic accident, the order of severity of traffic accident is carried out to modeling, the characteristic parameter that energy automatic rejection and accident are irrelevant, distinguishes the key factor that affects traffic safety, takes measure more targetedly while contributing to roading. This method is being predicted traffic accident, is analyzing cause of accident, is being instructed road construction improvement to have actual engineering using value.

Description

Traffic accident modeling and control method based on sparse multi-output regression
Technical field
The invention belongs to traffic intelligent management and control technology field, use the sparse multi-output regression model proposing, set up a kind of roadRoad traffic accident method, and can automatically detect the key factor that affects traffic safety, in prediction traffic accident, analysis thingTherefore reason, instruct road construction to improve to there is actual engineering using value.
Background technology
Along with the develop rapidly of China's road traffic cause, traffic accident is surged and is become the serious problem of society today receiving much concern nowOne of. Just because of the serious consequence that traffic accident causes, a large amount of human and material resources and wealth have poured into accident prevention and countermeasure in ChinaPower, formulated comparatively perfect control of traffic and road law, regulation and relevant policies, but the road traffic accident of China causesFactual Damage but still serious, and the trend in rising always, annual toll on traffic occupy first place, the world.Therefore, in order to ensure economic development and social stability, ensure people's safety and cargo transportation security, to Chinese traffic peaceHolotype gesture is made the prediction of science, becomes very necessary and extremely urgent.
At present, there is several different methods to be applied to accident modeling both at home and abroad. Existing research has been set up accident and has been occurred and traffic stream characteristics,Road environment condition, the size of population, the relation between the data such as vehicle guaranteeding organic quantity, but great majority are all to enter from macroscopical angleRow is analyzed, and the analysis result drawing is the prediction of the traffic accident number of times in a region, and this rises for prevention and minimizing traffic accidentLess than too large effect. Although also there is part Study impact on traffic accident from the factor of the angle analysis road of microcosmic,,Ignored that to affect the factor of traffic accident a lot, and some factors wherein exist as NF, if do not reject thisA little characteristic parameters and carry out traffic accident modeling, this may hide real feature and the inherent law of traffic accident.
Traffic accident is the complicated system being produced by many factors acting in conjunction, and relation between these factors is difficult to separatingThe method of analysing is described, and passes through the study to historical data based on sparse multi-output regression model, can build with degree of precisionThe model of vertical traffic accident.
For the accuracy of modeling, do not omit important information, we can consider as far as possible many factors, also must reject and thing simultaneouslyTherefore irrelevant characteristic parameter, therefore, the present invention has introduced sparse analytic approach, initial data is carried out to sparse analysis, thereby distinguishAffect the key factor of traffic safety, carry out quantitative analysis thereby effectively utilize a large amount of statistics, disclose between variableInternal relation, and the improvement that can be road construction provides targetedly and instructs.
Summary of the invention
The problem to be solved in the present invention is the above-mentioned deficiency that overcomes prior art, provides a kind of simple to operate, and accuracy is high, is suitable forThe method of traffic accident modeling in road safety checks.
Technical scheme of the present invention is:
A) obtain the historical traffic accident data sample N group { X in somewhere1,X2,...,XN, each XI(I=1,2,3,...,N)Comprise following population, vehicle, road and environmental information: total number of people x1, schooling x2, income level x3, car steeringMember's quantity x4, truck conspicuity marking quantity x5, cyclecar quantity x6, bicycle quantity x7, main section total length x8,Main road section traffic volume lamp quantity x9, average weather conditions x10
B) obtain corresponding traffic accident result sample { Y1,Y2,...,YN, each YI(I=1,2,3 ..., N) comprise accident frequencyy1, death toll y2, number of injured people y3And direct economic loss y4
C) traffic accident data pretreatment. By above-mentioned N group traffic accident data sample { X1,X2,...,XNNormalization, and leave square inIn battle array X. Above-mentioned N is organized to corresponding traffic accident result sample { Y1,Y2,...,YNLeave in matrix Y.
D) supposition XIWith YI(I=1,2,3 ..., N) meet and be related to YI=XIW+ ε, W is coefficient, ε is Gaussian noise. If λ isBalance parameters, can be taken as 0.01≤λ≤100, and W can obtain by solving following sparse multi-output regression model
min J ( W ) = Σ I = 1 N | | X I W - Y I | | 2 + λ | | W | | 21 = | | XW - Y | | 2 + λ | | W | | 21
E) choose coefficient matrix W, make J (W) reach minimum of a value, establishing D is unit matrix.
F) by J (W) to W differentiation, and make result equal 0, obtain
XT(XW-Y)+λDW=0
G) solve above formula and obtain W=(XTX+λD)-1XTY。
If D' is diagonal matrix, diagonal element is the inverse of the L2 norm of W corresponding row.
H) whether evaluation algorithm restrains, || and whether D-D'|| is less than given several ε (as ε=0.001), if so, exports W,Otherwise D' is copied in D, forward step F to).
I) identify the factor irrelevant to traffic accident modeling according to the row that in W is entirely zero, according to the positive and negative of the capable value of non-zero in W andThe impact of the different factors of size discrimination on traffic accident. As driver's quantity x4Coefficient for negative, road that should be afterwardsIn the traffic programme of road, increase driver's quantity to reduce following traffic accident.
J) according to the traffic accident characteristic parameter X'=(x in the above-mentioned area gathering later1,x2,...,x10), predict that the accident in above-mentioned area occursNumber of times, death toll, number of injured people and direct economic loss are WX'.
Technique effect of the present invention is as follows:
Based on sparse multi-output regression model. More existing research methods directly utilize initial data to instruct as input variablePractice, but these initial data are deposited problem both ways, the one, data volume is large, and dimension is many, and the one, between data, there is correlation,Information has redundancy, and this will cause the reduction of the not high and training speed of the precision of result. Initial data is carried out to sparse analysis,The characteristic parameter that automatic rejection and accident are irrelevant, carries out quantitative analysis thereby effectively utilize a large amount of statistics, disclose variable itBetween internal relation.
Choose influence factor from the angle of microcosmic and carry out modeling. Traditional Predictive Methods of Road Accidents is all to utilize some macroscopical shadowsThe factor of sound is analyzed as input variable, and what obtain is the data of an overall traffic accident. The present invention is micro-by roadSight factor, as input data, is analyzed the relation of these factors and traffic accident, thus obtain these factors to traffic accident respectivelyFrom impact, this conclusion can be applied in Highway traffic safety assessment, for design and the trimming of road and surrounding enviroment provideFoundation.
Brief description of the drawings
Fig. 1 is the traffic accident model FB(flow block) based on sparse multi-output regression of the present invention.
Detailed description of the invention
At society, traffic safety gets more and more people's extensive concerning, and research that road safety is evaluated is for traffic accidentReduce and have very large effect. The present invention is from the angle of microcosmic, and the modeling method of application based on sparse multi-output regression is to a certain areaRelation etc. multiple variable and traffic accident is analyzed and researched.
Introduce in detail workflow of the present invention below:
One, choosing of data
Road traffic is the dynamical system being made up of key elements such as people, car, road, environment. Traffic accident is in road traffic system,The event accidentally happening suddenly due to the cooperation imbalance of people, car, road, all key elements of environment. Therefore, choosing the impact of traffic accidentBecause of the period of the day from 11 p.m. to 1 a.m, analyze from above aspect.
1. human factor: according to the statistics of traffic accidents result of 2002, the death toll causing due to people's reason accounts for worked as year traffic88.98% of death by accident sum, what wherein automobile driver was main cause accounts for 78.56%, and bicycle driver account for4.20%, pedestrian and rider account for 6.22%. Visible, in general, the key of traffic accident is automobile driver.Therefore chosen a certain regional periphery total number of people and average schooling and income, four indexs of driver's quantity.
2. the factor of car: causing the second largest key element of traffic accident is vehicle. The vehicle travelling on road, existing motor vehicle,Also have bicycle, wherein motor vehicle (especially truck conspicuity marking) is the one vehicles fast, energy maximum, protectionProperty is also best. But this protective is only protected driver and occupant, therefore, with respect to bicycle and other bicycles,Motor vehicle is traffic powerhouse. But bicycle quantity increase also can to traffic accident number exert an influence, therefore chosen largeType vehicles number, cyclecar quantity, three indexs of bicycle quantity.
3. road and environmental factor: the environmental factor that affects traffic safety can be divided into natural environment and artificial environment, natural environmentMainly comprise geographical position, meteorological condition, the time etc., artificial environment comprises that situation of land use, trackside disturb, road barricadeThing etc. The present invention chosen the main section total length in above-mentioned area, main road section traffic volume lamp quantity and average weather conditions this threeIndividual index is as factor of influence.
Two, set up road accident risk model:
A) obtain the historical traffic accident data sample N group { X in somewhere1,X2,...,XN, each XI(I=1,2,3 ..., N) comprisePopulation, vehicle, road and environmental information: total number of people x1, schooling x2, income level x3, driver's quantity x4,Truck conspicuity marking quantity x5, cyclecar quantity x6, bicycle quantity x7, main section total length x8, hand in main sectionLogical lamp quantity x9, average weather conditions x10
B) obtain corresponding traffic accident result sample { Y1,Y2,...,YN, each YI(I=1,2,3 ..., N) comprise accident frequencyy1, death toll y2, number of injured people y3And direct economic loss y4
C) traffic accident data pretreatment. By above-mentioned N group traffic accident data sample { X1,X2,...,XNNormalization, and leave square inIn battle array X. Above-mentioned N is organized to corresponding traffic accident result sample { Y1,Y2,...,YNLeave in matrix Y.
D) supposition XIWith YI(I=1,2,3 ..., N) meet and be related to YI=XIW+ ε, W is coefficient, ε is Gaussian noise. If λ isBalance parameters, can be taken as 0.01≤λ≤100, and W can obtain by solving following sparse multi-output regression model
min J ( W ) = Σ I = 1 N | | X I W - Y I | | 2 + λ | | W | | 21 = | | XW - Y | | 2 + λ | | W | | 21
E) choose coefficient matrix W, make J (W) reach minimum of a value, establishing D is unit matrix.
F) by J (W) to W differentiation, and make result equal 0, obtain
XT(XW-Y)+λDW=0
G) solve above formula and obtain W=(XTX+λD)-1XTY。
If D' is diagonal matrix, diagonal element is the inverse of the L2 norm of W corresponding row.
H) whether evaluation algorithm restrains, || and whether D-D'|| is less than given several ε (as ε=0.001), if so, exports W,Otherwise D' is copied in D, forward step F to).
I) identify the factor irrelevant to traffic accident modeling according to the row that in W is entirely zero, according to the positive and negative of the capable value of non-zero in W andThe impact of the different factors of size discrimination on traffic accident. As driver's quantity x4Coefficient for negative, road that should be afterwardsIn the traffic programme of road, increase driver's quantity to reduce following traffic accident.
J) according to the traffic accident characteristic parameter X'=(x in the above-mentioned area gathering later1,x2,...,x10), predict that the accident in above-mentioned area occursNumber of times, death toll, number of injured people and direct economic loss are WX'.
Practice process is divided into model learning and model uses two processes.
Model learning: as Fig. 1, collect traffic accident data and corresponding traffic accident result data in a certain regional certain hour.The traffic accident function of setting up sparse multi-output regression according to aforementioned modeling procedure, calculates YI=XIThe sparse matrix W of W+ ε.
Model uses: identify the factor irrelevant to traffic accident modeling according to the row that in W is entirely zero, and capable according to non-zero in WThe impact of positive and negative different with the size discrimination factors of value on traffic accident. As driver's quantity x4Coefficient for negative, shouldIn traffic planninng afterwards, increase driver's quantity to reduce following traffic accident. And, according to same groundThe traffic accident feature X'=(x gathering in certain time of district1,x2,...,x10), utilize the accident frequency of W to this area, deadThe number of dying, number of injured people and direct economic loss are predicted X'W.

Claims (3)

1. the traffic accident modeling method based on sparse multi-output regression, is characterized in that, concrete steps are:
A) obtain the historical traffic accident data sample N group { X in somewhere1,X2,...,XN, each XIComprise following population, vehicle,Road and environmental information: total number of people x1, schooling x2, income level x3, driver's quantity x4, large-scale motor-drivenCar quantity x5, cyclecar quantity x6, bicycle quantity x7, main section total length x8, main road section traffic volume lamp numberAmount x9, average weather conditions x10; Wherein, I=1,2,3 ..., N;
B) obtain corresponding traffic accident result sample { Y1,Y2,...,YN, each YIComprise accident frequency y1, death tolly2, number of injured people y3And direct economic loss y4
C) traffic accident data pretreatment; By above-mentioned N group traffic accident data sample { X1,X2,...,XNNormalization, and leave square inIn battle array X; Above-mentioned N is organized to corresponding traffic accident result sample { Y1,Y2,...,YNLeave in matrix Y;
D) supposition XIWith YIMeet and be related to YI=XIW+ ε, W is coefficient matrix, ε is Gaussian noise; If λ is balance parameters,Be taken as 0.01≤λ≤100, W obtains by solving following sparse multi-output regression model
min J ( W ) = Σ I = 1 N | | X I W - Y I | | 2 + λ | | W | | 21 = | | X W - Y | | 2 + λ | | W | | 21 ;
E) choose coefficient matrix W, make J (W) reach minimum of a value, establishing D is unit matrix;
F) by J (W) to W differentiation, and make result equal 0, obtain
XT(XW-Y)+λDW=0;
G) solve above formula and obtain W=(XTX+λD)-1XTY;
If D' is diagonal matrix, diagonal element is the inverse of the L2 norm of W corresponding row;
H) whether evaluation algorithm restrains, || and whether D-D'|| is less than given several ε, if so, exports W, otherwise D' is copiedIn D, forward step F to);
I) identify the factor irrelevant to traffic accident modeling according to the row that in W is entirely zero, according to the positive and negative of the capable value of non-zero in W andThe impact of the different factors of size discrimination on traffic accident;
J) according to the traffic accident characteristic parameter X'=(x in the above-mentioned area gathering later1,x2,…,x10), predict that the accident in above-mentioned area is sent outRaw number of times, death toll, number of injured people and direct economic loss are WX'.
2. the traffic accident modeling method based on sparse multi-output regression according to claim 1, is characterized in that described stepH) in, given several ε: 0.001≤ε≤0.1.
3. the traffic accident modeling method based on sparse multi-output regression according to claim 1, is characterized in that described stepRapid I) in, the method for the impact according to the different factors of positive and negative and size discrimination of the capable value of non-zero in W on traffic accident is: ifThe value that in W, non-zero is capable is for just, and this correspondence factor has negative effect to accident forecast, should in the time of roading, consider to reduce thisFactor; If the value that in W, non-zero is capable is for negative, this correspondence factor has positive influences to accident forecast, should in the time of roading, examineConsider this factor that increases; As driver's quantity x4Coefficient for negative, should in traffic planninng afterwards, increase vapourCar driver's quantity is to reduce following traffic accident.
CN201310596073.7A 2013-11-22 2013-11-22 Traffic accident modeling and control method based on sparse multi-output regression Expired - Fee Related CN103646533B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310596073.7A CN103646533B (en) 2013-11-22 2013-11-22 Traffic accident modeling and control method based on sparse multi-output regression

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310596073.7A CN103646533B (en) 2013-11-22 2013-11-22 Traffic accident modeling and control method based on sparse multi-output regression

Publications (2)

Publication Number Publication Date
CN103646533A CN103646533A (en) 2014-03-19
CN103646533B true CN103646533B (en) 2016-05-25

Family

ID=50251740

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310596073.7A Expired - Fee Related CN103646533B (en) 2013-11-22 2013-11-22 Traffic accident modeling and control method based on sparse multi-output regression

Country Status (1)

Country Link
CN (1) CN103646533B (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106816020B (en) * 2015-12-02 2020-06-30 大陆投资(中国)有限公司 Traffic accident information processing method based on data analysis
CN106485919B (en) * 2016-10-11 2019-06-18 东南大学 A method of determining that through street fixed point tachymeter influences traffic accident quantity
CN106530712B (en) * 2016-12-16 2019-09-13 上海斐讯数据通信技术有限公司 The system and method for vehicular casualty rate is estimated by big data
CN108776650A (en) * 2018-06-05 2018-11-09 浙江大学 A kind of non-motorized lane service level evaluation method based on mixed Gauss model
CN109087505B (en) * 2018-07-11 2020-05-22 北京交通大学 Traffic safety regulation effectiveness analysis method for motorcycle
CN113887604A (en) * 2021-09-27 2022-01-04 同济大学 Key management area identification method for safety education improvement

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004126757A (en) * 2002-09-30 2004-04-22 Toshiba Corp Data predicting device, data predicting method and data prediction program
DE102006044086A1 (en) * 2006-09-20 2008-04-10 Audi Ag Traffic situation e.g. accident-critical driving situation, simulating system, has reference vehicle, where behavior of vehicle temporally takes place based on information detected by detection unit during simulated traffic situation
CN102103715A (en) * 2010-11-18 2011-06-22 上海海事大学 Negative binomial regression-based maritime traffic accident investigation analysis and prediction method
CN102306450A (en) * 2011-08-30 2012-01-04 同济大学 Layout method for traffic detectors of sparse road network
CN102360525A (en) * 2011-09-28 2012-02-22 东南大学 Discriminant analysis-based high road real-time traffic accident risk forecasting method
CN102855395A (en) * 2012-08-16 2013-01-02 长安大学 Method for distinguishing road black spot
CN103226818A (en) * 2013-04-25 2013-07-31 武汉大学 Single-frame image super-resolution reconstruction method based on manifold regularized sparse support regression
CN103400346A (en) * 2013-07-18 2013-11-20 天津大学 Video super resolution method for self-adaption-based superpixel-oriented autoregression model

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004126757A (en) * 2002-09-30 2004-04-22 Toshiba Corp Data predicting device, data predicting method and data prediction program
DE102006044086A1 (en) * 2006-09-20 2008-04-10 Audi Ag Traffic situation e.g. accident-critical driving situation, simulating system, has reference vehicle, where behavior of vehicle temporally takes place based on information detected by detection unit during simulated traffic situation
CN102103715A (en) * 2010-11-18 2011-06-22 上海海事大学 Negative binomial regression-based maritime traffic accident investigation analysis and prediction method
CN102306450A (en) * 2011-08-30 2012-01-04 同济大学 Layout method for traffic detectors of sparse road network
CN102360525A (en) * 2011-09-28 2012-02-22 东南大学 Discriminant analysis-based high road real-time traffic accident risk forecasting method
CN102855395A (en) * 2012-08-16 2013-01-02 长安大学 Method for distinguishing road black spot
CN103226818A (en) * 2013-04-25 2013-07-31 武汉大学 Single-frame image super-resolution reconstruction method based on manifold regularized sparse support regression
CN103400346A (en) * 2013-07-18 2013-11-20 天津大学 Video super resolution method for self-adaption-based superpixel-oriented autoregression model

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
一种基于贝叶斯网络模型的交通事故预测方法;秦小虎等;《计算机仿真》;20051130;第22卷(第11期);第230-231页 *
基于灰色系统的道路交通事故预测模型研究;韩文涛等;《西安建筑科技大学学报(自然科学版)》;20050930;第37卷(第03期);第375-377页 *
道路交通事故预测模型的构建与应用;姜华平;《山东交通学院学报》;20030331;第11卷(第01期);第35-37页 *

Also Published As

Publication number Publication date
CN103646533A (en) 2014-03-19

Similar Documents

Publication Publication Date Title
CN103646533B (en) Traffic accident modeling and control method based on sparse multi-output regression
CN108198415B (en) A kind of city expressway accident forecast method based on deep learning
CN108765235A (en) Automatic driving vehicle test scene construction method and test method based on the destructing of traffic accident case
CN103971523B (en) A kind of mountain road traffic safety dynamic early-warning system
CN103646534B (en) A kind of road real-time traffic accident risk control method
Chen et al. Motorcycle Ban and Traffic Safety: Evidence from a Quasi‐Experiment at Zhejiang, China
CN106383920B (en) A kind of particularly serious traffic accident causation recognition methods based on correlation rule
CN107845262A (en) A kind of Predictive Methods of Road Accidents and device
CN104318757B (en) Bus bus or train route section Forecasting Methodology working time on a kind of public transportation lane
CN113380033A (en) Urban traffic safety early warning method and system based on man-machine hybrid enhanced intelligence
CN108682149A (en) The linear causation analysis method in highway accident stain section based on binary logistic regression
CN110324336B (en) Internet of vehicles data situation sensing method and device based on network security
CN108460685A (en) Method and apparatus for excavating the correlation rule between vehicle insurance risks and assumptions
CN100481153C (en) Method for automatically inspecting highway traffic event based on offset minimum binary theory
CN103632535A (en) Judgment method for section pedestrian crossing signal lamp arrangement
Xu et al. Traffic accident hotspots identification based on clustering ensemble model
Dutta et al. Critical gaps at three-legged unsignalised intersections using microsimulation
Park et al. Using vehicle data as a surrogate for highway accident data
Persaud et al. Can crash modification factors be estimated from surrogate measures of safety?
CN106709595A (en) Accident delay time prediction method and system based on unformatted information
Al Jaman et al. Effects of automated traffic enforcement on driver behavior at a signalized intersection in Saudi Arabia
Deng et al. Research on the Active Safety Monitoring Evaluation System of Key Operating Vehicles
Tian et al. Understanding scenarios for cooperative V2P safety applications using connected vehicle datasets
Dong et al. Risk Analysis of Vehicle Rear‐End Collisions at Intersections
Wei et al. Predicting and factor analysis of rider injury severity in two-wheeled motorcycle and vehicle crash accidents based on an interpretable machine learning framework

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160525

Termination date: 20161122