CN102289932B - Dynamic OD (Origin Destination) matrix estimating method based on AVI (Automatic Vehicle Identification) device - Google Patents

Dynamic OD (Origin Destination) matrix estimating method based on AVI (Automatic Vehicle Identification) device Download PDF

Info

Publication number
CN102289932B
CN102289932B CN 201110163206 CN201110163206A CN102289932B CN 102289932 B CN102289932 B CN 102289932B CN 201110163206 CN201110163206 CN 201110163206 CN 201110163206 A CN201110163206 A CN 201110163206A CN 102289932 B CN102289932 B CN 102289932B
Authority
CN
China
Prior art keywords
path
vehicle
avi
information
mean
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.)
Active
Application number
CN 201110163206
Other languages
Chinese (zh)
Other versions
CN102289932A (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.)
Tongji University
Original Assignee
Tongji University
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 Tongji University filed Critical Tongji University
Priority to CN 201110163206 priority Critical patent/CN102289932B/en
Publication of CN102289932A publication Critical patent/CN102289932A/en
Application granted granted Critical
Publication of CN102289932B publication Critical patent/CN102289932B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Traffic Control Systems (AREA)

Abstract

The invention belongs to the field of traffic plan and management, in particular relating to a dynamic OD (Origin Destination) matrix estimating method based on an AVI (Automatic Vehicle Identification) device. The method comprises the following steps of: introducing a part of path information detected through AVI, dynamic travel time information and detector measurability criterion; performing lost path range determination and selection probability modification on the path-lost vehicle information according to a Bayes estimation algorithm firstly; simulating any vehicle by utilizing monte carlo random simulation and selecting lost paths so as to obtain an initially modified OD matrix based on a part of path of an individual vehicle; and finally correcting the initially modified OD matrixby utilizing the flow information of the AVI detection to obtain a final OD matrix estimation value. Through the method disclosed by the invention, the defects that the dependency on the priori information is high and aspects, such as detection precision and the like, are never considered because the road section flow and the travel time information are only considered in the traditional method can be overcome.

Description

Dynamic OD Matrix Estimation method based on automatic vehicle identification equipment
Technical field
The invention belongs to traffic programme and management domain, be specifically related to a kind of Dynamic OD Matrix Estimation method based on automatic vehicle identification equipment.
Background technology
Vehicle driving OD matrix is the key foundation information of traffic system planning, design and operation management.The precision of OD information directly affects urban traffic control person, and traffic programme engineering technical personnel the accurate judgement of situation may occur for current situation of traffic and future transportation, and therefore can directly affect the validity of traffic management measure, the rationality of traffic programme.Therefore accurate OD matrix no matter for urban traffic control personnel and Urban Traffic Planning personnel's final decision with according to being vital.
Vehicle driving OD matrix is to describe one group of math matrix that in transportation network, vehicle is reached home from starting point arbitrarily, and it has directly reacted between different districts set out vehicle number and arrival vehicle number.Yet obtaining of OD matrix is all a difficulties in traffic administration and control all the time.The acquisition methods of traditional OD matrix normally carries out large-scale vehicle driving sample survey, yet this method exists, technique for investigation difficulty, funds cost are huge, the shortcoming such as long that expends time in, and the traffic data obtained exists post-processed loaded down with trivial details and can not be applied to the dynamic management of urban transportation.Therefore for fear of the problem existed in traditional OD matrix acquisition methods, from 1978 s, Van Zuylen and Willumsen utilize the vehicle flow data of the magnetic test coil acquisition be laid in road to carry out the research of OD matrix estimation method.Up to now, the method for estimation of OD matrix has mainly comprised the types such as least square method, state-space method, Information Theory Model.Yet these methods are subject to restriction and the restriction of detection means, the Main Analysis data of the highway section vehicle flow data in the detecting device as research have only been utilized, therefore applying traditional " section type " detecting device carries out the OD Matrix Estimation to be subject to artificial assumed condition too much, the serious restriction of the conditions such as detection information is few, and the checkout equipment precision is low.
In recent years, along with automatic vehicle identification (Automatic Vehicle Identification, AVI) technology and equipment are at the propagation and employment in China one line city, and the traffic information collection technology turns to " wide area type " detection technique from traditional " section type " detection technique rapidly.The core of automatic vehicle identification technology is to detect vehicle ID(license plate number), by time and vehicle position information.Based on existing documents and materials are studied, the method for at present based on the automatic vehicle identification technology, carrying out the OD Matrix Estimation is only all a kind of improvement of traditional OD matrix estimation method.
These classic methods exist following problem and challenge:
(1) for the dynamic OD estimation under the AVI environment, current method remains by classical OD estimation model is improved, and adds new AVI detection information to improve the OD estimated accuracy.In fact, the precision that the routing information that AVI detects is estimated OD is most important.
(2) OD estimates not only and network topology, link flow is relevant, also with the precision of prior imformation, close relationship is arranged, and research in the past all supposes that having obtained reliable prior imformation calculates OD usually, and long due to the OD survey interval in reality, priori OD information often precision is not high.
(3) in the OD that detects information based on AVI estimates, do not consider AVI accuracy of detection problem.
Summary of the invention
The object of the invention is to the deficiency existed in the research and technology for conventional dynamic OD matrix estimation method, particularly for AVI information excavating degree of depth deficiency, the problem that the prior imformation dependence is excessively strong, proposed a kind of new Dynamic OD Matrix Estimation method based on automatic vehicle identification equipment.
This method has following four characteristics: one, broken through and only relied on the main information that link flow and journey time are estimated as OD in the classic method; Two, propose the measurability criterion of detecting device, be specially adapted to the accuracy of detection problem of AVI checkout facility outwardness; Three, do not require for the road network topology structure, go for any type of open road network; Four, solve the dependence of classic method to prior imformation, can under random prior imformation, extrapolate the OD of degree of precision.
For reaching above target, the present invention proposes the Dynamic OD Matrix Estimation method based on automatic vehicle identification equipment.At first each is estimated to the vehicle license data that obtain in each checkout equipment in the period, vehicle due in data, the detecting device numbering is extracted; Then will extract the data that obtain and take the vehicle license data as according to carrying out Data classification, and convert the vehicle sections routing information to according to vehicle license data and detecting device numbering, and carry out Data classification according to the building form of vehicle sections routing information; Then take priori OD information as basis, according to the installation position of detecting device, based on detecting device measurability criterion, reduce the scope of optional residual paths and the selection probability of residual paths is revised simultaneously; The process that further Actual path of the method simulating vehicle of employing random simulation is selected, to obtain the fullpath of vehicle; And the initial OD matrix of complete rear acquisition is repaired in the path that gathers all vehicles that obtain, the link flow finally obtained in conjunction with the AVI detection carries out the initial OD matrix and revises.
AVI measurability (Measurability) concept definition that said method is mentioned is: the characteristics that lack based on AVI detection information data and the dynamic travel time information that has detected vehicle, analyze the possibility of judgement vehicle through a certain path.According to the result of measurability criterion, can greatly dwindle the range of choice of vehicle residual paths, the reliability that raising vehicle residual paths judges and then the precision that significantly improves dynamic OD estimation.Concrete steps are as follows:
(1) extraction of AVI information data and classification
The information that AVI is obtained is usingd the vehicle license data as index data, by vehicle license information, identical data are mated combination, and according to the quantity of information of arbitrary same vehicle licence plate, information of vehicles is divided into to following 3 large classes: 1, comprise the vehicle starting point, the data of all path nodes of terminal and process are called path omniscient type vehicle data; 2, the data of two path nodes that at least comprise the process of vehicle are called path and partly know the type vehicle data; 3, the data that only comprise a path node of vehicle process are called path and singly know the type vehicle data.
(2) the expansion sample of prior imformation
The vehicle number that the AVI of take catches is basis, and prior imformation is expanded to the sample processing according to the vehicle number of actual acquisition.
(3) filtration path omniscient type vehicle data
Convert the vehicle data of path omniscient type to the OD matrix information, and this OD matrix information is subtracted each other with the prior imformation after the expansion sample (if in prior imformation, corresponding OD is less than 0 to flow, give a minimum flow 1), obtain the priori OD information after filtering.
(4) estimate the path journey time in the period
According to the laying information of different AVI, the journey time of the vehicle of two groups of AVI processes of arbitrary neighborhood is added up and obtained its average travel time; According to the different vehicle journey time, distribute in addition, the principle of Based on Probability opinion, set up AVI journey time probability distribution function, and by the possibility of all the other these two groups of AVI detecting devices of vehicles faileds process of this Functional Analysis; Then pass through highway section-path journey time relation function, the journey time of way to acquire, and using this as judging whether vehicle is estimating the criterion that enters road network in the period.
(5) path estimation of type vehicle is partly known in path
The part path of path partly being known to type vehicle disappearance is divided into upstream miss path and downstream miss path by direction of traffic.According to different miss path, with reference to detecting device measurability foundation, by the Bayesian Estimation algorithm, dwindle the optional miss path of any vehicle hunting zone.And on the basis of filtering priori OD information, the selection probability of possible path is revised; Finally by the Monte Carlo random simulation, any vehicle is made to routing in its feasible path, and carry out the path reparation with this.The routing information of type vehicle is partly known in acquired any one group of path, can pass through the path journey time, and the judgement vehicle enters the time of road network and leaves the time of road network.
(6) secondary filtration of priori OD information
Path after repairing is partly known to the type data-switching becomes the OD matrix information, and this OD matrix information is subtracted each other to (if in prior imformation, corresponding OD is less than 0 to flow, giving a minimum flow 1) with the prior imformation after the expansion sample.Then according to remaining path, singly know that the type vehicle number carries out secondary and expands sample, finally obtain the priori OD information after secondary filtration.
(7) path estimation of type vehicle is singly known in path
Path singly knows that the method that the type vehicle route is estimated partly knows that with path the method for type vehicle route estimation is similar, specifically can partly know type vehicle route method of estimation referring to path, and obtain singly knowing the OD matrix information of type information of vehicles with this.
(8) initially repairing the OD matrix obtains
The OD matrix information that front 3 classes are obtained adds up to obtain initial reparation OD matrix.
(9) initially repair the OD matrix information based on the flow adjustment
The data based highway section of link flow of OD matrix information and the AVI information interception of initially reparation-path flow funtcional relationship is dynamically strolled to matching.When path flow and highway section measured discharge relative error are 5%, stop strolling matching, and as final OD matrix information.
In the present invention, the path estimation of type vehicle and the path estimation that the type vehicle is singly known in the middle path of step (7) are partly known in path described in step (5), are specially:
Arrive next highway section detecting device when a certain vehicle has grace time, but the routing information of vehicle is not while comprising next highway section detecting device, has any in following situation:
(1), vehicle passes through next highway section, but is not recorded owing to detecting error;
(2), vehicle, through laying the highway section of detecting device, does not directly arrive the highway section that sensorless is laid;
To the possible residual paths of vehicle, set up initial routing probability set according to prior imformation,
Figure 15878DEST_PATH_IMAGE001
,then analyze vehicle according to following formula and can within this period, pass through next highway section, and dwindle the range of choice to residual paths with this;
Figure 577440DEST_PATH_IMAGE002
(1.1)
Figure 532233DEST_PATH_IMAGE003
(1.2)
Figure 583366DEST_PATH_IMAGE004
(1.3)
Wherein
Figure 383963DEST_PATH_IMAGE005
: mean the journey time between any two AVI;
Figure 545954DEST_PATH_IMAGE006
: mean the time of i car through arbitrary neighborhood AVI detecting device;
Figure 748396DEST_PATH_IMAGE007
: mean the distance between any two adjacent AVI detecting devices;
Figure 967500DEST_PATH_IMAGE008
: mean the distance between two nodes at place between two AVI detecting devices;
Figure 989814DEST_PATH_IMAGE009
: mean the stochastic error in journey time calculating.
Figure 955496DEST_PATH_IMAGE010
: mean that h is in the period, the selection probability in i the journey time interval by two AVI detecting devices;
Figure 12445DEST_PATH_IMAGE011
: mean that h is in the period, the vehicle number in i the journey time interval by two AVI detecting devices;
Figure 343063DEST_PATH_IMAGE012
: mean that h is in the period, by the vehicle fleet of two AVI detecting devices;
Figure 115323DEST_PATH_IMAGE013
: mean that h is in the period, be numbered the vehicle of j by the journey time of two AVI detecting devices;
Figure 556799DEST_PATH_IMAGE014
: mean that h, in the period, detects the vehicle average travel time according to two AVI detecting devices;
S: mean random chance numerical value;
N: mean time interval length, i.e. unit interval length;
Figure 530572DEST_PATH_IMAGE015
: mean the stochastic error of calculating;
For the residual paths that can select, according to the selection probability of following formula revised residual paths, obtain new routing probability set;
Figure 32091DEST_PATH_IMAGE016
(2.1)
Figure 91314DEST_PATH_IMAGE017
(2.2)
Figure 333552DEST_PATH_IMAGE018
: the prior imformation that means n group miss path in the travel direction of downstream;
Figure 365093DEST_PATH_IMAGE019
: the prior imformation that means n group miss path in the travel direction of upstream;
Figure 99831DEST_PATH_IMAGE020
: mean to meet the total j group of the AVI detecting device that can survey condition in the travel direction of n group miss path downstream;
Figure 584033DEST_PATH_IMAGE021
: mean to meet the total k group of the AVI detecting device that can survey condition in the travel direction of n group miss path upstream;
Figure 429629DEST_PATH_IMAGE022
: be illustrated in the moment that vehicle i is detected for the last time;
Figure 336185DEST_PATH_IMAGE023
: be illustrated in the moment that vehicle i is detected for the first time;
: the zero hour that means to estimate the period;
: the length that means to estimate the period.
Compared with prior art, the present invention has broken through usings the defect of link flow as dynamic OD estimation main information foundation on the traditional sense, repair the Main Means as dynamic OD estimation by usining the microcosmic path of vehicle, take link flow information as auxiliary correction means, make AVI information be fully used.The method is any open road network of the extraordinary adaptation of energy simultaneously, and can under low AVI coverage rate condition, improve the lower OD prior imformation of precision.This is not enough for domestic OD survey information, precision is not high actually has a considerable realistic meaning, China's traffic administration with control reality in boundless application space is arranged.
The accompanying drawing explanation
Fig. 1 is the Dynamic OD Matrix Estimation method flow diagram that the present invention proposes.
Schematic diagram is repaired in the path that Fig. 2 is the open road network based on the measurability criterion that proposes of the present invention.
Fig. 3 is the road network figure that the embodiment of the present invention 1 adopts.
Fig. 4 is the OD Matrix Estimation result schematic diagram under lower accuracy coverage rate and prior imformation prerequisite in the embodiment of the present invention 1.
Fig. 5 be in the embodiment of the present invention 1 to precision the improved effect schematic diagram of lower initial OD.
Embodiment
Below in conjunction with 3 pairs of embodiments of the invention of accompanying drawing, elaborate: the present embodiment is implemented take technical solution of the present invention under prerequisite, provided detailed embodiment and concrete operating process, but protection scope of the present invention is not limited to following embodiment.
Embodiment 1: this method of estimation is for certain the ground velocity road system shown in Fig. 3, this urban expressing system has the AVI equipment of 17 gateways and 9 sections, by the right number of the AVI facility detectable OD of actual covering, be wherein 216 groups, account for all OD to 74.74% of number.The equipment input needed: at road section or ring road gateway, lay the AVI video detector.The acquisition input message that requires: the vehicle license information after identification, vehicle due in, AVI detecting device numbering.
After obtaining above-mentioned input message, take certain path and partly know that repair in the path of type vehicle data is example, repair the driving path of any vehicle according to Fig. 2, step is as follows:
(1) when the known vehicle licence plate vehicle that is A respectively by being numbered AVI 1and AVI 4the time, can obtain respectively vehicle A at the T time of arrival of these two groups of detecting devices 1and T 4.And, according to the topological structure of road network, analyze the optional residual paths of vehicle A and also with the prior imformation in different paths, it is demarcated
Figure 800292DEST_PATH_IMAGE026
.
(2) calculate vehicle A by path journey time function and arrive AVI 3journey time.
(1)
Figure 880036DEST_PATH_IMAGE028
(2)
Figure 135568DEST_PATH_IMAGE029
(3)
Wherein
Figure 526229DEST_PATH_IMAGE030
: mean the journey time between any two AVI;
Figure 449186DEST_PATH_IMAGE031
: mean the time of i car through arbitrary neighborhood AVI detecting device;
Figure 634311DEST_PATH_IMAGE032
: mean the distance between any two adjacent AVI detecting devices;
Figure 374209DEST_PATH_IMAGE033
: mean the distance between two nodes at place between two AVI detecting devices;
Figure 302982DEST_PATH_IMAGE009
: mean the stochastic error in journey time calculating;
: mean that h is in the period, the selection probability in i the journey time interval by two AVI detecting devices;
Figure 436471DEST_PATH_IMAGE011
: mean that h is in the period, the vehicle number in i the journey time interval by two AVI detecting devices;
Figure 932175DEST_PATH_IMAGE012
: mean that h is in the period, by the vehicle fleet of two AVI detecting devices;
Figure 396130DEST_PATH_IMAGE013
: mean that h is in the period, be numbered the vehicle of j by the journey time of two AVI detecting devices;
: mean that h, in the period, detects the vehicle average travel time according to two AVI detecting devices;
S: mean random chance numerical value;
N: mean time interval length, i.e. unit interval length;
Figure 820606DEST_PATH_IMAGE015
: mean the stochastic error of calculating.
(3) according to detecting device measurability criterion, when having grace time, vehicle A arrives next highway section detecting device AVI3, but the routing information of vehicle A inclusion test device AVI not 3the time, there are two kinds of sights:
Sight 1: vehicle A passes through AVI 5, but be not recorded owing to detecting error.
Sight 2: vehicle is not through laying AVI 5highway section, and be short to most and reach D based on travel time 6.
Figure 538026DEST_PATH_IMAGE034
(4)
(5)
Figure 560657DEST_PATH_IMAGE018
: the prior imformation that means n group miss path in the travel direction of downstream;
: the prior imformation that means n group miss path in the travel direction of upstream;
: mean to meet the total j group of the AVI detecting device that can survey condition in the travel direction of n group miss path downstream;
Figure 472222DEST_PATH_IMAGE021
: mean to meet the total k group of the AVI detecting device that can survey condition in the travel direction of n group miss path upstream;
Figure 78784DEST_PATH_IMAGE022
: be illustrated in the moment that vehicle i is detected for the last time;
: be illustrated in the moment that vehicle i is detected for the first time;
Figure 636597DEST_PATH_IMAGE024
: the zero hour that means to estimate the period;
: the length that means to estimate the period.
Therefore can revise the probability of these optional residual paths
Figure 975622DEST_PATH_IMAGE035
.
(4) revised routing probability is approached to the process that actual vehicle route is selected by the method maximum of Monte Carlo random simulation, obtain the fullpath of vehicle.
(5) the incomplete vehicle data in all the other paths can carry out the reparation in path and convert the OD matrix of initial reparation to by above-mentioned method, and repairs initial matrix by highway section-path flow funtcional relationship.
Figure 15253DEST_PATH_IMAGE036
(6)
Wherein
Figure 129315DEST_PATH_IMAGE037
: be illustrated in the revised routing probability of h period n bar miss path;
: the priori that is illustrated in h period n bar miss path is selected probability;
Figure 860959DEST_PATH_IMAGE039
: the stochastic error that means the routing probability;
All the other alphabetical meanings are identical with above formula.
What wherein Fig. 4 meaned is hanging down under the coverage rate condition, the relative error of estimation OD and actual OD, and what Fig. 5 meaned is the schematic diagram to the improved effect of priori OD information.Even, when the AVI coverage rate is 60%, when priori OD precision is 40%, its relative error is only also 21.8%.

Claims (2)

1. the Dynamic OD Matrix Estimation method based on the AVI checkout equipment is characterized in that concrete steps are as follows:
(1) extraction of AVI information data and classification
The information that AVI is obtained is usingd the vehicle license data as main index data, and by vehicle license information, identical data are mated combination, and according to the quantity of information of arbitrary same vehicle licence plate, information of vehicles are divided into to following 3 large classes:
1., comprise the vehicle starting point, the data of all path nodes of terminal and process are called path omniscient type vehicle data;
2., the data of two path nodes that at least comprise the process of vehicle are called path and partly know the type vehicle data;
3., the data that only comprise a path node of vehicle process are called path and singly know the type vehicle data;
(2) the expansion sample of prior imformation
The vehicle number that the AVI of take catches is basis, and prior imformation is expanded to the sample processing according to the vehicle number of actual acquisition;
(3) filtration path omniscient type vehicle data
Convert path omniscient type vehicle data to the OD matrix information, and the prior imformation after this OD matrix information and expansion sample is subtracted each other, if OD corresponding in prior imformation is less than 0 to flow, gives a minimum flow 1, and then obtain the priori OD information after filtering;
(4) estimate the path journey time in the period
According to the laying information of different AVI, the journey time of the vehicle of two groups of AVI processes of arbitrary neighborhood is added up and obtained its average travel time; According to the different vehicle journey time, distribute, Based on Probability opinion model, set up AVI journey time probability distribution function, and by the possibility of all the other these two groups of AVI detecting devices of vehicles faileds process of this Functional Analysis; Then pass through highway section-path journey time relation function, the journey time of way to acquire, and this is judged to whether vehicle is estimating the criterion that enters road network in the period as one;
(5) path estimation of type vehicle is partly known in path
The part path of path partly being known to type vehicle disappearance is divided into upstream miss path and downstream miss path by direction of traffic; According to different miss path, with reference to detecting device measurability foundation, by the Bayesian Estimation algorithm, the optional miss path of any vehicle is dwindled rapidly to hunting zone, and on the basis of filtering priori OD information, the selection probability of possible path is revised; Finally by the Monte Carlo random simulation, any vehicle is made to routing in its feasible path, and carry out the path reparation with this; The routing information of type vehicle is partly known in acquired any one group of path, and by the path journey time, the judgement vehicle enters the time of road network and leaves the time of road network;
(6) secondary filtration of priori OD information
Path after repairing is partly known to the type data-switching becomes the OD matrix information, and the prior imformation after this OD matrix information and expansion sample is subtracted each other, if OD corresponding in prior imformation is less than 0 to flow, give a minimum flow 1, simultaneously according to remaining, singly know that the type vehicle number carries out secondary and expands sample, finally obtain the OD information after secondary filtration;
(7) path estimation of type vehicle is singly known in path
Path singly knows that method that the type vehicle route estimates and path partly know that the method for path estimation of type vehicle is similar, partly knows type vehicle route method of estimation according to the path of step (5), and obtains singly knowing the OD matrix information of type information of vehicles with this;
(8) initially repairing the OD matrix obtains
The OD matrix information that front 3 classes are obtained adds up to obtain initial reparation OD matrix;
(9) initially repair the OD matrix information based on the flow adjustment
The link flow information of the OD matrix information of initially reparation and AVI detection is dynamically strolled to matching according to highway section-path flow funtcional relationship; When path flow and highway section measured discharge relative error are 5%, stop strolling matching, and as final OD matrix information.
2. the Dynamic OD Matrix Estimation method based on the AVI checkout equipment according to claim 1, is characterized in that path described in step (5) partly knows that path in the path estimation of type vehicle and step (7) singly knows the path estimation of type vehicle, is specially:
Arrive next highway section detecting device when a certain vehicle has grace time, but the routing information of vehicle is not while comprising next highway section detecting device, has any in following situation:
(1), vehicle passes through next highway section, but is not recorded owing to detecting error;
(2), vehicle, through laying the highway section of detecting device, does not directly arrive the highway section that sensorless is laid;
To the possible residual paths of vehicle, set up initial routing probability set according to prior imformation,
Figure 783173DEST_PATH_IMAGE001
, then according to following formula, analyze vehicle and can within this period, pass through next highway section, and dwindle the range of choice to residual paths with this;
Figure 266107DEST_PATH_IMAGE002
(1.1)
Figure 863311DEST_PATH_IMAGE003
(1.2)
Figure 39077DEST_PATH_IMAGE004
(1.3)
Wherein
Figure 26625DEST_PATH_IMAGE005
: mean the journey time between any two AVI;
Figure 313250DEST_PATH_IMAGE006
: mean the time of i car through arbitrary neighborhood AVI detecting device;
: mean the distance between any two adjacent AVI detecting devices;
Figure 846048DEST_PATH_IMAGE008
: mean the distance between two nodes at place between two AVI detecting devices;
Figure 383209DEST_PATH_IMAGE009
: mean the stochastic error in journey time calculating;
Figure 473524DEST_PATH_IMAGE010
: mean that h is in the period, the selection probability in i the journey time interval by two AVI detecting devices;
Figure 723283DEST_PATH_IMAGE011
: mean that h is in the period, the vehicle number in i the journey time interval by two AVI detecting devices;
Figure 975273DEST_PATH_IMAGE012
: mean that h is in the period, by the vehicle fleet of two AVI detecting devices;
Figure 937413DEST_PATH_IMAGE013
: mean that h is in the period, be numbered the vehicle of j by the journey time of two AVI detecting devices;
Figure 565840DEST_PATH_IMAGE014
: mean that h, in the period, detects the vehicle average travel time according to two AVI detecting devices;
S: mean random chance numerical value;
N: mean time interval length, i.e. unit interval length;
Figure 929825DEST_PATH_IMAGE015
: mean the stochastic error of calculating;
For the residual paths that can select, according to the selection probability of following formula revised residual paths, obtain new routing probability set;
Figure 352717DEST_PATH_IMAGE016
(2.1)
Figure 802152DEST_PATH_IMAGE017
(2.2)
: the prior imformation that means n group miss path in the travel direction of downstream;
Figure 187183DEST_PATH_IMAGE019
: the prior imformation that means n group miss path in the travel direction of upstream;
Figure 46555DEST_PATH_IMAGE020
: mean to meet the total j group of the AVI detecting device that can survey condition in the travel direction of n group miss path downstream;
Figure 780025DEST_PATH_IMAGE021
: mean to meet the total k group of the AVI detecting device that can survey condition in the travel direction of n group miss path upstream;
Figure 15834DEST_PATH_IMAGE022
: be illustrated in the moment that vehicle i is detected for the last time;
Figure 88832DEST_PATH_IMAGE023
: be illustrated in the moment that vehicle i is detected for the first time;
Figure 119105DEST_PATH_IMAGE024
: the zero hour that means to estimate the period;
Figure 277554DEST_PATH_IMAGE025
: the length that means to estimate the period;
Figure 45615DEST_PATH_IMAGE027
: mean that vehicle i organizes the journey time of j the AVI detecting device in miss path downstream from current location to n;
Figure 2011101632062100001DEST_PATH_IMAGE029
: mean the journey time of vehicle i from k of miss path upstream of n group AVI detecting device to the vehicle current location.
CN 201110163206 2011-06-17 2011-06-17 Dynamic OD (Origin Destination) matrix estimating method based on AVI (Automatic Vehicle Identification) device Active CN102289932B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201110163206 CN102289932B (en) 2011-06-17 2011-06-17 Dynamic OD (Origin Destination) matrix estimating method based on AVI (Automatic Vehicle Identification) device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201110163206 CN102289932B (en) 2011-06-17 2011-06-17 Dynamic OD (Origin Destination) matrix estimating method based on AVI (Automatic Vehicle Identification) device

Publications (2)

Publication Number Publication Date
CN102289932A CN102289932A (en) 2011-12-21
CN102289932B true CN102289932B (en) 2013-10-23

Family

ID=45336311

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201110163206 Active CN102289932B (en) 2011-06-17 2011-06-17 Dynamic OD (Origin Destination) matrix estimating method based on AVI (Automatic Vehicle Identification) device

Country Status (1)

Country Link
CN (1) CN102289932B (en)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103440764B (en) * 2013-08-19 2015-10-28 同济大学 A kind of city road network vehicle driving path reconstruction method based on automatic vehicle identification data
CN103646187B (en) * 2013-12-27 2017-02-08 中国科学院自动化研究所 Method for obtaining vehicle travel path and OD (Origin-Destination) matrix in statistic period
CN103903437B (en) * 2014-02-27 2016-06-08 中国科学院自动化研究所 The Vehicle emission OD matrix acquisition methods of data is detected based on video traffic
CN105631945B (en) * 2014-10-30 2018-11-23 国际商业机器公司 The method and apparatus of start point/end point pairing for identification
CN104915731B (en) * 2015-06-11 2018-04-27 同济大学 A kind of grand microcosmic integration method of vehicle driving reconstructing path based on automatic vehicle identification data
CN105679037B (en) * 2016-03-29 2019-03-08 招商局重庆交通科研设计院有限公司 A kind of dynamic path planning method based on user's trip habit
CN109035784B (en) * 2018-09-17 2021-01-26 江苏智通交通科技有限公司 Dynamic traffic flow OD estimation method based on multi-source heterogeneous data
CN109461307B (en) * 2018-11-16 2021-09-24 中电科新型智慧城市研究院有限公司 Method for estimating average vehicle flow and OD (origin-destination) requirements of road section
CN110245423B (en) * 2019-06-14 2023-01-31 重庆大学 Method for analyzing flow relation between highway toll stations
CN110491157B (en) * 2019-07-23 2022-01-25 中山大学 Vehicle association method based on parking lot data and checkpoint data
CN110910647A (en) * 2019-12-17 2020-03-24 江苏智城慧宁交通科技有限公司 Electronic police based traffic travel OD matrix estimation system
CN110956808B (en) * 2019-12-23 2020-11-24 北京交通大学 Heavy truck traffic flow prediction method based on non-full-sample positioning data

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0902404A2 (en) * 1997-09-11 1999-03-17 Siemens Aktiengesellschaft Method for transmitting traffic information
CN101739824A (en) * 2009-11-12 2010-06-16 上海第二工业大学 Data fusion technology-based traffic condition estimation method

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0902404A2 (en) * 1997-09-11 1999-03-17 Siemens Aktiengesellschaft Method for transmitting traffic information
CN101739824A (en) * 2009-11-12 2010-06-16 上海第二工业大学 Data fusion technology-based traffic condition estimation method

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
OD矩阵估计中的AVI设施布局优化;孙剑等;《公路交通科技》;20100831;第27卷(第8期);70-75 *
基于视频牌照识别的动态交通OD估计仿真优化;孙剑等;《公路交通科技》;20090831;第26卷(第8期);130-134 *
孙剑等.OD矩阵估计中的AVI设施布局优化.《公路交通科技》.2010,第27卷(第8期),70-75.
孙剑等.基于视频牌照识别的动态交通OD估计仿真优化.《公路交通科技》.2009,第26卷(第8期),130-134.

Also Published As

Publication number Publication date
CN102289932A (en) 2011-12-21

Similar Documents

Publication Publication Date Title
CN102289932B (en) Dynamic OD (Origin Destination) matrix estimating method based on AVI (Automatic Vehicle Identification) device
US9240124B2 (en) Traffic-volume prediction device and method
CN102288193B (en) Motor vehicle travel route determination method based on historical data
CN103646187B (en) Method for obtaining vehicle travel path and OD (Origin-Destination) matrix in statistic period
CN101944288B (en) Method for setting stop stations of urban bus line
CN104751642B (en) A kind of advanced road real-time predictor method of traffic flow operation risk
US20180211350A1 (en) Urban road network asset valuation method, apparatus and system
CN103440764A (en) Urban road network vehicle travel path reconstruction method based on vehicle automatic identification data
CN103927872A (en) Method for predicting multi-period travel time distribution based on floating vehicle data
CN102663890B (en) Method for determining left-turning forbiddance traffic limitation intersection by using floating car data
CN102930718A (en) Intermittent flow path section travel time estimation method based on floating car data and coil flow fusion
CN110019569B (en) Method for acquiring urban rail transit operation state information
CN104021674A (en) Method for rapidly and accurately forecasting travel time of vehicles for passing through road sections
CN103093625B (en) City road traffic condition real-time estimation method based on reliability verification
Zhi-Peng et al. An improved adaptive exponential smoothing model for short-term travel time forecasting of urban arterial street
CN106097717A (en) The signalized intersections average transit time method of estimation merged based on two class floating car datas
CN108765985A (en) The signalized intersections entrance driveway delay estimation method reached based on first car
Marisamynathan et al. A new approach to estimate pedestrian delay at signalized intersections
CN102156909A (en) Method for identifying rail transit trip mode based on mobile phone signal data
CN105427394A (en) Congestion charging optimal toll rate determining method based on trial-and-error method and motor vehicle flow
CN102184636A (en) Floating vehicle technology based road widening identification method
CN112767690B (en) Regional road network heavy truck traffic estimation method based on section survey data
Hui et al. Estimation of time-varying OD demands incorporating FCD and RTMS data
CN112201041A (en) Trunk road path flow estimation method integrating electric alarm data and sampling trajectory data
Gholami et al. Innovative detector layout for automated traffic turning volume counting

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant