Prediction of Highway Lane Changes Based on Prototype Trajectories
David Augustin Marius Hofmann Ulrich Konigorski
cs.LG cs.RO stat.ML



Prediction of Highway Lane Changes Based on Prototype Trajectories  



David Augustin M.Sc., Marius Hofmann B.Sc., Opel Automobile GmbH; Prof. Dr.-Ing. Ulrich Konigorski, TU Darmstadt 
 






Abstract 

The vision of automated driving is to increase both road safety and efficiency, while offering passengers a convenient travel experience. This requires that autonomous systems correctly estimate the current traffic scene and its likely evolution. In highway scenarios early recognition of cut-in maneuvers is essential for risk-aware maneuver planning.
In this paper, a statistical approach is proposed, which advantageously utilizes a set of prototypical lane change trajectories to realize both early maneuver detection and uncertainty-aware trajectory prediction for traffic participants. Generation of prototype trajectories from real traffic data is accomplished by ahc. During clustering, the alignment of the cluster prototypes to each other is optimized and the cohesion of the resulting prototype is limited when two clusters merge. In the prediction stage, the similarity of observed vehicle motion and typical lane change patterns in the data base is evaluated to construct a set of significant features for maneuver classification via bdt. The future trajectory is predicted combining typical lane change realizations in a mixture model. B-splines based trajectory adaptations  guarantee continuity during transition from actually observed to predicted vehicle states. Quantitative evaluation results demonstrate the proposed concept's improved performance for both maneuver and trajectory prediction compared to a previously implemented reference approach.

Keywords maneuver detection, trajectory prediction, motion pattern recognition, highly automated driving, boosted decision trees



	





§ INTRODUCTION


The development of automated driving functions is a central activity of industry and science. Currently, first assistance systems combining both longitudinal and lateral guidance are introduced on the market.


In order to autonomously drive in public traffic these systems need to correctly evaluate the current traffic scene and its likely evolution. This includes estimating the driving intentions of traffic participants, detecting maneuvers being performed, and predicting future motion of surrounding traffic. All these tasks are non-deterministic and uncertainties in maneuver and motion estimation need to be considered.

In this contribution, highly automated driving (HAD) on highways is addressed with the focus on early lane change maneuver detection and uncertainty-aware trajectory prediction.

Before outlining the basic concepts of the proposed approach in Subsection <ref> an overview of related work in the field of motion prediction is provided.



 §.§ Related Work
 

Motion prediction is a research area since decades. The survey articles <cit.> and <cit.> provide overviews about a variety of published approaches and fields of application related to road transport. Additionally, both provide suggestions on how motion prediction approaches can be categorized.


Regarding driver intention recognition, a large number of approaches do assess the traffic environment on an abstract level while neglecting the prediction of physical motion which is crucial for criticality assessment and trajectory planning. The opposite concept relying purely on physically-based motion prediction is only valid for a limited time horizon of less than a second <cit.>. Approaches that realize both motion detection and trajectory prediction can be classified as methods of pattern recognition and methods fusing dynamic motion models with behavior descriptions <cit.>. Methods of pattern recognition do not only consider the current state of an object but also its history. These states are compared to a database of previously generated motion patterns, each linked to a specific maneuver type. In scenarios like intersections, in which the road topology predefines typical trajectories, the patterns can be directly created from map data <cit.>. In <cit.> trajectory prediction for lane changes is realized combining a physical motion model with artificially created motion patterns. At each time step a new set of motion patterns is generated based on the current vehicle state, the road parameters, and the detected maneuver. Another possibility for pattern creation is to analyze real traffic data. 
The approach <cit.> utilizes a set of recorded highway lane change trajectories as prototypes in a Bayesian network classifier for maneuver detection. Clustering methods are not applied <cit.>. Motion prediction for traffic participants is realized using multi-agent simulation.
In <cit.> a cluster-based technique that learns the typical motion patterns using pairwise clustering is designed for a visual pedestrian tracking system placed in an entry hall. For static environments like interior rooms of buildings or intersections it is fairly easy to define regions of start and end points of motion patterns compared to environments with varying topology like highways. 




 §.§ Proposed Approach
 

In this contribution we propose an extension to a previously published approach for lane change maneuver detection and uncertainty-aware trajectory prediction in highway scenarios <cit.>. The approach is centered around the idea that the variety of lane change realizations can be clustered into typical motion patterns which support early maneuver recognition and trajectory prediction. Similar lane change courses are aggregated into groups, each represented by a prototype trajectory. Figure <ref> shows a schematic diagram of the proposed approach.


The objectives of learning stage are to produce a set of representative lane change trajectories, which can be utilized to formulate descriptive features for maneuver recognition, and to learn a convenient classification algorithm. Prototype trajectories are generated from recorded highway footage of traffic participants applying ahc in combination with a preprocessing step optimizing the trajectory alignment. Additionally, constraints are introduced to improve cluster cohesion and define  the end condition for pairwise clustering. To reduce the diversity of lane change embodiments and thus limit the number of required representatives we categorize lane change trajectories solely by their lateral course. This enables maneuvers of different absolute velocity but similar lateral movement to be cumulated in the same cluster. Each lane change cluster is represented by lateral position and velocity prototypes. Details of typical motion pattern extraction from real traffic data are provided in Chapter <ref>.

In the estimation stage, which is covered in Section <ref>, we buffer the latest states for each voi. The vehicle state histories are referred to as partially observed trajectories. A distance value is computed for each pair of partial and prototype trajectory under optimal alignment for lateral position and velocity independently. For each lane change type (lcl or lcr) we select the pairs with minimal distances as best matches. The distance values of the best matches for lcl and lcr are utilized to construct features for maneuver detection. Additionally, the lateral offset of a vehicle within its lane and its lateral velocity is applied for maneuver inference. All these features are derived from motion states of each voi, individually. Vehicle interdependencies or infrastructure information are not taken into account. This limits the possible prediction horizon because potential causes for lane change maneuvers are not considered. Maneuvers can only be detected during execution. On the other hand, approaches with lower abstraction level limit the requirements on environment sensing system (Section <ref>), i.e. there is no need to detect the vehicle driving in front of a voi in order to reason about its influence on driving intentions. Due to the choice of features the proposed approach can be categorized as maneuver-based driver intention estimation.

The extensions compared to our previous publication <cit.> mainly lie in advances of prototype usage for both maneuver detection and maneuver-based trajectory prediction.
We train an ensemble of decision trees by adaptive boosting for maneuver classification and discuss its performance in comparison to the generative classifier gda. The trajectories of the voi are predicted according to the estimated maneuvers.
In the case of lane changes, the future lateral motion is expressed as a weighted sum of all maneuver-related prototypes forming a gmm. The weight of each prototype correlates with its proximity to the partially observed trajectory.
Initial offsets between the actual lateral position of a voi and the trajectory prediction are compensated by applying a B-spline approximation. Longitudinal motion is covered independently by a (nearly) constant acceleration model.
Both classification and motion prediction are evaluated with a test set of recorded highway data. The results are discussed in Section <ref>.









§ ENVIRONMENT SENSING
 

The test vehicle is equipped with various radars and cameras for all-round object recognition. As we are interested in cut-in maneuver detection, we focus on vehicles located in the front of the test car.











A tracking system produces, for each sensed object, a local estimate of its state vector X^(target):

    X^(target) = [x, y, θ, v]^T

where x and y are local Cartesian coordinates, θ the relative heading angle, and v the absolute velocity.
The camera system detects lane markings and approximates their course as cubic functions in the local coordinate system:

    y(x) = c_3 x^3+c_2 x^2 + c_1 x + c_0

with polynomial coefficients c_3, c_2, c_1, and c_0.

To be independent of the ego motion of the test vehicle, the states of all tracked objects are transformed into the Frenet frame. Centerlines of corresponding driving lanes of each traffic object are utilized as reference curves.	
Vehicle trajectories γ_i(t):[0, T_i] →ℂ are expressed in the lane-relative configuration space ℂ which is [s, ṡ, d, ḋ] with s being the traveled distance from the start of the maneuver, ṡ the velocity in longitudinal direction, while d and ḋ being the lateral offset to the lane's centerline and the lateral velocity, respectively.
To cut down diversity of maneuver executions, the configuration space ℂ is reduced to [d,  ḋ] for cluster learning and representation. The lateral offset, d, can be derived using the minimal distance of the vehicle to its left, d_l, and right lane markings, d_r:

    d = 1/2 (d_l + d_r)
 
Thereby, the lane width is taken into account inherently. To reduce computational effort for the calculation of d_l and d_r, lane marker polynomials are locally linearized. The maximal error of this approximation is below our measurement accuracy. The lateral velocity, ḋ, is derived by applying a low-pass filtered difference quotient. The longitudinal velocity, ṡ, is assumed to be equal to the vehicle's absolute velocity, v. The traveled distance, s, between two consecutive measurements is computed applying a constant velocity model <cit.>.









§ MOTION PATTERN RECOGNITION
 

Due to diversity of traffic situations and individual driving styles the courses of lane changes highly vary. Nevertheless, we assume that the variety of lane change realizations can be mapped to a finite number of typical motion patterns. Such a pattern represents a group of similar maneuver executions. As the number of clusters and the course of their representing trajectories are unknown, real highway traffic data is analyzed by a prototype-based clustering technique called ahc. Each maneuver type is analyzed separately. The structured highway environment limits vehicle maneuvers with respect to their lateral motion to lcl, lk, and lcr. Each data point in the data set is assigned to one maneuver type via automatic labeling.



 §.§ Labeling
 

Lane changes can easily be detected in a data set by searching for lane marker crossings. A more challenging aspect for reliable automatic labeling is to define start and end time of a maneuver. Before a vehicle passes the centerline of its starting lane we regard the motion as lane centering and assign it to the lane keeping maneuver. Also a possible overshoot on the target lane is treated as lane keeping. Additionally, it was determined that a maneuver bound is reached when the vehicle's absolute lateral velocity exceeds or undercuts a velocity threshold of 0.2  m/s. To improve the lateral velocity estimate, we approximate the vehicle's position course via B-spline curve and calculate lateral velocity, ḋ, via derivation.




 §.§ Preprocessing and Clustering
 

Objectives of traffic-related motion pattern recognition are to automatically assign similar vehicle trajectories to the same cluster and to compute a common representation for each group, which we denote as prototype trajectory. Due to formulation of goals prototype-based clustering methods, like ahc or k-Means, are suitable solution approaches. Both methods were investigated and compared in <cit.> with ahc performing superior for highway scenarios.

ahc is an algorithm iteratively merging the two closest clusters until only one single cluster exists or a termination condition is fulfilled <cit.>. Initially, each cluster is formed by a unique trajectory of the training data set. When two clusters are merged, they are replaced by a common prototype and their similarity to all other cluster representatives must be reevaluated. Cluster representations consist of time-dependent mean, μ(t), and variance, σ(t), values. The mean of one cluster, C_m, is calculated from its comprised trajectories by:

    μ_m(t) = 1/N_m∑_i=1^N_md_i(t)

with N_m denoting the total number of trajectories in the cluster. For prototype generation, ahc is only applied to position trajectories, d∈𝔻. 𝔻 denotes the set of position trajectories relative to the centerline of the driving lane. For maneuver classification, it turned out to be beneficial to utilize both position and velocity prototypes, so in this case d∈𝔻×𝕍 with 𝕍 denoting the set of velocity trajectories <cit.>.

The scatter of cluster trajectories around the prototype mean is expressed via the variance

    σ_m^2(t) = 1/N_m∑_i=1^N_m (d_i(t) - μ_m(t))^2

In order to express the relation between two trajectories a measure of proximity is required. We chose the average Euclidean distance between the two trajectories d_i, d_j ∈𝔻  as dissimilarity measure:

    δ(d_i, d_j)    = ( 1/T∫_t=t_min^t_max (d_i(t) - d_j(t))^2dt )^1/2
    
    t_min   = min(t_0,i, t_0,j) 
    
    t_max   = max(t_0,i + T_i , t_0,j + T_j) 
    
    T    = t_max - t_min

where t_0,i, t_0,j, T_i, and T_j are the starting times and the total maneuver durations of trajectories d_i and d_j, respectively.	
For comparison of trajectories of unequal length, the shorter ones are extended by keeping their initial and/or final value. The factor 1/T averages the Euclidean distance to negate the influence of duration of the two trajectories being compared on their dissimilarity value.A major influence factor on the dissimilarity between two lane change trajectories is their mutual alignment <cit.>, which in the case of highway trajectories is mainly determined by the applied labeling approach. In order to handle possible labeling inconsistencies we studied different approaches for active alignment of trajectories.	
Best results were achieved by an adjustment strategy minimizing the dissimilarity of each pair of trajectories by computing an individual optimal configuration <cit.>. In ahc the alignment of cluster representatives needs to be updated after each clustering iteration.



Figure <ref> shows prototypical lane change motion patterns created using the presented ahc approach.










§ MANEUVER CLASSIFICATION
 

Motion prediction for traffic participants can be split into two consecutive tasks. In the first step, the performed maneuvers need to be detected. Subsequently, the most probable future trajectories according to the maneuver estimate can be calculated. Maneuver detection is a typical classification problem. The highway topology limits the set of possible lateral maneuvers to C = {LCL, LK, LCR}. We implemented maneuver inference based on a set of features which can be derived from the observed motion of individual vehicles. Vehicle interaction is not taken into account. This approach reduces computational and sensory  requirements. Additionally, not relying on strong assumptions about vehicle interdependencies increases robustness. On the other hand, maneuvers can only be recognized during their execution, thus limiting the possible prediction horizon. The applied features for maneuver detection of a voi are lateral offset within its lane, lateral velocity and proximity of its motion path to prototype lane change trajectories in the database.
The proximity of a partially observed trajectory, d_p, and a cluster prototype, c_m, can be calculated via average Mahalanobis L_1 distance:

    Δ(d_p, c_m, τ_m)    = 1/T_m∫_-T_m^0( ( d_p(t) - μ_m(t+τ_m) )^2/σ_m^2(t+τ_m))^1/2 dt 
    
    T_m    = min(T_buffer, τ_m)

The evaluation period, T_m, is the minimum of buffered partial trajectory and prototype trajectory length. The cluster prototype is shifted by τ_m optimizing the trajectories' adjustment by minimizing their dissimilarity within the recent 0.5  s <cit.>. Compared to the Euclidean distance, utilized during clustering, the Mahalanobis distance measure takes the non-uniform uncertainty of the prototype course into account. Dissimilarities in regions of low variance have increased weight due to the term 1/σ_m^2.
Mahalanobis distance of position trajectories is denoted as Δ_p. The notation in case of velocity trajectories is Δ_v.

During feature generation a vehicle's partial trajectory is matched individually to all prototype trajectories in the database. For each pair {d_p,c_m} the optimal shifting coefficient, τ_m, is calculated minimizing Δ_p,[-0.5:0](d_p,c_m). The velocity-dependent distance Δ_v(d_p,c_m) related to the same cluster is calculated for an equal value of τ_m. As the best matching prototypes for lcl and lcr we select those minimizing the dissimilarity Δ_p for the whole evaluation period. Their dissimilarity values are denoted Δ_p,LCL, Δ_v,LCL, Δ_p,LCR and Δ_v,LCR.

In <cit.> a Quadratic Gaussian Discriminant Analysis (gda) classifier is applied for maneuver inference. In this contribution its performance is compared to bdt, an ensemble of shallow decision trees sequentially trained by a boosting method, described in Section <ref>.



















 §.§ Boosting
 

Boosting is a sequential method for consecutively fitting a set of basic learners to variations of the training data. After each iteration, the whole training set is evaluated <cit.>. Misclassified data points are given greater weights when used to train the next classifier. Basic learners, also referred to as "weak learners" <cit.>, typically have low classification performance when applied individually. In an ensemble of N basic learners, each unit provides a class vote dependent on observed features. Those votes are weighted by α_n and combined like:

    y(f)     =    ∑_n=1^N α_n h_n(f)

with h_n denoting the hypothesis of the n-th basic learner. We apply the Adaboost.M2 algorithm which is applicable in classification problems with multiple class instances <cit.>. In AdaBoost.M2 weighting coefficients α_n are derived from weighted pseudo-loss functions, providing more accurate basic classifiers with a greater weight.

As basic classifiers we use shallow decision trees limiting the maximum number of branch nodes to 15. The ensemble classifier comprises 90 basic decision trees. Treated individually, decisions are easily interpretable. They are not limited to a specific type of data. The flexibility of decision trees is extended by enlargement to ensemble methods at the expense of easy interpretability. AdaBoost application does not require hyperparameter training and showed best performance in an extensive empirical comparison <cit.>.









§ MANEUVER-BASED TRAJECTORY PREDICTION
 

After successful driver intention recognition, the future trajectory of a vehicle is predicted according to the estimated maneuver. In <cit.> the potential of prototype trajectories for both longitudinal and lateral motion prediction was demonstrated. In case of detected lane changes, the best matching prototype of the respective maneuver type, lcl or lcr, was selected as an estimate for future lateral motion. Here, we extend this idea by formulating a Gaussian mixture model combining and weighting the predictions of all prototype trajectories of the estimated maneuver type. The predicted lateral course can be calculated as a weighted sum of the prototype means:

    μ(t) = ∑_m w_m μ̂_m(t)

and the combined variance can be computed by:

    σ^2(t)     =    ∑_m w_m ( σ̂^2_m(t) + μ̂^2_m(t) ) - μ^2(t)

The final values of mean and variance are hold when individual prototypes are evaluated outside their definition range.

    μ̂_m(t)     =    μ_m(τ_m + t),    τ_m + t < T_m
    μ_m(T_m),    otherwise
    σ̂^2_m(t)     =    σ^2_m(τ_m + t),    τ_m + t < T_m
    σ^2_m(T_m),    otherwise.

The weighting factors, w_m, are derived from the ratio of inversed Mahalanobis distance of each prototype representing a motion pattern of the detected maneuver to the partial trajectory:

    w_m = 1/Δ(d_p, μ_m)/∑_m1/Δ(d_p, μ_m).

Figure <ref> depicts the mean and variance values of the prototype ensemble for different prediction horizons.

The longitudinal motion during a lane change maneuver can be described by a (nearly) constant acceleration model <cit.>, also referred to as discrete Wiener process acceleration model <cit.>. It is given by:

    [   s_k+1; v_s_k+1; a_s_k+1 ] 
    =
    [      1      T 1/2T^2;      0      1      T;      0      0      1 ][   s_k; v_s_k;     a ]
    +
    [ 1/2T^2;      T;      1 ]ω_a

with

    ω_a = 𝒩(0, σ^2_a )

In <cit.> the constant acceleration value is calculated averaging the longitudinal accelerations of all trajectories comprised in the cluster, which is the best fit to the vehicles' partial trajectory. σ_a^2 is the related variance. Using the  Gaussian mixture model, the new values of a and σ_a^2 can be computed via Equations <ref> and <ref>, respectively.






 §.§ Smooth Transition to Prototype-based Trajectory Predictions


The future lateral motion course, predicted based on the previously proposed model of mixed cluster representatives, does not necessarily pass through the currently observed state of the vehicle. Therefore, the transition from current lateral offset to future values is generally discontinuous. To solve this problem, we adapt all cluster prototypes that contribute to the predicted trajectory within the interval [τ_m, τ_m+1  s]. This ensures continuity in lateral position, velocity, and acceleration also in the near future. The adaptations are based on transformation of the lateral position course into B-splines. B-splines and their basis functions have favorable characteristics allowing adaptations in limited scope (local support characteristic)  and guaranteeing continuity of a function and their derivatives of certain degree (continuity characteristic) <cit.>.

B-splines are piecewise polynomial curves of degree d defined on an interval [a, b]. They are represented by a knot vector t = [t_0, t_1, ... , t_m] with m + 1 elements and a set of n + 1 control points P = {P_0,...,P_n}. The relation n + 1 = m + d is always valid. 

B-Spline curves are represented by 

    C(t) = ∑_i=0^nP_i N_i,d(t)

where the basis functions N_i,d(t) are recursively defined as

    N_i,0(t)     =    
    1,     t ∈ [t_i,t_i+1)
    
    0,    otherwise
    
    N_i,d(t)     =    t - t_i/t_i+d - t_i N_i,d-1(t) + t_i+d+1 - t/t_i+d+1 - t_i+1 N_i+1,d-1(t).

Initially, a cluster's polynomial representation f(t) = a_0 + a_1 t + ... + a_n t^n is converted into a B-spline curve by scaling the polynomial to the interval [0,1] and calculating the control points <cit.> by 

    P_i = ∑_j=0^iij/dja_j.

Due to the local support characteristic, the control point P_i only affects the interval [t_i,t_i+p+1). We use this property to modify the trajectory within the first second of prediction by adding three knots at the interval boundaries, respectively. C(t) is d-k = 2 times continuously differentiable at these breakpoints, with multiplicity k = 3 corresponding to the number of identical knots in the knot vector. The insertion of a new knot t̂∈ [t_s, t_s+1) causes a new control point. As stated in <cit.>, the elements of the control point set are updated by 

    P̂_i     =    P_i,     0 ≤ i ≤ s-d
    
    (1-α_i) P_i-1 + α_i P_i,     s-d+1 ≤ i ≤ s
    P_i-1,     s+1 ≤ i ≤ n+1
    
    α_i     =    t̂ - t_i/t_i+d - t_i.

Within the first prediction second only a limited number of basis functions are unequal to zero. Due to  knot insertion the non-zero basis functions are N_4, 5, N_5, 5 and N_6, 5. The control points of these basis functions need to be adjusted such that the B-spline meets the currently observed vehicle states which are position, velocity and acceleration.

According to <cit.>, the r-th derivative of a B-spline is defined as  

    C^(r)(t) = ∑_i=0^n-rP^(r)_i N^(r)_i,d-r(t)


where


    P^(r)_i = (d-r+1)P^(r-1)_i+1-P^(r-1)_i/t_i+d+1 - t_i+r


with P^0_i = P_i. The knot vector is reduced to t_r, ..., t_m-r.
This yields the following equation system with three undefined variables P_4, P_5, and P_6:

    x     =    P_4 N_4, 5(τ_m) + P_5 N_5, 5(τ_m) + P_6 N_6, 5(τ_m)
    
    
    x'     =    P^(1)_4 N^(1)_4, 4(τ_m) + P^(1)_5 N^(1)_5, 4(τ_m)
    
    
    x”    =    P^(2)_4 N^(2)_4, 3(τ_m).


The derivatives can be calculated via Equation <ref>.
The vehicle's lateral velocity and acceleration are estimated via derivation of a third-order polynomial approximation of the partial trajectory and inserted into Equations <ref>, <ref>, and <ref> to solve for the values of the required control points.









§ RESULTS
 

For evaluation and comparison of proposed approaches for maneuver estimation and subsequent trajectory prediction we use the data set proposed in <cit.>. Aggregating hours of real highway footage, the whole data set comprises a total number of 434 lane change maneuvers of traffic participants of which 156 are lane changes to the left and 278 to the right neighbor lane, respectively.
Data is divided with a ratio of 70/30 into training and test set. During data division it is specially requested that ratio of maneuver types is equal in both fractions of the data set.	
In terms of maneuver classification we compute performance measures juxtaposing gda and boosting decision trees. As the distribution of classes is highly unequal, with lane keeping clearly outnumbering lane change maneuvers, balanced versions of performance measures are applied.
The advances in trajectory prediction, namely combining prototype trajectories in a Gaussian mixture model and B-spline based trajectory smoothing, are compared to the previously implemented approach <cit.>, which utilized solely the best matching cluster representative for future motion prediction.



 §.§ Evaluation of Maneuver Detection
 

Due to the highly unequal class incidences, balanced measures <cit.> must be applied for classification performance evaluation, whenever a quality criterion combines actual positive and actual negative portions of the regarded class (e.g. precision):


    TPR    =    TP/TP+FN = 1 - FNR
    PRC    =    TPR/TPR+FPR
    F_1     =    2 ·PRC·TPR/PRC + TPR


The True Positive Rate TPR, also denoted recall, states how likely a maneuver will be predicted by the respective approach. The precision PRC indicates how likely a predicted lane change will actually happen. The F_1 score is the harmonic mean of both precision and recall. Those values are computed on the whole test data set. Table <ref> shows the classification performance for the proposed methods, Quadratic Gaussian Discriminant Analysis and Boosted Decision Trees, using two different sets of features, respectively. As a reference, to investigate the influence of features constructed from cluster prototypes, both methods are evaluated utilizing a two-dimensional feature vector f=[d,ḋ]^T comprising lateral offset, d, and lateral velocity, ḋ, of a vehicle in Frenet coordinates. In case of gda we constructed two features utilizing distance values to the best matching prototypes by:


    f_3     = Δ_p,LCR - Δ_p,LCL
    
    f_4     = Δ_v,LCR - Δ_v,LCL


Thereby, the dimension of the feature vector can be kept low while achieving good prediction results. For bdt classifier we expand the feature vector to six dimensions f=[d, ḋ, Δ_p,LCR, Δ_p,LCL, Δ_v,LCR, Δ_v,LCL]^T, allowing the ensemble model to learn mutual dependencies of features. Application of this feature vector for gda worsens prediction quality, most likely because model assumptions are violated.
For the sake of comparison all classifiers were trained allowing a miss rate for the lane keeping maneuver of MCR_LK = 0.11. This is realized via adaptation of prior probabilities for gda and initial weights of the data points for bdt.






 §.§ Discussion of Maneuver Detection Performance
 

First, the influence of features constructed from cluster prototypes on the classification performance is investigated for both learning methods.
gda strongly benefits from additional information provided by prototype trajectories. The average prediction time can be increased by 0.14  s to a total of 1.67  s. The performance measures show improvements, as well. The increase in average prediction time of 0.06  s for the bdt method by adding prototype-derived features is comparatively small, but again, the prediction quality measures denote performance improvements.


For the two-dimensional feature vector f = [d, ḋ]^T bdt performs significantly better reaching an improved prediction time and higher F_1-scores. Whereas the precision values are on equal level, the major improvements are reported by increased recall values for lane change maneuvers, which means that actual lane changes are more likely to be detected. Incorporating information provided by motion pattern recognition, both gda and bdt achieve equal average prediction horizons of roughly 1.7  s, while, again, bdt method shows better values for performance measures. The optimal average prediction horizon according to the automatically labeled data can be calculated to 2  s. 




 §.§ Evaluation of Motion Prediction
 

An advantage of prototype-based maneuver detection lies in the usability of cluster prototypes for the future motion prediction of traffic participants. Here, the extended approach for trajectory prediction, as described in Section <ref>, is evaluated and compared to the original version presented in <cit.>.


Figure <ref> shows the mean absolute lateral and longitudinal errors between traffic participants' predicted and actual positions for different prediction horizons. The evaluation results for the reference approach is illustrated on left, while the new motion prediction results are depicted on the right side of the figure.
The prediction horizon [0  s, 4 s], which is plotted on the horizontal axis, is subdivided into intervals of Δ t = 0.5  s.
For the n-th prediction interval, the mean absolute position error is calculated by


    ϵ_n = 1/N_n·∑_τ∈ [n Δ T, (n+1) Δ T) | x_t - x̂_(t|t-τ) |


where N_n is the total number of summands, x_t = [s_t, d_t] denotes the actual observation at time t and x̂_(t|t-τ) is the prediction of vehicle state at time t calculated at time t - τ. The maximum evaluation time is limited by the length of the prototype ensemble used for motion prediction and capped to a value of 4  s.
The evaluation starts when the correct lane change maneuver is detected continuously and ends at the point in time when the vehicle crosses the lane marking.





 §.§ Discussion of Motion Prediction Performance
 

For longitudinal trajectory prediction during lane change maneuvers both approaches perform equally well. Comparing both lane change types, the errors are generally larger in case of lcl maneuvers. The reason for this lies in the larger diversity of longitudinal accelerations during lane changes to the left.


The performance in lateral motion prediction is significantly improved by the extensions presented in Section <ref>. For small prediction horizons, the absolute position errors correspond approximately to the measurement accuracy. For a prediction of horizon of 2 - 4  s, the absolute lateral error is steady with an approximate value of 0.3  m. 		
In case of  the reference approach, error peaks for medium prediction horizons can be observed. These peaks result from incorrect selections of prototype trajectories in an early stage of the detected maneuver. At the beginning many cluster prototypes are similar so that wrong prototype selection is quite likely. The combination of prototypes in a Gaussian mixture model can compensate this issue.


Figure <ref> plots the average Mahalanobis distance between traffic participants' predicted and actual lateral offsets. The upper bound for all considered prediction horizons is 0.7  σ, which indicates that the  variance of the prototype ensemble, which is calculated utilizing Equation <ref>, is a good estimate of prediction uncertainty.

  








§ CONCLUSION AND OUTLOOK


In this paper, a statistical approach is proposed, which successfully utilizes a set of prototypical lane change trajectories to realize both early maneuver detection and uncertainty-aware trajectory prediction for traffic participants. Generation of prototype trajectories from real traffic data is accomplished by ahc, optimizing mutual cluster configurations and constraining cohesion during the clustering process <cit.>. The main contributions of this paper exist in targeted adaptations of prediction methods to improve utilization of motion patterns and thus enhancing prediction results. For maneuver recognition, we implemented bdt classification increasing the detection rate of lane change maneuvers compared to previously applied gda. The future trajectory is predicted according to typical realizations of the estimated maneuver. We introduced a mixture model of cluster prototypes and demonstrated increased accuracy for lateral motion prediction. Additionally, a B-splines based adaptation technique is described guaranteeing continuity during transition from actually observed to predicted vehicle states.
In future work, the prediction framework is extended by taking interaction between traffic participants into account to further increase the prediction horizon and lower misclassification rates. 



unsrt









