CN103761430A - Method for identifying peak periods of road networks on basis of floating cars - Google Patents

Method for identifying peak periods of road networks on basis of floating cars Download PDF

Info

Publication number
CN103761430A
CN103761430A CN201410012796.2A CN201410012796A CN103761430A CN 103761430 A CN103761430 A CN 103761430A CN 201410012796 A CN201410012796 A CN 201410012796A CN 103761430 A CN103761430 A CN 103761430A
Authority
CN
China
Prior art keywords
section
tci
road
sigma
floating car
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
Application number
CN201410012796.2A
Other languages
Chinese (zh)
Other versions
CN103761430B (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.)
Anhui Keli Information Industry Co Ltd
Original Assignee
Anhui Keli Information Industry Co Ltd
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 Anhui Keli Information Industry Co Ltd filed Critical Anhui Keli Information Industry Co Ltd
Priority to CN201410012796.2A priority Critical patent/CN103761430B/en
Publication of CN103761430A publication Critical patent/CN103761430A/en
Application granted granted Critical
Publication of CN103761430B publication Critical patent/CN103761430B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

The invention relates to a method for identifying peak periods of road networks on the basis of floating cars. The method includes steps of computing speeds of single vehicle samples of road sections by the aid of GPS (global positioning system) data of the floating cars; extracting average travel speeds of the road sections; computing periodic traffic congestion indexes TCI; extracting morning and evening peak hour starting and ending time points. Compared with the prior art, the method has the advantages that the shortcoming of deficiency of a method for identifying peak periods of road networks on the basis of a floating car technology at present can be overcome; time-varying traffic laws of the road networks can be analyzed via the existing floating car technology and an urban road congestion analysis system, so that the peak periods of the road networks can be extracted.

Description

A kind of road network recognition methods peak period based on Floating Car
Technical field
The present invention relates to traffic planninng technical field, is a kind of road network recognition methods peak period based on Floating Car specifically.
Background technology
Floating car technology is a kind of technology of obtaining road situation according to pavement of road operational vehicle dynamic location information, utilization can Real-time Collection vehicle with the Floating Car (taxi or bus) of GPS information displacement information, seasonal effect in time series vehicle location coordinate is mated with map, can obtain the speed data of unsteady vehicle.Floating car technology can store the data that gather a year in database into, utilizes cycle section velocity information to obtain cycle link flow information.
The annual day magnitude of traffic flow of AADT(road, annual average daily traffic) be the very important parameter of traffic model and management decision, in research fields such as traffic programme, highway layout, traffic safety, transport need analysis, traffic controls, there is crucial effect.Existing AADT calculates and does not recycle traditional manual counts, and also by the road-section average speed of utilizing floating car technology to obtain, the intellectuality that realizes road AADT by a series of models calculating is accurately estimated.The demand data that has met traffic programme, traffic design, traffic administration, has improved work efficiency.
Refer to the section in rush hour road traffic morning and evening causing due to the traffic of travelling frequently peak period (peak hours).The morning peak period is generally: 7:00-9:00, and evening peak is generally 17:00-19:00, and when concrete, discontinuous point there are differences with the difference in region, category of roads, section.Peak period, algorithm for estimating can be realized calculating peak period in city road network and region, for road conditions issuing service provides basic data source, in traffic administration and traffic-information service, played an important role.Utilizing floating car technology to calculate in the process of AADT, when calculating hour magnitude of traffic flow, the judgement of its peak period is according to the road traffic of setting forth in < < urban road traffic congestion assessment indicator system > > (exposure draft) peak period standard sooner or later.But the peak period in each city is all not identical, if utilize unified standard variant actual peak period obvious and each city, thereby affect the intelligent estimation precision of AADT.How developing a kind of recognition methods that can carry out based on floating car technology road network peak period for different cities has become and has been badly in need of the technical matters that solves.
Summary of the invention
To the object of the invention is the defect that there is no road network recognition methods peak period based on floating car technology in prior art in order solving, to provide a kind of road network recognition methods peak period based on Floating Car to solve the problems referred to above.
To achieve these goals, technical scheme of the present invention is as follows:
Road network recognition methods peak period based on Floating Car, comprises the following steps:
Utilize Floating Car gps data to calculate section bicycle sample speed;
Extract road-section average travel speed;
Computation period traffic congestion index TCI;
Extract beginning and ending time peak hour point sooner or later.
The described Floating Car gps data calculating section bicycle sample speed of utilizing comprises the following steps:
By Floating Car gps data obtain sample vehicle j the routing information { P of adjacent 2 of front and back of process i, i=1,2, L, n};
By path Δ d jwith mistiming Δ t jobtain the Average Travel Speed in this section of path
Figure BDA0000455314710000021
If approach section number only have one or kilometer/hour time, will be assigned to section P 1; Otherwise, the instantaneous velocity v by four kinds of traffic behavior principles in conjunction with starting point 1instantaneous velocity v with terminal 2, to each section speed of approach assignment respectively.
The determination methods of four kinds of described traffic behavior principles is as follows:
Deceleration regime, meets
Figure BDA0000455314710000024
time, initial section velocity amplitude is composed and is
Figure BDA0000455314710000025
other section velocity amplitude is total travel time Δ t jdeduct the travel time in initial section, then by distance, divided by this time, obtain speed;
Acceleration mode, meets
Figure BDA0000455314710000031
time, termination section velocity amplitude is composed and is
Figure BDA0000455314710000032
other section velocity amplitude is total travel time Δ t jdeduct the travel time in initial section, then by distance, divided by this time, obtain speed;
First slow down and accelerate afterwards, initial section velocity amplitude is composed as v 1, stop section velocity amplitude and compose as v 2, middle section velocity amplitude is total travel time Δ t jdeduct the travel time in initial section, then by distance, divided by this time, obtain speed;
First accelerate to slow down afterwards, approach section velocity amplitude is composed and is
Figure BDA0000455314710000033
Described extraction road-section average travel speed computing formula is
V i = V i = l i t i = l i ( &Sigma; j = 1 n l i v j ) / n = n &Sigma; j = 1 n 1 v j , if n i &NotEqual; 0 V i = V i &OverBar; , if n i = 0
Wherein, V ifor segmental arc P iaverage velocity, l ifor segmental arc P ilength, t ijbe j car segmental arc P in path ion travel time, n ifor segmental arc P ithe number of vehicles that upper participation is calculated.
Described computation period traffic congestion index TCI comprises the following steps:
Based on road-section average travel speed V icarry out congestion status identification, judge the section that blocks up;
Calculate the section mileage ratio RCR that blocks up, calculate respectively through street mileage ratio RCRf, trunk roads mileage ratio RCRa, secondary distributor road mileage ratio RCRm and the branch road mileage ratio RCRl that blocks up that blocks up that blocks up that blocks up, computing formula is as follows:
Figure BDA0000455314710000035
Figure BDA0000455314710000036
Figure BDA0000455314710000037
RCR=RCRf*ω 1+RCRa*ω 2+RCRm*ω 3+RCRl*ω 4
Wherein,
L (i) is the length of section i, the length that Lc (i) is the section i that gets congestion,
N f: the total number in section, through street,
N a: the total number in trunk roads section,
N m: the total number in secondary distributor road section,
N l: the total number in branch road section,
W1, w2, w3, w4 represents respectively the weight of each grade road,
Figure BDA0000455314710000042
Calculate the road grid traffic index TCI that blocks up, computing formula is as follows:
TCI = a 2 ( 0 &le; a &le; 4 ) 2 + a - 4 2 ( 4 < a &le; 8 ) 4 + ( a - 8 ) &times; 2 3 ( 8 < a &le; 11 ) 6 + ( a - 11 ) &times; 2 3 ( 11 < a &le; 14 ) 8 + a - 14 5 ( 14 < a &le; 24 ) 10 ( a &GreaterEqual; 24 )
Wherein: a=RCR*100.
Described extraction sooner or later beginning and ending time peak hour point comprises the following steps:
Judge in one day 24 hours whether Normal Distribution of TCI curve, if Normal Distribution enters next step, calculate, if disobey normal distribution, represent traffic abnormity on the same day, reject data and also reselect data;
Set confidence value c, c is the error range that estimated value and population parameter allow;
According to 24 hours TCI changing values, maximal value and the minimum value of getting TCI,
Take 0 o'clock to 12 o'clock as dividing, TCI maximal value is max_a, the number of cycles that wherein a is 1-288, and the cycle is 5 minutes; The anterior minimum value of TCI is min_t1, and wherein t1 is periodicity corresponding to minimum value; TCI rear portion minimum value min_t2, wherein t2 is periodicity corresponding to minimum value;
Take 12 o'clock to 24 o'clock as dividing, TCI maximal value max_p, wherein p is the number of cycles of 1-288, the cycle is 5 minutes; The anterior minimum value min_t3 of TCI, wherein t3 is periodicity corresponding to minimum value; TCI rear portion minimum value min_t4, wherein t4 is periodicity corresponding to minimum value;
Zoning total area S1, S2, S3, S4,
S 1 = &Sigma; t 1 a ( TCI _ i - min _ t 1 ) i &Element; [ t 1 , a ] ,
S 2 = &Sigma; a t 2 ( TCI _ i - min _ t 2 ) i &Element; [ a , t 2 ] ,
S 3 = &Sigma; t 3 p ( TCI _ i - min _ t 3 ) i &Element; [ t 3 , p ] ,
S 4 = &Sigma; p t 4 ( TCI _ i - min _ t 4 ) i &Element; [ p , t 4 ] ;
Calculate variance area S 1', S 2', S 3', S 4',
S ' 1 = &Sigma; j 1 a ( TCI _ j 1 - min _ t 1 ) j &Element; [ t 1 , a ] ,
S 2 ' = &Sigma; a j 2 ( TCI _ j 2 - min _ t 2 ) j &Element; [ a , t 2 ] ,
S 3 ' = &Sigma; j 3 p ( TCI _ j 3 - min _ t 3 ) j &Element; [ t 3 , p ] ,
S 4 ' = &Sigma; p j 4 ( TCI _ j 4 - min _ t 4 ) j &Element; [ p , t 4 ] ;
By S1, S2, S3, S4 and S 1', S 2', S 3', S 4' substitution formula c=S corresponding to difference i'/S isolve, by solving, obtain respectively j1, j2, j3, j4, wherein, i=1,2,3,4, S i' be variance area, S iit is interval area;
Determine that the morning peak period is T1 to T2, determine that the evening peak period is T3 to T4, wherein T1, T2, T3, the T4 difference start time in cycle of corresponding j1, j2, j3, j4 successively.
Described judge TCI curve whether the formula of Normal Distribution be:
0.9 < M x &OverBar; < 1.1 And x &OverBar; > 3 s ,
Wherein: be arithmetic mean, M is median, and s is standard deviation.
Beneficial effect
A kind of road network recognition methods peak period based on Floating Car of the present invention, compared with prior art can pass through existing floating car technology and urban road congestion analysis system, becomes rule while therefrom analyzing road grid traffic, extracts road network peak period.Obtain the most serious period of road grid traffic load, for traffic administration person and traffic programme person provide Data support, improve the intelligent estimation precision of AADT.
Accompanying drawing explanation
Fig. 1 is method flow diagram of the present invention
Fig. 2 is TCI24 hour curvilinear motion figure and relevant parameter mark figure
Embodiment
For making that architectural feature of the present invention and the effect reached are had a better understanding and awareness, in order to preferred embodiment and accompanying drawing, coordinate detailed explanation, be described as follows:
As shown in Figure 1, a kind of road network recognition methods peak period based on Floating Car of the present invention, comprises the following steps:
The first step, utilizes Floating Car gps data to calculate section bicycle sample speed.
Utilize Floating Car gps data to calculate bicycle sample mean travelling speed in the measurement period in section, first, by Floating Car gps data obtain sample vehicle j the routing information { P of adjacent 2 of front and back of process i, i=1,2, L, n}.Secondly, based on gps data, can pass through path Δ d jwith mistiming Δ t jobtain the Average Travel Speed in this section of path
Figure BDA0000455314710000071
again, when approach section number only have one represent do not cross over crossing or
Figure BDA0000455314710000072
while kilometer/hour being unimpeded state, will be assigned to section P 1; Otherwise, by the instantaneous velocity v in conjunction with starting point 1instantaneous velocity v with terminal 2, point four kinds of traffic behaviors are to each section speed of approach assignment respectively.
When Floating Car is in deceleration regime, meet
Figure BDA0000455314710000074
time, initial section velocity amplitude is composed and is
Figure BDA0000455314710000075
other section velocity amplitude is total travel time Δ t jdeduct the travel time in initial section, then by distance, divided by this time, obtain speed, then speed is assigned to section P 1.
When Floating Car is in acceleration mode, meet
Figure BDA0000455314710000076
time, termination section velocity amplitude is composed and is
Figure BDA0000455314710000077
other section velocity amplitude is total travel time Δ t jdeduct the travel time in initial section, then by distance, divided by this time, obtain speed, then speed is assigned to section P 1.
When Floating Car is in first slowing down and accelerate afterwards, initial section velocity amplitude is composed as v 1, stop section velocity amplitude and compose as v 2, middle section velocity amplitude is total travel time Δ t jdeduct the travel time in initial section, then by distance, divided by this time, obtain speed, then speed is assigned to section P 1.
When Floating Car deceleration after first accelerating, approach section velocity amplitude is composed and is
Figure BDA0000455314710000078
again will
Figure BDA0000455314710000079
be assigned to section P 1.
Second step, extracts road-section average travel speed.
Extracting road-section average travel speed computing formula is
V i = V i = l i t i = l i ( &Sigma; j = 1 n l i v j ) / n = n &Sigma; j = 1 n 1 v j , if n i &NotEqual; 0 V i = V i &OverBar; , if n i = 0
Wherein, V ifor segmental arc P iaverage velocity, l ifor segmental arc P ilength, t ijbe j car segmental arc P in path ion travel time, n ifor segmental arc P ithe number of vehicles that upper participation is calculated.Here, work as n iequal 0, while there is no data cover on this section, we supplement by the historical average speeds of one week different time sections of historical accumulation; Work as n ibe not equal at 0 o'clock, road-section average travel speed is the harmonic average speed of multiple samples.
The 3rd step, computation period traffic congestion index TCI.
Cycle traffic congestion index refers in a measurement period (normally 5 minutes), describes the level of blocking up of current region road network with the numerical value of one 0~10, is the quantizating index of describing the degree of blocking up.Calculation procedure is as follows:
First, based on road-section average travel speed V icarry out congestion status identification, judge the section that blocks up.
According to the speed interval table in table 1, judge current section and whether belong to congestion status, the section data for congestion status are put forward for next step processing.Threshold value table Data Source in table 1 is according to stipulating in the < < urban traffic management assessment indicator system > > of announcement in 2010, the corresponding threshold value of motor vehicle average travel speed in city thoroughfare, and according to Hefei City's actual traffic characteristic, carried out that fine setting obtains.While using in different cities, can adapt to fine setting according to local actual traffic characteristic.
The road condition of table 1 based on road-section average travelling speed divided speed interval table
Classification 1 grade: Free-flow 2 grades: unimpeded 3 grades: jogging 4 grades: crowded 5 grades: block up
Through street More than 65km/h 65-40km/h 40-30km/h 30-15km/h Below 15km/h
Major trunk roads More than 50km/h 50-30km/h 30-22km/h 22-12km/h Below 12km/h
Subsidiary road More than 40km/h 40-25km/h 25-18km/h 18-10km/h Below 10km/h
Branch road More than 35km/h 35-23km/h 23-16km/h 16-8km/h Below 8km/h
Secondly, calculate the section mileage ratio RCR that blocks up.
Calculate respectively through street mileage ratio RCRf, trunk roads mileage ratio RCRa, secondary distributor road mileage ratio RCRm and the branch road mileage ratio RCRl that blocks up that blocks up that blocks up that blocks up, computing formula is as follows:
Figure BDA0000455314710000081
Figure BDA0000455314710000082
Figure BDA0000455314710000083
Figure BDA0000455314710000084
RCR=RCRf*ω 1+RCRa*ω 2+RCRm*ω 3+RCRl*ω 4
Wherein,
L (i) is the length of section i, the length that Lc (i) is the section i that gets congestion,
N f: the total number in section, through street,
N a: the total number in trunk roads section,
N m: the total number in secondary distributor road section,
N l: the total number in branch road section,
W1, w2, w3, w4 represents respectively the weight of each grade road,
Figure BDA0000455314710000091
W1, w2, w3, w4 represent respectively the weight of each grade road, from the total current mileage historical data statistical study of each grade road of road network vehicle, draw.In table 2, have respectively working day weights recommendation tables and festivals or holidays weights recommendation tables, in practice, also can with reference to provincial standard carry out w1, w2, w3, w4 appointment calculate.
Table 2 < < urban road traffic congestion assessment indicator system > > Beijing provincial standard
Working day Through street Trunk roads Secondary distributor road Branch road Add up to
Peak period 0.19 0.43 0.15 0.23 1.00
Full-time average 0.20 0.41 0.16 0.23 1.00
Festivals or holidays Through street Trunk roads Secondary distributor road Branch road Add up to
Full-time average 0.20 0.41 0.16 0.23 1.00
Finally, calculate the road grid traffic index TCI that blocks up.
Computing formula is as follows:
TCI = a 2 ( 0 &le; a &le; 4 ) 2 + a - 4 2 ( 4 < a &le; 8 ) 4 + ( a - 8 ) &times; 2 3 ( 8 < a &le; 11 ) 6 + ( a - 11 ) &times; 2 3 ( 11 < a &le; 14 ) 8 + a - 14 5 ( 14 < a &le; 24 ) 10 ( a &GreaterEqual; 24 )
Wherein: a=RCR*100, according to the value difference of a, selects corresponding computing formula.
The 4th step, extracts beginning and ending time peak hour point sooner or later.
Comprise the steps:
(1) as shown in Figure 2, first judge in one day 24 hours whether Normal Distribution of TCI curve, if Normal Distribution enters next step, calculate, if disobey normal distribution, represent traffic abnormity on the same day, have emergency situations to occur.Data are difficult for, as with reference to computational data, therefore rejecting data and reselecting data.Differentiate the TCI curve method whether Normal Distribution adopts, it is one of existing statistics the inside normal distribution-test method, to judge by the ratio of sample median M and arithmetic mean and the relation of arithmetic mean and standard deviation, reflection peak shape and kurtosis, formula is as follows:
0.9 < M x &OverBar; < 1.1 And x &OverBar; > 3 s ,
Wherein:
Figure BDA0000455314710000103
be arithmetic mean, M is median, and s is standard deviation.
(2) set confidence value c, c is the error range that estimated value and population parameter allow.Confidence value is the estimated value of judgement, can specify according to city and decision maker's actual needs, in general, in order to guarantee larger confidence level, generally gets confidence value and is greater than 90.
(3) according to 24 hours TCI changing values, maximal value and the minimum value of getting TCI.
Take 0 o'clock to 12 o'clock as dividing, TCI maximal value is max_a, the number of cycles that wherein a is 1-288, and the cycle is 5 minutes, 288 is to divide and get take 5 minutes as a cycle according to 24 hours.The anterior minimum value of TCI is min_t1, and wherein t1 is periodicity corresponding to minimum value.TCI rear portion minimum value min_t2, wherein t2 is periodicity corresponding to minimum value.
Take 12 o'clock to 24 o'clock as dividing, TCI maximal value max_p, wherein p is the number of cycles of 1-288, the cycle is 5 minutes.The anterior minimum value min_t3 of TCI, wherein t3 is periodicity corresponding to minimum value.TCI rear portion minimum value min_t4, wherein t4 is periodicity corresponding to minimum value.
(4) zoning total area S1, S2, S3, S4, its computing formula is as follows:
S 1 = &Sigma; t 1 a ( TCI _ i - min _ t 1 ) i &Element; [ t 1 , a ] ,
S 2 = &Sigma; a t 2 ( TCI _ i - min _ t 2 ) i &Element; [ a , t 2 ] ,
S 3 = &Sigma; t 3 p ( TCI _ i - min _ t 3 ) i &Element; [ t 3 , p ] ,
S 4 = &Sigma; p t 4 ( TCI _ i - min _ t 4 ) i &Element; [ p , t 4 ] .
(5) calculate variance area S 1', S 2', S 3', S 4', its computing formula is as follows:
S ' 1 = &Sigma; j 1 a ( TCI _ j 1 - min _ t 1 ) j &Element; [ t 1 , a ] ,
S 2 ' = &Sigma; a j 2 ( TCI _ j 2 - min _ t 2 ) j &Element; [ a , t 2 ] ,
S 3 ' = &Sigma; j 3 p ( TCI _ j 3 - min _ t 3 ) j &Element; [ t 3 , p ] ,
S 4 ' = &Sigma; p j 4 ( TCI _ j 4 - min _ t 4 ) j &Element; [ p , t 4 ] .
(6) by S1, S2, S3, S4 and S 1', S 2', S 3', S 4' substitution formula c=S corresponding to difference i'/S isolve, by solving, obtain respectively j1, j2, j3, j4, wherein, i=1,2,3,4, S i' be variance area, S iit is interval area.
(7) determine that the morning peak period is T1 to T2, determine that the evening peak period is T3 to T4, wherein T1, T2, T3, the T4 difference start time in cycle of corresponding j1, j2, j3, j4 successively.Because j1, j2, j3, j4 this representative are, by 24 hours take 5 minutes 288 periodicities as the cycle, by j1, j2, j3, the concrete time point T1 of j4 representative, T2, T3, T4, thereby just can judge the morning peak period, be that T1-T2, evening peak period are T3-T4.
More than show and described ultimate principle of the present invention, principal character and advantage of the present invention.The technician of the industry should understand; the present invention is not restricted to the described embodiments; what in above-described embodiment and instructions, describe is principle of the present invention; the present invention also has various changes and modifications without departing from the spirit and scope of the present invention, and these changes and improvements all fall in claimed scope of the present invention.The protection domain that the present invention requires is defined by appending claims and equivalent thereof.

Claims (7)

1. road network recognition methods peak period based on Floating Car, is characterized in that, comprises the following steps:
1) utilize Floating Car gps data to calculate section bicycle sample speed;
2) extract road-section average travel speed;
3) computation period traffic congestion index TCI;
4) extract beginning and ending time peak hour point sooner or later.
2. a kind of road network recognition methods peak period based on Floating Car according to claim 1, is characterized in that: the described Floating Car gps data calculating section bicycle sample speed of utilizing comprises the following steps:
21) by Floating Car gps data obtain sample vehicle j the routing information { P of adjacent 2 of front and back of process i, i=1,2, L, n};
22) by path Δ d jwith mistiming Δ t jobtain the Average Travel Speed in this section of path v &OverBar; = &Delta; d j / &Delta; t j ;
23) if approach section number only have one or
Figure FDA0000455314700000012
kilometer/hour time, will
Figure FDA0000455314700000013
be assigned to section P 1; Otherwise, the instantaneous velocity v by four kinds of traffic behavior principles in conjunction with starting point 1instantaneous velocity v with terminal 2, to each section speed of approach assignment respectively.
3. a kind of road network recognition methods peak period based on Floating Car according to claim 2, is characterized in that, the determination methods of four kinds of described traffic behavior principles is as follows:
31) deceleration regime, meets time, initial section velocity amplitude is composed and is
Figure FDA0000455314700000015
other section velocity amplitude is total travel time Δ t jdeduct the travel time in initial section, then by distance, divided by this time, obtain speed;
32) acceleration mode, meets
Figure FDA0000455314700000016
time, termination section velocity amplitude is composed and is
Figure FDA0000455314700000017
other section velocity amplitude is total travel time Δ t jdeduct the travel time in initial section, then by distance, divided by this time, obtain speed;
33) first slow down and accelerate afterwards, initial section velocity amplitude is composed as v 1, stop section velocity amplitude and compose as v 2, middle section velocity amplitude is total travel time Δ t jdeduct the travel time in initial section, then by distance, divided by this time, obtain speed;
34) first accelerate to slow down afterwards, approach section velocity amplitude is composed and is
Figure FDA0000455314700000018
4. a kind of road network recognition methods peak period based on Floating Car according to claim 1, is characterized in that: described extraction road-section average travel speed computing formula is
V i = V i = l i t i = l i ( &Sigma; j = 1 n l i v j ) / n = n &Sigma; j = 1 n 1 v j , if n i &NotEqual; 0 V i = V i &OverBar; , if n i = 0
Wherein, V ifor segmental arc P iaverage velocity, l ifor segmental arc P ilength, t ijbe j car segmental arc P in path ion travel time, n ifor segmental arc P ithe number of vehicles that upper participation is calculated.
5. a kind of road network recognition methods peak period based on Floating Car according to claim 1, is characterized in that, described computation period traffic congestion index TCI comprises the following steps:
51) based on road-section average travel speed V icarry out congestion status identification, judge the section that blocks up;
52) calculate the section mileage ratio RCR that blocks up, calculate respectively through street mileage ratio RCRf, trunk roads mileage ratio RCRa, secondary distributor road mileage ratio RCRm and the branch road mileage ratio RCRl that blocks up that blocks up that blocks up that blocks up, computing formula is as follows:
Figure FDA0000455314700000022
Figure FDA0000455314700000023
Figure FDA0000455314700000025
RCR=RCRf*ω 1+RCRa*ω 2+RCRm*ω 3+RCRl*ω 4
Wherein,
L (i) is the length of section i, the length that Lc (i) is the section i that gets congestion,
N f: the total number in section, through street,
N a: the total number in trunk roads section,
N m: the total number in secondary distributor road section,
N l: the total number in branch road section,
W1, w2, w3, w4 represents respectively the weight of each grade road,
Figure FDA0000455314700000032
53) calculate the road grid traffic index TCI that blocks up, computing formula is as follows:
TCI = a 2 ( 0 &le; a &le; 4 ) 2 + a - 4 2 ( 4 < a &le; 8 ) 4 + ( a - 8 ) &times; 2 3 ( 8 < a &le; 11 ) 6 + ( a - 11 ) &times; 2 3 ( 11 < a &le; 14 ) 8 + a - 14 5 ( 14 < a &le; 24 ) 10 ( a &GreaterEqual; 24 )
Wherein: a=RCR*100.
6. a kind of road network recognition methods peak period based on Floating Car according to claim 1, is characterized in that, described extraction sooner or later beginning and ending time peak hour point comprises the following steps:
61) judge in one day 24 hours whether Normal Distribution of TCI curve, if Normal Distribution enters next step, calculate, if disobey normal distribution, represent traffic abnormity on the same day, reject data and also reselect data;
62) set confidence value c, c is the error range that estimated value and population parameter allow;
63) according to 24 hours TCI changing values, maximal value and the minimum value of getting TCI,
Take 0 o'clock to 12 o'clock as dividing, TCI maximal value is max_a, the number of cycles that wherein a is 1-288, and the cycle is 5 minutes; The anterior minimum value of TCI is min_t1, and wherein t1 is periodicity corresponding to minimum value; TCI rear portion minimum value min_t2, wherein t2 is periodicity corresponding to minimum value;
Take 12 o'clock to 24 o'clock as dividing, TCI maximal value max_p, wherein p is the number of cycles of 1-288, the cycle is 5 minutes; The anterior minimum value min_t3 of TCI, wherein t3 is periodicity corresponding to minimum value; TCI rear portion minimum value min_t4, wherein t4 is periodicity corresponding to minimum value;
64) zoning total area S1, S2, S3, S4,
S 1 = &Sigma; t 1 a ( TCI _ i - min _ t 1 ) i &Element; [ t 1 , a ] ,
S 2 = &Sigma; a t 2 ( TCI _ i - min _ t 2 ) i &Element; [ a , t 2 ] ,
S 3 = &Sigma; t 3 p ( TCI _ i - min _ t 3 ) i &Element; [ t 3 , p ] ,
S 4 = &Sigma; p t 4 ( TCI _ i - min _ t 4 ) i &Element; [ p , t 4 ] ;
65) calculate variance area S 1', S 2', S 3', S 4',
S ' 1 = &Sigma; j 1 a ( TCI _ j 1 - min _ t 1 ) j &Element; [ t 1 , a ] ,
S 2 ' = &Sigma; a j 2 ( TCI _ j 2 - min _ t 2 ) j &Element; [ a , t 2 ] ,
S 3 ' = &Sigma; j 3 p ( TCI _ j 3 - min _ t 3 ) j &Element; [ t 3 , p ] ,
S 4 ' = &Sigma; p j 4 ( TCI _ j 4 - min _ t 4 ) j &Element; [ p , t 4 ] ;
66) by S1, S2, S3, S4 and S 1', S 2', S 3', S 4' substitution formula c=S corresponding to difference i'/S isolve, by solving, obtain respectively j1, j2, j3, j4, wherein, i=1,2,3,4, S i' be variance area, S iit is interval area;
67) determine that the morning peak period is T1 to T2, determine that the evening peak period is T3 to T4, wherein T1, T2, T3, the T4 difference start time in cycle of corresponding j1, j2, j3, j4 successively.
7. a kind of road network recognition methods peak period based on Floating Car according to claim 6, is characterized in that, described judge TCI curve whether the formula of Normal Distribution be:
0.9 < M x &OverBar; < 1.1 And x &OverBar; > 3 s ,
Wherein:
Figure FDA0000455314700000053
be arithmetic mean, M is median, and s is standard deviation.
CN201410012796.2A 2014-01-10 2014-01-10 A kind of road network peak period recognition methods based on Floating Car Active CN103761430B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410012796.2A CN103761430B (en) 2014-01-10 2014-01-10 A kind of road network peak period recognition methods based on Floating Car

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410012796.2A CN103761430B (en) 2014-01-10 2014-01-10 A kind of road network peak period recognition methods based on Floating Car

Publications (2)

Publication Number Publication Date
CN103761430A true CN103761430A (en) 2014-04-30
CN103761430B CN103761430B (en) 2017-07-07

Family

ID=50528666

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410012796.2A Active CN103761430B (en) 2014-01-10 2014-01-10 A kind of road network peak period recognition methods based on Floating Car

Country Status (1)

Country Link
CN (1) CN103761430B (en)

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104574967A (en) * 2015-01-14 2015-04-29 合肥革绿信息科技有限公司 City large-area road network traffic sensing method based on plough satellite
CN105070056A (en) * 2015-07-23 2015-11-18 合肥革绿信息科技有限公司 Intersection traffic congestion index calculation method based on floating car
CN105139647A (en) * 2015-07-27 2015-12-09 福建工程学院 Real-time road congestion detection method
CN105139645A (en) * 2015-07-23 2015-12-09 合肥革绿信息科技有限公司 Urban regional road network operation index assessment method based on floating car technology
CN105261210A (en) * 2015-07-23 2016-01-20 合肥革绿信息科技有限公司 Beidou-equipment-based calculating method of traffic congestion index of road
CN105405294A (en) * 2015-12-30 2016-03-16 杭州中奥科技有限公司 Early warning method of traffic congestion roads
CN105489016A (en) * 2016-02-01 2016-04-13 北京交通发展研究中心 Urban road operation condition evaluation method
CN105608895A (en) * 2016-03-04 2016-05-25 大连理工大学 Local abnormity factor-based urban heavy-traffic road detection method
CN105788289A (en) * 2014-12-17 2016-07-20 上海宝康电子控制工程有限公司 Method and system for realizing traffic condition assessment and analysis based on computer software system
CN105931463A (en) * 2016-06-27 2016-09-07 安徽四创电子股份有限公司 Method for calculating road traffic performance index based on traffic surface radar
CN106448159A (en) * 2016-09-09 2017-02-22 蔡诚昊 Road traffic hierarchical early warning method based on dynamic traffic information
CN106600965A (en) * 2017-01-19 2017-04-26 上海理工大学 Sharpness-based automatic identification method for morning and evening peak periods of traffic flows
CN106781509A (en) * 2017-03-06 2017-05-31 长安大学 A kind of collaborative urban road congestion detection method based on V2V
CN106846806A (en) * 2017-03-07 2017-06-13 北京工业大学 Urban highway traffic method for detecting abnormality based on Isolation Forest
CN109697854A (en) * 2019-02-25 2019-04-30 公安部交通管理科学研究所 The urban road traffic state assessment method of various dimensions
CN114627642A (en) * 2022-02-25 2022-06-14 青岛海信网络科技股份有限公司 Traffic jam identification method and device

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6577946B2 (en) * 2001-07-10 2003-06-10 Makor Issues And Rights Ltd. Traffic information gathering via cellular phone networks for intelligent transportation systems
CN101794507A (en) * 2009-07-13 2010-08-04 北京工业大学 Method for evaluating macroscopic road network traffic state based on floating car data
CN103000027A (en) * 2012-12-19 2013-03-27 安徽科力信息产业有限责任公司 Intelligent traffic guidance method based on floating car under congestion condition

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6577946B2 (en) * 2001-07-10 2003-06-10 Makor Issues And Rights Ltd. Traffic information gathering via cellular phone networks for intelligent transportation systems
CN101794507A (en) * 2009-07-13 2010-08-04 北京工业大学 Method for evaluating macroscopic road network traffic state based on floating car data
CN103000027A (en) * 2012-12-19 2013-03-27 安徽科力信息产业有限责任公司 Intelligent traffic guidance method based on floating car under congestion condition

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
佟阳春: "基于浮动车和WSN的交通数据处理方法研究", 《中国优秀硕士学位论文全文数据库信息科技辑》, no. 10, 15 October 2012 (2012-10-15) *
郭继孚等: "北京市交通拥堵宏观评价指标体系开发及其应用", 《2007第三届中国智能交通年会论文集》, 1 December 2007 (2007-12-01) *
黄赟等: "基于交通拥堵指数城市高峰小时估计方法研究", 《科技资讯》, no. 29, 13 October 2012 (2012-10-13), pages 199 - 200 *

Cited By (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105788289A (en) * 2014-12-17 2016-07-20 上海宝康电子控制工程有限公司 Method and system for realizing traffic condition assessment and analysis based on computer software system
CN104574967A (en) * 2015-01-14 2015-04-29 合肥革绿信息科技有限公司 City large-area road network traffic sensing method based on plough satellite
CN104574967B (en) * 2015-01-14 2016-08-24 合肥革绿信息科技有限公司 A kind of city based on Big Dipper large area road grid traffic cognitive method
CN105070056A (en) * 2015-07-23 2015-11-18 合肥革绿信息科技有限公司 Intersection traffic congestion index calculation method based on floating car
CN105139645A (en) * 2015-07-23 2015-12-09 合肥革绿信息科技有限公司 Urban regional road network operation index assessment method based on floating car technology
CN105261210A (en) * 2015-07-23 2016-01-20 合肥革绿信息科技有限公司 Beidou-equipment-based calculating method of traffic congestion index of road
CN105139647A (en) * 2015-07-27 2015-12-09 福建工程学院 Real-time road congestion detection method
CN105139647B (en) * 2015-07-27 2017-12-08 福建工程学院 A kind of method that congestion in road detects in real time
CN105405294A (en) * 2015-12-30 2016-03-16 杭州中奥科技有限公司 Early warning method of traffic congestion roads
CN105489016A (en) * 2016-02-01 2016-04-13 北京交通发展研究中心 Urban road operation condition evaluation method
CN105489016B (en) * 2016-02-01 2018-07-10 北京交通发展研究中心 A kind of urban road operating condition appraisal procedure
CN105608895A (en) * 2016-03-04 2016-05-25 大连理工大学 Local abnormity factor-based urban heavy-traffic road detection method
CN105931463A (en) * 2016-06-27 2016-09-07 安徽四创电子股份有限公司 Method for calculating road traffic performance index based on traffic surface radar
CN106448159A (en) * 2016-09-09 2017-02-22 蔡诚昊 Road traffic hierarchical early warning method based on dynamic traffic information
CN106448159B (en) * 2016-09-09 2018-11-02 蔡诚昊 A kind of road traffic grading forewarning system method based on dynamic information
CN106600965A (en) * 2017-01-19 2017-04-26 上海理工大学 Sharpness-based automatic identification method for morning and evening peak periods of traffic flows
CN106781509A (en) * 2017-03-06 2017-05-31 长安大学 A kind of collaborative urban road congestion detection method based on V2V
CN106781509B (en) * 2017-03-06 2019-09-10 长安大学 A kind of collaborative urban road congestion detection method based on V2V
CN106846806A (en) * 2017-03-07 2017-06-13 北京工业大学 Urban highway traffic method for detecting abnormality based on Isolation Forest
CN109697854A (en) * 2019-02-25 2019-04-30 公安部交通管理科学研究所 The urban road traffic state assessment method of various dimensions
CN109697854B (en) * 2019-02-25 2021-07-16 公安部交通管理科学研究所 Multi-dimensional urban road traffic state evaluation method
CN114627642A (en) * 2022-02-25 2022-06-14 青岛海信网络科技股份有限公司 Traffic jam identification method and device

Also Published As

Publication number Publication date
CN103761430B (en) 2017-07-07

Similar Documents

Publication Publication Date Title
CN103761430A (en) Method for identifying peak periods of road networks on basis of floating cars
WO2019047905A1 (en) Road traffic analysis system, method and apparatus
CN102737504B (en) Method for estimating bus arrival time in real time based on drive characteristics
CN111091720B (en) Congestion road section identification method and device based on signaling data and floating car data
CN102855760B (en) On-line queuing length detection method based on floating vehicle data
CN108629973A (en) Road section traffic volume congestion index computational methods based on fixed test equipment
CN102419905B (en) Traffic-wave theory-based traffic influence area determining method of expressway accidents
CN110363985B (en) Traffic data analysis method, device, storage medium and equipment
CN105355049B (en) A kind of highway evaluation of running status method based on macroscopical parent map
CN106384509A (en) Urban road driving time distribution estimation method considering taxi operation states
CN105070056A (en) Intersection traffic congestion index calculation method based on floating car
CN105139645A (en) Urban regional road network operation index assessment method based on floating car technology
CN104732075A (en) Real-time prediction method for urban road traffic accident risk
CN105825669A (en) System and method for identifying urban expressway traffic bottlenecks
CN102968901A (en) Method for acquiring regional congestion information and regional congestion analyzing device
CN103903433A (en) Real-time dynamic judgment method and device for road traffic state
CN104778834A (en) Urban road traffic jam judging method based on vehicle GPS data
CN110633558B (en) Urban traffic system modeling system
CN112767694B (en) Traffic optimization method and device for relieving road congestion
CN112950946B (en) Method for judging traffic jam reason in construction period of municipal road
CN106803347A (en) Urban intersection traffic state judging method based on RFID data
CN103632541B (en) Traffic incident road chain detection and data filling method
CN103942952B (en) A kind of road network functional hierarchy state grade appraisal procedure
CN103632537B (en) A kind of urban road AADT method of estimation based on Floating Car
CN110827537B (en) Method, device and equipment for setting tidal lane

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