# Real-Time Corrected Traffic Correlation Model for Traffic Flow Forecasting

Authors: Hua-pu Lu; Zhi-yuan Sun; Wen-cong Qu; Ling Wang
Journal: Mathematical Problems in Engineering
 (2015)

---

## Abstract

This paper focuses on the problems of short-term traffic flow forecasting. The main goal is to put forward traffic correlation model and real-time correction algorithm for traffic flow forecasting. Traffic correlation model is established based on the temporal-spatial-historical correlation characteristic of traffic big data. In order to simplify the traffic correlation model, this paper presents correction coefficients optimization algorithm. Considering multistate characteristic of traffic big data, a dynamic part is added to traffic correlation model. Real-time correction algorithm based on Fuzzy Neural Network is presented to overcome the nonlinear mapping problems. A case study based on a real-world road network in Beijing, China, is implemented to test the efficiency and applicability of the proposed modeling methods.


---

## Body

## 1. Introduction

It is of practical significance to predict traffic flow quickly, precisely, and timely. Short-term traffic flow forecasting provides an important basis for traffic guidance and control. Existing studies of short-term traffic flow forecasting can be classified into six categories in transportation literature:(a)
linear system theory based models, such as Autoregressive Integrated Moving Average (ARIMA) [1] and Kalman Filtering model [2];
(b)
data mining based models, such as Neural Network [3], Nonparametric Regression [4], and Support Vector Machine [5];
(c)
nonlinear system theory based models, such as Wavelet Analysis [6], Catastrophe Theory [7], and Chaos Theory [8];
(d)
simulation based models [9];
(e)
combination model based models [10];
(f)
the other models.In the era of big data, it brings both opportunities and challenges to short-term traffic flow forecasting. During data processing, traffic big data meets the same difficulties with the general big data, such as capture, storage, search, sharing, analytics, and visualization. Therefore, short-term traffic flow forecasting method needs to have the capacity to deal with traffic big data. Traffic big data holds several characteristics, such as temporal correlation, spatial correlation, historical correlation, and multistate. Considering the advantages of traffic big data, data-driven based mathematical models can be set up. The physical meaning of these models can by described clearly. In addition, we can put forward real-time correction algorithm to improve the accuracy of traffic flow forecasting.However, taking into account all the present researches in this field, there is still a lack of consideration of traffic big data and real-time correction for traffic flow forecasting. Further researches remain to be conducted on the direction of traffic big data analysis. In this paper, the method of short-term traffic flow forecasting is proposed in detail. The remainder of this paper is organized as follows. Section2 presents basic mathematical model. In Section 3, real-time corrected traffic correlation model is established. A case study based on a real-world road network is carried out in Section 4 to demonstrate the performance and applicability of the proposed method. Finally, conclusions are drawn in Section 5.


## 2. Basic Mathematical Model

### 2.1. Big Data Driven Based Traffic Correlation Model

Traffic big data has a strong temporal-spatial-historical correlation as follows.(i)
In the temporal series, the traffic flow of last moment can be regarded as a continuation of current traffic flow. Dynamic traffic flow data continuously change over time with a certain trend.
(ii)
In the spatial series, the traffic flow of downstream sections can be seen as a continuation of the upstream traffic flow. There exists a spatial association between traffic flow data of neighboring junctions or sections and that of target junctions or sections.
(iii)
In the historical series, the traffic demand characteristics determine that traffic flow characteristics of the same day in the same period are similar. The law of traffic flow cycle is especially evident.Therefore, the basic form of traffic correlation model [11] is expressed as(1)
V
s
t
=
γ
T
V
T
t
+
γ
S
V
S
t
+
γ
H
V
H
t
,
s
.
t
.
γ
T
+
γ
S
+
γ
H
=
1
,
γ
T
,
γ
S
,
γ
H
≥
0
,where V
s
t is the traffic flow parameter of section s at time t, representing flow q
s
t, speed v
s
t, or occupancy ρ
s
t. V
T
t, V
S
t, and V
H
t are the estimated value of V
s
t. γ
T, γ
S, and γ
H are coefficients of these three variables.V
T
t is calculated by temporal correlation analysis, generally based on Regression Analysis Model [12]. V
S
t is calculated by spatial correlation analysis, generally based on Neighbor Regression Model [13]. V
H
t is calculated by historical correlation analysis, generally based on Discrete Fourier Transform Model [14].Thus, simplified equation of traffic correlation model is obtained:(2)
V
s
t
=
∑
∑
ω
j
d
V
j
t
-
d
,
s
.
t
.
∑
ω
=
1
,
ω
>
0
,
η
∈
N
*
,where ω
j
d is regression coefficient of V
j
t
-
d. η is the number of ω
j
d.


### 2.2. Correction Coefficients Optimization Algorithm

It is found that the speed and accuracy of data processing are both important for big data driven method. To improve the speed of traffic correlation model, the number of unknown variables in formulation (2) must be reduced. However, variables reduction may decrease the accuracy of the model.Therefore, this paper defines a threshold of computing speed and derives the maximum of acceptable number of variables. Thus, the number of unknown variableη can be achieved. The correlation coefficients (r
j
d) between the studied section s and the related section j are used to choose variables, the number of which is η:(3)
r
j
d
=
∑
t
=
1
n
-
d
V
s
t
V
j
t
-
d
-
1
n
-
d
∑
t
=
1
n
-
d
V
s
t
∑
t
=
1
n
-
d
V
j
t
-
d
r
j
d
W
·
∑
t
=
1
n
-
d
V
s
2
t
-
1
n
-
d
∑
t
=
1
n
-
d
V
s
t
2
W
r
j
d
W
×
∑
t
=
1
n
-
d
V
j
2
t
-
d
-
1
n
-
d
∑
t
=
1
n
-
d
V
j
t
-
d
2
-
1
/
2
,
s
.
t
.
d
≥
0
.In addition, forV
j
·, a unique value of V
j
t
-
d is determined. Therefore, a lot of variables are reduced. When the max of correlation coefficient (r
j) corresponding to each section j is calculated, corresponding time delay (d) and the unique value of V
j
t
-
d are obtained:(4)
r
j
=
max
⁡
r
j
d
.If the value ofr
j is large enough, the variable V
j
t
-
d is preserved. Otherwise, the variable V
j
t
-
d is reduced. After several data tests, it is found that the values of r
j are not very different with relatively high values when η is less than 4; the rapid decay of r
j is observed with relatively low values when η is more than 8. Therefore, η
∈
4,8.Variables reduction makesω
j
d meaningless. So, a new variable ω
j, which is normalized r
j, is selected to replace the variable ω
j
d:(5)
ω
j
=
r
j
∑
j
=
1
η
r
j
,
∑
j
=
1
η
ω
j
=
1
,
ω
j
>
0
.Since the alternative process will bring some errors, which are likely to be systematic, a linear correction algorithm is present. Two correction variablesφ and ε are introduced for calibration error. Simplified traffic correlation model is(6)
min
⁡
V
s
t
-
V
^
s
t
,
i
s
.
t
.
V
s
t
=
φ
∑
j
=
1
η
ω
j
V
j
t
-
d
+
ε
,where V
^
s
t is the actual value of V
s
t.


## 2.1. Big Data Driven Based Traffic Correlation Model

Traffic big data has a strong temporal-spatial-historical correlation as follows.(i)
In the temporal series, the traffic flow of last moment can be regarded as a continuation of current traffic flow. Dynamic traffic flow data continuously change over time with a certain trend.
(ii)
In the spatial series, the traffic flow of downstream sections can be seen as a continuation of the upstream traffic flow. There exists a spatial association between traffic flow data of neighboring junctions or sections and that of target junctions or sections.
(iii)
In the historical series, the traffic demand characteristics determine that traffic flow characteristics of the same day in the same period are similar. The law of traffic flow cycle is especially evident.Therefore, the basic form of traffic correlation model [11] is expressed as(1)
V
s
t
=
γ
T
V
T
t
+
γ
S
V
S
t
+
γ
H
V
H
t
,
s
.
t
.
γ
T
+
γ
S
+
γ
H
=
1
,
γ
T
,
γ
S
,
γ
H
≥
0
,where V
s
t is the traffic flow parameter of section s at time t, representing flow q
s
t, speed v
s
t, or occupancy ρ
s
t. V
T
t, V
S
t, and V
H
t are the estimated value of V
s
t. γ
T, γ
S, and γ
H are coefficients of these three variables.V
T
t is calculated by temporal correlation analysis, generally based on Regression Analysis Model [12]. V
S
t is calculated by spatial correlation analysis, generally based on Neighbor Regression Model [13]. V
H
t is calculated by historical correlation analysis, generally based on Discrete Fourier Transform Model [14].Thus, simplified equation of traffic correlation model is obtained:(2)
V
s
t
=
∑
∑
ω
j
d
V
j
t
-
d
,
s
.
t
.
∑
ω
=
1
,
ω
>
0
,
η
∈
N
*
,where ω
j
d is regression coefficient of V
j
t
-
d. η is the number of ω
j
d.


## 2.2. Correction Coefficients Optimization Algorithm

It is found that the speed and accuracy of data processing are both important for big data driven method. To improve the speed of traffic correlation model, the number of unknown variables in formulation (2) must be reduced. However, variables reduction may decrease the accuracy of the model.Therefore, this paper defines a threshold of computing speed and derives the maximum of acceptable number of variables. Thus, the number of unknown variableη can be achieved. The correlation coefficients (r
j
d) between the studied section s and the related section j are used to choose variables, the number of which is η:(3)
r
j
d
=
∑
t
=
1
n
-
d
V
s
t
V
j
t
-
d
-
1
n
-
d
∑
t
=
1
n
-
d
V
s
t
∑
t
=
1
n
-
d
V
j
t
-
d
r
j
d
W
·
∑
t
=
1
n
-
d
V
s
2
t
-
1
n
-
d
∑
t
=
1
n
-
d
V
s
t
2
W
r
j
d
W
×
∑
t
=
1
n
-
d
V
j
2
t
-
d
-
1
n
-
d
∑
t
=
1
n
-
d
V
j
t
-
d
2
-
1
/
2
,
s
.
t
.
d
≥
0
.In addition, forV
j
·, a unique value of V
j
t
-
d is determined. Therefore, a lot of variables are reduced. When the max of correlation coefficient (r
j) corresponding to each section j is calculated, corresponding time delay (d) and the unique value of V
j
t
-
d are obtained:(4)
r
j
=
max
⁡
r
j
d
.If the value ofr
j is large enough, the variable V
j
t
-
d is preserved. Otherwise, the variable V
j
t
-
d is reduced. After several data tests, it is found that the values of r
j are not very different with relatively high values when η is less than 4; the rapid decay of r
j is observed with relatively low values when η is more than 8. Therefore, η
∈
4,8.Variables reduction makesω
j
d meaningless. So, a new variable ω
j, which is normalized r
j, is selected to replace the variable ω
j
d:(5)
ω
j
=
r
j
∑
j
=
1
η
r
j
,
∑
j
=
1
η
ω
j
=
1
,
ω
j
>
0
.Since the alternative process will bring some errors, which are likely to be systematic, a linear correction algorithm is present. Two correction variablesφ and ε are introduced for calibration error. Simplified traffic correlation model is(6)
min
⁡
V
s
t
-
V
^
s
t
,
i
s
.
t
.
V
s
t
=
φ
∑
j
=
1
η
ω
j
V
j
t
-
d
+
ε
,where V
^
s
t is the actual value of V
s
t.


## 3. Model Improvement

### 3.1. Real-Time Correction Problem Statement

Basic mathematical model can be used for traffic flow forecasting. The error of traffic flow forecasting is written as(7)
Δ
t
=
V
^
s
t
-
V
s
t
=
V
^
s
t
-
φ
∑
j
=
1
η
ω
j
V
j
t
-
d
+
ε
.Therefore,V
^
s
t is achieved:(8)
V
^
s
t
=
φ
∑
j
=
1
η
ω
j
V
j
t
-
d
+
ε
+
Δ
t
.The variation range ofΔ
t reflects the accuracy of traffic flow forecasting. To improve the accuracy of traffic flow forecasting, Δ
t is made as a compensation variable for V
s
t. Then, V
^
s
t is replaced by V
s
t in formulation (8):(9)
V
s
t
=
φ
∑
j
=
1
η
ω
j
V
j
t
-
d
+
ε
+
Δ
t
.For different traffic state, the propagation of traffic congestion is different. So, the temporal-spatial-historical correlation variables are dynamic. As shown in Section2, basic mathematical model for traffic flow forecasting is put forward based on temporal-spatial-historical correlation characteristic of traffic big data. However, the multistate characteristic is largely ignored; the temporal-spatial-historical correlation variables are seen as static variables. Although a linear fitting method is used to improve accuracy, the error which is the main part of Δ
t still exists.Effective analysis of traffic correlation model is shown in Figure1. The error of traffic congestion stage is larger than that of traffic smooth stage. The error of last moment may affect the error of current moment.Figure 1
Effective analysis of traffic correlation model.Δ
t responds to changes with traffic flow state and temporal series. It is difficult for mathematical model to describe the characteristics of Δ
t. This paper tries to present real-time correction algorithm based on nonlinear mapping.It assumes thatΔ
t is influenced by level of service at time t
-
1 (LOS
t
-
1) and Δ
t
-
1. Then, the nonlinear mapping can be written as(10)
f
:
Δ
t
-
1
,
LOS
t
-
1
⟶
Δ
t
.


### 3.2. Real-Time Correction Algorithm

The main goal of real-time correction algorithm is to calculate the error term (Δ
t) in formulation (9). Because of nonlinear mapping, Fuzzy Neural Network (FNN) can be used to overcome this problem. The structure of Fuzzy Neural Network is shown in Figure 2.Figure 2
The structure of Fuzzy Neural Network.Every input unit in Fuzzy Neural Network is corresponding to certain fuzzy subset of the input variables, which areLOS
t
-
1 and Δ
t
-
1. The format of input signal is(11)
μ
A
LOS
,
μ
B
LOS
,
μ
C
LOS
,
μ
D
LOS
,
μ
E
LOS
,
E
μ
F
LOS
,
μ
V
S
Δ
,
μ
S
Δ
,
μ
M
Δ
,
μ
L
Δ
,
μ
V
L
Δ
.Every output unit in Fuzzy Neural Network is corresponding to certain fuzzy subset of the output variable, which isΔ
t. The format of output signal is(12)
μ
V
S
Δ
,
μ
S
Δ
,
μ
M
Δ
,
μ
L
Δ
,
μ
V
L
Δ
.To get membership degree model ofLOS and Δ, SAGA-FCM is used to get the clustering centers. LOS is divided into six types [15], as LOS
=
A
,
B
,
C
,
D
,
E
,
F. The membership degree of LOS is calculated by(13)
μ
i
LOS
t
=
x
t
-
a
i
x
t
-
a
A
+
x
t
-
a
i
x
t
-
a
B
μ
i
LOS
t
W
=
+
⋯
+
x
t
-
a
i
x
t
-
a
F
-
1
,
s
.
t
.
i
∈
A
,
B
,
C
,
D
,
E
,
F
,where μ
i
LOS
t is the membership degree of which sample data x
t that belongs to LOS
i, x
t
=
v
t
,
ρ
t. a is the set of clustering centers for LOS, a
=
a
A
,
a
B
,
a
C
,
a
D
,
a
E
,
a
F, and a
i
=
v
i
,
ρ
i.Δ is divided into five types, as Δ
=
V
S
,
S
,
M
,
L
,
V
L. The membership degree model of Δ is calculated by(14)
μ
k
Δ
t
=
Δ
t
-
b
k
Δ
t
-
b
V
S
+
Δ
t
-
b
k
Δ
t
-
b
S
μ
k
Δ
t
=
W
+
⋯
+
Δ
t
-
b
k
Δ
t
-
b
V
L
-
1
,
s
.
t
.
k
∈
V
S
,
S
,
M
,
L
,
V
L
,where μ
k
Δ
t is the membership degree of which sample data Δ
t belongs to Δ
k. b is the set of clustering centers for Δ, b
=
b
V
S
,
b
S
,
b
M
,
b
L
,
b
V
L.Historical date is applied to train the Fuzzy Neural Network. Input signal is corresponding to target output. After training, the network can be seen as a container of fuzzy relations and if we want to get other conclusions from the network, the only thing that needs to be done is to input the real value after defuzzification, as shown in(15)
Δ
t
=
∑
μ
k
Δ
t
·
b
k
∑
μ
k
Δ
t
.


### 3.3. Real-Time Corrected Traffic Correlation Model

Real-time corrected traffic correlation model, which is seen as the improved traffic correlation model, is composed of static part and dynamic part. The static part is “φ
∑
j
=
1
η
ω
j
V
j
t
-
d
+
ε,” which shows the physical meaning of traffic correlation. The dynamic part is “Δ
t,” which shows the physical meaning of traffic multistate characteristic. Traffic flow forecasting framework is outlined in Figure 3.Figure 3
Traffic flow forecasting framework.The steps of traffic flow forecasting are as follows.Step 1 (traffic correlation model).
Based on historical data, temporal data, and spatial data, basic traffic correlation model, as shown in formulation (1), is built.Step 2 (simplified traffic correlation model).
Simplified traffic correlation model, as shown in formulation (6), is proposed based on correction coefficients optimization algorithm. Thus, “φ
∑
j
=
1
η
ω
j
V
j
t
-
d
+
ε” is obtained.Step 3 (real-time correction algorithm).
Based on Fuzzy Neural Network, real-time correction algorithm is put forward. Thus, “Δ
t” is obtained.Step 4 (real-time corrected traffic correlation model).
Real-time corrected traffic correlation model, as shown in formulation (8), is stored in the database as system knowledge.Step 5 (short-term traffic flow forecasting).
Based on system knowledge, real-time data is processed to calculate the forecasting results.


## 3.1. Real-Time Correction Problem Statement

Basic mathematical model can be used for traffic flow forecasting. The error of traffic flow forecasting is written as(7)
Δ
t
=
V
^
s
t
-
V
s
t
=
V
^
s
t
-
φ
∑
j
=
1
η
ω
j
V
j
t
-
d
+
ε
.Therefore,V
^
s
t is achieved:(8)
V
^
s
t
=
φ
∑
j
=
1
η
ω
j
V
j
t
-
d
+
ε
+
Δ
t
.The variation range ofΔ
t reflects the accuracy of traffic flow forecasting. To improve the accuracy of traffic flow forecasting, Δ
t is made as a compensation variable for V
s
t. Then, V
^
s
t is replaced by V
s
t in formulation (8):(9)
V
s
t
=
φ
∑
j
=
1
η
ω
j
V
j
t
-
d
+
ε
+
Δ
t
.For different traffic state, the propagation of traffic congestion is different. So, the temporal-spatial-historical correlation variables are dynamic. As shown in Section2, basic mathematical model for traffic flow forecasting is put forward based on temporal-spatial-historical correlation characteristic of traffic big data. However, the multistate characteristic is largely ignored; the temporal-spatial-historical correlation variables are seen as static variables. Although a linear fitting method is used to improve accuracy, the error which is the main part of Δ
t still exists.Effective analysis of traffic correlation model is shown in Figure1. The error of traffic congestion stage is larger than that of traffic smooth stage. The error of last moment may affect the error of current moment.Figure 1
Effective analysis of traffic correlation model.Δ
t responds to changes with traffic flow state and temporal series. It is difficult for mathematical model to describe the characteristics of Δ
t. This paper tries to present real-time correction algorithm based on nonlinear mapping.It assumes thatΔ
t is influenced by level of service at time t
-
1 (LOS
t
-
1) and Δ
t
-
1. Then, the nonlinear mapping can be written as(10)
f
:
Δ
t
-
1
,
LOS
t
-
1
⟶
Δ
t
.


## 3.2. Real-Time Correction Algorithm

The main goal of real-time correction algorithm is to calculate the error term (Δ
t) in formulation (9). Because of nonlinear mapping, Fuzzy Neural Network (FNN) can be used to overcome this problem. The structure of Fuzzy Neural Network is shown in Figure 2.Figure 2
The structure of Fuzzy Neural Network.Every input unit in Fuzzy Neural Network is corresponding to certain fuzzy subset of the input variables, which areLOS
t
-
1 and Δ
t
-
1. The format of input signal is(11)
μ
A
LOS
,
μ
B
LOS
,
μ
C
LOS
,
μ
D
LOS
,
μ
E
LOS
,
E
μ
F
LOS
,
μ
V
S
Δ
,
μ
S
Δ
,
μ
M
Δ
,
μ
L
Δ
,
μ
V
L
Δ
.Every output unit in Fuzzy Neural Network is corresponding to certain fuzzy subset of the output variable, which isΔ
t. The format of output signal is(12)
μ
V
S
Δ
,
μ
S
Δ
,
μ
M
Δ
,
μ
L
Δ
,
μ
V
L
Δ
.To get membership degree model ofLOS and Δ, SAGA-FCM is used to get the clustering centers. LOS is divided into six types [15], as LOS
=
A
,
B
,
C
,
D
,
E
,
F. The membership degree of LOS is calculated by(13)
μ
i
LOS
t
=
x
t
-
a
i
x
t
-
a
A
+
x
t
-
a
i
x
t
-
a
B
μ
i
LOS
t
W
=
+
⋯
+
x
t
-
a
i
x
t
-
a
F
-
1
,
s
.
t
.
i
∈
A
,
B
,
C
,
D
,
E
,
F
,where μ
i
LOS
t is the membership degree of which sample data x
t that belongs to LOS
i, x
t
=
v
t
,
ρ
t. a is the set of clustering centers for LOS, a
=
a
A
,
a
B
,
a
C
,
a
D
,
a
E
,
a
F, and a
i
=
v
i
,
ρ
i.Δ is divided into five types, as Δ
=
V
S
,
S
,
M
,
L
,
V
L. The membership degree model of Δ is calculated by(14)
μ
k
Δ
t
=
Δ
t
-
b
k
Δ
t
-
b
V
S
+
Δ
t
-
b
k
Δ
t
-
b
S
μ
k
Δ
t
=
W
+
⋯
+
Δ
t
-
b
k
Δ
t
-
b
V
L
-
1
,
s
.
t
.
k
∈
V
S
,
S
,
M
,
L
,
V
L
,where μ
k
Δ
t is the membership degree of which sample data Δ
t belongs to Δ
k. b is the set of clustering centers for Δ, b
=
b
V
S
,
b
S
,
b
M
,
b
L
,
b
V
L.Historical date is applied to train the Fuzzy Neural Network. Input signal is corresponding to target output. After training, the network can be seen as a container of fuzzy relations and if we want to get other conclusions from the network, the only thing that needs to be done is to input the real value after defuzzification, as shown in(15)
Δ
t
=
∑
μ
k
Δ
t
·
b
k
∑
μ
k
Δ
t
.


## 3.3. Real-Time Corrected Traffic Correlation Model

Real-time corrected traffic correlation model, which is seen as the improved traffic correlation model, is composed of static part and dynamic part. The static part is “φ
∑
j
=
1
η
ω
j
V
j
t
-
d
+
ε,” which shows the physical meaning of traffic correlation. The dynamic part is “Δ
t,” which shows the physical meaning of traffic multistate characteristic. Traffic flow forecasting framework is outlined in Figure 3.Figure 3
Traffic flow forecasting framework.The steps of traffic flow forecasting are as follows.Step 1 (traffic correlation model).
Based on historical data, temporal data, and spatial data, basic traffic correlation model, as shown in formulation (1), is built.Step 2 (simplified traffic correlation model).
Simplified traffic correlation model, as shown in formulation (6), is proposed based on correction coefficients optimization algorithm. Thus, “φ
∑
j
=
1
η
ω
j
V
j
t
-
d
+
ε” is obtained.Step 3 (real-time correction algorithm).
Based on Fuzzy Neural Network, real-time correction algorithm is put forward. Thus, “Δ
t” is obtained.Step 4 (real-time corrected traffic correlation model).
Real-time corrected traffic correlation model, as shown in formulation (8), is stored in the database as system knowledge.Step 5 (short-term traffic flow forecasting).
Based on system knowledge, real-time data is processed to calculate the forecasting results.


## 4. Case Study

### 4.1. Data Characteristics

Taking a section of the Second Ring Road (Section  1, as shown in Figure4) and its surrounding trunk road in Beijing, China, as the object of study, it verifies the effectiveness and feasibility of the proposed method. Basic traffic flow data are detected by microwave detectors. Interval time of traffic flow data is 5 min.Figure 4
Spatial location of research object.


### 4.2. Effective Analysis

#### 4.2.1. Basic Mathematical Model

Based on historical data, basic mathematical model, including traffic correlation model and simplified traffic correlation model, is built. The static part of real-time corrected traffic correlation model, which is “φ
∑
j
=
1
η
ω
j
V
j
t
-
d
+
ε,” is achieved, as shown in (16)
q
1
t
=
1.0629
×
0.2246
·
q
9
t
+
0.2151
·
q
1
t
-
5
q
1
t
=
1.0629
×
w
+
0.2145
·
q
1
H
t
-
5
+
0.1787
·
q
2
t
q
1
t
=
1.0629
×
w
+
0.1670
·
q
6
t
-
10
+
22.652
,
v
1
t
=
1.029
×
0.2677
·
v
9
t
+
0.2630
·
v
1
t
-
5
v
1
t
=
1.029
×
w
+
0.2461
·
v
1
H
t
+
0.2232
·
v
2
t
v
1
t
=
+
1.8173
,
ρ
1
t
=
1.0148
×
0.2791
·
ρ
9
t
+
0.2747
·
ρ
1
t
-
5
ρ
1
t
=
1.0148
×
w
+
0.2294
·
ρ
2
t
+
0.2168
·
ρ
1
H
t
ρ
1
t
=
+
0.4793
.Goodness of Fit is shown in Table1, in which the value of R
2 shows the high effectiveness of formulation (16).Table 1
Goodness of Fit.
Parameter of traffic flow
Flow
Speed
Occupancy
R
2
0.8885
0.9655
0.9664


#### 4.2.2. Model Improvement

Real-time correction algorithm is presented to obtain the dynamic part of real-time corrected traffic correlation model, which is “Δ
t.” The steps are as follows.Step 1 (clustering center search).
Based on SAGA-FCM,LOS is divided into six types; Δ is divided into five types. Clustering Centers are shown in Tables 2 and 3.Table 2
Clustering centers of LOS.
i
Speed (km/h)
Occupancy (%)
A
68.1
2.1
B
58.3
3.1
C
51.2
11.4
D
46.0
14.9
E
38.0
20.5
F
23.4
35.1Table 3
Clustering centers ofΔ.
k
Flow (pcu/h)
Speed (km/h)
Occupancy (%)
V
S
−1043
−9.7998
−7.4513
S
−498
−2.5634
−2.2416
M
−2.9
0.2030
−0.2448
L
291.6
2.3692
1.6348
V
L
812.2
9.9043
7.0243Step 2 (calculation of membership degree).
Formulations are shown is Section3.2.Step 3 (training Fuzzy Neural Network).
Input layer has 11 nodes; the transmission function of the nerve cells in the hidden layer is transig; the output layer has 5 nodes and the transmission function of the nerve cells in the output layer is logsig, while the training function is traingdx. Historical data is used to train Fuzzy Neural Network.


#### 4.2.3. Short-Term Traffic Flow Forecasting

Making one day as an example, the result is shown in Figure5.Figure 5
Effective analysis of real-time correction algorithm.Mean Absolute Percentage Error (MAPE) is used to evaluate the accuracy of the proposed model:(17)
MAPE
=
1
n
∑
t
=
1
n
V
^
s
t
-
V
s
t
V
^
s
t
×
100
%
.Evaluation result is shown in Table4. Comparing with these two kinds of model, real-time corrected traffic correlation model is better than the basic mathematical model. That is, the static part of real-time corrected traffic correlation model reflects the temporal-spatial-historical correlation by way of mathematical model; the dynamic part of real-time corrected traffic correlation model modifies the mathematical model based on multistate characteristics of traffic big data.Table 4
Evaluation result of proposed models.
Model
MAPE of flow (%)
MAPE of speed (%)
MAPE of occupancy (%)
Basic mathematical model
13.14%
9.18%
14.48%
Model improvement
10.93%
7.17%
10.54%


## 4.1. Data Characteristics

Taking a section of the Second Ring Road (Section  1, as shown in Figure4) and its surrounding trunk road in Beijing, China, as the object of study, it verifies the effectiveness and feasibility of the proposed method. Basic traffic flow data are detected by microwave detectors. Interval time of traffic flow data is 5 min.Figure 4
Spatial location of research object.


## 4.2. Effective Analysis

### 4.2.1. Basic Mathematical Model

Based on historical data, basic mathematical model, including traffic correlation model and simplified traffic correlation model, is built. The static part of real-time corrected traffic correlation model, which is “φ
∑
j
=
1
η
ω
j
V
j
t
-
d
+
ε,” is achieved, as shown in (16)
q
1
t
=
1.0629
×
0.2246
·
q
9
t
+
0.2151
·
q
1
t
-
5
q
1
t
=
1.0629
×
w
+
0.2145
·
q
1
H
t
-
5
+
0.1787
·
q
2
t
q
1
t
=
1.0629
×
w
+
0.1670
·
q
6
t
-
10
+
22.652
,
v
1
t
=
1.029
×
0.2677
·
v
9
t
+
0.2630
·
v
1
t
-
5
v
1
t
=
1.029
×
w
+
0.2461
·
v
1
H
t
+
0.2232
·
v
2
t
v
1
t
=
+
1.8173
,
ρ
1
t
=
1.0148
×
0.2791
·
ρ
9
t
+
0.2747
·
ρ
1
t
-
5
ρ
1
t
=
1.0148
×
w
+
0.2294
·
ρ
2
t
+
0.2168
·
ρ
1
H
t
ρ
1
t
=
+
0.4793
.Goodness of Fit is shown in Table1, in which the value of R
2 shows the high effectiveness of formulation (16).Table 1
Goodness of Fit.
Parameter of traffic flow
Flow
Speed
Occupancy
R
2
0.8885
0.9655
0.9664


### 4.2.2. Model Improvement

Real-time correction algorithm is presented to obtain the dynamic part of real-time corrected traffic correlation model, which is “Δ
t.” The steps are as follows.Step 1 (clustering center search).
Based on SAGA-FCM,LOS is divided into six types; Δ is divided into five types. Clustering Centers are shown in Tables 2 and 3.Table 2
Clustering centers of LOS.
i
Speed (km/h)
Occupancy (%)
A
68.1
2.1
B
58.3
3.1
C
51.2
11.4
D
46.0
14.9
E
38.0
20.5
F
23.4
35.1Table 3
Clustering centers ofΔ.
k
Flow (pcu/h)
Speed (km/h)
Occupancy (%)
V
S
−1043
−9.7998
−7.4513
S
−498
−2.5634
−2.2416
M
−2.9
0.2030
−0.2448
L
291.6
2.3692
1.6348
V
L
812.2
9.9043
7.0243Step 2 (calculation of membership degree).
Formulations are shown is Section3.2.Step 3 (training Fuzzy Neural Network).
Input layer has 11 nodes; the transmission function of the nerve cells in the hidden layer is transig; the output layer has 5 nodes and the transmission function of the nerve cells in the output layer is logsig, while the training function is traingdx. Historical data is used to train Fuzzy Neural Network.


### 4.2.3. Short-Term Traffic Flow Forecasting

Making one day as an example, the result is shown in Figure5.Figure 5
Effective analysis of real-time correction algorithm.Mean Absolute Percentage Error (MAPE) is used to evaluate the accuracy of the proposed model:(17)
MAPE
=
1
n
∑
t
=
1
n
V
^
s
t
-
V
s
t
V
^
s
t
×
100
%
.Evaluation result is shown in Table4. Comparing with these two kinds of model, real-time corrected traffic correlation model is better than the basic mathematical model. That is, the static part of real-time corrected traffic correlation model reflects the temporal-spatial-historical correlation by way of mathematical model; the dynamic part of real-time corrected traffic correlation model modifies the mathematical model based on multistate characteristics of traffic big data.Table 4
Evaluation result of proposed models.
Model
MAPE of flow (%)
MAPE of speed (%)
MAPE of occupancy (%)
Basic mathematical model
13.14%
9.18%
14.48%
Model improvement
10.93%
7.17%
10.54%


## 4.2.1. Basic Mathematical Model

Based on historical data, basic mathematical model, including traffic correlation model and simplified traffic correlation model, is built. The static part of real-time corrected traffic correlation model, which is “φ
∑
j
=
1
η
ω
j
V
j
t
-
d
+
ε,” is achieved, as shown in (16)
q
1
t
=
1.0629
×
0.2246
·
q
9
t
+
0.2151
·
q
1
t
-
5
q
1
t
=
1.0629
×
w
+
0.2145
·
q
1
H
t
-
5
+
0.1787
·
q
2
t
q
1
t
=
1.0629
×
w
+
0.1670
·
q
6
t
-
10
+
22.652
,
v
1
t
=
1.029
×
0.2677
·
v
9
t
+
0.2630
·
v
1
t
-
5
v
1
t
=
1.029
×
w
+
0.2461
·
v
1
H
t
+
0.2232
·
v
2
t
v
1
t
=
+
1.8173
,
ρ
1
t
=
1.0148
×
0.2791
·
ρ
9
t
+
0.2747
·
ρ
1
t
-
5
ρ
1
t
=
1.0148
×
w
+
0.2294
·
ρ
2
t
+
0.2168
·
ρ
1
H
t
ρ
1
t
=
+
0.4793
.Goodness of Fit is shown in Table1, in which the value of R
2 shows the high effectiveness of formulation (16).Table 1
Goodness of Fit.
Parameter of traffic flow
Flow
Speed
Occupancy
R
2
0.8885
0.9655
0.9664


## 4.2.2. Model Improvement

Real-time correction algorithm is presented to obtain the dynamic part of real-time corrected traffic correlation model, which is “Δ
t.” The steps are as follows.Step 1 (clustering center search).
Based on SAGA-FCM,LOS is divided into six types; Δ is divided into five types. Clustering Centers are shown in Tables 2 and 3.Table 2
Clustering centers of LOS.
i
Speed (km/h)
Occupancy (%)
A
68.1
2.1
B
58.3
3.1
C
51.2
11.4
D
46.0
14.9
E
38.0
20.5
F
23.4
35.1Table 3
Clustering centers ofΔ.
k
Flow (pcu/h)
Speed (km/h)
Occupancy (%)
V
S
−1043
−9.7998
−7.4513
S
−498
−2.5634
−2.2416
M
−2.9
0.2030
−0.2448
L
291.6
2.3692
1.6348
V
L
812.2
9.9043
7.0243Step 2 (calculation of membership degree).
Formulations are shown is Section3.2.Step 3 (training Fuzzy Neural Network).
Input layer has 11 nodes; the transmission function of the nerve cells in the hidden layer is transig; the output layer has 5 nodes and the transmission function of the nerve cells in the output layer is logsig, while the training function is traingdx. Historical data is used to train Fuzzy Neural Network.


## 4.2.3. Short-Term Traffic Flow Forecasting

Making one day as an example, the result is shown in Figure5.Figure 5
Effective analysis of real-time correction algorithm.Mean Absolute Percentage Error (MAPE) is used to evaluate the accuracy of the proposed model:(17)
MAPE
=
1
n
∑
t
=
1
n
V
^
s
t
-
V
s
t
V
^
s
t
×
100
%
.Evaluation result is shown in Table4. Comparing with these two kinds of model, real-time corrected traffic correlation model is better than the basic mathematical model. That is, the static part of real-time corrected traffic correlation model reflects the temporal-spatial-historical correlation by way of mathematical model; the dynamic part of real-time corrected traffic correlation model modifies the mathematical model based on multistate characteristics of traffic big data.Table 4
Evaluation result of proposed models.
Model
MAPE of flow (%)
MAPE of speed (%)
MAPE of occupancy (%)
Basic mathematical model
13.14%
9.18%
14.48%
Model improvement
10.93%
7.17%
10.54%


## 5. Conclusions

Traffic big data strongly shows temporal-spatial-historical correlation and multistate characteristic. Traffic correlation model is established based on temporal-spatial-historical correlation. Correction coefficients optimization algorithm is put forward to reduce parameters and ensure calculation accuracy. In order to improve the effectiveness of short-term traffic flow forecasting, real-time correction algorithm is presented based on multistate characteristic. Fuzzy Neural Network is used to overcome the problem of nonlinear mapping. Case study shows that real-time correction algorithm can improve the effectiveness of traffic correlation model.The core of this paper is to present a short-term traffic flow forecasting based on traffic big data analysis. The advantages of real-time corrected traffic correlation model for traffic flow forecasting are as follows.(1)
The temporal-spatial-historical correlation, which is considered in the static part of model, explains the physical meaning of traffic flow forecasting by mathematical model.
(2)
The multistate, which is considered in the dynamic part of model, explains the dynamic characteristic of traffic flow.
(3)
Real-time correction algorithm improves the accuracy of traffic flow forecasting. Case study shows the high efficiency and applicability of the proposed methods.Moreover, the proposed methods can be extended to 15- or 30-minute-ahead forecasting. The next steps of this work are to study traffic incident and its influence for short-term traffic flow forecasting. In addition, how to deal with the long period traffic flow forecasting like one hour or even longer can also be focused on.



---

*Source: 348036-2015-09-09.xml*