CN112201035A - Signal control intersection periodic flow estimation method driven by track data - Google Patents

Signal control intersection periodic flow estimation method driven by track data Download PDF

Info

Publication number
CN112201035A
CN112201035A CN202011020112.5A CN202011020112A CN112201035A CN 112201035 A CN112201035 A CN 112201035A CN 202011020112 A CN202011020112 A CN 202011020112A CN 112201035 A CN112201035 A CN 112201035A
Authority
CN
China
Prior art keywords
flow
period
tensor
unknown
periodic
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.)
Pending
Application number
CN202011020112.5A
Other languages
Chinese (zh)
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 CN202011020112.5A priority Critical patent/CN112201035A/en
Publication of CN112201035A publication Critical patent/CN112201035A/en
Pending legal-status Critical Current

Links

Images

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
    • 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/0129Traffic data processing for creating historical data or processing based on historical data
    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/065Traffic control systems for road vehicles by counting the vehicles in a section of the road or in a parking area, i.e. comparing incoming count with outgoing count
    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/07Controlling traffic signals
    • G08G1/081Plural intersections under common control
    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/07Controlling traffic signals
    • G08G1/085Controlling traffic signals using a free-running cyclic timer

Abstract

The invention relates to a trajectory data-driven signal-controlled intersection periodic flow estimation method, which comprises the following steps: 1) dividing the periodic flow into a known part and an unknown part according to the parking position and the expected arrival time of the last queuing sampling track; 2) extracting parameters related to flow of the known part and the unknown part in each period, constructing a three-dimensional tensor comprising the flow of the unknown part in each period, and converting the periodic flow estimation problem into a data repair problem; 3) and repairing the constructed tensor based on a tensor decomposition method to obtain the flow of each period. Compared with the prior art, the method has the advantages of track data driving, high effectiveness, wide applicability and the like.

Description

Signal control intersection periodic flow estimation method driven by track data
Technical Field
The invention relates to the field of traffic information, in particular to a trajectory data-driven signal control intersection periodic flow estimation method.
Background
Traffic flow is one of important inputs of urban road planning, operation and management, periodic traffic flow estimation is very important for dynamic evaluation and optimization of urban road signal control intersections, existing traffic flow detection modes mainly rely on fixed-point detectors, such as fixed-point coils, microwave radars, electric alarms and the like, however, the coverage range of the fixed-point detectors in reality is extremely limited due to factors such as installation and maintenance costs, and in addition, practical application of the fixed-point detectors is also influenced by faults and data quality.
In recent years, due to the rapid development of technologies such as mobile internet, intelligent internet and the like, the acquisition of massive high-frequency vehicle track data is possible, generally, vehicles which are provided with a GPS, a Beidou and the like and can upload position information in real time are called internet vehicles or sampling vehicles, at present, such sampling vehicles can upload the instantaneous positions and speeds of the vehicles at intervals of 1-5 s, so that abundant traffic flow information is provided, compared with a traditional fixed-point detector, the track of the sampling vehicle can provide an economic and efficient data source for monitoring the traffic running condition of a road network, and meanwhile, the track of the sampling vehicle has great potential in the aspects of dynamic evaluation and optimization of signal control intersection running, and it needs to be explained that the positioning accuracy of the sampling vehicle is still not enough to realize lane-level positioning of the vehicles at present.
In the last few years, much research has focused on the use of sampled vehicle trajectories to evaluate the behavior of signalized intersections, where a lot of research has focused on estimating the length of queues at signalized intersections, using deterministic or stochastic methods based mainly on traffic wave theory (such as Cheng et al, 2012; Ramezani and gerrolimminis, 2015; Li et al, 2017; Yin et al, 2018), which achieves the estimation of the length of queues per cycle by determining the critical motion points of the sampled vehicle trajectories, i.e. joining or leaving queues, and reconstructing the signalized intersections to control the formation and dissipation processes of the queues at intersections, but when reconstructing the traffic, it is generally necessary to assume that the remaining vehicles between two consecutive sampled vehicles arrive at a constant flow rate, which may generate large errors in low permeability conditions (< 10%), and stochastic methods by assuming that the arrival of vehicles obeys a time-varying distribution, therefore, the periodic queuing length is constructed as a random variable based on methods such as probability theory and the like so as to estimate the maximum possible queuing length or the expected value of the queuing length in the period (such as Commert and Cetin, 2009, 2011; Hao et al, 2013; Hao et al, 2014; Tan et al, 2019; Zhang et al, 2019), the method is generally more robust than a deterministic method, but generally needs other inputs such as permeability, signal timing and the like, in addition, the estimation precision of the method under low permeability still cannot meet the practical application, recently, for different traffic states, Zhang et al (2019) assume different arrival distributions, estimate the arrival rate of each period by adopting an expected value maximization algorithm, then estimate the maximum queuing length and the residual queuing length based on the traffic wave theory, the method can still generate a more accurate estimation result even under a sparse environment, and in addition, tan et al (2019) fully utilize historical trajectory data to obtain arrival distribution in a period, then convert a queuing length estimation problem into a parameter estimation problem based on a probability theory, and finally realize the estimation of the maximum queuing length and the red light tail queuing length.
In addition to queue length estimation, some scholars have attempted to estimate traffic flow based on sample trajectories, Zheng and Liu (2017) et al have translated the flow estimation problem into a maximum likelihood estimation problem based on the assumption that vehicle arrivals obey time-varying poisson, and finally solved using an expectation maximization algorithm, which is the first attempt to estimate traffic flow based on sample trajectories, but have limitations in that given that the arrivals of each cycle within each time interval are consistent and do not achieve estimation of cycle-level flow, Wong et al (2019) and Zhao et al (2019) achieve estimation of signal control intersection vehicle permeability by analyzing the parking positions of the sample vehicles queued within the time interval, and then Zhao et al (2019) have proposed various extended methods to further estimate the number of queued vehicles and the time interval flow based on the estimated permeability, since the proposed methods do not require information such as signal timing, the method can be easily expanded to large-scale road network application, and recently Yao et al (2019) propose a probability model and traffic wave theory hybrid estimation method, the method realizes periodic flow estimation for the first time, the method firstly estimates the queuing length, namely parking flow, based on the traffic wave theory, then converts the parking flow estimation into a parameter estimation problem based on the assumption that the vehicle arrives and obeys bounded time-varying Poisson distribution, and finally realizes the periodic queuing length estimation.
In summary, the existing methods for estimating the queuing length and the flow based on the sampling trajectory have obvious advantages and disadvantages, and the main limitation is that the existing methods generally rely on the assumption of the vehicle arrival process, and such assumption usually requires that the arrival of the vehicles with different flow directions on different lanes are consistent and the queuing principle of first-in first-out is followed, which causes that a plurality of lanes exist in each flow direction of a signal control intersection, and when the queuing trajectory of first-in first-out is not suitable, such methods generate large errors, and furthermore, the more important factor is that the arrival distribution of the vehicles is rather uncertain and unknown in practice, the existing methods assume that a specific arrival distribution is difficult to accurately describe the arrival situation of real vehicles, and the existing methods are limited in that they are all traffic flow model-driven methods, which directly causes the problem that their estimation accuracy under low permeability is insufficient, therefore, a trajectory data driven method is needed to realize accurate estimation of the periodic traffic flow of the signalized intersection.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provide a signal control intersection periodic flow estimation method driven by track data.
The purpose of the invention can be realized by the following technical scheme:
a signal control intersection periodic flow estimation method driven by track data comprises the following steps:
1) dividing the periodic flow into a known part and an unknown part according to the parking position and the expected arrival time of the last queuing sampling track;
2) extracting parameters related to flow of the known part and the unknown part in each period, constructing a three-dimensional tensor comprising the flow of the unknown part in each period, and converting the periodic flow estimation problem into a data repair problem;
3) and repairing the constructed tensor based on a tensor decomposition method to obtain the flow of each period.
In step 1), the flow rate in each period is divided into a known part and an unknown part by the last queued sampling trajectory, and then:
Figure BDA0002700335970000031
Figure BDA0002700335970000032
wherein v iskIs the flow rate of the k-th cycle,
Figure BDA0002700335970000033
the number of vehicles arriving before the last queued sample trajectory, i.e. the known flow,
Figure BDA0002700335970000034
the number of vehicles arriving after the last queued sample trajectory, i.e. the unknown flow,
Figure BDA0002700335970000035
parking position of the sampling trajectory for the last queue, d0The average parking headway distance is shown, and omega is the number of lanes in the lane group.
The step 2) specifically comprises the following steps:
21) for the period of observing the queued sampling track, extracting two known and unknown parameters related to the flow in each period, including the number of the observed sampling tracks and the length of an observation period, and constructing a period characteristic matrix;
22) for the period of the sampling track which is not observed in the queue, replacing the periodic feature matrix of the previous observed track;
23) and finally, constructing a three-dimensional tensor by using the feature matrix of each period in the study period.
In the step 21), the periodic feature matrix is specifically:
Figure BDA0002700335970000041
wherein, XkA feature vector for the kth observed queued sampling trajectory period;
Figure BDA0002700335970000042
and
Figure BDA0002700335970000043
corresponding to the number of observed sample traces for the known and unknown portions respectively,
Figure BDA0002700335970000044
and
Figure BDA0002700335970000045
the lengths of the observation periods corresponding to the known and unknown portions, respectively.
In said step 22), for the period in which no queued sampling trace is observed, the period feature matrix of the previous observed trace is replaced, so that:
Xi=Xi-1
wherein, XiIs the eigenvector of the ith cycle, which does not observe the queued sample traces, Xi-1The feature vector for the i-1 th cycle, which has observed the queued sample traces.
In the step 23), the constructed three-dimensional tensor expression is as follows:
Figure BDA0002700335970000046
wherein the content of the first and second substances,
Figure BDA0002700335970000047
for an initial tensor composed of K cycles of eigenvectors over the study period, vbefAnd vaftVector consisting of known and unknown flow for K cycles in the study period, nbefAnd naftRespectively, a vector consisting of the number of sampling trajectories observed for the K periods of known and unknown flow components within the study period, tbefAnd taftRespectively, vectors consisting of the lengths of observation periods corresponding to the K periods of known and unknown flow portions within the study period.
The step 3) specifically comprises the following steps:
31) initial tensor containing unknown partial flow based on Tucker decomposition
Figure BDA0002700335970000048
Decomposing the data into a form of multiplying a kernel tensor by a factor matrix in three dimension directions, and converting a Tucker decomposition process into an optimization problem;
32) solving an optimization problem based on a gradient descent method, and realizing Tucker decomposition of an initial tensor;
33) and restoring the initial tensor based on the kernel tensor and the factor matrix, realizing tensor repair, obtaining the flow of the unknown part of each period, and finally finishing the estimation of the periodic flow of the intersection.
In the step 31), the initial tensor
Figure BDA0002700335970000051
The expression of the Tucker decomposition of (1) is:
Figure BDA0002700335970000052
wherein the content of the first and second substances,
Figure BDA0002700335970000053
is the decomposed nuclear tensor, an
Figure BDA0002700335970000054
Figure BDA0002700335970000055
I.e. a non-negative tensor of size P x Q x R,
Figure BDA0002700335970000056
and
Figure BDA0002700335970000057
respectively factor matrices with corresponding dimensions of 3 XP, 2 XQ and KXR, respectively, a reference bookξAnd xi is 1,2 and 3, which is the direction product calculation and represents the inner product of the tensor and the matrix in the dimension xi.
In the step 32), the Tucker decomposition process is converted into an optimization problem, and the following steps are performed:
Figure BDA0002700335970000058
wherein λ is a regularization strength parameter, | | · |. non-calculation2The norm of L2 is shown,
Figure BDA0002700335970000059
is a regularization term.
In step 33), the initial tensor is restored based on the kernel tensor and the factor matrix, and then:
Figure BDA00027003359700000510
the flow rate estimates for each cycle are as follows:
Figure BDA00027003359700000511
wherein the content of the first and second substances,
Figure BDA00027003359700000512
in order to restore the initial tensor in the original state,
Figure BDA00027003359700000513
to investigate the flow estimate for period k over the period,
Figure BDA00027003359700000514
is an unknown partial flow estimated based on a tensor decomposition method.
Compared with the prior art, the invention has the following advantages:
firstly, driving track data: the method is different from the existing model-driven flow estimation method based on the traffic flow theory, converts the periodic flow estimation problem into the pure track-driven problem of tensor repair, does not need any prior traffic flow hypothesis, and can be used by non-traffic background professionals.
Secondly, the effectiveness is high: the flow estimation method provided by the invention has higher precision than the existing method in the aspects of cycle-level estimation and time-period-level estimation.
Thirdly, the applicability is wide: the sampling track adopted by the method comprises a network car booking track, data of various geomap navigation companies and the like, the coverage of the whole road network can be basically realized at the present stage, and the method has no requirement on the signal control type, so the method can be applied to the urban road network in a large scale.
Drawings
FIG. 1 is a schematic diagram of the periodic flow division in the present invention.
FIG. 2 is a schematic diagram of the periodic feature extraction in the present invention.
Fig. 3 is a schematic diagram of the construction of the three-dimensional tensor according to the present invention.
Fig. 4 is a decomposition diagram of tensor Tucker in the present invention.
FIG. 5 is a diagram illustrating an exemplary verification scenario.
Fig. 6 is a comparison diagram of results of example verification and the conventional method, wherein fig. 6a is a comparison diagram of cycle-level results, fig. 6b is a comparison diagram of 10-minute results, fig. 6c is a comparison diagram of 30-minute results, and fig. 6d is a comparison diagram of 60-minute results.
Detailed Description
The invention is described in detail below with reference to the figures and specific embodiments.
Examples
The invention provides a track data-driven signal-controlled intersection periodic flow estimation method, which adopts a data driving method to estimate the passing flow of a signal-controlled intersection in each signal period based on sampled vehicle track data and comprises the following steps:
1) the periodic flow is divided into two parts, known and unknown, according to the parking position and the expected arrival time of the last queued sample trajectory, as shown in fig. 1:
Figure BDA0002700335970000061
Figure BDA0002700335970000062
wherein v iskA flow rate of period k;
Figure BDA0002700335970000063
the number of vehicles arriving before the last queued sample trajectory is the known flow;
Figure BDA0002700335970000064
the number of vehicles arriving after the last queued sampling trajectory is the unknown flow;
Figure BDA0002700335970000065
parking position of the sampling trajectory queued for the last; d0The average parking vehicle head distance is obtained; ω is the number of lanes in the study lane group.
2) And extracting parameters related to the flow, such as the number of sampling tracks and the time interval length of each part of each period, constructing a three-dimensional tensor comprising the flow of the unknown part of each period, and converting the periodic flow estimation problem into a data repair problem. The method specifically comprises the following steps:
21) for the period of observing the queued sampling track, extracting two parts of parameters related to the flow, namely the number of the observed sampling tracks and the length of an observation period, of each period, and constructing a period feature matrix, as shown in fig. 2:
Figure BDA0002700335970000066
wherein, XkA feature vector for the kth observed queued sampling trajectory period;
Figure BDA0002700335970000071
and
Figure BDA0002700335970000072
respectively corresponding to the observed sampling track numbers of the known flow part and the unknown flow part;
Figure BDA0002700335970000073
and
Figure BDA0002700335970000074
the lengths of the observation periods corresponding to the known and unknown flow portions, respectively.
22) For the period of the sample trace not observed in the queue, the periodic feature matrix of the previous observed trace is replaced by:
Xi=Xi-1
wherein, XiNo queued sampling trajectory is observed for the feature vector of the ith cycle; xi-1The queued sample traces are observed for the i-1 th cycle of the feature vector.
23) Finally, a three-dimensional tensor is constructed by the feature matrix of each period in the study period, as shown in fig. 3:
Figure BDA0002700335970000075
wherein the content of the first and second substances,
Figure BDA0002700335970000076
for K cycles in the study periodAn initial tensor composed of eigenvectors; v. ofbefAnd vaftOne-dimensional vectors consisting of known and unknown flows of K periods in a research period respectively; n isbefAnd naftRespectively forming one-dimensional vectors by the number of sampling tracks correspondingly observed by K periods of known and unknown flow parts in a research period; t is tbefAnd taftAnd the one-dimensional vectors are respectively formed by the lengths of observation periods corresponding to the known and unknown flow parts of K periods in the research period.
3) And (4) repairing the constructed tensor based on a tensor decomposition method, so that the flow of each period can be obtained. The method specifically comprises the following steps:
31) the problem is converted into an optimization problem by decomposing an initial tensor containing unknown flow into a kernel tensor multiplied by a factor matrix of three dimensional directions based on the Tucker decomposition.
As shown in fig. 4, the Tucker decomposition of the initial tensor is as follows:
Figure BDA0002700335970000077
wherein the content of the first and second substances,
Figure BDA0002700335970000078
to be the decomposed nuclear tensor,
Figure BDA0002700335970000079
i.e., a non-negative tensor of size P × Q × R;
Figure BDA00027003359700000710
and
Figure BDA00027003359700000711
the dimensions of the factor matrixes are three, and the sizes of the factor matrixes are 3 multiplied by P, 2 multiplied by Q and K multiplied by R respectively; is prepared fromξ(xi ═ 1,2,3) is the directional product, representing the inner product of the tensor over dimension xi with the matrix.
This Tucker decomposition process can be translated into an optimization problem as follows:
Figure BDA00027003359700000712
wherein λ is a regularization strength parameter; the second term of the objective function is a regularization term.
32) And solving the optimization problem based on a gradient descent method to realize tensor Tucker decomposition.
Inputting: initial tensor
Figure BDA0002700335970000081
An error threshold epsilon;
initializing a kernel tensor with random values
Figure BDA0002700335970000082
And factor matrix A, B, C, initializing learning rate
Figure BDA0002700335970000083
Marking theta as a target function value;
if Δ θ > ε, cycle:
β←β+1
Figure BDA0002700335970000084
for the initial tensor
Figure BDA0002700335970000085
Each of which is not null
Figure BDA0002700335970000086
Figure BDA0002700335970000087
Figure BDA0002700335970000088
Figure BDA0002700335970000089
Figure BDA00027003359700000810
Figure BDA00027003359700000811
Figure BDA00027003359700000812
Figure BDA00027003359700000813
When the cycle is over, return to
Figure BDA00027003359700000814
A、B、C。
Wherein the content of the first and second substances,
Figure BDA00027003359700000815
representing the outer product; Δ θ is the difference between the two iterations; β is a parameter for initializing the learning rate.
33) And restoring the initial tensor based on the kernel tensor and the factor matrix, and realizing tensor repair, namely estimating the known flow of each period. And restoring the initial tensor based on the nuclear tensor and the factor matrix:
Figure BDA00027003359700000816
finally, the flow rate estimates for each cycle are as follows:
Figure BDA00027003359700000817
wherein the content of the first and second substances,
Figure BDA00027003359700000818
is the flow estimation value of period k in the research period;
Figure BDA00027003359700000819
is an unknown partial flow estimated based on a tensor decomposition method.
4) And taking the Shenzhen emperor Lufu middle road intersection as a case to verify the precision of the flow estimation method provided by the invention.
The case verification scene of the invention is 4 straight lanes of the north approach at the Fuzhonglu-Huang post road intersection in Shenzhen city, as shown in FIG. 5. The sampling vehicle is a drip operating vehicle, and the track uploading frequency is 3 seconds. The data acquisition period is 2017, 4 and 13, and 10: 30 to 14: 30. the real traffic flow is recorded by a high-definition camera erected on the spot. The average permeability is 8.6%, and the cart ratio is 4.2%. The intersection is controlled by a SMOOTH self-adaptive control system widely applied in Shenzhen city, and the time length of each period and the time length of green light fluctuate.
The Mean Absolute Error (MAE) and Mean relative Error (MAPE) were used to evaluate the estimation accuracy, calculated as follows:
Figure BDA0002700335970000091
Figure BDA0002700335970000092
the values of the algorithm-related parameters are shown in the following table:
TABLE 1 Algorithm-related parameters
Parameter(s) Value taking
P 3
Q 2
R 5
λ 0.01
ε 0.001
β 108
Two existing representative methods are selected for comparison, namely the periodic-level flow estimation method of Yao et al (2019) and the time-interval-level flow estimation method of Zheng et al (2017), and the results are shown in fig. 6 and table 2. Therefore, the method provided by the invention is obviously higher than the existing method in the flow estimation precision of the periodic level and the time interval level. The periodic error is only 15.3%, and the hourly estimation error is only 5.1%.
TABLE 2 comparison of the estimation results of the present invention with the prior art methods
Figure BDA0002700335970000093

Claims (10)

1. A signal control intersection periodic flow estimation method driven by track data is characterized by comprising the following steps:
1) dividing the periodic flow into a known part and an unknown part according to the parking position and the expected arrival time of the last queuing sampling track;
2) extracting parameters related to flow of the known part and the unknown part in each period, constructing a three-dimensional tensor comprising the flow of the unknown part in each period, and converting the periodic flow estimation problem into a data repair problem;
3) and repairing the constructed tensor based on a tensor decomposition method to obtain the flow of each period.
2. The method for estimating the periodic flow of the signal-controlled intersection driven by the track data according to claim 1, wherein in the step 1), the flow in each period is divided into a known part and an unknown part by a last queued sampling track, and the method comprises the following steps:
Figure FDA0002700335960000011
Figure FDA0002700335960000012
wherein v iskIs the flow rate of the k-th cycle,
Figure FDA0002700335960000013
the number of vehicles arriving before the last queued sample trajectory, i.e. the known flow,
Figure FDA0002700335960000014
the number of vehicles arriving after the last queued sample trajectory, i.e. the unknown flow,
Figure FDA0002700335960000015
parking position of the sampling trajectory for the last queue, d0The average parking headway distance is shown, and omega is the number of lanes in the lane group.
3. The trajectory data-driven signal-controlled intersection periodic flow estimation method according to claim 2, wherein the step 2) specifically comprises the following steps:
21) for the period of observing the queued sampling track, extracting two known and unknown parameters related to the flow in each period, including the number of the observed sampling tracks and the length of an observation period, and constructing a period characteristic matrix;
22) for the period of the sampling track which is not observed in the queue, replacing the periodic feature matrix of the previous observed track;
23) and finally, constructing a three-dimensional tensor by using the feature matrix of each period in the study period.
4. The method for estimating the periodic traffic of the signal-controlled intersection driven by the trajectory data according to claim 3, wherein in the step 21), the periodic feature matrix is specifically:
Figure FDA0002700335960000021
wherein, XkA feature vector for the kth observed queued sampling trajectory period;
Figure FDA0002700335960000022
and
Figure FDA0002700335960000023
corresponding to the number of observed sample traces for the known and unknown portions respectively,
Figure FDA0002700335960000024
and
Figure FDA0002700335960000025
the lengths of the observation periods corresponding to the known and unknown portions, respectively.
5. The method for estimating the periodic traffic of the signalized intersection driven by the track data as claimed in claim 4, wherein in the step 22), for the period in which the queued sampling track is not observed, the period feature matrix of the previous observed track is replaced by the period feature matrix of the previous observed track, and the method comprises the following steps:
Xi=Xi-1
wherein, XiIs the eigenvector of the ith cycle, which does not observe the queued sample traces, Xi-1The feature vector for the i-1 th cycle, which has observed the queued sample traces.
6. The method for estimating the periodic traffic of the signal-controlled intersection driven by the track data according to claim 4, wherein in the step 23), the constructed three-dimensional tensor expression is as follows:
Figure FDA0002700335960000026
wherein the content of the first and second substances,
Figure FDA0002700335960000027
for an initial tensor composed of K cycles of eigenvectors over the study period, vbefAnd vaftVector consisting of known and unknown flow for K cycles in the study period, nbefAnd naftRespectively, a vector consisting of the number of sampling trajectories observed for the K periods of known and unknown flow components within the study period, tbefAnd taftRespectively, vectors consisting of the lengths of observation periods corresponding to the K periods of known and unknown flow portions within the study period.
7. The trajectory data-driven signal-controlled intersection periodic flow estimation method according to claim 6, wherein the step 3) specifically comprises the following steps:
31) initial tensor containing unknown partial flow based on Tucker decomposition
Figure FDA0002700335960000028
Decomposing the data into a form of multiplying a kernel tensor by a factor matrix in three dimension directions, and converting a Tucker decomposition process into an optimization problem;
32) solving an optimization problem based on a gradient descent method, and realizing Tucker decomposition of an initial tensor;
33) and restoring the initial tensor based on the kernel tensor and the factor matrix, realizing tensor repair, obtaining the flow of the unknown part of each period, and finally finishing the estimation of the periodic flow of the intersection.
8. The method for estimating the periodic traffic of the signal-controlled intersection driven by the track data as claimed in claim 7, wherein in the step 31), an initial tensor is used
Figure FDA0002700335960000031
The expression of the Tucker decomposition of (1) is:
Figure FDA0002700335960000032
wherein the content of the first and second substances,
Figure FDA0002700335960000033
is the decomposed nuclear tensor, an
Figure FDA0002700335960000034
Figure FDA0002700335960000035
I.e. a non-negative tensor of size P x Q x R,
Figure FDA0002700335960000036
and
Figure FDA0002700335960000037
respectively factor matrices with corresponding dimensions of 3 XP, 2 XQ and KXR, respectively, a reference bookξXi is 1,2,3 is directionThe product computation represents the inner product of the tensor over dimension ξ with the matrix.
9. The trajectory data-driven signal-controlled intersection periodic flow estimation method according to claim 8, wherein in the step 32), if a Tucker decomposition process is converted into an optimization problem, the method comprises the following steps:
Figure FDA0002700335960000038
wherein λ is a regularization strength parameter, | | · |. non-calculation2The norm of L2 is shown,
Figure FDA0002700335960000039
is a regularization term.
10. The method for estimating the periodic traffic of the signal-controlled intersection driven by the track data according to claim 9, wherein in the step 33), the initial tensor is restored based on the kernel tensor and the factor matrix, and the method includes:
Figure FDA00027003359600000310
the flow rate estimates for each cycle are as follows:
Figure FDA00027003359600000311
wherein the content of the first and second substances,
Figure FDA00027003359600000312
in order to restore the initial tensor in the original state,
Figure FDA00027003359600000313
to investigate the flow estimate for period k over the period,
Figure FDA00027003359600000314
is an unknown partial flow estimated based on a tensor decomposition method.
CN202011020112.5A 2020-09-24 2020-09-24 Signal control intersection periodic flow estimation method driven by track data Pending CN112201035A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011020112.5A CN112201035A (en) 2020-09-24 2020-09-24 Signal control intersection periodic flow estimation method driven by track data

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011020112.5A CN112201035A (en) 2020-09-24 2020-09-24 Signal control intersection periodic flow estimation method driven by track data

Publications (1)

Publication Number Publication Date
CN112201035A true CN112201035A (en) 2021-01-08

Family

ID=74007226

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011020112.5A Pending CN112201035A (en) 2020-09-24 2020-09-24 Signal control intersection periodic flow estimation method driven by track data

Country Status (1)

Country Link
CN (1) CN112201035A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113240902A (en) * 2021-03-25 2021-08-10 同济大学 Signal control road network path flow estimation method based on sampled vehicle trajectory data
CN114141013A (en) * 2021-11-25 2022-03-04 东南大学 Intersection signal timing method based on low-permeability network connection track data

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107123276A (en) * 2016-08-25 2017-09-01 苏州华川交通科技有限公司 Utilize the intersection vehicles queue length evaluation method of low sampling rate gps data
CN108805348A (en) * 2018-06-05 2018-11-13 北京京东金融科技控股有限公司 A kind of method and apparatus of intersection signal timing control optimization
CN108898829A (en) * 2018-06-07 2018-11-27 重庆邮电大学 The dynamic short-time traffic flow forecast system with Sparse is divided for the indifference opposite sex
CN109377760A (en) * 2018-11-29 2019-02-22 北京航空航天大学 The detection of loss traffic data and restorative procedure based on iteration tensor algorithm
CN109544915A (en) * 2018-11-09 2019-03-29 同济大学 A kind of queue length distribution estimation method based on sample path data
CN109830102A (en) * 2019-02-14 2019-05-31 重庆邮电大学 A kind of short-term traffic flow forecast method towards complicated urban traffic network
CN110149228A (en) * 2019-05-20 2019-08-20 湖南友道信息技术有限公司 A kind of top-k elephant stream prediction technique and system based on the filling of discretization tensor

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107123276A (en) * 2016-08-25 2017-09-01 苏州华川交通科技有限公司 Utilize the intersection vehicles queue length evaluation method of low sampling rate gps data
CN108805348A (en) * 2018-06-05 2018-11-13 北京京东金融科技控股有限公司 A kind of method and apparatus of intersection signal timing control optimization
CN108898829A (en) * 2018-06-07 2018-11-27 重庆邮电大学 The dynamic short-time traffic flow forecast system with Sparse is divided for the indifference opposite sex
CN109544915A (en) * 2018-11-09 2019-03-29 同济大学 A kind of queue length distribution estimation method based on sample path data
CN109377760A (en) * 2018-11-29 2019-02-22 北京航空航天大学 The detection of loss traffic data and restorative procedure based on iteration tensor algorithm
CN109830102A (en) * 2019-02-14 2019-05-31 重庆邮电大学 A kind of short-term traffic flow forecast method towards complicated urban traffic network
CN110149228A (en) * 2019-05-20 2019-08-20 湖南友道信息技术有限公司 A kind of top-k elephant stream prediction technique and system based on the filling of discretization tensor

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
KESHUANG TANG 等: "A tensor decomposition method for cycle-based traffic volume estimation using sampled", 《TRANSPORTATION RESEARCH PART C: EMERGING TECHNOLOGIES》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113240902A (en) * 2021-03-25 2021-08-10 同济大学 Signal control road network path flow estimation method based on sampled vehicle trajectory data
CN114141013A (en) * 2021-11-25 2022-03-04 东南大学 Intersection signal timing method based on low-permeability network connection track data
CN114141013B (en) * 2021-11-25 2022-12-13 东南大学 Intersection signal timing method based on low-permeability network connection track data

Similar Documents

Publication Publication Date Title
CN107180530B (en) A kind of road network trend prediction method based on depth space-time convolution loop network
CN106981198B (en) Deep learning network model for travel time prediction and establishing method thereof
CN108345666B (en) Vehicle abnormal track detection method based on time-space isolated points
CN109272157A (en) A kind of freeway traffic flow parameter prediction method and system based on gate neural network
CN109872533B (en) Lane-level real-time traffic information processing method based on spatial data
CN113591380B (en) Traffic flow prediction method, medium and equipment based on graph Gaussian process
CN111709549A (en) Navigation reminding method for short-time traffic flow prediction based on SVD-PSO-LSTM
CN109118787B (en) Vehicle speed prediction method based on deep neural network
CN112201035A (en) Signal control intersection periodic flow estimation method driven by track data
AU2020316538B2 (en) Meteorological parameter-based high-speed train positioning method and system in navigation blind zone
CN112862171B (en) Flight arrival time prediction method based on space-time neural network
CN107730889B (en) Target vehicle retrieval method based on traffic video
CN113129605B (en) Electronic police data-based intersection lane queuing length estimation method
CN112633602B (en) Traffic congestion index prediction method and device based on GIS map information
CN112201037B (en) Intersection arrival rate estimation method based on sampling trajectory data
CN114202120A (en) Urban traffic travel time prediction method aiming at multi-source heterogeneous data
CN111383273B (en) High-speed rail contact net part positioning method based on improved structure reasoning network
CN116307152A (en) Traffic prediction method for space-time interactive dynamic graph attention network
CN115985104A (en) Traffic flow prediction device, prediction method and prediction model construction method
CN115099328A (en) Traffic flow prediction method, system, device and storage medium based on countermeasure network
CN114860700A (en) Missing track filling method and system
CN110458370A (en) A kind of urban roads safety prediction technique and system based on open data
CN116543603A (en) Flight path completion prediction method and device considering airspace situation and local optimization
CN114996659B (en) Vehicle track prediction method and system with continuous learning capability
Yang et al. Arterial roadway travel time distribution estimation and vehicle movement classification using a modified Gaussian Mixture Model

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20210108

RJ01 Rejection of invention patent application after publication