The majority of studies on road traffic flow prediction have focused on the flow of passenger cars or the flow of traffic as a whole, which ignore the significant impact of trucks with different sizes and operational characteristics on traffic flow efficiency. Therefore, in this paper, we focus on truck traffic flow and propose a Multifeatures Spatial-Temporal-Based Neural Network model (M-BiCNNGRU) to improve its prediction. The proposed model not only comprises conventional temporal characteristics and spatial relationships but also includes a range of multifeatures. These multifeatures include policy limit, optimal time delay, road resistance, and traffic congestion state. The impacts of upstream and downstream road sections are considered on the spatial relationship by using a Convolutional Neural Network (CNN). A Bidirectional Gated Recurrent Unit (Bi-GRU) is employed to account for the temporal characteristics. To evaluate the proposed model, traffic flow data were collected from a major expressway in Beijing and the results were compared with those derived from existing models. The results show that the prediction accuracy of the BiCNNGRU model, with spatial-temporal characteristics, and the M-BiGRU model, with multifeatures and temporal, is, respectively, 4.13% and 2.15% greater than that of the Bi-GRU model, with temporal characteristics. The prediction accuracy of the proposed M-BiCNNGRU model is 92.86%, which is 7.12% greater than that of the Bi-GRU model and 13.83% greater than that of the Support Vector Regression (SVR) model. In general, therefore, the proposed M-BiCNNGRU model, which combines multifeatures, temporal characteristics, and spatial relationships, can significantly improve accuracy in predicting truck traffic flow.
It is generally acknowledged that the efficient operation of modern supply chains relies heavily on the use of trucks. This is particularly the case in China which has seen a growth of over 40% in truck road traffic during the last decade [
Despite the increasing congestion problems caused by trucks, currently, most research related to congestion is associated with cars [
In terms of traffic flow prediction, the models currently employed can be classified as either parametric models or nonparametric models [
In recent years, not least with improvements in computer performance, the field of deep learning has developed rapidly. Of particular relevance, deep learning methods can capture the characteristics of nonlinear traffic data and achieve excellent results. However, each deep learning model has its fields and characteristics. For example, the Convolutional Neural Network (CNN) has an excellent performance in processing spatial data, and the Recurrent Neural Network (RNN) has advantages in capturing temporal features [
Compared with passenger cars, trucks are required to adhere to special operational policies, mainly in terms of driving time and space. For example, the transportation department stipulates that trucks are prohibited from passing on certain urban roads within 24 hours or trucks are prohibited from passing through certain areas from 7:00–8:30 and 16:30–18:30 during passenger car rush hours. Also, traffic condition, the impedance of the road, and the mixed flow rate of trucks all can affect traffic efficiency but there are relatively few researches into the influence of these factors on truck traffic flow prediction; most rely instead on the use of single-source data from historical traffic flow records.
In general, to improve the prediction accuracy of truck traffic flow, additional factors that affect such flow will be analyzed from the perspective of policy, traffic status, and road impedance. Within this approach, multiple factors are defined as variables to form multiple features, and the input variables are, therefore, changed from one-dimensional traffic flow data to multidimensional feature data. Finally, the spatial relationship between upstream and downstream road sections is added to form a multifeatures spatial-temporal-based joint deep learning model to predict truck traffic flow and evaluate the prediction results. The overall aim and main contributions of this study are summarized as follows: The primary objective of the study is to improve the prediction accuracy of truck traffic flow. Although passenger car traffic flow forecasting has been widely studied, the characteristics of trucks and passenger cars are different, it is not feasible to apply the forecasting method developed for passenger cars to trucks. For this reason, it is necessary to incorporate truck characteristics into the method for predicting truck traffic flow. A Multifeatures Spatial-Temporal-Based Neural Network model (M-BiCNNGRU) is proposed in this study to improve the prediction accuracy of truck traffic flow. In the proposed model, three types of factors are considered, namely, “spatial relationship,” “temporal series,” and “road operation.” Moreover, the impacts of upstream and downstream road sections are considered on the spatial relationship using a CNN. A Bidirectional Gated Recurrent Unit (Bi-GRU) is employed to account for the temporal characteristics. The road operation factor includes a range of multifeatures. These multifeatures include “policy limit,” “optimal time delay,” “road resistance,” and “traffic congestion state.” These multifeatures are also through the CNN model to explore the relationship between the truck flow and multifeatures. The prediction performance of the M-BiCNNGRU model is explored. Several baseline models are also trained and predicted, including (i) the BiCNNGRU model, which considers spatial relationships and temporal characteristics; (ii) the M-BiGRU model, which considers multifeatures and temporal characteristics; (iii) the Bi-GRU model, which considers only temporal characteristics; (iv) the SVR model, which again, only considers temporal characteristics; (v) finally, the ARIMA model, which likewise only considers temporal characteristics.”
This paper is organized as follows. Section
In previous models, traffic flow predictions associated with passenger cars were only inferred from historical traffic flow data. However, other than historical traffic flow, there are multiple other features that affect truck traffic flow, such as current road conditions, road resistance, traffic congestion index, and traffic restriction policies. To consider these additional features, in this paper, we propose a Multifeatures Spatial-Temporal-Based Neural Network prediction method for truck traffic flow, which is termed M-BiCNNGRU. The CNN model is used for multifeatures and spatial feature learning and Bi-GRU model is employed for temporal characteristics learning. The proposed M-BiCNNGRU model structure is shown in Figure
The proposed M-BiCNNGRU model structure.
Five features have been analyzed and incorporated into the model to improve the prediction performance of truck traffic flow. The interpretation and definition of the considered features are as follows: Truck traffic flow: Truck traffic flow has a self-correlation, which means that historical truck traffic flow will have an impact on future truck traffic. The truck traffic flow Policy limit: For different road sections and different periods, truck operation will be restricted by policy, which will cause changes in truck traffic flow. The policy limit variable is defined as a discrete variable. If trucks are allowed to operate during a period, it is defined as 1; if not, it is defined as 0. Optimal time delay: In this paper, the road section where traffic flow is to be predicted is called the target section. For situations when truck flow “upstream” of the target section changes, a period of time is incurred before the traffic information can be disseminated to the target section. Therefore, it is necessary to find a time delay difference In this paper, the Pearson correlation coefficient (PCC) is used to calculate the correlation of truck flow series under a different time delay Road resistance: In previous models, the BPR function was often used to calculate the road resistance of passenger cars. However, for truck traffic, the traditional BPR function is inadequate and led to the development of an improved truck road resistance function ( where Congestion coefficient: At present, the standards for measuring congestion indicators are not uniform. Due to the complex and diverse traffic operating environment of urban roads, the classification of traffic operating status is often not accurate enough and has a certain degree of ambiguity [
The classification of truck type and the Passenger Car Equivalent (PCE) conversion factor.
Truck type | Conversion factor | Load and passenger capacities |
---|---|---|
Small truck | 1 | Load capacity ≤ 2.0 tons |
Medium truck | 1.5 | 2.0 tons < load capacity ≤ 7.0 tons |
Large truck | 2 | 7.0 tons < load capacity ≤ 14.0 tons |
Extra large truck | 3 | Load capacity > 14.0 ton |
Schematic diagram of time delay.
When
For the fuzzy comprehensive evaluation method, the threshold range according to the congestion degree needs to be defined first. For the CC, the larger its value means the more congested the road section is. The CC of a road network is ranked in a descending order and divided into five levels according to a certain proportion. In this paper, the definition of congestion is obtained with reference to Beijing’s “Evaluation Index System for Urban Road Traffic Operation” (EISU) published on April 28, 2011. EISU proposed the recommended conversion relationship between the Road Network Operation Level (RNOL) and Traffic Performance Index (TPI). When the TPI is between the values [0, 2], [2, 4], [4, 6], [6, 8], and [8, 10], the RNOL is defined as “unblocked,” “basically smooth,” “light congestion,” “moderate congestion,” and “severe congestion,” respectively. Moreover, the respective variables associated with these RNOL definitions are 1, 2, 3, 4, and 5. Further details of the above variable definitions can be found in Table
The details of the variable definition.
Symbol | Feature | Type | Description |
---|---|---|---|
Policy limit | Discrete variable | 1: open to traffic; 0: closed to traffic | |
Optimal time delay | Discrete variable | 1: optimal time delay is 1; 2: optimal time delay is 2; 3: optimal time delay is 3 | |
Road resistance | Continue variable | Value of truck road resistance function | |
Truck traffic flow | Continue variable | Value of truck traffic flow | |
Congestion coefficient | Discrete variable | 1: unblocked; 2: basically smooth; 3: light congestion; 4: moderate congestion; 5: the severe congestion |
Studies have shown that the traffic conditions of the upstream and downstream sections have an impact on the target section flow. However, the target road section may have multiple upstream road sections as well as multiple downstream road sections, as shown in Figure
Multiple scenarios linking the target section with the upstream and downstream. (a) One upstream and one downstream section. (b) Multiple upstream and multiple downstream sections.
The feature variables of the upstream section, the downstream section, and the target section are defined as
If the target road section has only one upstream and one downstream road section, as shown in Figure
Therefore, when there are multiple upstream road sections, this paper proposes aggregating the feature variables of multiple upstream sections into one upstream section. The aggregated rule is as follows: For the truck traffic flow variable, it is formed by summing up the multiple upstream sections. For the policy limit, it is set as the policy for the road section with the highest levels of truck traffic flow. For the optimal time delay, CC, and road resistance, they are recalculated according to the sum of the truck traffic flow in the multiple upstream sections.
The rule for the downstream road section is the same as the above upstream section. Through the above processing, multiple scenarios linking the target section with the upstream and downstream section are aggregated into a unified dimension of the input feature vector.
To mine the spatial relationship between the upstream and downstream sections and the target section, the CNN was chosen. This is because CNN has advantages in mining spatial local related information and has achieved significant results in other fields such as natural language processing. The conventional 1D CNN is applied to learn the spatial relationship with the locality structure of key features. The CNN contains a convolutional layer, activation layer, and a pooling layer. By performing 1D convolution on the input feature vector
Truck traffic flow has a long-term temporal dependency. In this paper, a Bi-GRU is introduced to learn the temporal characteristics associated with truck traffic flow. This is because the Bi-GRU has advantages in capturing time characteristics. Recently, an LSTM method has been applied in numerous problems as it, in general, performs well in terms of prediction accuracy, but compared with LSTM, the Bi-GRU contains fewer parameters and has achieved a better performance in the field of natural languages. However, for truck traffic flow, few scholars have studied its potential benefits. Specifically, in the Bi-GRU model, the reset gate and update gate were applied to control the update status of the time series. The reset gate
In summary, the Multifeatures Spatial-Temporal-Based Neural Network prediction method, namely, the M-BiCNNGRU model, is a new joint neural network model. Compared with the CNN model, the proposed M-BiCNNGRU model adds the structure of mining temporal characteristics from the Bi-GRU neural network. As far as the authors are aware, this study is the first to combine the Bi-GRU and CNN neural networks for truck traffic flow tasks. The proposed model has three advantages; firstly, the impact of the truck traffic flow in the upstream and downstream sections on the target road section can be explored using the spatial relationship. Therefore, the road conditions of the upstream and downstream sections can be well represented in the proposed model. Secondly, the Bi-GRU neural network is well capable of mining the temporal characteristics. Through the full connectivity layer, the output result of the spatial relationship is placed in the Bi-GRU model which can then continue to explore the time relationship. Thirdly, the proposed model is more capable of expressing the characteristics of truck flow from the traffic perspective.
The data for the period June 1st to June 5th, 2019, were provided by the traffic survey in Beijing, China, and contains truck traffic flow as well as passenger car traffic flow information. The data were collected by manual counting methods. The data were provided by Beijing Transport Institute, which invited professional survey companies to conduct data collection work, and a total of 30 investigators participated in this work. The time interval of truck traffic flow is one hour, and the size of the data in the CSV format is 3.48 GB. As an important part of Beijing's northwest freight corridor, the Sixth Loop Expressway is selected as the research object and is divided into 1,293 road sections and numbered accordingly. To consider different scenarios, six road sections (S1–S6) were randomly selected as target road sections, as shown in Figure
Truck traffic flow locations.
From Figure
The average daily truck traffic flow.
The properties of the six selected target sections are shown in Table
The properties of the six selected sections.
Section | Mean | Std | Mixed ratio | Peak mean | Rush std | Low peak mean | Low peak std |
---|---|---|---|---|---|---|---|
S1 | 400.84 | 90.22 | 0.49 | 454.79 | 43.93 | 251.25 | 13.88 |
S2 | 216.11 | 38.02 | 0.42 | 232.64 | 25.95 | 150.00 | 12.04 |
S3 | 124.84 | 37.91 | 0.42 | 143.21 | 29.17 | 77.50 | 11.19 |
S4 | 46.21 | 14.90 | 0.21 | 52.21 | 9.31 | 20.75 | 4.02 |
S5 | 277.89 | 75.98 | 0.43 | 327.64 | 22.08 | 163.25 | 26.02 |
S6 | 182.47 | 46.67 | 0.32 | 210.21 | 24.04 | 113.00 | 15.60 |
To show that the upstream and downstream sections have a spatial relationship, PCC is usually applied to study the interaction between sections. The PCC is defined as the quotient of the covariance and standard deviation between the truck traffic flows of two road sections. The value of PCC is between −1 and 1. When it is closer to 1, it indicates that there is a positive correlation between road sections. The PCC with the upstream and downstream sections of the six selected target sections is shown in Figure
The Pearson correction coefficient with the upstream and downstream sections.
As can be seen from Figure
Using the fuzzy evaluation method, the TPI of six road sections is obtained as shown in Figure
The TPI for the six road sections. (a) S1 (b) S2. (c) S3. (d) S4. (e) S5. (f) S6.
The congestion situation in the upstream and downstream will affect the flow of the target road section. However, by calculating the CC, the operating status of the road section can become well understood, which is helpful for the prediction of truck flow. Similarly, the CCs of the upstream and downstream sections of the six sections are also calculated.
To analyze the performance of the proposed model, model training and evaluation of the truck traffic in the six selected road sections were carried out. The first four days of data were used for training purposes, and data from the last day were used for verification, within which the Mean Absolute Error (MAE), Mean Absolute Percentage Error (MAPE), and Root Mean Square Error (RMSE) were applied as parameters for evaluation.
To explore the prediction performance of the M-BiCNNGRU model, several baseline models were also trained and used in predicting truck traffic flows. These model include the following: BiCNNGRU, which considers spatial relationships and temporal characteristics. M-BiGRU, which considers multifeatures and temporal characteristics. Bi-GRU, which considers temporal characteristics. SVR, which considers temporal characteristics. ARIMA, which considers temporal characteristics. For the deep learning algorithm, several superparameters that affect the performance of the proposed model need to be discussed; these superparameters include batch size, hidden unit size, hidden layer number, time steps, and epochs. When the values of the evaluation parameter of MAPE, MAE, and RMSE reach their minimum, the superparameters are selected as the optimal parameters. The learning rate in the study is set to 0.001; the batch size is set to 20. Epochs are 100, 230, 130, 240, 160, and 210, respectively, for the six sections of S1, S2, S3, S4, S5, and S6, as determined from multiple iterations of training to find the optimal parameters. For the training environment, Keras and TensorFlow provided the deep learning packages, and Python 3.6 was used as a general-purpose programming language. These were hosted on a desktop computer, comprising an Intel Xeon ( For choosing the optimal hidden units, the evaluation value with different hidden units was compared. Firstly, the hidden units from [5, 10, 15, 30, 50, 100] were tested, respectively. The results indicated that the best parameter occurs for a unit of 5. Then the hidden units from [1, 2, 3, 4, 5] were used in further testing, where it was found that the respective MAPE value for each hidden unit for the road section S1 is 8.32%, 9.31%, 8.59%, 7.62%, and 9.65%. The value of hidden unit 4 represented by the lowest value of MAPE (7.62%) is, therefore, the optimal parameter of the hidden unit. For different time steps, the proposed M-BiCNNGRU model has different evaluation performances. The evaluation values of MAPE, MAE, and RMSE with different time steps [4, 8, 12, 16, 20, 24] were tested, and the corresponding evaluation distributions are shown in Figure
From Figure
The evaluation values of MAPE, MAE, and RMSE with different time steps.
Forecasting effect of peak and low peak.
Period | Section | S1 | S2 | S3 | S4 | S5 | S6 | Mean |
---|---|---|---|---|---|---|---|---|
Peak | RMSE | 28.88 | 13.68 | 12.84 | 3.30 | 37.40 | 16.88 | 18.83 |
MAE | 24.73 | 11.82 | 10.14 | 2.63 | 32.77 | 11.62 | 15.62 | |
MAPE | 5.35 | 5.10 | 6.73 | 4.78 | 9.37 | 5.26 | 6.10 | |
Low peak | RMSE | 39.23 | 6.69 | 12.79 | 3.90 | 22.10 | 10.08 | 15.80 |
MAE | 32.12 | 5.20 | 8.16 | 3.17 | 19.43 | 8.89 | 12.83 | |
MAPE | 11.53 | 2.79 | 8.14 | 10.37 | 10.31 | 7.07 | 8.37 | |
Difference | RMSE | −10.35 | 6.99 | 0.05 | -0.60 | 15.29 | 6.80 | 3.03 |
MAE | −7.39 | 6.63 | 1.98 | −0.54 | 13.34 | 2.73 | 2.79 | |
MAPE | −6.18 | 2.31 | −1.42 | −5.59 | −0.94 | −1.81 | −2.27 |
The RMSEs, MAPEs, and MAEs of various models.
Model | Evaluation | S1 | S2 | S3 | S4 | S5 | S6 | Mean | |
---|---|---|---|---|---|---|---|---|---|
M-BiCNNGRU | Multifeatures/spatial/temporal | RMSE | 26.46 | 18.71 | 12.82 | 3.53 | 26.09 | 14.74 | 17.06 |
MAPE | 7.62 | 5.50 | 7.25 | 6.84 | 9.71 | 5.93 | 7.14 | ||
MAE | 21.96 | 12.74 | 9.41 | 2.83 | 22.28 | 10.62 | 13.31 | ||
BiCNNGRU | Spatial/temporal | RMSE | 42.04 | 27.05 | 13.88 | 8.79 | 28.82 | 24.02 | 24.10 |
MAPE | 8.42 | 9.95 | 9.92 | 12.61 | 9.16 | 10.72 | 10.13 | ||
MAE | 33.58 | 21.19 | 11.69 | 7.25 | 24.51 | 19.71 | 19.65 | ||
M-BiGRU | Multifeatures/temporal | RMSE | 36.66 | 26.70 | 16.87 | 6.43 | 23.20 | 22.88 | 22.12 |
MAPE | 9.08 | 12.17 | 15.03 | 13.40 | 9.59 | 13.41 | 12.11 | ||
MAE | 28.43 | 18.24 | 13.84 | 5.08 | 18.30 | 17.72 | 16.94 | ||
Bi-GRU | Temporal | RMSE | 46.60 | 23.76 | 14.54 | 9.03 | 30.13 | 24.55 | 24.77 |
MAPE | 11.47 | 15.84 | 15.51 | 18.67 | 10.19 | 13.89 | 14.26 | ||
MAE | 38.20 | 18.80 | 11.25 | 7.85 | 23.11 | 18.39 | 19.60 | ||
SVR | Temporal | RMSE | 58.59 | 49.99 | 27.18 | 25.91 | 51.23 | 34.04 | 41.16 |
MAPE | 14.88 | 21.18 | 22.14 | 28.92 | 15.55 | 23.15 | 20.97 | ||
MAE | 25.02 | 22.54 | 19.19 | 28.10 | 33.14 | 19.29 | 24.55 | ||
ARIMA | Temporal | RMSE | 66.43 | 48.89 | 30.00 | 15.38 | 50.03 | 41.74 | 42.08 |
MAPE | 15.43 | 20.78 | 21.46 | 29.28 | 16.39 | 18.33 | 20.28 | ||
MAE | 53.20 | 34.88 | 23.22 | 11.32 | 40.96 | 29.28 | 32.14 |
Comparison of prediction performance of various models.
It can be seen in Table
The BiCNNGRU model was designed by considering the factors of spatial relationships and temporal characteristics; its purpose is to explore whether the road operation factors of multifeatures play a positive role in improving the prediction accuracy of truck traffic flow. The average MAPE value of the BiCNNGRU truck traffic flow prediction model considering spatial and temporal factors is 10.13, which is 2.99 higher than the proposed M-BiCNNGRU model. However, without the road operation factors of multifeatures, the prediction accuracy reduces by 2.99%; hence, the results show that the operation factors of multifeatures is effective in improving the prediction accuracy of truck flow.
Concerning the M-BiGRU model, it was designed by considering the factors of temporal characteristics and the road operation factors of multifeatures; its purpose is to explore whether spatial relationships of upstream and downstream sections play a positive role in improving prediction accuracy. The average MAPE value of the BiCNNGRU model is 12.11, which is 4.97 higher than the proposed M-BiCNNGRU model. The result shows that, without the factor of spatial relationships, the prediction accuracy is reduced by 4.97%. The reduced accuracy of 4.97% by not considering spatial relationships is higher than the reduced accuracy of 2.99% by not considering multifeatures. This result indicates that, compared with multifeatures, the consideration of spatial relationships has a higher accuracy contribution.
The Bi-GRU model was designed by only considering the factor of temporal characteristics. It is a deep learning model which has been widely studied and applied as a result of its good performance in the mining of temporal characteristics. However, does this good performance also apply in truck traffic flow prediction? How much of an increase in truck flow prediction accuracy can be achieved using the Bi-GRU model to mine the temporal characteristics compared to the parameter-type model ARIMA and machine learning model SVR?
The average MAPE of the Bi-GRU model is 14.26, which is 6.71 and 6.02 lower than that of SVR model and ARIMA model, respectively, showing that the Bi-GRU model selected in this paper to capture time characteristics has a good prediction performance and perform better than the parameter-type models and machine learning models. It also indicates that the Bi-GRU model can play a fundamental role in obtaining a high prediction accuracy in the constructed M-BiCNNGRU model.
Furthermore, the average MAPE of the M-BiGRU model considering multifeatures and temporal characteristics is 12.11, which is 2.14 lower than that of Bi-GRU model only considering temporal characteristics. It also indicates that the operation factor of multifeatures is helpful in improving truck flow prediction accuracy, and the improved prediction accuracy value is 2.14%. The average MAPE of the BiCNNGRU model considering the spatial relationship and temporal characteristics is 4.13 lower than the Bi-GRU model only considering temporal characteristics. It also indicates that the factor of spatial relationship is conducive to improving the truck flow prediction accuracy, and the improved prediction accuracy value is 4.13%, which is higher than the improved prediction accuracy value 2.14% of multifeatures. The result shows that the contribution to improving the prediction accuracy of spatial relationship factor is greater than the operation factor of multifeatures.
In summary, the proposed M-BiCNNGRU model, considering three type factors of spatial relationship, temporal characteristics, and road operation multifeatures, outperforms other models in predicting truck traffic flow. Besides, for the contribution to improving the prediction accuracy, the temporal characteristics are the biggest, followed by the spatial relationship and then the road operation factor of multifeatures.
To further discuss the performance of the M-BiCNNGRU model, the original truck traffic flow is compared with the predicted truck traffic flow values. At the same time, this paper divides the day into two periods, low peak period (4:00–10:00) and peak period (11:00–22:00), and compares the prediction results of truck traffic flow in the two periods. The comparison results are shown in Figure
Comparison of the predicted truck flow of peak and low peak.
For an improved quantitative analysis of the truck traffic flow prediction model in both peak and low peak periods, MAE, MAPE, and RMSE indexes were also used for evaluation. The results are shown in Table
It can be seen from Figure
Prediction performance of peak and low peak.
The paper proposes a Multifeatures Spatial-Temporal-Based model (M-BiCNNGRU) for the accurate forecasting of truck traffic flow. The multifeatures include factors such as truck operational policy, road resistance, truck traffic flow rate, traffic congestion state, and optimal time delay which were all incorporated into the model. The novelty of the model is that it includes the impact factors of multifeatures, as well as temporal characteristics and spatial relationships. The model was trained using data from six road sections from a major expressway network in China collected over a five-day period, where the final day’s data were excluded from the training but were instead employed in verifying the performance of the model. To gain further insight into the prediction results of the M-BiCNNGRU model, five baseline models, namely, BiCNNGRU, M-BiGRU, Bi-GRU, SVR, and ARIMA were also applied. The results show that the prediction accuracy of the BiCNNGRU model is 89.87%, which is 4.13%, 10.84%, and 10.15% higher than that of the Bi-GRU, SVR, and ARIMA models. This indicates that the prediction performance of truck traffic flow can be improved by increasing the spatial relationship between upstream and downstream sections. The prediction accuracy of the M-BiGRU model is 87.89%, which is 2.15%, 8.86%, and 8.17% higher than that of the Bi-GRU, SVR, and ARIMA models. This indicates that the incorporation of multifeatures can help to improve prediction accuracy. The prediction accuracy of the M-BiCNNGRU is 92.86%, which is 13.83% and 13.14% higher than that of the SVR and ARIMA models with single temporal characteristics. In summary, the proposed model based on multifeatures, temporal characteristics, and spatial relationships outperforms baseline algorithms as demonstrated in this paper. However, in this paper, the contributions of each individual feature to the overall prediction effectiveness were not specifically analyzed but it is suggested that this could form part of a future study to determine the influence of each feature in terms of weighting factors.
Some data and code used during the study are available in a repository in accordance with funder data retention policies (
The authors declare that they have no conflicts of interest.
The authors confirm contribution to the paper as follows: Chunfu Shao and Shengyou Wang conceptualized and designed the study; Yan Zheng performed data collection; Shengyou Wang, Song Xue, and Yajiao Zhai performed analysis and interpretation of results; Shengyou Wang and Yajiao Zhai prepared the original manuscript. All authors reviewed the results and approved the final version of the manuscript.
The research was supported by funding provided by the Fundamental Research Funds for the Central Universities (Grant no. 2019RC027) and the National Key R&D Program of China (Grant no. 2019YFF0301403). The authors would like to thank Professor Chunjiaodong from Beijing Jiaotong University and Dr. Chaoying Yin from Nanjing Forestry University for their assistance in editing the paper and providing data.