Summary of the invention
The object of the present invention is to provide a kind of urban road traffic state method of discrimination based on section detecting device.The basic thought of the method is with the magnitude of traffic flow, speed, these three the traffic flow parameter structure evaluation index of comprehensively blocking up---traffic congestion indexes of time occupancy, utilizes traffic congestion index to differentiate section traffic behavior.For achieving the above object, the traffic state judging method that the present invention proposes comprises that information issues section traffic flow data obtains in interval level and smooth step, the speed of step traffic flow parameter block up exponential sum occupation rate block up step, critical velocity that index calculates the block up step that index calculates, the step of traffic congestion index calculating, the step of road section traffic volume condition discrimination.
The traffic state judging method that the present invention proposes, the laying situation that has comprised three kinds of section detecting devices: there is one group of section detecting device traffic state judging, have many group detecting device traffic state judgings, sensorless traffic state judging.
Having one group of section detecting device traffic state judging is a kind of urban road traffic state method of discrimination, differentiate section and only have one group of detecting device, each issue obtains traffic flow data by detecting device in interval, determine traffic congestion index, according to predefined traffic congestion index ranking interval, judge traffic behavior.
Having many group section detecting device traffic state judgings is to differentiate section to have many group detecting devices, according to the each group of detector location COMPREHENSIVE CALCULATING road section traffic volume index that blocks up, according to predefined traffic congestion index ranking interval, carries out traffic behavior judgement.
Sensorless traffic state judging is the one that urban road traffic state is differentiated, and differentiating section does not have detecting device, according to the traffic behavior of upstream and downstream line, realizes the traffic state judging in this section by setting upstream and downstream traffic behavior weight.
The basic step of this method is as follows:
C1, from each track, in each section detecting device, obtain these three traffic flow parameters of the magnitude of traffic flow, speed and time occupancy in this this track of detection section according to the pre-determined sampling interval time; Then obtain the magnitude of traffic flow, speed and time occupancy that information issue interval characterizes this this track traffic stream characteristics of detection section; Traffic flow parameter is carried out to pre-service, and the information that obtains is issued interval and is characterized the magnitude of traffic flow, speed and the time occupancy of this detection section traffic stream characteristics.
C2, traffic flow parameter level and smooth.
C3, according to the pretreated traffic flow parameter computing velocity exponential sum occupation rate index that blocks up that blocks up.
C4, according to dividing the critical velocity corresponding to the critical speed calculation of the category of roads index that blocks up.
C5, calculate according to the speed exponential sum occupation rate index that blocks up that blocks up index---the traffic congestion index that comprehensively blocks up.
C6, the critical velocity obtaining according to the step c4 traffic congestion index that exponential sum c5 obtains that blocks up judges section traffic behavior.
The process of obtaining arithmetic for real-time traffic flow parameter in step c1 comprises:
C11, determine the section of required detection and section section Loop detector layout situation.
C12, specified data sampling interval: choosing sampling interval is 1 minute.
C13, by the magnitude of traffic flow, speed and time occupancy data on every track in each sampling interval of ring section detector acquisition.
C14, the magnitude of traffic flow, speed and time occupancy data by every track in each sampling interval calculate each information and issue the magnitude of traffic flow, speed and the time occupancy data on Nei Meitiao track, interval.
C15, the each track arithmetic for real-time traffic flow parameter obtaining in step c14 is carried out to pre-service, obtain characterizing the traffic flow parameter of this detection section.
In step c13, for the each section detecting device on every track, obtain the magnitude of traffic flow and the time occupancy in corresponding track, and the speed of each car;
From section detecting device, obtaining traffic flow parameter specifically comprises:
C131, obtain traffic flow parameter;
In each sampling interval, the magnitude of traffic flow is passed through the vehicle number of this section detecting device in one minute, and unit is veh;
In formula:i---of this sectionibar track;
k---thekindividual sampling interval;
n---in sampling interval by the of this section detecting devicencar;
---thekindividual sampling intervalion bar track, pass through the vehicle number of this section detecting device;
---the
kindividual sampling interval
ithe magnitude of traffic flow in bar track;
C132, acquisition speed parameter;
In each sampling interval, average velocity passes through the average velocity of all vehicles of this section detecting device in one minute, and unit is km/h;
In formula:
---the
kindividual sampling interval
ithe average velocity of bar track vehicle.
---the
nthe car speed during by section detecting device;
C133, acquisition time occupation rate parameter;
In each sampling interval, time occupancy takies the T.T. of detecting device and the number percent of 1 minute by all vehicles of this section detecting device in one minute.
In formula:
---the
ncar occupies the time of detecting device during by section detecting device;
---the
kindividual sampling interval
ithe occupation rate in bar track;
The abnormal data that needs to reject each track section detecting device in step c14, can adopt threshold value screening method, rejects the magnitude of traffic flow, speed and the time occupancy data that exceed certain threshold value; Then also need that the qualified data information of carrying out is issued to interval data synthetic, obtain each information and issue the magnitude of traffic flow, speed and the time occupancy data on Nei Meitiao track, interval, detailed step is as follows:
C141, each information are issued the traffic flow parameter processing on Nei Meitiao track, interval;
In issue interval, the magnitude of traffic flow is issued in interval and is passed through the vehicle number of this section detecting device, and unit is veh, that is:
In formula:
k---the
kindividual sampling interval,
;
m---information is issued interval,m=1,2,3,4,5,10,15 ... Deng;
---theibar lane information is issued the flow in interval;
C142, each information are issued the speed parameter processing on Nei Meitiao track, interval;
Issuing average velocity in interval is, the weighted mean value of each sampling interval average velocity, and taking the flow of each sampling interval as weight, unit is km/h, that is:
In formula:
---the
ibar lane information is issued the average velocity in interval.
C143, each information are issued the time occupancy parameter processing on Nei Meitiao track, interval
In issue interval, time occupancy is got the mean value of each sampling interval time occupancy, that is:
In formula:
---the
ibar lane information is issued the occupation rate in interval;
5. the traffic flow parameter of in step c15, each information being issued on Nei Meitiao track, interval synthesizes, and obtains the corresponding section magnitude of traffic flow, time occupancy and speed;
C151, the processing of section traffic flow parameter;
Section flow is to issue in interval to pass through the vehicle number of this section, i.e. each track flow sum, and unit is veh, that is:
In formula:j---thejindividual information is issued interval;
---the
jindividual information is issued the section flow at interval;
C152, the processing of section speed parameter
Section speed is the weighted mean value of track average velocity, and taking track flow as weight, unit is km/h, that is:
In formula:
---the
jindividual information is issued the velocity amplitude at interval.
C152, the processing of section time occupancy parameter;
Section time occupancy is the mean value of Ratio of driveway occupancy time, that is:
In formula:---thejindividual information is issued interlude occupation rate.
Step c2 detects to reality the profile data obtaining and carries out smoothing processing, and the stability of guarantee system operation, reduces the interference of enchancement factor, and the continuity and the stability that keep traffic behavior to issue, be calculated as follows in detail:
In formula:---the after level and smoothjindividual information is issued time occupancy in interval;
---the after level and smooth
jindividual information is issued velocity amplitude in interval;
---the
jindividual information is issued the section flow in interval;
---the
jindividual information is issued the time occupancy in interval;
---the
jindividual information is issued the speed in interval;
β---smoothing factor, hasβ1,β2,β3.
Step c3 is normalized traffic flow data, the speed of the obtaining exponential sum time occupancy index that blocks up that blocks up, and detailed step is as follows.
The c31 speed index that blocks up;
Suppose that the block up relation of index of speed and speed is linear, and the minimum value of speed (0) and maximal value (
) the corresponding speed index that blocks up is respectively 1,0.The block up computing formula of index of speed is:
In formula:jv---the speed index that blocks up;
---free stream velocity;
---the section speed after level and smooth.
The c32 time occupancy index that blocks up;
Time occupancy and the time occupancy index that blocks up is linear, the minimum value (0) of time occupancy and maximal value (
) the corresponding speed index that blocks up is respectively 0,1.The block up computing formula of index of time occupancy is as follows:
In formula:jo---the time occupancy index that blocks up;
---the historical maximal value of time occupancy;
---the section time occupancy after level and smooth.
Step c4 is according to the critical velocity value of urban road grade classification
,
urban road traffic state is divided into block up, slow and unimpeded three grades, according to critical velocity value
,
calculate the corresponding critical velocity index that blocks up:
In formula:j1,j2---the critical value that traffic behavior changes;
,
---the critical velocity that traffic behavior is divided;
vf---free stream velocity.
The step c5 fusion speed exponential sum time occupancy index that blocks up that blocks up, sets up the comprehensive index of blocking up---traffic congestion index, and detailed step is as follows:
Mono-group of detecting device section of c51
In formula,j---the composite target of traffic state judging, is called traffic congestion index;
jv---the speed index that blocks up;
jo---the time occupancy index that blocks up;
η---the weight coefficient of speed index and time occupancy index, value is 0-1, system can be defaulted as 0.5, and adjusts according to actual conditions.
C52 organizes detecting device section more
While there is many group section detecting devices, need to carry out comprehensive distinguishing and calculate according to the position of section detecting device the traffic congestion index of line.
In formula:jdi---detecting deviceitraffic congestion index;
ldi---detecting deviceito the distance of its downstream detector;
ld1---the distance apart from the nearest detecting device of downstream stop line apart from stop line.
C53 is without detecting section
When upstream and downstream section equal sensorless, this road section traffic volume state is grey, unknown.
Upstream or downstream are provided with detecting device, and this road section traffic volume state above downstream road section traffic behavior is foundation, utilize upstream and downstream traffic congestion index to calculate the traffic congestion index in this section.
In formula:jup---the traffic congestion index of upstream detector;
jdown---the traffic congestion index of downstream detector;
---the weight coefficient of upstream and downstream detecting device traffic congestion index, in the time of downstream road section sensorless,
; In the time of the sensorless of section, upstream,
; When upstream and downstream all has detecting device,
.
Step c6 has considered the problem of processing at the critical localisation of traffic behavior variation, the credibility interval that adopts traffic behavior to change.Definition ± Δjfor the normal fluctuation interval of state variation, the detailed step of traffic state judging is as follows.
If c61jin individual issue interval, when traffic behavior is red, judgej+1issue interval traffic state judging according to being:
1. work as
time, traffic behavior is red;
2. work as
time, traffic behavior is yellow;
3. work astime, traffic behavior is green.
If c62jin individual issue interval, when traffic behavior is yellow, judge thej+1issue interval traffic state judging according to being:
1. work as
time, traffic behavior is red;
2. at that time, traffic behavior was yellow;
3. work astime, traffic behavior is green.
If c63jin individual issue interval, when traffic behavior is green, judgej+1issue interval traffic state judging according to being:
1. work as
time, traffic behavior is red;
2. work as
time, traffic behavior is yellow;
3. work as
time, traffic behavior is green.
In above-mentioned rule, Δj1with Δj2value need to be according to actual conditions setting.
Beneficial effect of the present invention: the transport information that the present invention is based on some detection sections on section, consider multiple traffic flow parameter, this section of automatic discrimination traffic behavior of living in, the method adopts and tries one's best few discrimination threshold and make full use of existing resource simultaneously, is easy to Project Realization.
Embodiment
Below in conjunction with accompanying drawing, the present invention will be described in detail, and as shown in Figure 1, concrete steps of the present invention are:
Step 1 is obtained track sampling interval data:
Certain sectioniin the sampling interval of bar track, the computing formula of flow, average velocity and time occupancy is as follows:
(2)
(3)
In formula:i---of this sectionibar track;
k---thekindividual sampling interval;
n---in sampling interval by the of this section detecting devicencar;
---the
kindividual sampling interval
ion bar track, pass through the vehicle number of this section detecting device;
---the
ncar occupies the time of detecting device during by section detecting device;
---the
nthe car speed during by section detecting device;
---the
kindividual sampling interval
ithe flow in bar track;
---thekindividual sampling intervalithe occupation rate in bar track;
---the
kindividual sampling interval
ithe average velocity of bar track vehicle.
Step 2 is calculated lane information and is issued interval data:
Suppose that the time interval that information is issued is
mminute,
m=1,2,3,4,5,10,15 ... Deng,
ibar lane information is issued the flow in interval
for
mthe algebraic sum of individual sampling interval, time occupancy
for
mthe mean value of individual sampling interval, speed
for
mthe weighted mean value of individual sampling interval, computing formula is as follows.
In formula:
k---the
kindividual sampling interval,
;
m---information is issued interval,m=1,2,3,4,5,10,15 ... Deng;
---theibar lane information is issued the flow in interval;
---the
ibar lane information is issued the occupation rate in interval;
---the
ibar lane information is issued the average velocity in interval.
Step 3 is synthesized section traffic flow data:
Suppose that this detection section haslbar track, the computing formula of section traffic flow data is as follows:
In formula:j---thejindividual information is issued interval;
---the
jindividual information is issued the section flow at interval;
---the
jindividual information is issued interlude occupation rate;
---the
jindividual information is issued the velocity amplitude at interval.
Step 4 is carried out the level and smooth of traffic flow modes parameter:
Section time occupancy and speed adopt following methods to carry out smoothing processing:
In formula:
---the after level and smooth
jindividual information is issued time occupancy in interval;
---the after level and smoothjindividual information is issued velocity amplitude in interval;
---thejindividual information is issued the section flow in interval;
---the
jindividual information is issued the time occupancy in interval;
---the
jindividual information is issued the speed in interval;
β---smoothing factor, hasβ1,β2,β3.
Step 5 is determined block up index and interval:
(1) speed index
Suppose that the block up relation of index of speed and speed is linear, as Fig. 2 speed and speed are blocked up as shown in several relation curves.The block up computing formula of index of speed is:
In formula:jv---the speed index that blocks up;
---the section speed after level and smooth.
(2) the time occupancy index that blocks up
The block up relation curve of index of time occupancy and time occupancy blocks up as shown in exponential relationship curve as Fig. 3 time occupancy and time occupancy, and the block up computing formula of index of time occupancy is as follows:
(13)
In formula:jo---the time occupancy index that blocks up;
---the historical maximal value of time occupancy;
---the section time occupancy after level and smooth.
(3) speed block up index ranking divide
According to practical experience both domestic and external, the division of urban road traffic state should be the travel speed in section according to major parameter.Based on this consideration, the critical velocity value of dividing according to urban road grade and traffic behavior is divided into urban road traffic state and blocks up, slow and unimpeded three grades, as shown in the table.
Table 1 urban road speed divided rank
Can obtain the block up divided rank of index of speed according to the block up normalization formula of index and table 1 of speed, in table 2.
In formula:j1,j2---the critical value that traffic behavior changes;
,
---the critical velocity that traffic behavior is divided;
vf---free stream velocity.
The table 2 urban road speed index divided rank of blocking up
Step 6 is set up the comprehensive index of blocking up:
(1) calculate single group detecting device road section traffic volume index that blocks up
Consider the impact of speed and time occupancy, the composite target of setting up traffic state judging is as follows:
In formula,j---the composite target of traffic state judging, is called traffic congestion index;
jv---the speed index that blocks up;
jo---the time occupancy index that blocks up;
η---the weight coefficient of speed index and time occupancy index, value is 0-1, system can be defaulted as 0.5, and adjusts according to actual conditions.
(2) calculate many group detecting device road section traffic volumes index that blocks up
While there is many group section detecting devices, need to carry out comprehensive distinguishing and calculate according to the position of section detecting device the traffic congestion index of linej.Line Loop detector layout schematic diagram is as shown in Fig. 4 line Loop detector layout schematic diagram.
In figure,xupwithxdownrepresent respectively the coordinate position (with the coordinate replacement of intersection parking line) of upstream and downstream crossing;xd1,xd2,xdi...,xdnrepresent respectively from 1 tonthe coordinate position of individual section detecting device.
The traffic congestion index of line can obtain with the weighting processing of all detecting device traffic congestion indexes in this line, and its computing formula is as follows:
In formula:jdi---detecting deviceitraffic congestion index;
ldi---detecting deviceito the distance of its downstream detector;
ld1---the distance apart from the nearest detecting device of downstream stop line apart from stop line.
(3) determine without detecting the road section traffic volume index that blocks up
When upstream and downstream section equal sensorless, this road section traffic volume state is grey, unknown.
Upstream or downstream are provided with detecting device, and this road section traffic volume state above downstream road section traffic behavior is foundation, and traffic congestion formula of index is as follows:
In formula:jup---the traffic congestion index of upstream detector;
jdown---the traffic congestion index of downstream detector;
---the weight coefficient of upstream and downstream detecting device traffic congestion index, in the time of downstream road section sensorless,
; In the time of the sensorless of section, upstream,
; When upstream and downstream all has detecting device,
.
Step 7 is carried out traffic state judging:
In the time carrying out traffic state judging, the problem that the critical localisation that needs consideration to change at traffic behavior is especially processed, with the continuous and stable that ensures that traffic behavior changes.Therefore,, in the time carrying out traffic state judging, need to consider the credibility interval of state variation.Definition ± Δjfor the normal fluctuation interval of state variation,jdiscriminate interval can be expressed as traffic state judging figure (Fig. 5).
(1) ifjin individual issue interval, when traffic behavior is red, judgej+1issue interval traffic state judging according to being:
1. work as
time, traffic behavior is red;
2. work as
time, traffic behavior is yellow;
3. work as
time, traffic behavior is green.
(2) ifjin individual issue interval, when traffic behavior is yellow, judge thej+1issue interval traffic state judging according to being:
1. work as
time, traffic behavior is red;
2. work astime, traffic behavior is yellow;
3. work as
time, traffic behavior is green.
(3) ifjin individual issue interval, when traffic behavior is green, judgej+1issue interval traffic state judging according to being:
1. work as
time, traffic behavior is red;
2. work as
time, traffic behavior is yellow;
3. work astime, traffic behavior is green.
In above-mentioned rule, Δj1with Δj2value need to be according to actual conditions setting.