Stochastic Parameter Identification Method for Driving Trajectory Simulation Processes Based on Mobile Edge Computing and Self-Organizing Feature Mapping

Shenyang Institute of Automation Chinese (Guangzhou) Academy of Sciences, Guangzhou 511458, China Shenyang Institute of Automation Chinese Academy of Sciences, Shenyang 110016, China School of Electronics and Communication Engineering, Sun Yat-Sen University, Guangzhou 510006, China Huangpu Customs District Technical Center, Guangzhou 510730, China School of Communications and Information Engineering, Xi’an University of Posts and Telecommunications, Xi’an 710061, China Guangdong Zhongke Zhenheng Information Technology Co., Ltd, Foshan 528225, China South China Agricultural University, Guangzhou 510642, China Yaz Technology Co., Ltd., Guangzhou 510630, China


Introduction
With the development of automated driving technologies, especially the onset of 5G technology, the demand for traffic data collection in the field of intelligent networking has increased significantly. Conventional driving environment detection technologies use road monitoring equipment with different traffic sensing technologies including video, GPS, geomagnetism, and radar. ese devices are installed on roads by local public security agencies to ensure a safe and smooth traffic flow, thus strengthening the law enforcement against traffic violations and the intelligent control of citylevel transportation. However, these driving environment detection technologies cannot satisfy automated driving communication requirements of short time delay, high reliability, wide coverage, and vehicle-to-vehicle communication. Under this backdrop, the application of electronic registration identification of motor vehicles has been proposed and promoted on a national level through such measures as establishing national standards. Electronic registration identification of motor vehicles have been promoted in a few cities and applied to several areas such as smart parking [1,2], intelligent signal control [3,4], intelligent connected vehicles [5], special vehicle traffic management and special vehicle identification monitoring [6], and traffic operation supervision and environmental protection restrictions [7]. However, as a means of law enforcement, the application of vehicle interactions is still lacking due to different technical problems; thus, it is difficult to achieve the required accuracy for vehicle driving trajectory identification. Hence, researchers have proposed different identification methods to solve this problem; for example, grid-based methods have been successfully used in multiobjective optimization algorithms [8][9][10], such as the positioning of driverless cars [11][12][13]. Grewe et al. [14] explored a series of MEC-enabled, high-quality, and reliable vehicle-borne services (such as electronic horizon, which assist vehicle movement), summarized the challenging problems encountered in the application of MEC technology to the network of vehicles, and proposed some potential solutions. Truong et al. [15] advocated the use of a combination of software-defined networking and MEC technology to address the problems of existing vehicle-borne self-organizing networks, such as inadequate coverage of communication range, unscalable network communication capacity, and a management that lacks intelligence and flexibility.
ey designed an edge computing network framework (dubbed FSDN) for the Internet of vehicles and MEC server that acts as a distributed zone controller to provide different local services, such as streaming of media content distribution and vehicle lane-change prompts, aimed at providing excellent decision-making capabilities. Real-time and accurate traffic information is the prerequisite for an improved urban traffic efficiency; Dikaiakos et al. [16] proposed the vehicle information transmission protocol (VITP) and implemented dynamic route planning using the forward traffic status and roadside points of interest obtained from query requests to downstream vehicles through the geographic routing protocol. Dombush and Joshi [17] studied the automatic discovery mechanism of traffic congestion and the distributed clustering of abnormal traffic flow and achieved the aggregation of perceived data through interactive sensing between vehicles. Zhang and Zhao [18] constructed a mathematical model to automatically collect, aggregate, process, and transmit traffic information and dynamically update the traffic coverage of the entire road network. Gramaglia et al. [19] used external data sources such as vehicle beacon messages and weather conditions to detect the degree of traffic congestion through complex event processing. Terroso-Saenz et al. [20] and Li et al. [21] detected traffic congestion from intervehicle collaboration where each vehicle used the collected beacon messages to estimate the surrounding traffic congestion using fuzzy logic and then corrected the individual estimation error through consultation. Bauza et al. [22] conducted traffic condition sensing through VANET collaboration and made short-term traffic condition forecasts using a linear least squares method. However, these studies have only solved the problems of cooperative traffic condition sensing and neither consider implementation issues in automated driving scenario nor analyze important driving trajectories in vehicle-to-vehicle communication. In the present paper, based on electronic registration identification of motor vehicles and readers, a closed loop from perception to service was formed using mobile edge computing and self-organizing feature mapping to identify stochastic parameters for driving trajectory simulation processes and solve different application problem of lane-change scenarios in automated driving.

Driving Behavior, Speed, and Trajectory
Calculation Based on Electronic Registration Identification of Motor Vehicles e vehicle data collection method integrated with electronic registration identification of motor vehicles reads the vehicle information of the vehicle from the on-board unit with the electronic registration identification of motor vehicles reader installed on the gantry frame. e collected information is then uploaded to the comprehensive sensing base station on the test section of the road for edge computing. e overall deployment of the comprehensive sensing base station can complete data collection, verification, transmission, and processing in actual applications.
e overall layout and actual installation are shown in Figure 1: e main method of acquiring the driving speed of a vehicle with an electronic registration identification of motor vehicles reader uses the UHF radio frequency identification technology. e directional horizontally-polarized UHF antenna of the registration reader interacts with the unit installed on the front windshield of the vehicle and calculates the speed of the vehicle by measuring the time difference for the vehicle to pass through a fixed distance in the identification zone. Specifically, the positions of the vehicle at different times are calculated from RSSI values returned by the electronic registration identification of motor vehicles reader at different times, and the travel speed can then be calculated from the time differences. By analyzing RSSI values returned by the vehicle after entering the identification cross-section, it can be determined whether the location of the vehicle is at the boundary between the direct illumination zone and the blind zone or at the farthest point in the reflection zone; subsequently, the vehicle speed is calculated by

Mobile Edge Computing and Self-Organizing
Feature Mapping

Mobile Edge
Computing. e fundamental concept of mobile edge computing is to move the cloud computing platform to the edge of the mobile access network in an attempt to deeply integrate a cellular network of conventional telecommunication with Internet services and also to reduce the end-to-end delay of mobile service delivery. It changes the state of a separate network in conventional wireless communication systems by exploiting the inherent capabilities of wireless networks and improves the user experience. By adding mobile edge computing platform network elements to a conventional wireless network, a business platform including contents, services, and applications, can be lowered to the edge of a mobile network to provide computing and data storage services to users. e basic characteristics of mobile edge computing include service localization, short-distance, low time-delay service delivery, user location awareness, and other network service capabilities. It can, therefore, bring new changes to the mode of operation for telecom operators and establish a new ecosystem of industrial chains and networks. e edge computing architecture relies on the communication infrastructure and services provided by Edge Cloud Collaboration and LET/5G. e edge side mainly includes vehicle-borne edge computing units and RSU or MEC servers. Vehicle units depend on the cooperation of RSU or MEC servers. On the vehicle side, the control systems of intelligent networked vehicles are becoming increasingly complex with the use of AI and V2X applications. Load integration and simplified control systems are used to integrate different systems into HMI, and complex vehicleborne information services including ADAD, IVI, digital instruments, and HUD are made to run on the same hardware platform through simulation technology. Moreover, load integration based on virtualization and the hardware abstraction layer HAL are more amenable to flexible business orchestration, deep learning model update, and software and firmware upgrade for a cloud-based vehicle driving system.
Due to several hardware constraints, the generation of a large amount of V2X data, and the imbalance of task processing, real-time scheduling is necessary for computing tasks. In the present work, a hybrid critical computing task scheduling technique was employed. e proposed method is linearly complex but is low cost and highly schedulable and can lower the time consumption of a system by reducing the number of repeated calculations in the passive delay priority updating. By establishing a more accurate upper bound of the peak period based on the quantitative relationship between the task number and the peak period, the spatial overhead of the system can be reduced and the scheduling performance can be improved.

Stochastic Parameter Identification Method Based on Self-
Organizing Feature Mapping. According to the current practice of data optimization through neural networks, a self-organizing feature map (SOFM) was used in the current work to perform data optimization. Self-organizing feature map (SOFM) is a neural network with a clustering function [23]. Self-organizing mapping is a multidimensional scaling method, which can map any dimensional data in the input space to an output space of lower dimension. A self-organizing map (SOM) is an array of neurons at a regular location. Neurons can be placed on grid nodes in one, two, or higher dimensions, however usually in a one-dimensional or two-dimensional grid. e operation of a SOFM network can be divided into two stages: training and work [24,25]. In the training stage, samples from the training set are randomly fed into the network. For a specific input mode, the output layer generally has a neuron that produces the largest response to win. At the beginning of the training phase, it is uncertain which neuron in the output layer will respond most to the input mode. When the input mode changes to a different category, the winning neuron in the two-dimensional plane also changes. Neurons surrounding the winner neuron have a greater response due to the lateral mutual excitation. erefore, weight vectors connected to the winner neuron and its surrounding neurons adjust in different degrees to the direction of the input vector. e degree of adjustment gradually diminishes as the distance from the winner neuron to its surrounding neurons increases. e network adjusts weights in a selforganizing manner by a large number of training samples. Finally, each neuron in the output layer becomes a sensitive neuron for a specific mode, and the corresponding inner weight vector becomes the center vector of each input mode. When the features of two input modes approach each other, neurons representing these two mode categories also approach each other, resulting in the formation of an ordered feature map reflecting the distribution of sample modes in the output layer.
After the completion of SOFM training, the specific relationship between each neuron in the output layer and each input mode category is completely determined, to use  Complexity the network as a mode classifier. At each input mode, the specific neuron in the output layer of the network generates the maximum response and automatically classifies the input. It should be noted that when the input mode of the network does not belong to any of the modes of the training set, SOFM can only assign it as belonging to the closest mode category. e classification process is performed by the following steps [24,25].
In the first step, the winner neuron is found, the dot product of the input mode and the weight vector are calculated, and the winner neuron with the largest dot product is selected.
In the second step, the winning neighborhood is defined and the weight adjustment domain at time t centered around the winner neuron is determine. Generally, the initial neighborhood N is larger and gradually shrinks with the increase of the training time during the training process.
In the third step, the weights of all neurons in the winning neighborhood are adjusted by the following formula: where η(t, N) is a function of the training time t and the topographical distance N between neuron j and the winner neuron j * in the neighborhood.
Many functions can meet the above rules, for example, the following function can be constructed: (2) e fourth step is the termination of the inspection. Unlike a BP network, the concept of output errors does not exist during the training of a SOFM network. As it is an unsupervised training, the training ends when the rate of learning η(t)decays to zero or a certain preset positive small number. If the condition is not met, the process returns to step 1.
Kalman filter, based on the model of a state space, estimates the recursive relationship of an algorithm for a given process of filtering the state vector and constructs a measurement function (structural parameter or state) as an independent variable for an unknown quantity. It therefore derives an estimate of the unknown quantity from the measurement data. Kalman filter is widely applied to system identification. e conventional linear Kalman filter generally suffers from several shortcomings, such as low accuracy, poor stability, and slow response to target maneuver. In the application, the nonlinear motion equation of motion of the estimated dynamic system can be linearized and applied to the convergence of Kalman filter estimation. Furthermore, considering the nonlinear characteristics of the dynamic system near the reference trajectory and the comparison between strategically estimated algorithms, the real system can be described by a linearized equation.
Kalman filter can be divided into two parts: state prediction equation (state prediction) and state correction equation (observation update). e state prediction equation is responsible to predict current state variables and estimate the error covariance in a timely manner to construct a priori estimate for the next time state. e state correction equation is responsible for giving feedback and combines the a priori estimate and the new observed variable to construct an improved a posteriori estimate.
Let the following equations describe an n-dimensional linear dynamic system and an m-dimensional linear observation system [26][27][28]: where x i+1 is ann-dimensional vector that represents the state of the system at the time instance i − 1 andA is a n × n ordered matrix, which becomes the state transition matrix of the system, and reflects the state transition of the system from the ith sampling time to the i + 1st sampling time. w i is an n-dimensional vector and represents the random interference acting on the system at the time instance i + 1, and it is assumed that w i is a Gaussian white noise sequence with a known zero mean and a covariance matrix Q i , z i is an m-dimensional observation vector, and H is a m × n ordered observation matrix and represents the transition from the state x i+1 to the observation vector z i . For the m-dimensional observation noise v i , it is also assumed that v i is a Gaussian white noise sequence with a known zero mean and a covariance matrix R i . erefore, the following recursive formula of the filter can be obtained: where Q i is an n × n-ordered covariance matrix of the model noise w i , R i is the m × m-ordered covariance matrix of the observed noise v i , K i is the n × m-ordered gain matrix, x i+1 is an n-dimensional vector and represents the estimated value after filtering at the time instancei + 1, and P i is the estimated error covariance matrix of order n × n.
According to the calculation method described above, one can start with x 1 and a given P 0 and then recursively calculate the estimate of each time state by the known matrices Q i , R i , H, and A and the observation value z i at the time i. If the linear system is stationary, then A and H are both constant matrices. If the model noise w i and the observation noise v i are stationary random sequences, then Q i and R i are constant matrices. Under such conditions, the constant gain discrete Kalman filter becomes asymptotically stable.

Tests and Results
In the present work, a single-vehicle changing lane was used as the main application scenario to establish a motion model and perform system identification. e schematic diagram of a lane-change process is shown in Figure 2.
Let t 0 be the running speed of the vehicle at time v 0 ; then, its lateral displacement and longitudinal position are, respectively, L h � v 0 t 0 cosθ and Lz � v 0 t 0 sinθ. Starting from the initial moment of the trajectory, a search could find two points with the smallest time intervals on both sides of the lane line. Now, setting them to t c1 and t c2 , respectively, the time for lane change can be obtained by the following formula:  e results are compared with the actual data, as shown in Figure 5.
In comparison to the lane-change test results with no vehicles in the target lane, the longitudinal distance for completing the lane-change operation slightly increased slightlyin this scenario. e average longitudinal distance for completing a left lane change increased by 1.42 meters, whereas the average lateral distance increased significantly (nearly 93.89 meters). It happened because the obtained results were affected by an increased lane changing distance, the speed of vehicles in the target lane, and a longer distance to ensure safety.
Under standard operating conditions, the standard deviations between the actual value and the simulated trajectory for the three left lane-change tests were 1.12, 1.21, and 1.19, respectively, thus demonstrating a normal distribution.
e lateral distances in right lane-change operations were greater than those in left lane-change operations, and it can be ascribed to the seat position of the driver. e standard deviations, for the three right lane-change tests were 1.36, 1.53, and 1.47, respectively (slightly larger than those of left lane-change operations), thereby still conforming to a normal distribution.

Conclusion
e conventional automobile motion trajectory is mainly constrained by the kinematic and dynamic characteristics and real-time issues of a vehicle. In recent years, the research on smart car movement trajectories has paid more attention to the application of V2X communication among conventional automobiles. Transportation vehicles generally face a highly complex and stochastic driving environment. In order to strengthen the driving safety management, the present paper considered the single-vehicle lane-change process as a research topic. Based on electronic registration identification of motor vehicles and the associated reader/writer, mobile edge computing, and self-organizing feature mapping algorithms, a stochastic parameter recognition method was proposed for the driving trajectory simulation process, and its feasibility was verified through testing. A vehicle trajectory data acquisition method was first proposed based on electronic registration identification of motor vehicles, and the detailed description from the wiring and installation of the readers were then provided to the acquisition, calibration, and application of data. e lane-change process was mathematically modeled by mobile edge computing and self-organizing feature mapping, and the obtained simulation results were compared with actual test data. It was found that the standard deviation between the actual value of the lane-change operation and the simulated trajectory conformed to a normal distribution. erefore, the proposed method can effectively improve the accuracy of intelligent vehicle trajectory planning and the driving safety during lane-change operation.

Data Availability
e original data used to support the findings of this study are restricted by the relevant law-enforcement departments in order to protect vehicle information privacy and law enforcement basis. Data are available from relevant law-enforcement departments for researchers who meet the criteria for access to confidential data.

Conflicts of Interest
e authors declare that they have no conflicts of interest.

Authors' Contributions
Jingfeng Yang conceptualized the study, developed methodology, carried out formal analysis, and wrote the original draft. Zhiyong Luo developed methodology and supervised the study. Nanfeng Zhang validated the study and carried out data curation. Jinchao Xiao carried out data curation, formal analysis, and visualization. Honggang Wang validated the study and contributed to resources. Shengpei Zhou wrote, reviewed, and edited the article. Xiaosong Liu contributed to software and resources. Ming Li investigated the study, contributec to funding acquisition, and administered project.