Probabilistic Forecasting of Traffic Flow Using Multikernel Based Extreme Learning Machine

Real-time and accurate prediction of traffic flow is the key to intelligent transportation systems (ITS). However, due to the nonstationarity of traffic flow data, traditional point forecasting can hardly be accurate, so probabilistic forecasting methods are essential for quantification of the potential risks and uncertainties for traffic management. A probabilistic forecasting model of traffic flow based on a multikernel extreme learning machine (MKELM) is proposed. Moreover, the optimal output weights of MKELM are obtained by utilizing Quantum-behaved particle swarm optimization (QPSO) algorithm. To verify its effectiveness, traffic flow probabilistic prediction using QPSO-MKELM was compared with other learning methods. Experimental results show that QPSO-MKELM is more effective for practical applications. And it will help traffic managers to make right decisions.


Introduction
Recently, the traffic flows maintain a steady growth in both urban and rural traffic leading to pollution, accidents, and congestion.To solve the problems, the intelligent transportation systems (ITS) are developed by many countries.The effectiveness of ITS is improved by using a lot of modern information technologies.According to the prediction period of time, traffic flow prediction can be divided into long-term, mid-term, and short-term prediction.Short-term traffic flow forecasting has become one of the main research areas of ITS.Prediction of real-time and accurate traffic flow becomes extremely important for effective traffic management systems, including traffic control, traffic induction, and vehicle routing.The problem induced by the randomness, nonlinearity, and complexity of traffic flow has compelled us to search for more reliable techniques to forecast traffic flow.
A lot of short-term traffic prediction algorithms are proposed in the literature [1][2][3][4][5][6][7].Conventionally, a majority of study focus on developing accurate point-prediction method structures and learning algorithms for short-term traffic flow prediction, but these methods cannot be used for the quantitative analysis of the uncertainty of the prediction.In fact, there are lots of traffic variables influencing the results of traffic flow prediction such as weather, date and time, types, and flow parameters.The aim of the traffic flow forecasting model is to utilize these variables to predict the traffic flow.
Because of the chaotic property of traffic flow, mistakes in prediction are simply inevitable.In other words, forecast with exact-point value becomes inadequate to describe the real world information.To deal with the problem, prediction intervals (PIs) are vital for quantifying the underlying risks and uncertainties.PIs are a list of ranges including upper and lower bounds where the targets will lie.On the basis of the PIs with corresponding confidence level, the quantitative uncertainties of traffic flow predictions can supply valuable information to traffic managers for good preparation for the most severe and excellent situations in advance.

Literature Review
The exact point-prediction methods only supply the point forecasting value, while PIs work as intervals that consist of upper and lower bounds as well as pointing out the probability of correct forecast.PIs not only indicate the prediction accuracy, but also provide a range that targets will locate 2 Scientific Programming [8,9].Some investigations in regression problems along with interval outputs have been conducted and they are able to be sorted into two classifications based on the amount of models they utilized.Studies in the first category use only one model, whereas they bring in other approaches to acquire interval outputs.These PIs construction approaches are usually set after a point-prediction model with particular preceding assumptions.Bayesian [10], mean-variance estimation [11], and bootstrap [12] are frequently used to obtain PIs.The main disadvantage of these existing methods is the high computational requirement.The other category focuses on learning a twin model; the PIs construction method named lower upper bound estimation (LUBE) approach is brought forward in [13].Nevertheless, traditional neural networks (NNs) applied in the LUBE approach have the issue of high computational cost and overtraining.
ELM [14] is brought forward to train single-hidden layer feedforward neural networks (SLFNs).Although ELM extensively enhances the training effectiveness, fluctuation caused by the random input and hidden layer weights influences the steadiness of ELM in the situation of identical training data as well as model parameter [15].The method that replaces ELM hidden layer with a kernel function makes ELM avoid choosing input and hidden layer weights randomly since the computation for hidden input is conducted by kernel function.The problem in ELM caused by random input and hidden layer parameters is solved by Kernel-ELM which gains higher stability by sacrificing the training rate [15].Moreover, one kernel function is usually used in the standard kernel learning algorithms [16].In our last paper [17], we proposed a single kernel extreme learning machine-(KELM-) based probabilistic forecasting method of traffic flow.On the basis of multikernel learning thought [18], the composites of two kernels may combine the good characteristics of them and have better performance than any other single kernel.The polynomial kernel function and the Gaussian kernel function are mixed as a kernel function which includes both kernels' advantages in this study.
A novel probabilistic traffic flow prediction approach is brought forward based on the multikernel ELM in this paper, which is applied to construct PIs for traffic flow.Then the output weights from MKELM models are optimized by Quantum-behaved particle swarm optimization (QPSO) [19].The proposed approach has been examined by the practical traffic flow data.Accurate prediction results have represented the good performance of the QPSO-MKELM approach for traffic flow.
The rest of this paper is organized as follows.In Section 3 the algorithms of framework of MKELM and QPSO are presented.Section 4 proposes a novel MKELM model to construct optimal PIs for traffic flow with QPSO.The experiments to verify the effectiveness of the model are carried out in Section 5. Finally, Section 6 draws the conclusion.[20], is a novel learning algorithm for SLFNs that randomly selects hidden nodes parameters as well as ascertaining the output layer parameters of SLFNs analytically.

Methodology
In a specified training dataset including  samples, (, ) is a training sample.The formula of the SLFN with  ℎ hidden nodes is able to be represented as The output equation of ELM is able to be represented as ( According to ELM theory [14], input as well as hidden layer parameters are able to be randomly distributed as long as the activation function becomes infinitely differentiable.As to fixed input weights   and the hidden layer biases   , to train an SLFN is simply equal to figuring out a least-squares solution β of the linear system  = .
The smallest norm least-squares solution of the above linear system is where  + represents the Moore-Penrose generalized inverse of matrix .
With a user-defined cost coefficient , Huang et al. [15] optimized the computation for the output weights .Among them, when the amount of the hidden nodes is larger compared to that of train data, the result of  is The hidden layer output of every sample ℎ(  ) is able to be considered as an nonlinear mapping of samples   .Then ] *  . ( When the hidden layer characteristic mapping ℎ() could not be unidentified, Huang et al. [15] suggested applying a kernel function.According to theory of kernel function, the kernel matrix Ω for ELM is able to be described as follows, where (  ,   ) is kernel function: So ( 5) can be deduced: Then the output function of Kernel-ELM can be represented as Any kernel function which meets Mercer's theorem [21] is able to be a kernel function of KELM.The kernels have their own advantages and shortcomings, which can be divided into two types, global and local kernels [16].In local kernels, only the data that are close to or in the proximity of each other have obvious effects on kernel values.The Gaussian kernel function is an example of a typical local kernel.In contrast, a global kernel allows data points that are far away from each other to have obvious effects on the kernel value as well.A typical example of a global kernel is the polynomial kernel.
In this work, we proposed a multikernel function which is taking the advantages of the polynomial kernel function and the Gaussian kernel function.The multikernel function is described as In multikernel function,  is an adjusting coefficient between   and  poly .When  approximates to 0, the value of  approaches  poly .In the same way, when  approximates to 1, the value of  approaches   .Moreover, the width  and the degree  are the kernel parameters of the Gaussian kernel and the polynomial kernel functions, separately.
For achieving better generalization performance, the four variables , , and  of MKELM and the punish parameter  need to be chosen appropriately.

QPSO.
As an advanced method, QPSO [18] is an intelligence searching method; in the meanwhile, the status of every particle is illustrated by using wave function when all the particles are transferred into the quantum space, rather than the position as well as velocity in traditional PSO [19].The particles of QPSO could appear much far away from current locations according to characteristics of the wave function, which enhances the possibility of escaping from the local optimal value.Considering that the QPSO algorithms include -sized particle population in the -dimensional searching space, the th particle   is renewed based on where  = 1, 2, . . ., ,  = 1, 2, . . ., ,  is a random number distributed equally on (0, 1),  is a positive number as well as ranges from  min to  max called the contraction expansion (CE) coefficient, which balances the local and the global optimum,   is known as the local attractor of every particle at -dimension on the basis of the trajectory analyses in [22],  best is the mean best position,   is the best previous position of particle ,   is the position of global best particle, and  1 and  2 are two different random vectors.

PIs Model Construction by QPSO-KELM
According to the theory of PIs, when targets lie in the PI nominal confidence (PINC), the prediction values need to be in the constructed PIs at the possibility of PINC equals 100(1 − )%.It can be expressed as where L  (  ) and Ũ  (  ) are the lower borders and upper borders of the prediction PIs at the nominal confidence level of the input   .
To construct PIs for the traffic flows, the KELM-based PIs establishing approach can be employed, which is illustrated in Figure 1.The KELM-based probabilistic prediction model targets are shown to directly create the upper borders and lower borders of two outputs.In Figure 1, () and () are the two output vectors of SLFN with respect to the input sample   ,   = [ 1 ,  2 , . . .,   ]  is the weight vector connecting the th hidden node and the input nodes, and   = [ 1 ,  2 , . . .,   ]  is the weight vector connecting the th hidden node and the output nodes.Besides, the determination for reliable PIs of the traffic flow via straight optimization for both sharpness and dependability will be demonstrated in the following part.

Construction of Optimal PIs for Traffic Flow.
To assess the properties of the PIs acquired, the PI coverage probability (PICP) as well as the PI normalized average width (PINAW), which represents the dependability and the sharpness of PIs, is utilized in [11].PICP is a crucial indicator regarding the dependability for the established PIs, which can be presented by where  test refers to the number for testing data.If the target   locates within the th lower border and the th upper border, then   = 1, or else   = 0.
In the procedure of interval prediction, the future targets are anticipated to locate within borders of constructed PIs at the PINC level.Nevertheless, it is able to be simply achieved by enlarging PIs from upper or lower bound.Actually, such PIs are meaningless for making a decision.Among the publications, PINAW is defined to indicate the PIs average width quantitatively, which is presented by where  refers to the targets' range and  is applied to standardize the average width of PIs in terms of percentage.
To guarantee the high properties of created traffic flow's PIs, the output layer parameters  of the MKELM model are improved to lead to higher PICP as well as lower PINAW value.The extensively utilized coverage width-based criterion (CWC) [11] in (11) can be considered for training the KELM model: In (11),  denotes the confidence probability 100(1 − )%. works to magnify any tiny discrepancy between PICP and , and  PICP is a function of PICP and defined as When PICP is more than the given confidence probability ,  PICP = 0 and CWC inclines to be identical to PINAW.
Otherwise,  PICP = 1 and the function is going to be calculated by CWC.However, it is difficult to decide the value of  when the CWC function is used.
In this paper, a width deviation rule (WDR) is brought forward to form a new objective function where the exponential punishment term in CWC will be substituted by the equation called WD as follows: The PIs width deviation for each sample can be expressed in (19). is set to 1/ 2 as a penalty coefficient in the present work.For PICP > , () will be equal to PINAW which is the same as CWC.Otherwise, the width deviation information of all samples is considered as a penalty to describe the sharpness more comprehensively.

QPSO Algorithm for PIs Optimization.
To obtain the effective interval forecasts of traffic flow, the optimal output weights  of MKELM need to be calculated by using QPSO algorithm.The properties of acquired PIs are mirrored by the fitness values of every particle, which is obtained by computing the aim function in (11).The main steps of PIs optimization by using QPSO are demonstrated as follows.
(2) Initialization.The interval prediction model based on MKELM with two outputs needs to be constructed.The upper bounds are set 30% higher than the targets.In the same way, the lower bounds are set 30% lower than the targets.Calculate the initial output matrix of hidden layer  int .The particles position [ 1 ,  2 , . . .,   ] is initialized by the initial output weights  int . (

3) Construction of PIs and Evaluation of Cost Function.
With the initialization parameters, the PIs model based on MKELM is constructed and the fitness and sharpness are calculated for each particle using ( 15) and ( 16).
(4) Renew the Position of Each Particle.All particles' locations are updated in the light of ( 11)-( 13).
(5) Renew   and   .With the renewed , establish a novel model and assess the fitness values of fresh particles.If present fitness value is better compared to that of   , afterwards   will be renewed.Moreover, if the fitness is better compared to   , then   will be renewed.
(6) Loop.On condition that the largest amount of repetition has not been achieved, go back to step (4).Or else, the process ends and the optimal output weights are acquired to establish MKELM-based model.To get an effective database, the original traffic data need to be processed in advance.For example, when the speed is more than 100 km/h, the sample is invalid and needs to be eliminated.When flow, speed, and occupancy are zeros at the same time, the sample also will be removed.Each simulation runs 4-fold cross validation apart from particular account.The training samples are randomly divided into 4 same scale subsamples, where a single one is kept to be the validation data to test the model, and the rest are utilized to be training data every time.

Application Studies
The detailed sets on May 15 are illustrated in Table 1.
The 258 attributes include the most important seven kinds of features such as date of week, time of day, traffic density, speed, flow, weather condition, and road logical region in Table 4. Generally, attributes either point out measurements for certain continuous scale such as traffic flow for the last time under this circumstance or demonstrate information about certain definite or separate features such as date type within the above characteristics.The characteristics utilized in this paper blend definite characteristics with real-valued characteristics and are supposed to be converted to all definite or real-valued characteristics.In our conditions, all definite characteristics possess only single value, which is set to be either 0 or 1.For example, characteristic date type from Monday to Friday is going to be converted to 7 new characteristics, each of which refers to either 0 or 1 for pointing out whether it occurs.Therefore, all the characteristics are real-valued and are able to be equally evaluated.For traffic flow prediction, feature extraction is regarded as a significant part and will impact the forecasting outcomes directly.It is essential to decide appropriate input variables for precise predictions.Apparently, some factors affect the traffic condition of roads, for example, whether there are traffic lights or not, where the road is, and if a school or a shopping center is located nearby.However, it is difficult to obtain these characteristics automatically.Consequently, these characteristics will not be taken into consideration.We chose the following characteristics.
Determining how to forecast the traffic flow accurately within a 5 min time interval is critical for ITS because the cycle of transportation induction is usually 5 minutes.
(iv) Traffic flow of the last time interval: continuous, normalization data varying from 0 to 1.
where Volume max, represents the maximum flow volume of the th day and Volume min, is the minimum flow volume of the th day.In addition,  donates the amount of the sample days.
Physical region, which partly represents the circumstance of roads, can play an important role in road congestion.For example, roads nearby the school could usually stay in a good situation apart from the time when students head for school on mornings and when the students go back home on afternoons.At the same time, traffic situation possesses strong space-time periodicity and area characteristic is supposed to be taken into consideration.Regrettably, the information of physical region cannot be automatically obtained by ITS, and it is easy to transform, particularly in the situation where a new road is in construction when ITS runs.In this paper, a logical region is utilized rather than physical region to demonstrate roads' circumstance characteristics because of the inaccuracy of the area labeled by individuals.
In this paper, the -means algorithm is used to cluster the logical regions into  clusters.The value of  is set to be 50 in this paper.
(vi) Road type: categorical, 1, 2, 3, 3 in total.This eigenvalue is defined with the numerical value.From branch way to major road the numerical values are 1, 2, and 3 in order.The higher the road grade, the higher the road's standard speed and maximum speed.
(vii) Number of lanes.
The number of lanes represents the road capacity.Under the circumstances of same traffic flow, the more lanes there are, the higher the speed is.
(viii) The average speed of all cars on the road section.Firstly, we calculate the average speed of every car in the interval of Δ on the road section.Then the average speed Therefore, the floating car 's average speed can be represented as   =   /(  −  0 ).
If the number of floating cars on the road section at the moment is , the average speed of all cars on the road section can be written as Ũ = (∑  1   )/.(ix) The speed distribution of all cars on the road section.The car speed is distinguished into different levels and the histogram is used to represent the distribution of speed data of all floating cars on the read section.The standard of division is based on the speed distribution of floating cars in the city (Figure 3).Car speed data is mainly no more than 75 km/h.Therefore, the car speed is distributed into 5 grades as shown in Table 3.
(x) The average stopping time of all cars on the road section.
When a car's speed is below 5 km/h, it is identified as a stopped car.The floating car 's stopping time can be represented as   = ∑ −1 =0 ( +1 −   )(  < 5).Thus the average stopping time of all cars on the road section can be written as T = (∑  1   )/.

Input Dimension Reduction.
Determining the most relevant features is an issue of research in itself.Taking all the factors into consideration will lead to computation complexity, dimensionality course, and overtraining.There are three sorts of methodologies for decreasing the input dimensionality including NNs sensitivity analysis, transformation, and correlation among features [23].The correlation approaches lead to their high computational complexity.
The problem with the transformation approaches is that the features of the original input no longer exist.To address the problems, we selected the sensitivity analysis method as the feature dimension reduction tool in this paper.Sensitivity analysis concentrates on how a network's output is influenced by its input perturbations.Irrelevant inputs are found and eliminated to reduce data collection cost and improve the network's performance.The underlying relationships between the inputs and outputs are found after sensitivity analysis.Sensitivity of an output   corresponding to its input   is especially Inputs which affect outputs most significantly are determined by examination of the mean square average (MSA) sensitivity matrix  ,avg .A minor value of  ,avg compared with others suggests that, as to the special th output for the network, the th input cannot remarkably donate per average to the output , and consequently could be probably negligible.The one input varies up or below the mean while all other inputs were kept fixed at their respective means.This process was the repeated for each input variable.The abbreviations of the chosen characteristics are shown in Table 4.And outcomes of sensitivity analysis are presented in Figure 4 in which seven most significant input parameters are date of week (DW), time of day (TD), traffic density (DEN), speed (SP), flow (FL), weather condition (WC), and road logical region (RLR).
When either one or more inputs possess comparably minor sensitivity compared to others, the NNs' dimension could be decreased by eliminating them; meanwhile a smaller-scale NNs are able to be retrained under most circumstances successfully.

Determination of Model's Parameters and Its Application to Traffic Congestion Prediction
(a) Optimization of Multikernel Model's Parameters.In (7), it can be observed that the prediction precision of MKELM is mainly affected by the four variables , , and  of MKELM and the punish parameter .The four variables need to be selected with optimizing method.Then the problem of how to choose proper parameters of MKELM can be converted into optimizing four variables:  = [, , , ].The vector  which has the best fitness is chosen in MKELM.In our study, QPSO is used to obtain optimized parameters.The function for fitness calculation is where (  ) is the predicted value,   is the target, and  is the amount of training data.The main procedures are described as follows: (1) Preprocesses the original samples and divide them into training and testing samples.(3) Calculate every particle's fitness of every particle and renew.
(4) Compare and update the top fitness value.
(5) If the ending conditions are met, the optimal parameters of MKELM are obtained.

(b) Traffic Congestion Prediction Using Optimized MKELM
Model.Driving speed refers to a significant indicator for judging the road situation.It is simple to obtain average speed of vehicles via gathering 5-minute speed.We primarily utilize speed for determining the Traffic Congestion Indicator (TCI).Equation ( 28) is the function applied to assess transport situation while its plot is illustrated in Figure 5.
) , where  refers to the speed of road segment (unit: km/h), whereas () refers to TCI.It can be seen that road situation becomes worst if the average speed turns to be zero, and consequently TCI is equal to 100.Otherwise, traffic situation becomes better as the speed increases.Thus, when speed inclines to be infinite, the value of TCI inclines to be 0, where  = 0 is the horizontal asymptote of ().
The output of assessment refers to a continuous value which changes between 0 and 100.The feasible result indicates that the speed-TCI model turns out to be rational and be able to represent traffic status appropriately.The assessment outcomes satisfy the true condition as illustrated in Figure 6.In the map, green represents smooth traffic, yellow shows average condition, and red means the road is congested.Seen from the image taken by surveillance cameras, the traffic evaluation accurately reflects the road traffic congestion at that time.The system creates assessment output every 5 minutes so as to portray the present urban traffic situation for every road segment.These data gradually become historical.The historical assessment data are able to be applied to train forecasting model, and it is clear that a large number of training data are generated every day.In practice, almost 1,500,000 valuable data are created every day.Large data call for a proper model with fast speed.It is known that MKELM  algorithm inclines to offer good generalization feature at a fast learning rate.
With the optimal parameters of MKELM, we just draw target values of 50 testing data selected randomly, single Kernel-ELM predictions, and MKELM, respectively.In standard single kernel model, the radial basis function (RBF) is used of which the parameter  and the cost coefficient  were selected in {20, 21, . . ., 210}. Figure 7 shows that the multikernel model of KELM has more accurate prediction results on traffic training data than standard one-kernel ELM.The difference between the true values and the ones predicted by multikernel model are very small for most of the cases.

Study Results and Discussions.
In this experiment, the kernel function of KELM is multikernel described in (10).For the parameters of QPSO, the population scale as well as the amount of repetitions is 200; the values of  min and  max are set to be 0.5 and 0.9.
The main target of the proposed QPSO-KELM approach is to originate dependable PIs.Moreover, traffic management needs valuable information along with higher confidence levels.Consequently, it is supposed to be more virtually significant to acquire high-confidence-level PIs for satisfying the requirements of traffic management.Diverse levels of PINC 100(1 − )% from 80% to 95% are taken into account in this research.
(1) Comparison of MKELM-Based Model Using CWC and WDR.To construct the PIs for traffic flow, whether the rule CWC or the rule WDR should be used in the MKELM-QPSO model needs to be decided firstly.
For the traffic flow, the dependability indicator PICP as well as the sharpness indicator PINAW of the constructed PIs at three different confidence levels 80%, 90%, and 95% is summarized in Table 5.Also the results using CWC and WDR are, respectively, summarized in Table 5.It is shown that, at every confidence level, the PICP values of WDR are closer to the confidence levels than the values of CWC.At the same time, CWC has a much higher sharpness than WDR which provides much more accurate PIs.The comparisons of CWC and WDR in Table 5 show that WDR generates the higher dependability and sharpness of the PIs.In the rest of the paper, the criterion WDR is applied to get the optimal PIs because of its better performance shown in Table 5.
To achieve a better visual effect, we draw the real values and forecasting PIs, respectively, of 100 testing samples which are selected randomly.The constructed optimal PIs at different confidence levels 80%, 90%, and 95% are shown in Figure 8.Most of the targets have been limited within the upper and the lower bound.For different confidence levels, most of the measured traffic flows are limited in the PIs established by the proposed method, which means that the remarkable performance will satisfy the requirement of traffic management.
Also the nonstationary and nonlinearity characteristics of traffic flow are clearly demonstrated in these graphs.
(2) Comparison of Different Methods for Constructing the PIs.For the sake of further demonstrating the effectiveness of the QPSO-MKELM method, different intelligence search and learning algorithms are compared using the same dataset.
As to the PIs dependability test, corresponding PICPs and PINAWs are demonstrated in Table 6.The corresponding PICPs and PINAWs outcomes of diverse methods are described in Table 6.At different confidence levels, the dependability indicator PICP and the sharpness indicator PINAW are summarized in Table 6, which are acquired by QPSO-MKELM, PSO-MKELM, and GA-MKELM.From the outcomes of PICP as well as PINAW indices, it is shown that QPSO-KELM methods offers far more precise PIs combined dependability (indicated by higher values of PICP) as well as sharpness (reflected by lower values of PINAW) on the three various confidence levels.
The SVM, ELM, standard KELM, and MKELM are compared on the basis of the training time, PICP, and PINAW at the same confidence level 90% in Table 7.Although the MKELM cannot compete with the ELM in terms of the training time, the MKELM provides much more accurate PIs than the ELM.We can also see that the SVM almost provides the same PICP and PINAW as the MKELM, but the training time for the SVM is more than 20000% of that for the KELM.Overall, the MKELM outperforms the SVM, the ELM, and the standard KELM.Thus, it is rational to draw a conclusion that the proposed method is an efficient probabilistic forecasting approach for traffic flow.

Conclusions
In order to develop an effective probabilistic forecasting method for traffic flow, a novel method on the basis of QPSO-MKELM has been described to establish the reliable PIs.MKELM has been developed to establish the reliable PIs, and the parameters of MKELM are optimized by using QPSO.The seven features including date of week, time of day, traffic density, speed, flow, weather condition, and road logical region are selected as inputs of KELM by sensitivity analysis.The experimental results have shown that QPSO-MKELM is an effective method to establish the optimal PIs.Moreover, the proposed method can offer far more precise PIs that combined higher dependability and sharpness at different confidence levels than other methods.Additionally, successful utilization in practical traffic flow prediction shows that QPSO-MKELM is an effective probabilistic forecasting method.In the current paper, only limited traffic flows are taken into consideration; other parameters such as accidents, traffic jams, or seasonal variation have not been discussed.In future study, more possible conditions will be included in a longer period of time.

Figure 2 :
Figure 2: Roads are split into many small pieces.Road pieces are the basic unit to evaluate or predict.

Figure 2 (
a) illustrates the road segments in Nanning city, and Figure 2(b) refers to the enlarged vision of road segments.

Figure 3 :
Figure 3: The speed distribution of floating cars in the city on 23rd January 2015.

Figure 6 :
Figure 6: The comparison of camera results and evaluation results shows that the evaluation model is reasonable.

Figure 7 :
Figure 7: Part of regression prediction shows that mixed kernel model has a great fitting result.

Table 1 :
Information of data on May 15.The traffic features such as flow parameters, date types, and environmental situations will be described in detail in the following part.Also the most relevant features for high prediction precision are selected by sensitivity analysis.
5.2.1.Feature Description.Actually, the roads are divided into minor road segments as fundamental components for assessment and forecast.The road in Nanning city consists of 18, 041 minor segments.

Table 2 :
Quantificational disposal of weather factors.

Table 3 :
Five levels of the speed of floating cars.In the interval of Δ, the velocity measurement sites of the floating car  on the terminal road section are distributed as shown in the figure below.Sequence { 0 ,  1 , . . .  } and sequence { 0 ,  1 , . . .  } are floating car's time sequence and speed sequence on the road section.The floating car's driving distance   is defined by

Table 4 :
List of input variables.