CN105869398A - Road traffic open degree judging method based on K-means cluster - Google Patents

Road traffic open degree judging method based on K-means cluster Download PDF

Info

Publication number
CN105869398A
CN105869398A CN201610304571.3A CN201610304571A CN105869398A CN 105869398 A CN105869398 A CN 105869398A CN 201610304571 A CN201610304571 A CN 201610304571A CN 105869398 A CN105869398 A CN 105869398A
Authority
CN
China
Prior art keywords
data
road
speed
value
traffic
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
CN201610304571.3A
Other languages
Chinese (zh)
Other versions
CN105869398B (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.)
Zhejiang University of Technology ZJUT
Original Assignee
Zhejiang University of Technology ZJUT
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 Zhejiang University of Technology ZJUT filed Critical Zhejiang University of Technology ZJUT
Priority to CN201610304571.3A priority Critical patent/CN105869398B/en
Publication of CN105869398A publication Critical patent/CN105869398A/en
Application granted granted Critical
Publication of CN105869398B publication Critical patent/CN105869398B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/01Detecting movement of traffic to be counted or controlled
    • G08G1/0104Measuring and analyzing of parameters relative to traffic conditions
    • G08G1/0108Measuring and analyzing of parameters relative to traffic conditions based on the source of data
    • G08G1/0112Measuring and analyzing of parameters relative to traffic conditions based on the source of data from the vehicle, e.g. floating car data [FCD]
    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/01Detecting movement of traffic to be counted or controlled
    • G08G1/0104Measuring and analyzing of parameters relative to traffic conditions
    • G08G1/0125Traffic data processing
    • G08G1/0133Traffic data processing for classifying traffic situation
    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/01Detecting movement of traffic to be counted or controlled
    • G08G1/052Detecting movement of traffic to be counted or controlled with provision for determining speed or overspeed

Abstract

A road traffic open degree judging method based on a K-means cluster includes the steps of clearing away invalid data through a basic data threshold, so data distribution is more concentrated after data is primarily cleared away, further processing of data is easy, efficiency of data analysis is improved and result precision is enhanced; obtaining initial K values according to the highest limit speed of a road to correspond to different traffic states respectively; clustering different vehicle-mounted GPS data to different K value points through the K-means cluster, analyzing the number of data contained by each K value point, and making the maximum number serve as the index of a current road traffic state. By clustering and analyzing the vehicle-mounted GPS speed data, road network traffic state information is collected, implementation cost is low and technological difficulty is small.

Description

A kind of road traffic unimpeded degree determination methods based on K-means cluster
Technical field
The present invention relates to a kind of the coast is clear degree judge, computer application field, in particular, Yi Zhongji Road traffic unimpeded degree determination methods in K-means cluster.
Background technology
Recently as the quickening of Urbanization in China, the surge of vehicles number in city, urban environment shape Condition is the severeest, cannot alleviate current traffic pressure by urban road infrastructure construction merely.Intelligence is handed over It is connected with and helps ensure the unimpeded of intersection, reach to alleviate the purpose of traffic congestion, and then solve to a certain extent Determine the social problems such as urban traffic accident, traffic jam, environmental pollution and energy resource consumption, wherein road traffic shape The collection of state information is the basis realizing intelligent transportation.
City real vehicles driving road-condition is complicated, vehicle travel process can comprise acceleration, deceleration, at the uniform velocity, quiet Only four kinds of states, the data gathered are then the set of these three data, in some data in particular cases also Cannot function as judging when data exception that the fault of the foundation of road conditions, such as collecting device causes, personnel get off Stationary vehicle data (speed is zero), vehicle peccancy are driven over the speed limit (speed is much larger than speed limit), vehicle is promptly made Speed data (such as running into pedestrian to bring to a halt) during Dong, because the change of these in particular cases car speeds Change is not that the odjective cause of traffic causes, it is therefore necessary to some the most invalid data rejected.
Same road vehicle quantity is more, and the velocity information speed gathered differs, can not direct reaction Go out the traffic status of road, the vehicle speed information of same road can be carried out by K-means clustering algorithm Sort out, thus judge this road current state in certain time period.Represent that the description of traffic behavior generally includes Traffic density, passage rate, occupation rate etc., although these factors can reflect the current shape of a road State, but traffic density and occupation rate etc. are difficult to measure.
Summary of the invention
In order to overcome by the method such as traditional video monitoring obtain road traffic state information implementation cost higher and The deficiency that technical difficulty is bigger, invention herein provides a kind of unobstructed degree of road traffic based on K-means cluster to sentence Disconnected method, by carrying out cluster analysis to vehicle GPS speed data, it is achieved road grid traffic passes through status information Gathering, implementation cost is relatively low and technical difficulty is less.
The technical solution adopted for the present invention to solve the technical problems is:
A kind of road traffic unimpeded degree determination methods based on K-means cluster, it is characterised in that: described detection Method comprises the following steps::
1) cleaning invalid data:
1.1) if car speed is that the data of zero are less than the 10% of total data on certain moment road, then it is assumed that vehicle Stopping is not that traffic congestion causes, it is judged that this data invalid;
1.2) if car speed is much larger than road speed limit, then it is assumed that this rule-breaking vehicle is driven, or collecting device Breaking down and cause gathering data exception, this type of data judging is invalid;
2) K value is selected:
Assuming that road exists N kind state, wherein the representation speed of every kind of state is Vn, by this N kind state Representation speed be inserted into initial data and concentrate that to carry out K-means cluster calculation as initial K value point quick Data set is gathered N number of apoplexy due to endogenous wind, then the n-th K value KnChoose as follows;
Kn=Vn, n=1,2,3 N (1)
Wherein KnValue along with the difference of road may be different, with the Maximum speed limit V of roadmaxFor upper Limit, carries out speed limit decile calculating, and takes speed V of correspondence positionnAs KnValue, speed VnAccording to Equation below obtains:
V n = V m a x × ( n - 1 ) N , n = 1 , 2 , 3 ... N - - - ( 2 )
3) cluster calculation:
3.1) N number of central point K is first selected1,K2,K3···KN, the most each K is worth exploitation such as formula (1) Shown in, they represent the corresponding current state of road respectively, and arranging iterations is g;
3.2) velocity information is traveled through successively, and the difference size of versus speed and central point, it is classified as difference minimum Apoplexy due to endogenous wind, it is judged that formula is as follows:
Si-Kn=min{ | Si-K1|,|Si-K2|,|Si-K3|···|Si-KN|, n=1,2,3 N (3)
Wherein SiRepresent i-th speed data, KnTable the n-th K value, if meeting above-mentioned condition, then S is describediBelong to the n-th class;
3.3) after iterative computation g time, adding up the sample size that each apoplexy due to endogenous wind comprises successively, the maximum class of quantity is wrapped
The K value point contained represents the current state of road.
The design thinking of the present invention is: same road vehicle quantity is more, and the velocity information speed gathered Differ, the traffic status of road can not be gone out by direct reaction, by K-means clustering algorithm can by with along with The vehicle speed information on road is sorted out, thus judges that the traffic in certain time period of this road is passed through state.
Beneficial effects of the present invention is mainly manifested in: carry out K-means cluster by gps data vehicle-mounted to road Analyze, vehicle speed information is converted into road conditions information more intuitively, solves road vehicle GPS number According to state dispersion can not the problem of direct reaction road conditions, the simplest K value selection has higher meter Calculate efficiency, the requirement to real-time in actual application can be met.
Accompanying drawing explanation
Fig. 1 is whole speed scatterplot.
Fig. 2 is the scatterplot after data scrubbing.
Fig. 3 is the flow chart of speed cluster.
Fig. 4 is cluster result schematic diagram.
Fig. 5 is the flow chart of a kind of road traffic unimpeded degree determination methods based on K-means cluster.
Detailed description of the invention
The present invention will be further described below in conjunction with the accompanying drawings.
With reference to Fig. 1~Fig. 5, a kind of road traffic unimpeded degree determination methods based on K-means cluster, including with Lower step:
1) cleaning invalid data:
Certain road is vehicle speed information data instance in 15 minutes, speed scatterplot distributed intelligence as it is shown in figure 1, Data that car speed on certain moment road is zero data less than the 10% of total data are rejected, by car Speed is rejected more than the data of road speed limit;
Result after cleaning as in figure 2 it is shown, wherein figure orbicular spot be invalid data, be placed in the figure lower left corner Because the data that speed is zero and speed is zero are judged as invalid data less than the 10% of total data, and It is placed in the speed in the figure upper right corner because being considered as hypervelocity or equipment event more than the speed limit of common road The invalid data that causes of barrier, after data are tentatively cleaned, the distribution of data is more concentrated, beneficially data enter one Step processes, and improves the efficiency of data analysis and strengthens result precision;
2) K value is selected:
Assuming that road exists N kind state, wherein the representation speed of every kind of state is Vn, the n-th K value Choose as follows:
Kn=Vn, n=1,2,3 N (1)
Maximum speed limit V with roadmaxFor the upper limit, speed limit is carried out decile calculating, and takes the speed of correspondence position Degree VnAs KnValue, speed VnCan obtain according to equation below:
V n = V m a x × ( n - 1 ) N , n = 1 , 2 , 3 ... N - - - ( 2 )
3) cluster calculation:
Cluster calculation process is as shown in Figure 3;
3.1) N number of central point K is first selected1,K2,K3···Kn, the most each K is worth exploitation such as formula (1) Shown in, they represent the corresponding current state of road respectively, and arranging iterations is g;
3.2) velocity information is traveled through successively, and the difference size of versus speed and central point, it is classified as difference minimum Apoplexy due to endogenous wind, it is judged that formula is as follows:
Si-Kn=min{ | Si-K1|,|Si-K2|,|Si-K3|···|Si-KN|, n=1,2,3 N (3)
Wherein SiRepresent i-th speed data, KnTable the n-th K value, if meeting above-mentioned condition, then S is describediBelong to the n-th class;
3.4) after iterative computation g time, adding up the sample size that each apoplexy due to endogenous wind comprises successively, the maximum class of quantity is wrapped The K value point contained represents the current state of road;
Cluster result is as shown in Figure 4;
56 vehicle GPS information in certain road 15min are carried out cluster analysis by this example, a kind of based on The road traffic unimpeded degree determination methods of K-means cluster, comprises the following steps:
1) cleaning invalid data:
The data that car speed on this road is zero are less than the 10% of total data by certain road speed limit 50km/h Data reject, by car speed more than road speed limit data reject, after cleaning, valid data are 44. 2) K value is selected:
Certain road speed limit is 50km/h, selects K1,K2,K3,K4Four values carry out cluster fortune as central point Calculate, represent respectively actual traffic road conditions be impassable, block up, normal, unobstructed four kinds of situations, pass through Judge four kinds of sizes clustered to judge the current state of road, four K values are the most as follows:
K 1 = 0 k m / h K 2 = 50 4 = 12.5 k m / h K 3 = V m a x 2 = 25 k m / h K 4 = 3 × V max 4 37.5 k m / h - - - ( 1 )
3) cluster calculation:
3.1) 4 central point K are first selected1,K2,K3,K4, the most each K is worth exploitation such as formula (1) Shown in, they represent the corresponding current state of road respectively, and taking iterations g is 100;
3.2) velocity information is traveled through successively, and the difference size of versus speed and central point, it is classified as difference minimum Apoplexy due to endogenous wind, it is judged that formula is as follows:
Si-Kn=min{ | Si-K1|,|Si-K2|,|Si-K3|···|Si-KN|, n=1,2,3 N (2)
Wherein SiRepresent i-th speed data, KnTable the n-th K value, if the condition of meeting a cassation, then S is describediBelong to the n-th class;
3.3) iteration 100 times, complete cluster calculation;
Cluster result data are as shown in table 1 below, the most altogether comprise vehicle speed information 44, Qi Zhongju Class is respectively 1,2,38 and 3 to the quantity of corresponding K value point, therefore can draw this road Road within this time period of 15 minutes the speed of most of vehicles all within normal range, therefore this road Road conditions are normal.
Table 1
The excellent effect of optimization that the embodiment that the present invention is given that described above is shows, it is clear that this Bright it is not only suitable for above-described embodiment, without departing from essence spirit of the present invention and without departing from involved by flesh and blood of the present invention And it can be done many variations on the premise of content and be carried out.

Claims (1)

1. a road traffic unimpeded degree determination methods based on K-means cluster, it is characterised in that: described detection Method comprises the following steps::
1) cleaning invalid data:
If on 1.1 certain moment road, car speed is that the data of zero are less than the 10% of total data, then it is assumed that vehicle Stopping is not that traffic congestion causes, it is judged that this data invalid;
If 1.2 car speeds are much larger than road speed limit, then it is assumed that this rule-breaking vehicle is driven, or collecting device Breaking down and cause gathering data exception, this type of data judging is invalid;
2) K value is selected:
Assuming that road exists N kind state, wherein the representation speed of every kind of state is Vn, by this N kind state Representation speed be inserted into initial data and concentrate that to carry out K-means cluster calculation as initial K value point quick Data set is gathered N number of apoplexy due to endogenous wind, then the n-th K value KnChoose as follows;
Kn=Vn, n=1,2,3 N (1)
Wherein KnValue along with the difference of road may be different, with the Maximum speed limit V of roadmaxFor upper Limit, carries out speed limit decile calculating, and takes speed V of correspondence positionnAs KnValue, speed VnAccording to Equation below obtains:
V n = V m a x × ( n - 1 ) N , n = 1 , 2 , 3 ... N - - - ( 2 )
3) cluster calculation:
3.1) N number of central point K is first selected1,K2,K3···KN, the most each K is worth exploitation such as formula (1) Shown in, they represent the corresponding current state of road respectively, and arranging iterations is g;
3.2) velocity information is traveled through successively, and the difference size of versus speed and central point, it is classified as difference minimum Apoplexy due to endogenous wind, it is judged that formula is as follows;
Si-Kn=min{ | Si-K1|,|Si-K2|,|Si-K3|···|Si-KN|, n=1,2,3 N (3)
Wherein SiRepresent i-th speed data, KnTable the n-th K value, if meeting above-mentioned condition, then S is describediBelong to the n-th class;
3.3) after iterative computation g time, adding up the sample size that each apoplexy due to endogenous wind comprises successively, the maximum class of quantity is wrapped The K value point contained represents the current state of road.
CN201610304571.3A 2016-05-10 2016-05-10 A kind of unimpeded degree judgment method of road traffic based on K-means clusters Active CN105869398B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610304571.3A CN105869398B (en) 2016-05-10 2016-05-10 A kind of unimpeded degree judgment method of road traffic based on K-means clusters

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610304571.3A CN105869398B (en) 2016-05-10 2016-05-10 A kind of unimpeded degree judgment method of road traffic based on K-means clusters

Publications (2)

Publication Number Publication Date
CN105869398A true CN105869398A (en) 2016-08-17
CN105869398B CN105869398B (en) 2018-10-23

Family

ID=56630624

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610304571.3A Active CN105869398B (en) 2016-05-10 2016-05-10 A kind of unimpeded degree judgment method of road traffic based on K-means clusters

Country Status (1)

Country Link
CN (1) CN105869398B (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106740829A (en) * 2017-03-23 2017-05-31 吉林大学 Based on the double semi-dragging truck riding stability automatic identifications of cluster analysis and early warning system
CN107424410A (en) * 2017-07-14 2017-12-01 中南大学 A kind of accident detection method calculated based on route travel time
CN108446432A (en) * 2018-02-06 2018-08-24 浙江工业大学 A kind of virtual bicycle rider based on model rides the computational methods of speed
CN110288824A (en) * 2019-05-20 2019-09-27 浙江工业大学 Based on Granger causality road network morning evening peak congestion and mechanism of transmission analysis method
TWI676397B (en) * 2018-10-18 2019-11-01 中華電信股份有限公司 Artificial intelligence traffic estimation system using mobile network signaling data and method thereof
CN111402581A (en) * 2020-03-11 2020-07-10 浙江大华技术股份有限公司 Traffic state detection method and detection device

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102360526A (en) * 2011-09-28 2012-02-22 东南大学 Real-time monitoring method for road section state of high road
US20120148094A1 (en) * 2010-12-09 2012-06-14 Chung-Hsien Huang Image based detecting system and method for traffic parameters and computer program product thereof
CN103150900A (en) * 2013-02-04 2013-06-12 合肥革绿信息科技有限公司 Traffic jam event automatic detecting method based on videos
US8970701B2 (en) * 2011-10-21 2015-03-03 Mesa Engineering, Inc. System and method for predicting vehicle location
CN105261217A (en) * 2015-10-03 2016-01-20 上海大学 Method for detecting urban traffic congestion state by using density-based clustering algorithm

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120148094A1 (en) * 2010-12-09 2012-06-14 Chung-Hsien Huang Image based detecting system and method for traffic parameters and computer program product thereof
CN102360526A (en) * 2011-09-28 2012-02-22 东南大学 Real-time monitoring method for road section state of high road
US8970701B2 (en) * 2011-10-21 2015-03-03 Mesa Engineering, Inc. System and method for predicting vehicle location
CN103150900A (en) * 2013-02-04 2013-06-12 合肥革绿信息科技有限公司 Traffic jam event automatic detecting method based on videos
CN105261217A (en) * 2015-10-03 2016-01-20 上海大学 Method for detecting urban traffic congestion state by using density-based clustering algorithm

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
倪升华: "基于数据的交通拥堵评价与预测方法", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106740829A (en) * 2017-03-23 2017-05-31 吉林大学 Based on the double semi-dragging truck riding stability automatic identifications of cluster analysis and early warning system
CN106740829B (en) * 2017-03-23 2018-04-10 吉林大学 Based on the double semi-dragging truck riding stability automatic identifications of cluster analysis and early warning system
CN107424410A (en) * 2017-07-14 2017-12-01 中南大学 A kind of accident detection method calculated based on route travel time
CN108446432A (en) * 2018-02-06 2018-08-24 浙江工业大学 A kind of virtual bicycle rider based on model rides the computational methods of speed
CN108446432B (en) * 2018-02-06 2021-12-17 浙江工业大学 Virtual bicycle rider riding speed calculation method based on model
TWI676397B (en) * 2018-10-18 2019-11-01 中華電信股份有限公司 Artificial intelligence traffic estimation system using mobile network signaling data and method thereof
CN111161529A (en) * 2018-10-18 2020-05-15 中华电信股份有限公司 Artificial intelligent traffic flow estimation system and method using mobile network signaling data
CN111161529B (en) * 2018-10-18 2021-06-08 中华电信股份有限公司 Artificial intelligent traffic flow estimation system and method using mobile network signaling data
CN110288824A (en) * 2019-05-20 2019-09-27 浙江工业大学 Based on Granger causality road network morning evening peak congestion and mechanism of transmission analysis method
CN111402581A (en) * 2020-03-11 2020-07-10 浙江大华技术股份有限公司 Traffic state detection method and detection device

Also Published As

Publication number Publication date
CN105869398B (en) 2018-10-23

Similar Documents

Publication Publication Date Title
CN105869398A (en) Road traffic open degree judging method based on K-means cluster
CN105913661B (en) A kind of express highway section traffic state judging method based on charge data
CN105868787A (en) Electric car driving range evaluation method based on working condition identification and fuzzy energy consumption
CN103839409B (en) Based on the traffic flow modes method of discrimination of multibreak facial vision sensing cluster analysis
CN103258430B (en) Road traveling time calculating and traffic road condition judging method and road traveling time calculating and traffic road condition judging device
CN103413263B (en) Travel time index entropy traffic circulation evaluation method
WO2019047905A1 (en) Road traffic analysis system, method and apparatus
CN104408925B (en) Crossing evaluation of running status method based on display radar
CN104484994B (en) Urban road network traffic based on array radar runs index number evaluation method
CN105355049B (en) A kind of highway evaluation of running status method based on macroscopical parent map
CN103903433A (en) Real-time dynamic judgment method and device for road traffic state
CN107240264B (en) A kind of non-effective driving trace recognition methods of vehicle and urban road facility planing method
CN103646542A (en) Forecasting method and device for traffic impact ranges
CN102360525A (en) Discriminant analysis-based high road real-time traffic accident risk forecasting method
CN101739825A (en) GPS floating vehicle-based traffic data fault identification and recovery method
CN101783074A (en) Method and system for real-time distinguishing traffic flow state of urban road
CN106529608A (en) Comprehensive stationing system for exhaust gas remote-measuring equipment of motor vehicle
CN103903441A (en) Road traffic state distinguishing method based on semi-supervised learning
CN103413435A (en) Road network running condition analysis method based on fuel consumption acquisition
CN107180534B (en) The express highway section average speed estimation method of support vector regression fusion
CN103413436A (en) Road network running condition analysis system based on fuel consumption acquisition
CN106803347A (en) Urban intersection traffic state judging method based on RFID data
CN114463868B (en) Traffic flow combination prediction method and system for toll station facing traffic flow control
CN107833459A (en) A kind of city bus operation conditions evaluation method based on gps data
CN111507879A (en) Urban space structure identification method

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