Summary of the invention
The objective of the invention is to overcome the deficiency that existing automatic accident detection algorithm exists, propose a kind of city expressway automatic detection method for traffic accident based on 1 minute occupation rate sequence of discontinuous slip, its verification and measurement ratio, detection time and rate of false alarm all are better than existing algorithm.
For reaching above purpose, solution of the present invention is:
A kind of city expressway automatic detection method for traffic accident based on non-continuous sliding sequence, it may further comprise the steps:
1) detection system is gathered the telecommunication flow information of a fixed step size in real time by earth coil;
2) judge according to the real-time traffic stream information of gathering whether have traffic hazard take place, if judgement has an accident then give the alarm, notice vehicle supervision department takes corresponding measure if detecting the highway section; Otherwise continue to read real time data, judge next time;
The collection of described real-time traffic stream information generates 1 minute occupation rate data sequence of discontinuous slip as the basic data of differentiating algorithm.
Structure serves as the 1 minute occupation rate sequence of discontinuous slip at interval of sliding with minimum period 20s, and its construction method is:
else,
Wherein: O
i---i value of 1 minute occupation rate sequence of discontinuous slip;
o
i 1m---i value of average 1 minute occupation rate sequence;
o
i Mov1m---with 20s is at interval, i value of the 1 minute occupation rate sequence of sliding;
o
i 20s---i value of 20s occupation rate sequence;
o
i Mov5m---with 20s is at interval, i value of the 5 minutes occupation rate sequences of sliding.
Judge the method whether highway section has an accident that detects:
1) handles the detection road section traffic volume stream information of gathering in real time, generate 1 minute occupation rate data sequence of discontinuous slip;
2) determine the threshold value of judge index;
Whether 3) the correlated judgment index of calculating 1 minute occupation rate data sequence of discontinuous slip is compared with decision threshold, judge to have an accident.
The correlated judgment index of 1 minute occupation rate data sequence of described discontinuous slip is respectively occupation rate O
T, i, 1 minute occupation rate of i section t discontinuous slip constantly; Occupation rate time difference d
T, i t, i section t constantly with the difference of 1 minute occupation rate of t-1 discontinuous slip constantly; Occupation rate space difference d
T, i s, the difference of t moment upstream i section and 1 minute occupation rate of the discontinuous slip of adjacent downstream i+1 section; Occupation rate space time difference d
T, i s→ t, t be occupation rate space difference and t-1 moment occupation rate space difference poor constantly.
Threshold value determination method is: based on the off-line calibration of historical casualty data or based on predetermined threshold value and in one day accident and non-casualty data online review one's lessons by oneself just definite.
Described occupation rate O
T, i, occupation rate time difference d
T, i t, occupation rate space difference d
T, i sWith occupation rate space time difference d
T, i S → tWhen all being not less than threshold value, the judgement accident takes place.
According to the needs that accident detects automatically, the basic data that is used for the judgement accident should satisfy following two basic demands: under accident situation, have good mutability, be convenient to carry out accident identification; Under non-accident situation, have good stationarity, reduce the false alarm under the normal traffic situation.Adopt 1 minute occupation rate sequence of discontinuous slip to carry out data processing, can realize this two targets effectively.
Owing to adopted such scheme, the present invention to have following characteristics: accident verification and measurement ratio height of the present invention, false alarm rate is low, and detection time is short, and index for selection has clear and definite meaning, is fit to Project Realization.
Embodiment
The present invention is further illustrated below in conjunction with the accompanying drawing illustrated embodiment.
Studies show that, for 20 seconds be the time average occupation rate data of step-length, its undulatory property is bigger, is difficult to directly apply to automatic accident detection algorithm.If employing was the time average occupation rate data sequence of step-length with 1 minute, have mutability under the best accident, but sense cycle was extended for 1 minute from 20 seconds; And 20 seconds serving as that the 1 minute occupation rate sequence of slip at interval of sliding has stationarity under the best non-accident, but the sudden change under the accident is obvious inadequately.The minimum period 20s that 1 minute occupation rate sequence of discontinuous slip is obtained with data slides at interval, combines the advantage of above two kinds of data sequences, and its construction method is as follows:
else,
In the formula: O
i---i value of 1 minute occupation rate sequence of discontinuous slip;
o
i 1m---i value of average 1 minute occupation rate sequence;
o
i Mov1m---with 20s is at interval, i value of the 1 minute occupation rate sequence of sliding;
o
i 20s---i value of 20s occupation rate sequence.
o
i Mov5m---with 20s is at interval, i value of the 5 minutes occupation rate sequences of sliding;
Before algorithm begins to make up 1 minute occupation rate sequence of continuous slip, should at least have 5 minutes 20s at interval raw data exist.Be exemplified below:
20s original occupation rate data sequence at interval:
18 19 17 18 19 20 21 19 19 20 21 19 21 17 22 24 21 28 35 42 43 45 48 41 36 20 17 19 18 20
1 minute occupation rate sequence (since the 5th minute end) of the corresponding discontinuous slip that makes up:
20.0 20.0 20.0 24.3 24.3 40.0 40.0 40.0 44.7 44.7 44.7 44.7 44.7 24.3 24.3
With the algorithm that 1 minute occupation rate data sequence of discontinuous slip is a fundamental construction, it is consistent with the shortest interval 20s of raw data that sense cycle can reach, and guaranteed the higher detection frequency, helps the timely discovery of accident.
As the basic data that the judgement accident takes place, determine four indexs of the automatic accident detection that this algorithm adopts with 1 minute occupation rate data sequence of discontinuous slip, adopt " serial judgement " to carry out accident and judge as many indexs accident associating determination methods:
(1) occupation rate (O
T, i), 1 minute occupation rate of i section t discontinuous slip constantly.
(2) occupation rate time difference (d
T, i t), i section t constantly with the difference of 1 minute occupation rate of (t-1) discontinuous slip constantly.
(3) occupation rate space difference (d
T, i s), the difference of t moment upstream i section and 1 minute occupation rate of the discontinuous slip of adjacent downstream i+1 section.
(4) occupation rate space time difference (d
T, i S → t), t is occupation rate space difference and t-1 moment occupation rate space difference poor constantly.
Order is judged above-mentioned 4 indexs, is judged as accident condition when satisfying condition simultaneously, and wherein any one index is ineligible then jumps out judgement, thinks that accident free takes place.
Implementation method of the present invention is as follows, with reference to shown in Figure 3:
1. arrange checkout equipment (earth coil) in city expressway surveyed area upstream and downstream, gather occupation rate parameter averaging time with certain sample unit (as 20 seconds), as the basic data of algorithm.
2. determine the threshold value of judge index.To some sections, the elementary section in the identical track in same detection highway section, adopt the same threshold group, promptly same threshold vector.To determining of metrics-thresholds, the present invention proposes the off-line calibration method of two kinds of method: A. based on historical casualty data, and B. is based on predetermined threshold value and accident and non-casualty data online from modification method in a day.
A. based on the off-line calibration method of historical casualty data
At first, delimit highway section and track to be calibrated,, transfer the data in respective stretch track in the historical data base according to the accident record of confirming with artificial supervision, telephone call or other modes;
Subsequently, calculate the maximal value that obtains each key parameter.Adopt limited exhaustive thought, avoid doing too much unnecessary computing.Promptly earlier determine the scope that it is exhaustive, needn't carry out computing surpassing its peaked scope according to each key parameter maximal value; According to its situation of change each key parameter is determined corresponding step-size in search respectively again, thereby obtain the threshold vector of some various combinations.
At last, these threshold vectors are moved detection algorithm one by one, obtain the detection effect of each threshold value grouping, choose and detect one group of best threshold value of effect as recommending threshold value to historical accident.
Be exemplified below table:
Table 1 is based on historical data off-line calibration key parameter
Parameter |
Occupation rate (%) |
Occupation rate time difference (%) |
Occupation rate space difference (%) |
Occupation rate space time difference (%) |
Maximal value |
71 |
44 |
44 |
34 |
Exhaustive scope |
0-75 |
0-45 |
0-45 |
0-35 |
Step-size in search |
5 |
5 |
5 |
5 |
The final threshold value of determining |
30 |
20 |
25 |
5 |
B. based on predetermined threshold value and accident and non-casualty data online in a day, judge relatively by three grades mainly to constitute that the algorithm detailed process is with reference to figure 6 from modification method.
At first, according to historical data analysis,, obtain predetermined threshold value vector T=[T for the default one group of threshold value in each track, highway section
1, T
2, T
3, T
4], T wherein
iIt is the predetermined threshold value of i index; Carry out one-day operation according to T,, acquire the set of bare maximum separately of following four indexs of normal traffic situation in a day 24 hours by confirming to detecting actual the having an accident in highway section
And the set of the relative minimum of following four indexs of accident situation, promptly can distinguish the minimum value set of accident condition and normal condition
Wherein
Max
j(t
i) represent that j plays accident i key parameter t takes place in back 3 minutes
iObtained maximal value.
Subsequently, with T and t
Max, t
MinCarry out vector respectively relatively.Wherein i, j, m, n, k, l ∈ 1,2,3,4}:
(1) when the predetermined threshold value of all indexs
The time, show that then preset threshold value might be judged as accident with normal condition:
1. if t
i MinBe sky, then accident free takes place in 24 hours same day.
ε=1;
2. if
T then
iBe provided with still higherly, can not detect whole accidents, need turn down to T '=t
Min
3. if
T then
iCan detect whole accidents, at this moment:
If
Show that then the i index occurs all having under each accident greater than the situation of mxm. under the normal condition, can heighten T
iExtremely
ε=1;
If
Show that then the i index just equals maximal value under the normal condition at the relative minimum under each accident, get this moment
If
Show that then relative minimum under this emergency conditions less than the maximal value under the normal condition, will cause the appearance of false alarm.Get this moment
To guarantee that whole accidents can detect;
If t
MaxIn the k index
The l index
Then get respectively
Promptly get T '=t
Min
4. if t
MinIn the m index
Refer to the n index
Then order respectively
Promptly get T '=t
Min
(2) when the predetermined threshold value of all indexs
The time, show that predetermined threshold value can not be judged as accident with normal condition:
1. if t
i MinBe sky, then accident free takes place in 24 hours same day.Get T '
i=T
i, i.e. T '=T;
2. if
Then all accidents can both be detected, and get T '
i=T
i, i.e. T '=T;
3. if t
MinIn the m index
And the n index
Then order respectively
Promptly get T '=t
Min
4. if
Then there is accident to be failed to report.At this moment:
If
Show that then the i index occurs all having under each accident greater than the situation of mxm. under the normal condition, can heighten T
iExtremely
ε=1;
If
Show that then the i index just equals maximal value under the normal condition at the relative minimum under each accident, get this moment
If
Show that then relative minimum under this emergency conditions less than the maximal value under the normal condition, will cause the appearance of false alarm.Reduce T this moment
iValue, order
To guarantee that whole accidents can detect;
If t
MaxIn the k index
The l index
Then get respectively
Promptly get T '=t
Min
(3) when
Then show part predetermined threshold value (T
i) normal condition might be judged as accident:
1. if t
MinBe sky, then accident free takes place in 24 hours same day.Get T '
k=T
k, i.e. T '=T;
2. if
Then all accidents can both be detected, and get T ' this moment
k=T
k, i.e. T '=T;
3. if
Then there is accident to be failed to report.This season
To guarantee that whole accidents can detect;
4. if t
MinIn the m index
And the n index
Then the m index can detect all accidents, and the n index has accident and fails to report.At this moment:
If
Then can improve
ε=1;
If
Then
Be T '=t
Max
If
Then have normal condition by reported by mistake for accident may, be detected for guaranteeing whole accidents this moment, gets
Be T '=t
Min
If t
MaxIn the k index
The l index
Then get respectively
Promptly get T '=t
Min
At last, with adjusted metrics-thresholds vector T '=[T '
1, T '
2, T '
3, T '
4] as new metrics-thresholds T, the accident that is used for next day is differentiated.
3. calculate four index O that the serial of 1 minute occupation rate data sequence of discontinuous slip is judged
T, i, d
T, i t, d
T, i sAnd d
T, i S → t
4. order is judged successively, relatively O
T, iWith T
1, d
T, i tWith T
2, d
T, i sWith T
3And d
T, i S → tWith T
4If all indexs have all been passed through serial arithmetic, satisfy simultaneously: O
T, i〉=T
1,
With
Then thinking has accident to take place, and sends warning.
Sample calculation analysis
The observation accident of Shanghai City city expressway is chosen in this experiment, and many indexs accident detection method that this research is set up and existing basic California algorithm and California 7# algorithm carry out parallel testing.Data are the Shanghai Urban through street that derives from through street, Shanghai supervisory system, are the coil section occupation rate averaging time data of step-length with 20 seconds.1 minute occupation rate data sequence of discontinuous slip algorithm adopts preceding 1 heavenly calendar historical event in advance so data are carried out the threshold value demarcation.
Test result is as shown in table 2.
Discontinuous slip 1 minute occupation rate data sequence algorithm of table 2 and the contrast of California test of heuristics
(detection zone 1 day data comprises 3 accidents)
Testing result shows:
(1) do not report by mistake in 13477 times are detected, rate of false alarm is 0; 2 wrong reports have taken place in basic California algorithm in 7200 times are detected, rate of false alarm is 0.03%; 9 wrong reports have taken place in California 7# algorithm in 7200 times are detected, rate of false alarm is 0.12%.The many indexs detection method that shows this method is better than basic California algorithm and California 7# algorithm on the performance of avoiding reporting by mistake;
(2) this algorithm all can effectively detect 3 test accidents with basic California algorithm, and California 7# algorithm can only detect 2 test accidents;
(3) this algorithm detection time and two kinds of California algorithms are basic identical, all find constantly early than artificial accident.
The above-mentioned description to embodiment is can understand and apply the invention for ease of those skilled in the art.The person skilled in the art obviously can easily make various modifications to these embodiment, and needn't pass through performing creative labour being applied in the General Principle of this explanation among other embodiment.Therefore, the invention is not restricted to the embodiment here, those skilled in the art should be within protection scope of the present invention for improvement and modification that the present invention makes according to announcement of the present invention.