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 PDFInfo
- 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
Links
Classifications
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/01—Detecting movement of traffic to be counted or controlled
- G08G1/0104—Measuring and analyzing of parameters relative to traffic conditions
- G08G1/0108—Measuring and analyzing of parameters relative to traffic conditions based on the source of data
- G08G1/0112—Measuring and analyzing of parameters relative to traffic conditions based on the source of data from the vehicle, e.g. floating car data [FCD]
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/01—Detecting movement of traffic to be counted or controlled
- G08G1/0104—Measuring and analyzing of parameters relative to traffic conditions
- G08G1/0125—Traffic data processing
- G08G1/0133—Traffic data processing for classifying traffic situation
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/01—Detecting movement of traffic to be counted or controlled
- G08G1/052—Detecting 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
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:
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:
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:
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:
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.
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)
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)
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 |
-
2016
- 2016-05-10 CN CN201610304571.3A patent/CN105869398B/en active Active
Patent Citations (5)
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)
Title |
---|
倪升华: "基于数据的交通拥堵评价与预测方法", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》 * |
Cited By (10)
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 |