# Identification of Aircraft Wake Vortex Based on SVM

Authors: Weijun Pan; Zhengyuan Wu; Xiaolei Zhang
Journal: Mathematical Problems in Engineering
 (2020)

---

## Abstract

The aircraft wake vortex has important influence on the operation of the airspace utilization ratio. Particularly, the identification of aircraft wake vortex using the pulsed Doppler lidar characteristics provides a new knowledge of wake turbulence separation standards. This paper develops an efficient pattern recognition-based method for identifying the aircraft wake vortex measured with the pulsed Doppler lidar. The proposed method is outlined in two stages. (i) First, a classification model based on support vector machine (SVM) is introduced to extract the radial velocity features in the wind fields by combining the environmental parameters. (ii) Then, grid search and cross-validation based on soft margin SVM with kernel tricks are employed to identify the aircraft wake vortex, using the test dataset. The dataset includes wake vortices of various aircrafts collected at the Chengdu Shuangliu International Airport from Aug 16, 2018, to Oct 10, 2018. The experimental results on dataset show that the proposed method can identify the aircraft wake vortex with only a small loss, which ensures the satisfactory robustness in detection performance.


---

## Body

## 1. Introduction

Aircraft wake turbulence is a function of an aircraft producing lift, which is invisible and harmful to the trailing aircraft, especially, hazardous in the takeoff or landing phases of the flight [1]. The study on the aircraft wake vortex can help in regulating the wake turbulence separation standards to ensure flight safety and avoid wake encounters. Thus, the conservative wake turbulence separation is the main factor that influences the increase of airspace capacity in the limited capacity of airspace.Computational fluid dynamics (CFD) and field measurements are two important methods to study the behavior of the aircraft wake vortex. The CFD supports the consistent investigation of wake vortex behavior under various environmental conditions but has not yet a significant impact on aircraft spacing for airport capacity enhancement [2]. The high accuracy and resolution of the lidar make it the most effective device at field measurement. Particularly, the new wake turbulence separation standards mainly rely on lidar measurements and technologies, which have emerged as an efficient tool for predicting and characterizing the aircraft wake vortex with different environmental parameters in different flight phases [3–9].However, the field measurement dataset collected by the lidar cannot be used directly because it contains nonwake wind field data. In addition, the dataset does not match with the correct flight because the relative position of the aircraft and the lidar is difficult to measure. On the contrary, the volume of data collected by the lidar is very large and it takes a long time to be processed manually even with the automatic dependent surveillance-broadcast (ADS-B) or quick access recorder (QAR) data. Particularly, the lower quality data caused by the sensor scanning and environmental parameters make the traditional manual methods incapable of working satisfactorily for the identification of large number of wake vortex data with high accuracy. Therefore, a fast and efficient method is required to automatically annotate the lidar data.Machine learning maybe a powerful method to classify the lidar data into nonwake and wake data by generating images from lidar data, but its performance relies on the visualization of lidar data and the configurations of color maps for different velocity. Unsuitable color maps may lead to low performance and loss of information in the process of data to images.In this paper, a classification model based on support vector machine (SVM) is proposed for automatically identifying the lidar data. The remainder of this paper is organized as follows. Section2 details our dataset, which was collected at the Chengdu Shuangliu International Airport (CSIA) from Aug 16, 2018, to Oct 10, 2018. The proposed method is described in Section 3. The results are presented in Section 4 followed by conclusion in Section 5.


## 2. Data Sources

The aircraft wake vortex behaves like large rolling air masses behind the flight path of the aircraft, and the main motion direction is perpendicular to the flight path of the aircraft in the plane. Thus, the lidar observational transverse mode is the range height indicator (RHI), as shown in Figure1.Figure 1
The observation mode of lidar, RHI.The mode of the lidar used for the experiment is Wind3D 6000, which is manufactured by Leica-Lidar Transient Technology Ltd\enleadertwodots. In addition, this lidar is based on all-fiber laser technology and the heterodyne detection technology, as shown in Figure2(a). The volume of the lidar is small, and its spatial and temporal resolutions are high enough for the detection of the aircraft wake vortex. The specifications of Wind3d 6000 are listed in Table 1.Figure 2
(a) Wind3D 6000. (b) Sketch map of field measurement.
(a)
(b)Table 1
The specifications of Wind3D 6000.
Qualification
Specifications
Wavelength (μm)
1.55
Pulse repetition (kHz)
10
Sampling rate (GHz)
1
Data update rate (Hz)
4
Pulse width (ns)
100–400
Pulse energy (μJ)
150
Measurement range (m)
45–6000
Range resolution (m)
15–50The field measurement is achieved at the Chengdu Shuangliu International Airport from Aug 16, 2018, to Oct 10, 2018. The position of the lidar is at pointA, and the threshold of runway 02R is at point B, as shown in Figure 2(b). The elongation line of the runway centerline and the detection plane of the RHI mode intersect at point C. The distance from the lidar to the cross-point C is 503 meters, and the distance between the threshold and cross-point C is 1468 meters.The specific configuration of field measurement is summarized in Table2. The runway 02R is for landing. For the aircraft following the standard approaching procedure, the altitude of the landing aircraft at point C is 78 meters, which can be covered by the scan range of the lidar. The scanning rate is 1°/s, and the elevation ranges from 0° to 10°. Thus, the scan duration is about 10 seconds, meeting the demands for wake vortex investigation.Table 2
Experimental configuration for the RHI scanning mode.
Parameters (unit)
Value
Azimuth angle (°)
112
Scanning rate (°/s)
1
Elevation range (°)
0–10
Elevation angle resolution (°)
0.2± 0.03
Detection radial range (m)
45–885
Longitudinal resolution (m)
15
Distance between pointA and B (m)
503
Distance between pointB and C (m)
1468A day with good weather was selected to avoid the influence of rain and fog. The dataset contained 344 samples in one day.


## 3. Methodology

The lidar can only measure the velocity at the direction of the line-of-sight (LOS) because of the limitation of the principle of the Doppler effect. The LOS velocity also called radial velocity is discrete, so the set of the radial velocity of the wind field is defined asVr:(1)Vr=vrri,θj,i∈1,n,j∈1,m,where ri is the i-th range which means the distance from the lidar, θj is j-th angle which denotes the pitch angle of the lidar, and vrri,θj indicates the radial velocity of the point at coordinate ri,θj.The sample data of Airbus 320 are visualized in Figure3, where the color bar represents the velocity. It should be noted that the area of the wake vortex colored with a pair red and green is only a small part of the entire detected area.Figure 3
Visualization of lidar data for Airbus-A320.
(a)
(b)
(c)


### 3.1. Feature Extraction

#### 3.1.1. Velocity Envelop at Different Angle or Range

The envelope of radial velocityvr is the main feature when the aircraft produces the wake vortices at the light wind. Drri and Daθj are introduced as follows:(2)Drri=maxvrri,θj−minvrri,θj,∀j∈1,m,Daθj=maxvrri,θj−minvrri,θj,∀i∈1,n,where Drri and Daθj are the max difference radial velocities at each range ri and each angle θj, respectively. Then, the core position of the vortices can be determined in terms of the difference between the maxima (for the positive envelope) and the minima (for the negative envelope).


#### 3.1.2. Velocity Difference of Rectangular Area

The envelopes in Section3.1.1 only consider the difference of maximum positive and negative velocities which may not be the true area of the aircraft wake vortices due to the existence of the nonuniform background wind field. As shown in Figure 4, the area colored red as the maximum velocity difference at the range about 675 meters does not contain the wake vortices. Actually, the most significant feature of the aircraft wake vortices is the pair of positive and negative velocities, which is located in a rectangular area. Thus, Drec is defined to express this feature:(3)Drecri,θj=∑vr,UL+vr,LR−∑vr,UR+vr,LL,where vr,UL, vr,UR, vr,LL, and vr,LR are the average radial velocities of upper-left, upper-right, lower-left, and lower-right areas, respectively, as shown in Figure 5. The specific size of the total area is 120 meters in width and 60 meters in height, and the area is divided into four same size subareas, which can cover most of the wake vortices.Figure 4
Example of the velocity envelopes in Airbus-A320 for the vortex circulation.Figure 5
Wake vortex regions in the detection area.


### 3.2. Environmental Parameters

The aircraft wake vortices are affected by various factors, such as temperature, humidity, wind shear, and turbulence. In order to increase the accuracy and robustness of the model, the temperatureT, the average velocity vwv, and the background wind field direction θwd and the air pressure P are considered. The data are collected by the automatic meteorological station of CAAC Southwest Air Traffic Management Bureau, Chengdu, with a data acquisition rate of 1 Hz.


### 3.3. Support Vector Machine

Support vector machine (SVM) is a powerful machine learning method that finds the hyperplane to classify the given vectors. The vector defined in Section3.1 and 3.2 is(4)xi=Dr,Da,Drec,T,vwv,θwd,P,yi∈−1,+1,where i∈1,N, yi=−1 denotes no wake vortex, and yi=+1 is the wake vortex.It becomes a typical classification problem with SVM: finding the maximum-margin hyperplane [10, 11]:(5)argminw,b,ξ12wTw+C∑i=1lξis.t.yiwTxi+b≥1−ξi,ξi≥0.Because feature vectorx1,y1,x2,y2,…,xn,yn is not linearly separable, the vectors are mapped into a higher dimensional space by kernel function ϕ:(6)argminw,b,ξ12wTw+C∑i=1lξis.t.yiwTϕxi+b≥1−ξi,ξi≥0,where C is the penalty parameter of the error term, ϕ is the function that maps the training vector xi into a higher dimensional space, ξ is the linear constraint, and wϕxi+b=0 is the hyperplane.The number of features is 7, which is less than the number of vectors in the dataset. Thus, the radial basis function (RBF) is chosen as the kernel function to reduce the complexity:(7)Kxi,xj=e−γx−y,where γ is the parameter of RBF kernel and is equal to 1/2σ2.The dimension of feature increases with the number of samples in the data set after choosing the RBF kernel andϕxj:(8)ϕxi=e−γxi−x1,e−γxi−x2,…,e−γxi−xm.


### 3.4. Algorithm for Wake Vortex Identification

The entire procedure for wake vortex identification is illustrated by the flowchart in Figure6.Figure 6
Flowchart of the proposed wake vortex finding algorithm.


## 3.1. Feature Extraction

### 3.1.1. Velocity Envelop at Different Angle or Range

The envelope of radial velocityvr is the main feature when the aircraft produces the wake vortices at the light wind. Drri and Daθj are introduced as follows:(2)Drri=maxvrri,θj−minvrri,θj,∀j∈1,m,Daθj=maxvrri,θj−minvrri,θj,∀i∈1,n,where Drri and Daθj are the max difference radial velocities at each range ri and each angle θj, respectively. Then, the core position of the vortices can be determined in terms of the difference between the maxima (for the positive envelope) and the minima (for the negative envelope).


### 3.1.2. Velocity Difference of Rectangular Area

The envelopes in Section3.1.1 only consider the difference of maximum positive and negative velocities which may not be the true area of the aircraft wake vortices due to the existence of the nonuniform background wind field. As shown in Figure 4, the area colored red as the maximum velocity difference at the range about 675 meters does not contain the wake vortices. Actually, the most significant feature of the aircraft wake vortices is the pair of positive and negative velocities, which is located in a rectangular area. Thus, Drec is defined to express this feature:(3)Drecri,θj=∑vr,UL+vr,LR−∑vr,UR+vr,LL,where vr,UL, vr,UR, vr,LL, and vr,LR are the average radial velocities of upper-left, upper-right, lower-left, and lower-right areas, respectively, as shown in Figure 5. The specific size of the total area is 120 meters in width and 60 meters in height, and the area is divided into four same size subareas, which can cover most of the wake vortices.Figure 4
Example of the velocity envelopes in Airbus-A320 for the vortex circulation.Figure 5
Wake vortex regions in the detection area.


## 3.1.1. Velocity Envelop at Different Angle or Range

The envelope of radial velocityvr is the main feature when the aircraft produces the wake vortices at the light wind. Drri and Daθj are introduced as follows:(2)Drri=maxvrri,θj−minvrri,θj,∀j∈1,m,Daθj=maxvrri,θj−minvrri,θj,∀i∈1,n,where Drri and Daθj are the max difference radial velocities at each range ri and each angle θj, respectively. Then, the core position of the vortices can be determined in terms of the difference between the maxima (for the positive envelope) and the minima (for the negative envelope).


## 3.1.2. Velocity Difference of Rectangular Area

The envelopes in Section3.1.1 only consider the difference of maximum positive and negative velocities which may not be the true area of the aircraft wake vortices due to the existence of the nonuniform background wind field. As shown in Figure 4, the area colored red as the maximum velocity difference at the range about 675 meters does not contain the wake vortices. Actually, the most significant feature of the aircraft wake vortices is the pair of positive and negative velocities, which is located in a rectangular area. Thus, Drec is defined to express this feature:(3)Drecri,θj=∑vr,UL+vr,LR−∑vr,UR+vr,LL,where vr,UL, vr,UR, vr,LL, and vr,LR are the average radial velocities of upper-left, upper-right, lower-left, and lower-right areas, respectively, as shown in Figure 5. The specific size of the total area is 120 meters in width and 60 meters in height, and the area is divided into four same size subareas, which can cover most of the wake vortices.Figure 4
Example of the velocity envelopes in Airbus-A320 for the vortex circulation.Figure 5
Wake vortex regions in the detection area.


## 3.2. Environmental Parameters

The aircraft wake vortices are affected by various factors, such as temperature, humidity, wind shear, and turbulence. In order to increase the accuracy and robustness of the model, the temperatureT, the average velocity vwv, and the background wind field direction θwd and the air pressure P are considered. The data are collected by the automatic meteorological station of CAAC Southwest Air Traffic Management Bureau, Chengdu, with a data acquisition rate of 1 Hz.


## 3.3. Support Vector Machine

Support vector machine (SVM) is a powerful machine learning method that finds the hyperplane to classify the given vectors. The vector defined in Section3.1 and 3.2 is(4)xi=Dr,Da,Drec,T,vwv,θwd,P,yi∈−1,+1,where i∈1,N, yi=−1 denotes no wake vortex, and yi=+1 is the wake vortex.It becomes a typical classification problem with SVM: finding the maximum-margin hyperplane [10, 11]:(5)argminw,b,ξ12wTw+C∑i=1lξis.t.yiwTxi+b≥1−ξi,ξi≥0.Because feature vectorx1,y1,x2,y2,…,xn,yn is not linearly separable, the vectors are mapped into a higher dimensional space by kernel function ϕ:(6)argminw,b,ξ12wTw+C∑i=1lξis.t.yiwTϕxi+b≥1−ξi,ξi≥0,where C is the penalty parameter of the error term, ϕ is the function that maps the training vector xi into a higher dimensional space, ξ is the linear constraint, and wϕxi+b=0 is the hyperplane.The number of features is 7, which is less than the number of vectors in the dataset. Thus, the radial basis function (RBF) is chosen as the kernel function to reduce the complexity:(7)Kxi,xj=e−γx−y,where γ is the parameter of RBF kernel and is equal to 1/2σ2.The dimension of feature increases with the number of samples in the data set after choosing the RBF kernel andϕxj:(8)ϕxi=e−γxi−x1,e−γxi−x2,…,e−γxi−xm.


## 3.4. Algorithm for Wake Vortex Identification

The entire procedure for wake vortex identification is illustrated by the flowchart in Figure6.Figure 6
Flowchart of the proposed wake vortex finding algorithm.


## 4. Experiment and Results

### 4.1. General Setting

In order to avoid the range of different features, the attributes are linearly scaled in the range [0, 1]. The values of features in the vector are represented as real numbers in SVM because the performance is significantly influenced if the attributes in greater numeric ranges dominate those in smaller numeric ranges [12]. Therefore, after the features are extracted from the radial velocity of the wind field, all the attributes are linearly scaling into 0,1, as shown in Table 3.Table 3
Range of attributes.
Attributes (unit)
Range
Average
D
r (m/s)
1.51–47.21
4.72
D
a (m/s)
0.00–31.66
4.59
D
rec (m/s)
0.00–13.30
1.44
T (°C)
25.10–29.10
27.32
v
wv (m/s)
0.14–4.13
2.63
θ
wd (°)
57.00–359.00
—
P (hPa)
1005.51–1008.37
1006.74


### 4.2. Model Training

The two parametersC and γ of RBF kernel affect the performance of SVM. In order to achieve the best performance, grid search is performed on C and γ, which is the most used method in SVM. The tool used for solving the SVM is the scikit-learn [13, 14] package in Python based on LIBSVM. The cross-validation is the method used to evaluate the performance of the model with different C and γ. The train dataset is divided into 5 folds. The first range of the grid search is C=2−5,2−4,…,214,215 and γ=2−15,2−14,…,22,23, the cross-validation score of the 1st time grid search is shown in Figure 7(a), and C=4096.00 and γ=1.00 are the best parameters with the score of 0.8296. The grid search is performed again in the range C=210,210.25,…,213.75,214 and γ=2−2,2−1.75,…,21.75,22, and the score increases to 0.8333 when C=2896.31 and γ=1.41 at this time. The results are shown in Figure 7(b).Figure 7
(a) Score of the 1st time grid search. (b) Score of the 2nd time grid search.
(a)
(b)


### 4.3. Performance Assessment

In order to test the feasibility of the proposed method in aircraft wake vortex identification, the performance was evaluated by measures of F1-score, receiver operating characteristic (ROC) curves [15], average conditional probability (ACP), and approximate correlation (AC). The F1-score conveys the balance between the precision and the recall. The ROC curves explore the effects on the true positive rate (TPR) and the false positive rate (FPR) as the position of an arbitrary decision threshold is varied, and the area under the ROC curve (AUC) with a larger value indicates more accurate identification. The ACP which ranges from 0 to 1 and the AC which ranges from −1 to 1 are used as the measures of overall identification accuracy.


### 4.4. Results

The performance of the final model withC=2896.31 and γ=1.41 is shown in Table 4. The number of support vectors is 42 among the 240 samples in the training set. To the best of the authors’ knowledge, this may be the first attempt to use machine leaning strategies in the field of wake vortex identification. In Pan et al.’s [16] method, the authors have used the k-nearest neighbor (KNN) to characterize the different features between wake vortices and background through the distance metric. To evaluate the validity of the proposed identification model, a comparison with the KNN method is conducted using the same dataset.Table 4
Confusion matrix of the SVM classifier withC=2896.31 and γ=1.41.
Dataset
True positive
False positive
False negative
True negative
Test set
29
6
17
51The performance of the test set predicted by KNN and SVM, which are trained using the same training set, is shown in Table5. It can be observed that the proposed method provides a consistently higher identification rate on the dataset than the KNN. Figure 8(a) is the plot of the ACP and AC at different classification thresholds. Using a threshold value of 44 percent, a maximum ACP of 0.805 and a maximum AC of 0.609 are obtained. Similarly, as shown in Figure 8(a), the KNN classifier with a k value of 38 provides the maximum classification accuracy of ACP = 0.803 and AC = 0.607. Clearly, the best performance is obtained with the SVM-based method. Considering the accuracy of the pulsed Doppler lidar measurements, the results of the proposed SVM-based method are sufficiently high. An additional classification experiment was designed to assess the proposed method. Figures 8(b) and 8(d) show the ROC curves of the proposed method by using the dataset. It can be observed that the two methods perform well (AUC > 0.8) in the classification of wake vortices. The method based on SVM achieves a low classification performance (AUC = 0.81). Consider that the wake vortex data measured by using pulsed the Doppler lidar provide signal uncertainty and eccentricity caused by sensor scanning patterns, background winds, and a receiver’s adjusted gain response. The results obtained using the proposed method are acceptable.Table 5
Performance of trained KNN and SVM classifiers on the same test set [16].
Model
Accuracy
Precision
Recall
F1-score
KNN
0.77
0.81
0.63
0.79
SVM
0.78
0.83
0.63
0.80Figure 8
AC and ROC Curves. (a) KNN AC curve. (b) KNN ROC curve. (c) SVM AC curve. (d) SVM ROC curve.
(a)
(b)
(c)
(d)


## 4.1. General Setting

In order to avoid the range of different features, the attributes are linearly scaled in the range [0, 1]. The values of features in the vector are represented as real numbers in SVM because the performance is significantly influenced if the attributes in greater numeric ranges dominate those in smaller numeric ranges [12]. Therefore, after the features are extracted from the radial velocity of the wind field, all the attributes are linearly scaling into 0,1, as shown in Table 3.Table 3
Range of attributes.
Attributes (unit)
Range
Average
D
r (m/s)
1.51–47.21
4.72
D
a (m/s)
0.00–31.66
4.59
D
rec (m/s)
0.00–13.30
1.44
T (°C)
25.10–29.10
27.32
v
wv (m/s)
0.14–4.13
2.63
θ
wd (°)
57.00–359.00
—
P (hPa)
1005.51–1008.37
1006.74


## 4.2. Model Training

The two parametersC and γ of RBF kernel affect the performance of SVM. In order to achieve the best performance, grid search is performed on C and γ, which is the most used method in SVM. The tool used for solving the SVM is the scikit-learn [13, 14] package in Python based on LIBSVM. The cross-validation is the method used to evaluate the performance of the model with different C and γ. The train dataset is divided into 5 folds. The first range of the grid search is C=2−5,2−4,…,214,215 and γ=2−15,2−14,…,22,23, the cross-validation score of the 1st time grid search is shown in Figure 7(a), and C=4096.00 and γ=1.00 are the best parameters with the score of 0.8296. The grid search is performed again in the range C=210,210.25,…,213.75,214 and γ=2−2,2−1.75,…,21.75,22, and the score increases to 0.8333 when C=2896.31 and γ=1.41 at this time. The results are shown in Figure 7(b).Figure 7
(a) Score of the 1st time grid search. (b) Score of the 2nd time grid search.
(a)
(b)


## 4.3. Performance Assessment

In order to test the feasibility of the proposed method in aircraft wake vortex identification, the performance was evaluated by measures of F1-score, receiver operating characteristic (ROC) curves [15], average conditional probability (ACP), and approximate correlation (AC). The F1-score conveys the balance between the precision and the recall. The ROC curves explore the effects on the true positive rate (TPR) and the false positive rate (FPR) as the position of an arbitrary decision threshold is varied, and the area under the ROC curve (AUC) with a larger value indicates more accurate identification. The ACP which ranges from 0 to 1 and the AC which ranges from −1 to 1 are used as the measures of overall identification accuracy.


## 4.4. Results

The performance of the final model withC=2896.31 and γ=1.41 is shown in Table 4. The number of support vectors is 42 among the 240 samples in the training set. To the best of the authors’ knowledge, this may be the first attempt to use machine leaning strategies in the field of wake vortex identification. In Pan et al.’s [16] method, the authors have used the k-nearest neighbor (KNN) to characterize the different features between wake vortices and background through the distance metric. To evaluate the validity of the proposed identification model, a comparison with the KNN method is conducted using the same dataset.Table 4
Confusion matrix of the SVM classifier withC=2896.31 and γ=1.41.
Dataset
True positive
False positive
False negative
True negative
Test set
29
6
17
51The performance of the test set predicted by KNN and SVM, which are trained using the same training set, is shown in Table5. It can be observed that the proposed method provides a consistently higher identification rate on the dataset than the KNN. Figure 8(a) is the plot of the ACP and AC at different classification thresholds. Using a threshold value of 44 percent, a maximum ACP of 0.805 and a maximum AC of 0.609 are obtained. Similarly, as shown in Figure 8(a), the KNN classifier with a k value of 38 provides the maximum classification accuracy of ACP = 0.803 and AC = 0.607. Clearly, the best performance is obtained with the SVM-based method. Considering the accuracy of the pulsed Doppler lidar measurements, the results of the proposed SVM-based method are sufficiently high. An additional classification experiment was designed to assess the proposed method. Figures 8(b) and 8(d) show the ROC curves of the proposed method by using the dataset. It can be observed that the two methods perform well (AUC > 0.8) in the classification of wake vortices. The method based on SVM achieves a low classification performance (AUC = 0.81). Consider that the wake vortex data measured by using pulsed the Doppler lidar provide signal uncertainty and eccentricity caused by sensor scanning patterns, background winds, and a receiver’s adjusted gain response. The results obtained using the proposed method are acceptable.Table 5
Performance of trained KNN and SVM classifiers on the same test set [16].
Model
Accuracy
Precision
Recall
F1-score
KNN
0.77
0.81
0.63
0.79
SVM
0.78
0.83
0.63
0.80Figure 8
AC and ROC Curves. (a) KNN AC curve. (b) KNN ROC curve. (c) SVM AC curve. (d) SVM ROC curve.
(a)
(b)
(c)
(d)


## 5. Conclusion

In this paper, an efficient method based on the soft margin SVM is proposed to identify the aircraft wake vortex. The method is applied to recognize the aircraft wake vortex by using the pulsed Doppler lidar characteristics acquired at the Chengdu Shuangliu International Airport from Aug 16, 2018, to Oct 10, 2018. From the measured data, the velocity envelops are first used to measure the significant feature of the aircraft wake vortex. Then, the grid search and cross-validation based on the soft margin SVM with the RBF kernel are applied to identify the characteristics of the aircraft wake vortex. The experimental results demonstrate that the final classifier withC=2896.31 and γ=1.41 has the best performance in cross-validation, and the precision is more than 70% and the AUC of ROC curve is 0.81. The results also validate that the proposed method has higher accuracy with satisfactory robustness in identifying the aircraft wake vortex. The accuracy of the aircraft wake vortex achieved by the method proposed in this paper is quite sufficient for the analysis of wake vortex behavior under wind field conditions.



---

*Source: 9314164-2020-05-12.xml*