CN107992536B - Urban traffic missing data filling method based on tensor decomposition - Google Patents
Urban traffic missing data filling method based on tensor decomposition Download PDFInfo
- Publication number
- CN107992536B CN107992536B CN201711184489.2A CN201711184489A CN107992536B CN 107992536 B CN107992536 B CN 107992536B CN 201711184489 A CN201711184489 A CN 201711184489A CN 107992536 B CN107992536 B CN 107992536B
- Authority
- CN
- China
- Prior art keywords
- missing data
- svd
- tensor
- data
- filling
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 46
- 238000000354 decomposition reaction Methods 0.000 title claims abstract description 35
- 239000013598 vector Substances 0.000 claims abstract description 23
- 238000005457 optimization Methods 0.000 claims abstract description 15
- 238000005065 mining Methods 0.000 claims abstract description 5
- 238000012549 training Methods 0.000 claims abstract description 3
- 239000011159 matrix material Substances 0.000 claims description 16
- 238000004364 calculation method Methods 0.000 claims description 11
- 230000008859 change Effects 0.000 claims description 10
- 230000008569 process Effects 0.000 claims description 9
- 230000009467 reduction Effects 0.000 claims description 4
- 238000011161 development Methods 0.000 claims description 2
- 238000002922 simulated annealing Methods 0.000 claims description 2
- 230000000694 effects Effects 0.000 description 12
- 230000006870 function Effects 0.000 description 10
- 238000013528 artificial neural network Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 3
- 238000012217 deletion Methods 0.000 description 3
- 230000037430 deletion Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000003086 colorant Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 239000013256 coordination polymer Substances 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000005429 filling process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/067—Enterprise or organisation modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/21—Design, administration or maintenance of databases
- G06F16/211—Schema design and management
-
- G06Q50/40—
Abstract
The invention relates to an urban traffic missing data filling method based on tensor decomposition, which comprises the following steps of: s1, constructing tensors of urban traffic data based on dimensions of road sections, dates and time periods; s2, pre-filling missing data to finish initialization of the missing data; s3, performing truncation singular value decomposition on the missing data obtained by pre-filling, and mining to obtain left singular vectors of the missing data in dimensions of road sections, dates and time periods; s4, calculating to obtain a core tensor by using the left singular vectors of the road section, the date and the time interval dimension; and S5, constructing a missing data filling model, inputting left singular vectors of road section, date and time interval dimensions and a core tensor training missing data filling model, continuously optimizing the missing data filling model by combining an optimization algorithm, and filling missing data through the missing data filling model after optimization is finished.
Description
Technical Field
The invention relates to the technical field of intelligent traffic, in particular to a tensor decomposition-based urban traffic missing data filling method.
Background
In the big data era, data-driven analytical methods have brought countless dividends to people. Taking the traffic industry as an example, by collecting and analyzing vehicle speed data, a city manager can master the running condition of a road section. However, in the data acquisition process, the inevitable problem is missing data, which directly hinders the urban managers from mastering the road section operation conditions, and thus effective management decisions cannot be made. Fortunately, urban traffic data often has strong space-time laws, such as morning and evening peaks and the like which often appear on workdays. Therefore, rules can be mined from historical data through a certain method, and missing data can be filled.
How to fully mine the existing data rules through a scientific method and realize the filling of large-scale traffic missing data so as to help city managers master traffic operation conditions and make effective management decisions, which is the key point and the purpose of the research of the invention.
The current domestic and foreign research and technical current situations are mainly as follows:
first, filling scattered missing data of urban traffic. Scholars at home and abroad make a great deal of research on filling methods for scattered urban traffic data loss, and classically, data points closely related to the scattered urban data loss are found out by utilizing an autocorrelation analysis theory, and then the scattered urban data loss is filled by adopting a weighted average method. The algorithm model mainly considers the characteristics among data, the characteristics of the data are less mined, and the algorithm model has the defects of complex calculation and high time consumption and is difficult to popularize on a large scale.
Secondly, filling up urban traffic continuous missing data. The existing continuous missing data filling method mainly takes machine learning related algorithms such as neural networks, support vector regression and the like as main steps, utilizes the existing historical data to train and obtain a regression model, and fills the continuous missing data. Although the filling effect of the algorithm is good, the algorithm also has the defects of complex calculation, more filling model parameters and poor interpretability, and is not beneficial to other users to understand the model and reproduce the model, so that the popularization of related application is limited to a certain extent.
Disclosure of Invention
The method aims to organize data from the time-space characteristics of urban traffic, excavate the internal rules of the data by a tensor decomposition method, establish a missing data filling model with strong interpretability by using the excavated rules, and finally realize effective filling of scattered or continuous missing data.
In order to realize the purpose, the technical scheme is as follows:
a tensor decomposition-based urban traffic missing data filling method comprises the following steps:
s1, constructing tensors of urban traffic data based on dimensions of road sections, dates and time periods;
s2, pre-filling missing data to finish initialization of the missing data;
s3, performing truncation singular value decomposition on the tensor obtained by pre-filling, and mining to obtain a left singular vector of the tensor in dimensions of a road section, a date and a time period;
s4, calculating to obtain a core tensor by using the left singular vectors of the road section, the date and the time interval dimension;
and S5, constructing a missing data filling model, inputting left singular vectors of road section, date and time interval dimensions and a core tensor training missing data filling model, continuously optimizing the missing data filling model by combining an optimization algorithm, and filling missing data through the missing data filling model after optimization is finished.
Preferably, the urban traffic data is vehicle speed data.
Preferably, the step S2 is to perform a pre-filling process on the missing data as follows:
s11, calculating an average value of observed vehicle speed data;
s12, initializing offset vectors in dimensions of road sections, dates and time periods;
s13, calculating a data estimation value by combining the offset vectors on the road section, date and time interval dimensions and the average speed:
whereinMu is the average of the observed vehicle speed data as an estimated value,for the biasing of the vehicle speed data in the road segment dimension,for the bias of the vehicle speed data in the date dimension,a bias in a time period dimension for vehicle speed data;
s14, calculating an objective function J:
wherein xijkIs an accurate value; a is an optimization function regular term parameter used for controlling the complexity of the model;
s15, judging whether the numerical value of the target function is minimum, if so, outputting a pre-filling value of the data to be filledIf so, outputting, finishing the initialization of missing data, otherwise, adjusting to a more optimal directionAnd then step S13 is performed.
Preferably, in step S3, the specific process of performing truncated singular value decomposition on the missing data obtained by pre-padding is as follows:
s21, performing modal expansion on missing data obtained by pre-filling to obtain 3 matrixes, and recording the matrixes as matrixes A, B and C;
s22, setting a singular value proportion threshold value p;
s23, singular value decomposition is carried out on the matrixes A, B and C, and left singular vectors U are obtained through calculation respectively(SVD)、V(SVD)And W(SVD)Wherein U is(SVD)、V(SVD)And W(SVD)Left singular vectors representing missing data in road segment, date and time segment dimensions, respectively.
Preferably, the modal 1 development of the core tensor is as follows:
Preferably, the specific implementation procedure of step S5 is as follows:
s31, constructing a missing data filling model;
s32, loadingU(SVD)、V(SVD)And W(SVD)As an initial value for the iterative update; inputting the data into a constructed missing data filling model to train the missing data filling model;
the data loss in the tensor is recorded for the binary tensor, and if some data in the tensor is lost, the tensor is formedThe record of the middle corresponding position is 0; if some data in the tensor is not missing, the tensorThe record of the middle corresponding position is 1;
s35, judging whether the error function is converged, if so, outputting a missing data filling model and a complete tensor x, and otherwise, performing optimization algorithm pairU(SVD)、V(SVD)And W(SVD)The optimization is performed and then step S32 is performed.
Preferably, the specific process of obtaining the tensor x by restoring the missing data filling model is as follows:
. Preferably, the step S35 outputs the factor matrix U when the error function converges(SVD)、V(SVD)And W(SVD)To U, to U(SVD)、V(SVD)And W(SVD)And analyzing to obtain the change modes of the vehicle speed data in the road section dimension, the date dimension and the time interval dimension, and improving the explanatory power of the missing data filling model.
Preferably, the step S35 adopts a pair of gradient optimization algorithm, newton method or simulated annealing algorithmU(SVD)、V(SVD)And W(SVD)And optimizing.
Compared with the prior art, the invention has the beneficial effects that:
(1) the method is different from the traditional data filling algorithm, can effectively fill scattered and continuous missing data, and has wide applicability;
(2) a filling model is established based on the data internal rule, the filling model has good interpretability and can be reproduced and applied in a large scale;
(3) the calculation method is simple, and the missing data filling efficiency is high.
Drawings
Fig. 1 is a frame diagram of a large-scale traffic missing data filling method based on tensor decomposition.
Fig. 2 is a graph showing the average vehicle speed variation in different road sections.
Fig. 3 is a graph showing the change in average vehicle speed on different days.
Fig. 4 is a graph showing the change in average vehicle speed at different times.
FIG. 5 is a flow chart of missing data initialization.
FIG. 6 is a flow diagram of truncated singular value decomposition.
FIG. 7 is a schematic diagram of a Tucker decomposition.
Fig. 8 is a flow chart of constructing a missing data padding model based on tensor decomposition.
FIG. 9 is a comparison graph of the effect of the neural network-based method of the present invention on filling up scattered data under different miss rates.
FIG. 10 is a comparison graph of the effect of the neural network-based method of the present invention on filling up continuous data loss at different loss rates.
FIG. 11 is a graph of a factor matrix filling the date dimension in the model.
Detailed Description
The drawings are for illustrative purposes only and are not to be construed as limiting the patent;
the invention is further illustrated below with reference to the figures and examples.
Example 1
As shown in fig. 1, the method provided by the present invention comprises the following steps:
one, constructing tensors
As shown in table 1 (where the time interval refers to the division of the time of day at intervals of 10 minutes), the vehicle speed data collected by the sensor is visualized, the time-space characteristic change rule is preliminarily explored, and as seen from the observation results in fig. 2 to fig. 4, the average vehicle speed has more obvious changes in different road sections, different dates and different time intervals: the average vehicle speed difference of the congested easy-to-send road section and the congested difficult-to-send road section is large; the average speed difference of working days, weekends and holidays is large; the average speed difference of the early peak, the late peak and the flat peak is large; fourthly, the change of the average speed has no obvious periodicity. Therefore, according to the result of the visual preliminary exploration, the tensor of the vehicle speed data based on the dimensions of the road section, the date and the time period is constructed.
TABLE 1 map matching results data sheet
Second, missing data initialization
From the analysis results of fig. 2-4, the vehicle speed data is affected by different road sections, different dates and different time periods, so the vehicle speed data can be regarded as the average vehicle speed plus the offset of the vehicle speed in the dimension of road section, day and time, as shown in equation (1).
Wherein the content of the first and second substances,mu is the average of the observed vehicle speed data as an estimated value,for the biasing of the vehicle speed data in the road segment dimension,for the bias of the vehicle speed data in the date dimension,is the offset of the vehicle speed data in the time period dimension.
Then, solving the magnitude of each offset value, and setting an objective function J, as shown in formula (2), which mainly comprises two parts, wherein one part is a measurement estimation valueAnd the exact value xijkAn error term for the error; the other part is a regularization term that measures the complexity of the model. By minimizing the objective function J, the optimum size of each offset value can be found.
The overall flow chart of this section is shown in fig. 5.
Three, truncated singular value decomposition
After the initialization step of missing data, although a complete tensor is obtained, the data filling effect is limited and the filling method is not strong enough in explanation, so that rules contained in the tensor are further mined. And performing modal expansion on the pre-filled tensor, namely matrixing the high-order tensor, and converting the abstract high-order tensor into a matrix of a two-dimensional space so as to facilitate subsequent matrix operation.
The invention constructs 3-order tensor for vehicle speed data, and 3 matrixes can be obtained after modal expansion and are recorded as matrixes A, B and C. Then, singular value decomposition is carried out on the matrixes A, B and C to obtain a left singular vector U(SVD)、V(SVD)And W(SVD)Wherein U is(SVD)、V(SVD)And W(SVD)And respectively representing left singular vectors of the missing data in the dimensions of the road section, the date and the time interval so as to obtain main characteristics of each matrix and remove data noise in the matrix. The singular values and the main left singular vectors of the matrix are controlled by setting a singular value proportion threshold p. If the singular value proportion threshold value p is set to be too small, the regular mining in the matrix is insufficient; if the setting is too large, noise of data in the matrix cannot be removed, and a method of repeated experiments is needed to determine a proper singular value proportion threshold value p in formal application.
After the step, left singular vectors of all dimensions are obtained. The overall flow chart of this section is shown in fig. 6.
Calculation of the four, core tensor
The concept of the core tensor comes from tensor decomposition, and common tensor decomposition includes CP decomposition, Tucker decomposition, etc., and here, a simple description is given by taking the Tucker decomposition as an example, as shown in fig. 7, and the calculation method is shown in formula (3). The tensor is used as a high-dimensional data structure, a complex data rule is hidden in the tensor, the tensor is decomposed into a plurality of factor matrixes and a core tensor through Tucker decomposition, the factor matrixes can reveal features of all dimensions of the tensor to guarantee interpretability of a final model, and the core tensor describes connection among all dimensions.
For the next step capable of decomposing the tensorThe structure is trained, and the core tensor needs to be obtained by calculationCore tensorThe mode 1 expansion form can be obtained by calculating the left singular vector and the matrix A obtained in the last step, and the calculation method is shown as a formula (4).
Missing data filling model based on tensor decomposition
Obtaining a factor matrix U through the steps of four and five(SVD)、W(SVD)And V(SVD)And core tensorThe tensor x can be obtained by reduction by the formula (3):
therefore, in this step, only the factor matrix U is required(SVD)、W(SVD)And V(SVD)And core tensorThe value of (A) is optimized, and the filled complete tensor x can be obtained. And in the optimizing process, a gradient optimizing algorithm is adopted to find the optimal tensor decomposition structure, so that an optimal missing data filling model is constructed.
EstablishingAn error function as shown in equation (5), whereThe data loss in the tensor is recorded for the binary tensor, and if some data in the tensor is lost, the tensor is formedThe record of the middle corresponding position is 0; if some data in the tensor is not missing, the tensorThe record of the corresponding position in (1). TensorThe existence of the model is beneficial to better judging the real filling effect of the model.
The process of model construction using the gradient optimization algorithm is shown in fig. 8.
By analysing the tensor decomposition structure U under the minimum error function(SVD)、V(SVD)And W(SVD)The change modes of the vehicle speed data in the road section dimension, the date dimension and the time interval dimension can be mined, and the explanatory power of the filling model is improved.
Example 2
In order to test the missing data filling effect of the method provided by the invention, 1855589 pieces of vehicle speed detection data of 214 road sections in Guangzhou city, 8/month 1/2016 to 9/month 30/2016 (61 days in total) are selected for experiment, and the following experiment scenes are set: testing the filling effect of scattered data deletion and continuous data deletion under different data deletion rates, and comparing the effect with the effect of a common neural network filling method; testing the interpretability of the filling model, namely the mining effect of the time-space change mode of the vehicle speed data.
And (3) introducing an average absolute error MAE as a filling effect evaluation index, wherein the calculation method is shown as a formula (6).
Where N is the number of padding data, yiTo fill in value, ziAre true values.
The error calculation results are shown in fig. 9-10, and it can be found that (i) under the condition of low loss rate or high loss rate, the method provided by the invention can better fill scattered data loss and continuous data loss, and the average absolute error of filling is between 2.5km/h and 3 km/h; secondly, the method provided by the invention is more stable. Whether filling scattered data loss or continuous data loss, the filling effect change fluctuation under different loss rates is small, and compared with a data filling method of a neural network, filling is more stable.
Taking the factor matrix of the date dimension as an example, the factor matrix of the date dimension in the final filling model is shown in fig. 11, and different colors represent different values. Observing the second column of the factor matrix, the occurrence of the dark blue stripe is found to have certain regularity, and combining with an actual calendar, several places where the dark blue stripe occurs are found to correspond to weekends (such as 8 months 6 days, 8 months 7 days, 8 months 13 days, 8 months 14 days and the like) or holidays (middle autumn festival from 9 months 15 days to 9 months 17 days), which shows that the vehicle speed change mode under different date conditions is successfully excavated by the method provided by the invention, and the method has strong explanatory property.
It should be understood that the above-described embodiments of the present invention are merely examples for clearly illustrating the present invention, and are not intended to limit the embodiments of the present invention. Other variations and modifications will be apparent to persons skilled in the art in light of the above description. And are neither required nor exhaustive of all embodiments. Any modification, equivalent replacement, and improvement made within the spirit and principle of the present invention should be included in the protection scope of the claims of the present invention.
Claims (7)
1. A tensor decomposition-based urban traffic missing data filling method is characterized by comprising the following steps: the method comprises the following steps:
s1, constructing tensors of urban traffic data based on dimensions of road sections, dates and time periods;
the urban traffic data is vehicle speed data;
the specific process of pre-filling the missing data in step S2 is as follows:
s11, calculating an average value of observed vehicle speed data;
s12, initializing offset vectors in dimensions of road sections, dates and time periods;
s13, calculating a data estimation value by combining the offset vectors on the road section, date and time interval dimensions and the average speed:
whereinMu is the average of the observed vehicle speed data as an estimated value,for the biasing of the vehicle speed data in the road segment dimension,for the bias of the vehicle speed data in the date dimension,a bias in a time period dimension for vehicle speed data;
s14, calculating an objective function J:
wherein xijkIs an accurate value; alpha is an optimization function regular term parameter and is used for controlling the complexity of the model;
s15, judging whether the numerical value of the target function is minimum, if so, outputting a pre-filling value of the data to be filledFinishing the initialization of missing data; otherwise, the direction is adjusted to a more optimal directionThen step S13 is performed;
s2, pre-filling missing data to finish initialization of the missing data;
s3, performing truncation singular value decomposition on the tensor obtained by pre-filling, and mining to obtain a left singular vector of the tensor in dimensions of a road section, a date and a time period;
s4, calculating to obtain a core tensor by using the left singular vectors of the road section, the date and the time interval dimension;
and S5, constructing a missing data filling model, inputting left singular vectors of road section, date and time interval dimensions and a core tensor training missing data filling model, continuously optimizing the missing data filling model by combining an optimization algorithm, and filling missing data through the missing data filling model after optimization is finished.
2. The tensor decomposition-based urban traffic missing data filling method according to claim 1, characterized in that: the specific process of performing truncated singular value decomposition on the missing data obtained by pre-padding in step S3 is as follows:
s21, performing modal expansion on missing data obtained by pre-filling to obtain 3 matrixes, and recording the matrixes as matrixes A, B and C;
s22, setting a singular value proportion threshold value p;
s23, singular value decomposition is carried out on the matrixes A, B and C, and left singular vectors U are obtained through calculation respectively(SVD)、V(SVD)And W(SVD)Wherein U is(SVD)、V(SVD)And W(SVD)Left singular vectors representing missing data in road segment, date and time segment dimensions, respectively.
4. The tensor decomposition-based urban traffic missing data filling method according to claim 3, wherein: the specific implementation process of step S5 is as follows:
s31, constructing a missing data filling model;
s32, loadingU(SVD)、V(SVD)And W(SVD)As an initial value for the iterative update; inputting the data into a constructed missing data filling model to train the missing data filling model;
the data loss in the tensor is recorded for the binary tensor, and if some data in the tensor is lost, the tensor is formedThe record of the middle corresponding position is 0; if some data in the tensor is not missing, the tensorThe record of the middle corresponding position is 1;
6. the tensor decomposition-based urban traffic missing data filling method according to claim 4, wherein: when the error function converges, the step S35 outputs the factor matrix U(SVD)、V(SVD)And W(SVD)To U, to U(SVD)、V(SVD)And W(SVD)And analyzing to obtain the change modes of the vehicle speed data in the road section dimension, the date dimension and the time interval dimension, and improving the explanatory power of the missing data filling model.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711184489.2A CN107992536B (en) | 2017-11-23 | 2017-11-23 | Urban traffic missing data filling method based on tensor decomposition |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711184489.2A CN107992536B (en) | 2017-11-23 | 2017-11-23 | Urban traffic missing data filling method based on tensor decomposition |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107992536A CN107992536A (en) | 2018-05-04 |
CN107992536B true CN107992536B (en) | 2020-10-30 |
Family
ID=62032593
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711184489.2A Active CN107992536B (en) | 2017-11-23 | 2017-11-23 | Urban traffic missing data filling method based on tensor decomposition |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107992536B (en) |
Families Citing this family (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108898829B (en) * | 2018-06-07 | 2021-02-09 | 重庆邮电大学 | Dynamic short-time traffic flow prediction system aiming at non-difference division and data sparseness |
CN109101638B (en) * | 2018-08-20 | 2019-05-28 | 河海大学 | A kind of Dam Deformation Monitoring continuity missing data complementing method |
CN109377760A (en) * | 2018-11-29 | 2019-02-22 | 北京航空航天大学 | The detection of loss traffic data and restorative procedure based on iteration tensor algorithm |
CN109684604B (en) * | 2018-12-06 | 2020-06-30 | 北京航空航天大学 | City dynamic analysis method based on context-aware nonnegative tensor decomposition |
CN109584557B (en) * | 2018-12-14 | 2021-02-26 | 北京工业大学 | Traffic flow prediction method based on dynamic decomposition mode and matrix filling |
CN109683029A (en) * | 2018-12-29 | 2019-04-26 | 云南电网有限责任公司电力科学研究院 | A kind of outlier identification method of metering automation system |
JP2020154828A (en) * | 2019-03-20 | 2020-09-24 | 富士通株式会社 | Data complement program, data complement method, and data complement device |
CN110162744B (en) * | 2019-05-21 | 2023-01-17 | 天津理工大学 | Tensor-based internet of vehicles data loss multiple estimation method |
CN110322693B (en) * | 2019-07-10 | 2020-09-29 | 山东大学 | Traffic data completion method, system, device and medium |
CN110428614B (en) * | 2019-07-11 | 2021-02-05 | 银江股份有限公司 | Traffic jam heat degree space-time prediction method based on non-negative tensor decomposition |
CN110689181A (en) * | 2019-09-18 | 2020-01-14 | 北京航空航天大学 | Travel time judgment method based on collaborative tensor decomposition |
CN110837888A (en) * | 2019-11-13 | 2020-02-25 | 大连理工大学 | Traffic missing data completion method based on bidirectional cyclic neural network |
CN111640298A (en) * | 2020-05-11 | 2020-09-08 | 同济大学 | Traffic data filling method, system, storage medium and terminal |
CN112382090A (en) * | 2020-11-11 | 2021-02-19 | 北京百度网讯科技有限公司 | Method, apparatus, device and storage medium for outputting information |
CN112488312B (en) * | 2020-12-07 | 2022-02-18 | 江苏自动化研究所 | Construction method of tensor-based automatic coding machine for network exchange data |
CN112836000A (en) * | 2021-01-28 | 2021-05-25 | 北京工业大学 | Abnormal resident travel mode mining method based on taxi OD data |
CN114120635B (en) * | 2021-11-05 | 2023-03-03 | 兴民智通(武汉)汽车技术有限公司 | Tensor decomposition-based urban road network linear missing flow estimation method and system |
CN116013087B (en) * | 2023-03-27 | 2023-05-30 | 东莞市城建规划设计院 | Traffic flow statistical method based on urban moving vehicle detection |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101794508A (en) * | 2009-12-30 | 2010-08-04 | 北京世纪高通科技有限公司 | Traffic information filling method, device and system |
US8818919B2 (en) * | 2011-08-05 | 2014-08-26 | International Business Machines Corporation | Multiple imputation of missing data in multi-dimensional retail sales data sets via tensor factorization |
CN105224507A (en) * | 2015-09-29 | 2016-01-06 | 杭州天宽科技有限公司 | A kind of disappearance association rule mining method based on tensor resolution |
CN106485212A (en) * | 2016-09-26 | 2017-03-08 | 天津大学 | A kind of resolution of tensor method for non-isometric video gesture identification |
CN107016649A (en) * | 2017-02-24 | 2017-08-04 | 同济大学 | A kind of vision data complementing method estimated based on local low-rank tensor |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103136239B (en) * | 2011-11-29 | 2015-03-25 | 北京理工大学 | Transportation data loss recovery method based on tensor reconstruction |
CN104217002B (en) * | 2014-09-14 | 2017-08-25 | 北京航空航天大学 | A kind of road condition information filling method obtained based on quality data |
CN104751169B (en) * | 2015-01-10 | 2018-02-02 | 哈尔滨工业大学(威海) | High ferro rail defects and failures sorting technique |
CN106646595B (en) * | 2016-10-09 | 2018-05-29 | 电子科技大学 | A kind of seismic data compression method that adaptive order based on tensor is blocked |
-
2017
- 2017-11-23 CN CN201711184489.2A patent/CN107992536B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101794508A (en) * | 2009-12-30 | 2010-08-04 | 北京世纪高通科技有限公司 | Traffic information filling method, device and system |
US8818919B2 (en) * | 2011-08-05 | 2014-08-26 | International Business Machines Corporation | Multiple imputation of missing data in multi-dimensional retail sales data sets via tensor factorization |
CN105224507A (en) * | 2015-09-29 | 2016-01-06 | 杭州天宽科技有限公司 | A kind of disappearance association rule mining method based on tensor resolution |
CN106485212A (en) * | 2016-09-26 | 2017-03-08 | 天津大学 | A kind of resolution of tensor method for non-isometric video gesture identification |
CN107016649A (en) * | 2017-02-24 | 2017-08-04 | 同济大学 | A kind of vision data complementing method estimated based on local low-rank tensor |
Non-Patent Citations (2)
Title |
---|
Novel methods for multilinear data completion and de-noising based on Tensor-SVD;Zemin Zhang et al;《2014 IEEE Conference on Computer Vision and Pattern Recognition》;20140628;第3842-3849页 * |
交通流量数据缺失值的插补方法;韩卫国等;《交通与计算机》;20050228;第23卷(第1期);第39-42页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107992536A (en) | 2018-05-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107992536B (en) | Urban traffic missing data filling method based on tensor decomposition | |
Zeng et al. | Application of a new grey prediction model and grey average weakening buffer operator to forecast China’s shale gas output | |
CA3014415C (en) | Reliability evaluation method for hoist main shaft of kilometer-deep mine considering multiple failure modes | |
CN112906982A (en) | GNN-LSTM combination-based network flow prediction method | |
US20190230913A1 (en) | Method of Predicting Central Fishing Ground of Flying Squid Family Ommastrephidae | |
CN107423814A (en) | A kind of method that dynamic network model is established using depth convolutional neural networks | |
Zou et al. | Use of skew-normal and skew-t distributions for mixture modeling of freeway speed data | |
CN110428175A (en) | A kind of Hot Strip Crown Prediction of Media method based on deep learning | |
CN111182564A (en) | Wireless link quality prediction method based on LSTM neural network | |
CN111898318A (en) | Industrial building concrete carbonization depth prediction method | |
CN114693064B (en) | Building group scheme generation performance evaluation method | |
CN108595803A (en) | Shale gas well liquid loading pressure prediction method based on recurrent neural network | |
CN112966871A (en) | Traffic jam prediction method and system based on convolution long-short term memory neural network | |
CN115730684A (en) | Air quality detection system based on LSTM-CNN model | |
CN113393488A (en) | Behavior track sequence multi-feature simulation method based on quantum migration | |
CN103049669B (en) | Wawter bloom based on the analysis of Multivariate Periodic Stationary Time Series and gray theory is predicted and factor approach | |
Gao et al. | An agent-based simulation system for evaluating gridding urban management strategies | |
CN111199298A (en) | Flood forecasting method and system based on neural network | |
CN112541256A (en) | Deep learning dimensionality reduction reconstruction-based strong heterogeneous reservoir history fitting method | |
CN114925152A (en) | Soil nutrient space interpolation method based on reverse distance weighting method and kriging method | |
CN110457295A (en) | Time series data completion system based on social relationships perception | |
CN111444516A (en) | Sensitivity-based deep reinforcement learning intelligent agent attack method | |
Shen et al. | Online teaching course recommendation based on autoencoder | |
CN114676887A (en) | River water quality prediction method based on graph convolution STG-LSTM | |
CN112270123A (en) | Basin reservoir group runoff random generation method based on convolution generation countermeasure network |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |