<figure><img src="image_1.png" /></figure>

HAL
open science

# Powered Two-Wheeler Riding Profile Clustering for an In-Depth Study of Bend-Taking Practices

Mohamed Diop, Abderrahmane Boubezoul, Latifa Oukhellou, Stéphane Espie

► To cite this version:

Mohamed Diop, Abderrahmane Boubezoul, Latifa Oukhellou, Stéphane Espie. Powered Two-Wheeler Riding Profile Clustering for an In-Depth Study of Bend-Taking Practices. Sensors, 2020, 20 (22), 19p. 10.3390/s20226696. [https://hal.science/hal-03609529](https://hal.science/hal-03609529)

HAL Id: hal-03609529

https://hal.science/hal-03609529

Submitted on 15 Mar 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

Article

# Powered Two-Wheeler Riding Profile Clustering for an In-Depth Study of Bend-Taking Practices

Mohamed Diop<fnref n="1" />, Abderrahmane Boubezoul <fnref n="1" />,* and Latifa Oukhellou <fnref n="2" />

Received: 17 October 2020; Accepted: 20 November 2020; Published: 23 November 2020

<figure><img src="image_4.png" /></figure>

check for
updates

**Abstract:** The understanding of rider/vehicle interaction modalities remains an issue, specifically in the case of bend-taking. This difficulty results both from the lack of adequate instrumentation to conduct this type of study and from the variety of practices of this population of road users. Riders have numerous explanations of strategies for controlling their motorcycles when taking bends. The objective of this paper is to develop a data-driven methodology in order to identify typical riding behaviors in bends by using clustering methods. The real dataset used for the experiments is collected within the VIROLO++ collaborative project to improve the knowledge of actual PTW riding practices, especially during bend taking, by collecting real data on this riding situation, including data on PTW dynamics (velocity, normal acceleration, and jerk), position on the road (road curvature), and handlebar actions (handlebar steering angle). A detailed analysis of the results is provided for both the Anderson-Darling test and clustering steps. Moreover, the clustering results are compared with the subjective data of subjects to highlight and contextualize typical riding tendencies. Finally, we perform an in-depth analysis of the bend-taking practices of one subject to highlight the differences between different methods of controlling the motorcycle (steering handlebar vs. rider's lean) using the rider action measurements made by pressure sensors.

**Keywords:** powered two-wheeler; riding profiles; riding activity analysis; data mining; wearable sensors

## 1. Introduction

In recent years, increasingly numerous cities in France and elsewhere in Europe have been experimenting with or setting up car-free zones. The objective is to limit the usage of private cars. Consequently, powered two-wheeled vehicles (PTWs) are becoming increasingly popular. In France, PTW users are particularly exposed to road accidents and fatalities. According to the French Observatory of Road Safety (ONISR) [1], PTWs barely constitute 2% of road traffic, but they account for 23% of road deaths and 44% of serious accidents, and PTW users are 18 times more likely to be killed than a car driver when normalizing for the number of kilometers traveled. In 2017, in France, 32% and 39% of fatal accidents of mopeds and motorcycles, respectively, occurred without an identified third party. In particular, loss of control in bends accounts for more than 42% of single-vehicle accidents. An understanding of the bend-taking behavior plays a significant role in scientific knowledge improvement and road safety. The in-depth analysis of rider/PTW interaction modalities remains an open question.

The literature on motorcycle riding study can be categorized into two categories: the first category addresses the rider's behavior analysis, and the second addresses rider modeling. A general overview

<footnote n="1">¹ TS2-SIMU&MOTO, Université Gustave Eiffel, IFSTTAR, F-77447 Marne-la-Vallée, France; mohamed.diop@univ-eiffel.fr (M.D.); stephane.espie@univ-eiffel.fr (S.E.)</footnote>

<footnote n="2">² COSYS-GRETTIA, Université Gustave Eiffel, IFSTTAR, F-77447 Marne-la-Vallée, France; latifa.oukhellou@univ-eiffel.fr</footnote>

<footnote>* Correspondence: abderrahmane.boubezoul@univ-eiffel.fr</footnote>
of this categorization is given in [2,3]. In the first category, the most important topics addressed are the handlebar control (steering torque versus steering angle), dominance of different types of control (handlebars versus rider lean) [4,5], and differences between experienced and novice riders [6]. More recently, in [7,8] the authors addressed the problem of motorcycle rider model identification based on experimental data recorded during road tests with a fully instrumented motorcycle.

The main scientific challenge to overcome is to design a model of steering control most similar to that achieved by PTW riders. To achieve this aim, it is essential to represent the sensorimotor processes used by a rider to control the vehicle trajectory through different modalities of actions on the vehicle at an appropriate level of complexity [9]. From a technical standpoint, the question relates to the embedded measurement of the rider/PTW interactions and how to capture these interactions. Accordingly, an experimental study has been carried out by using both pressure sensors and pressure gauges. Pressures have been recorded at the following contact points: left and right handlebar, saddle, tank, and footrests. These measures have been added to those related to vehicle dynamics.

After collecting all of these data, we address the following scientific question in this ongoing research work: Can the rider's behavior be reduced to a rider "profile" or categorized within classes? We postulate that the rider can be "reduced" to a psychometric problem and that these different strategies should be consequently identified using discriminative behavioral markers.

Various research works in the scientific literature focus on driving behavior study and driving style characterization. The approaches based on machine learning techniques are widely used for this purpose. The driver's behavior is analyzed using the signals issued from real experiments. Based on the features extracted from signals, the proposed approaches aim at clustering the subjects into different behavioral groups [10–12]. In this paper, we follow this line of research to study the global behavior of motorcyclists.

In comparison with the existing literature about driver profiles [13,14], research into PTW rider behavior based on real bend-taking manoeuvres is limited. Some research has been undertaken on simulators; in this regard, we can cite the works in [15,16].

In this paper, as a first step in riding profile classification “recognition”, exploratory data analysis is chosen, which thus permits us to analyze, discover, and potentially prove certain hypotheses about real riding practices. This approach will help us provide solutions to the following problems; how to visualize the rapid and robust classification based on the quantitative variables, how to combine the classification produced based on quantitative data with the remaining qualitative variables to interpret the classification, and how to provide evidence for the most important explanatory variables for each identified class.

The approach presented in this paper is an easy-to-use methodology based on data mining approaches suitable for the study of the riders' behavior and for the understanding of the way they interact with their vehicle when negotiating bends. In this study, we are trying to fill the knowledge gap on how a range of riders approach bends in real life. The intra- and inter-variability of riders, the particular dynamics of these vehicles, and the tight coupling between the rider and the vehicle make the problem particularly complex. The main contributions of this paper are summarized in the following.

* We propose an easy-to-use methodology for providing data analysis tools to road safety researchers to help them in their study of the motorcyclists' behavior during bend negotiation maneuver.

* We successfully apply the proposed framework to a real dataset of sensor data collected during experimentation that involves eight subjects with different profiles.

* We present a multi-sensors architecture to capture the rider's actions during the bend taking maneuver, and we analyze with a certain level of detail the collected data, which makes this manuscript a valuable reference for practitioners interested in such topic of research.
The human riding behavior is a complex concept, and its characterization may lead to a better understanding of the rider's decisions when encountering different situations. This characterization will prevent collisions and design the riding models, which is one of the core algorithms that might make the future of self-riding motorbikes possible. Autonomous vehicles have to interact with other vehicles (even non-autonomous ones), and understanding their driving style can provide valuable information to avoid traffic collisions.

The paper is organized as follows. Section 2 addresses the experimental protocol. Section 3 describes the proposed methodology for PTW rider clustering. The results and discussion of the clustering step are presented in Section 4, and an in-depth analysis of riding control of a selected subject is given in Section 5. The conclusion and perspectives are provided in Section 6.

## 2. Experimental Protocol: Data Collection And Description

The data used in this paper come from a real-life cornering experiment using a heavily instrumented motorcycle (cf. Figure 1a) on the La Ferte Gaucher track (cf. Figure 1b) [17].

<figure><img src="image_5.png" /><figcaption>**Figure 1.** The heavily instrumented motorcycle and the La Ferte Gaucher track</figcaption></figure>

### 2.1. *Experimental Protocol*

The La Ferte Gaucher circuit track in France was chosen for the experimental protocol. The main objective of this experimental protocol is to use sensors to collect data relating to the natural, handlebar, and body riding during bend-taking. Sensors and their description are summarized in Table 1. It is hypothesized that some riders approach a bend using the handlebars more than the body, while others do the opposite. Eight subjects participated in the experiment by undertaking three trials with each of the following instructions; free riding (FR), riding principally using body movements (BR), and riding principally using the handlebar (HR). In BR, the subjects use their body movements (foot, buttocks, knees) as much as possible, as opposed to HR, where they use more of the handlebar while negotiating bends. In FR, the subjects ride in their usual manner. Because of
the velocity, the duration of the trials may differ for the same subject for a given riding instruction. Similarly, this difference is observed among the subjects for all riding instructions (cf. Figure 2).

<table><caption>**Table 1.** Sensors and their description.</caption><thead><tr><th colspan="5">PTW Dynamic Measurements</th></tr><tr><th>N°</th><th>Sensor</th><th>Measurement</th><th>Symbol</th><th>Description</th></tr></thead><tbody><tr><td>1</td><td>Two Hall effect sensors</td><td>Wheel Speed</td><td>v</td><td>Velocity along the longitudinal axis</td></tr><tr><td>2</td><td>SICK DT35 Laser [18]</td><td>Roll angle measurement</td><td>φ</td><td>Two laser sensors are placed on both sides (right and left) of the motorbike to measure the roll angle of the motorcycle</td></tr><tr><td>3</td><td>Magnetic sensor AS5047P of AMS [19]</td><td>Steering angle</td><td>δ</td><td>To acquire the handlebar steering angle</td></tr><tr><td>4</td><td>MTi Xsens [20]</td><td>Three-dimensional (accelerometers, magnetometers, and gyroscopes)</td><td>—</td><td>To acquire inertial movements: longitudinal, lateral, vertical accelerations, and rotational velocities and angles (pitch, yaw, and roll)</td></tr><tr><th colspan="5">Rider Action Measurements</th></tr><tr><td>5</td><td>Strain gauges [21]</td><td>Applied forces on the handlebar</td><td>FrX and FIX</td><td>Strain gauges are placed on the half-handlebars (right and left) of the motorbike to measure the forces applied by the rider on each half-handlebar (right and left)</td></tr><tr><td>6</td><td>Mesurex D2 piezoelectric force button [21]</td><td>Applied forces on the foot-pegs</td><td>Fr<sub>f-p</sub> and F<sub>f-p</sub></td><td>Strain gauges are placed on the (right and left) foot-pegs of the motorbike to measure the forces applied by the rider on each foot-peg</td></tr><tr><td>7</td><td>XSENSOR LX100 and PX100 pressure matrix pads [22]</td><td>Left and right pressure of the buttock</td><td>LbP and RbP</td><td>To acquire pressure forces of the rider's buttocks</td></tr><tr><th colspan="5">Rider Motion Measurements</th></tr><tr><td>8</td><td>Tea Ergo CAPTIV Motion IMU [23]</td><td>Roll angle measurement</td><td>LbRa</td><td>To measure a lower-body roll angle</td></tr><tr><td>9</td><td>Tea Ergo CAPTIV Motion IMU [23]</td><td>Roll angle measurement.</td><td>HbRa</td><td>To measure a higher-body roll angle.</td></tr><tr><td>10</td><td>Tea Ergo CAPTIV Motion IMU [23]</td><td>Roll angle measurement</td><td>— — —</td><td>To measure a head roll angle.</td></tr><tr><th colspan="5">Context Information</th></tr><tr><td>11</td><td>RTK-GPS Septentrio Altus APS3G [24]</td><td>Latitude and longitude positions</td><td>GPS</td><td>To acquire a precise real-time kinematic positioning of the motorcycle.</td></tr><tr><td>12</td><td>Video camera</td><td>Context videos</td><td>—</td><td>Action camera on the top case, looking to the front (field of view including the back of the rider).</td></tr></tbody></table>

<figure><img src="image_6.png" /><figcaption>**Figure 2.** Three riding trials by eight subjects with each riding instruction: FR, HR, and BR.</figcaption></figure>
## 2.2. Data Collection

To collect data related to motorcycle riding, several sensors have been placed at strategic points on the motorcycle. As the data come from several sensors, an on-board computer is used for data centralization. For each riding instruction, the data relating to the motorcycle's dynamics, rider's actions, and their interactions are collected for all bends of the track. However, in this paper, we only use the data relative to the bend highlighted in red (cf. Figure 1b). Information about the eight subjects in the experiment is provided in Table 2. LA indicates the age of licenses, and km is the number of kilometers traveled on PTWs during the year preceding the experiment. The preference and instruction order designate the marks (out of ten) given to HR, the preference between BR and HR, and the instruction order.

<table><caption>**Table 2.** Subjects' declarative data.</caption><thead><tr><td>Subjects</td><td>LA</td><td>km</td><td>Preference</td><td>Instruction Order</td></tr></thead><tbody><tr><td>S1</td><td>2</td><td>8000</td><td>handlebar</td><td>FR, HR, BR</td></tr><tr><td>S2</td><td>7</td><td>3000</td><td>handlebar</td><td>FR, BR, HR</td></tr><tr><td>S3</td><td>0</td><td>0</td><td>body</td><td>FR, HR, BR</td></tr><tr><td>S4</td><td>11</td><td>0</td><td>body</td><td>FR, HR, BR</td></tr><tr><td>S5</td><td>2</td><td>20,000</td><td>handlebar</td><td>FR, BR, HR</td></tr><tr><td>S6</td><td>16</td><td>0</td><td>handlebar</td><td>FR, HR, BR</td></tr><tr><td>S7</td><td>18</td><td>2000</td><td>handlebar</td><td>FR, BR, HR</td></tr><tr><td>S8</td><td>12</td><td>3000</td><td>handlebar</td><td>FR, HR, BR</td></tr></tbody></table>

Among the multiple signals collected on the bend of interest, we only focus on the measurements of the handlebar steering angle ($\delta$), velocity ($v$), and GPS positions (longitude and latitude). In addition, derived signals such as the track curvature ($c$), normal acceleration ($an$), and jerk ($J$) are calculated using the velocity and GPS positions, see Figure 3.

<figure><img src="image_7.png" /><figcaption>**Figure 3.** Data collection process.</figcaption></figure>
## 2.3. Data Representation

The final dataset that we obtain at the end of the experimentation is composed of $\delta, v, an, J, C$. For this dataset, we use the following notations,

$$ \mathcal{D} = \langle \mathcal{D}^l \rangle_{l \in \{FR, HR, BR\}}, $$

with $\mathcal{D}^l = (x_k)_{k=1,...,N}$, where $N$ is the total number of sequences of riding instruction $l$ in $\{FR, HR, BR\}$. Each sequence $x_k = (x_{k,m})_{m=1,...,M}$ is a set of values observed over a period of time $(t_m)_{m=1,...,M}$, M (an integer between 82 and 124) is the total number of observations, and $x_{k,m}$ is given by

$$ x_{k,m} = \{\delta_{k,m}, v_{k,m}, an_{k,m}, J_{k,m}, C_{k,m}\}, $$

where $\delta_{k,m}, v_{k,m}, an_{k,m}, J_{k,m}$ and $C_{k,m}$ are the handlebar steering angle in degrees; velocity in $m/s$; normal acceleration in $m/s^2$; jerk, which is the third derivative of velocity in $m/s^3$; and curvature in $m^{-1}$ measured at time $t_m$ for sequence $k$, respectively.

# 3. Clustering Methodology

In this section, the clustering methodology for 8 subjects with respect to the parameters $\delta, v, an, J, C$ is described. A GPS-based approach for the behavior of PTW riders in bends [25] is used to select the curvature parameter, while the remaining parameters are chosen based on [26].

Three steps characterize our proposed methodology (see Figure 4): (1) Dimensionality reduction to extract a reduced number of descriptors, (2) Anderson-Darling test for the homogeneity study of the three trials, and (3) clustering of subjects.

<figure><img src="image_8.png" /><figcaption>Figure 4. Overview of the clustering methodology.</figcaption></figure>

## 3.1. Dimensionality Reduction Using the Piecewise Aggregate Approximation (PAA) Algorithm

The hierarchical agglomerative clustering (HAC) algorithm that we use in this paper is parameterized with metrics measuring the similarity between two samples. In this situation, as mentioned in [27], the most promising solutions involve first performing a dimensionality reduction on the data. In [27], the authors theoretically and empirically show that the PAA algorithm is superior to other dimensionality reduction techniques such as the singular value decomposition (SVD), discrete Fourier transform (DFT), and discrete wavelet transform (DWT) algorithms. Then, we use the PAA algorithm on each z-normalized $x_{k,m}$ as suggested in [27]. The main idea of the PAA algorithm is to reduce the input z-normalized $x_{k,m}$ dimensionality by splitting it into equally
sized segments (PAA-size) and averaging the values of points in each segment [28]. In our case, the PAA-size parameter is set to 10 so that we obtain approximately 10 samples in each segment per second, which coincides with the frequency of our data (10 Hz). After applying the PAA algorithm to each riding variable of a given set $x_{k,m}$, we obtain the reduced data $y_{k,l}$ with $l = 1, ..., L$ with $L = 10$. An example of the PAA algorithm applied on our data is given in Figure 5.

<figure><img src="image_9.png" /><figcaption>Figure 5. An example of the PAA algorithm applied on $x_{k,m}$.</figcaption></figure>

## 3.2. Homogeneity Study for the Three Attempts of Each Instruction Using the Anderson-Darling Test

Because the subjects were not used to riding the instrumented motorcycle, $\delta$, $v$, $an$, $J$, and $C$ may differ from one trial to another. Statistical tests such as the Kolmogorov-Smirnov [29] or Anderson-Darling [30] tests can be used to check this assumption. However, for small samples, the Anderson-Darling test is more powerful than the Kolmogorov-Smirnov test [31,32]. Thus, we use Anderson-Darling (AD) test, which states under its null hypothesis that the samples are drawn from the same distribution.

In our case, the AD test is performed on $y_{k,m}$ obtained after performing the PAA dimensionality reduction. Therefore, for each instruction and a given riding variable, we can test whether there is a significant difference among the distributions of the three trials. In the Anderson-Darling test, we reject the null hypothesis if the probability value ($p$-value) is less than the reference threshold $\alpha$, which is usually fixed at 5%. In other words, the $p$-value is the greatest threshold for which the null hypothesis is accepted.

Table 3 provides the Anderson-Darling test $p$-value for all riding instructions and variables of 8 subjects. In this table, the minimum value, which is shown in bold, is 0.24, which is greater than $\alpha$ (0.05). Therefore, for all riding instructions and variables, we cannot reject the null hypothesis. Thus, the existing difference among the three trials is not significant for a given subject in all riding instructions and variables.

<table><caption>Table 3. Anderson-Darling test $p$-value for all riding instructions and variables.</caption><thead><tr><th rowspan="2">Subjects</th><th colspan="3">δ</th><th colspan="3">v</th><th colspan="3">an</th><th colspan="3">J</th><th colspan="3">C</th></tr><tr><th>FR</th><th>HR</th><th>BR</th><th>FR</th><th>HR</th><th>BR</th><th>FR</th><th>HR</th><th>BR</th><th>FR</th><th>HR</th><th>BR</th><th>FR</th><th>HR</th><th>BR</th></tr></thead><tbody><tr><td>S1</td><td>0.69</td><td>0.51</td><td>0.96</td><td>0.99</td><td>0.99</td><td>0.91</td><td>0.99</td><td>0.99</td><td>0.99</td><td>0.70</td><td>0.30</td><td>0.74</td><td>0.99</td><td>0.99</td><td>0.99</td></tr><tr><td>S2</td><td>0.97</td><td>0.90</td><td>0.95</td><td>0.99</td><td>0.99</td><td>0.99</td><td>0.99</td><td>1</td><td>1</td><td><b>0.24</b></td><td>0.87</td><td>0.99</td><td>1</td><td>1</td><td>0.99</td></tr><tr><td>S3</td><td>0.96</td><td>0.66</td><td>0.70</td><td>0.99</td><td>0.99</td><td>0.99</td><td>1</td><td>0.99</td><td>0.99</td><td>0.84</td><td>0.90</td><td>0.92</td><td>1</td><td>0.99</td><td>0.99</td></tr><tr><td>S4</td><td>0.99</td><td>0.98</td><td>0.99</td><td>0.99</td><td>0.99</td><td>0.91</td><td>0.99</td><td>1</td><td>0.99</td><td>0.85</td><td>0.94</td><td>0.99</td><td>0.99</td><td>0.99</td><td>0.99</td></tr><tr><td>S5</td><td>0.74</td><td>0.92</td><td>0.99</td><td>0.90</td><td>0.99</td><td>0.99</td><td>1</td><td>0.99</td><td>0.99</td><td>0.76</td><td>0.92</td><td>0.77</td><td>1</td><td>0.99</td><td>1</td></tr><tr><td>S6</td><td>0.55</td><td>0.74</td><td>0.81</td><td>0.88</td><td>0.99</td><td>0.99</td><td>0.99</td><td>1</td><td>0.99</td><td>0.99</td><td>0.64</td><td>0.98</td><td>0.99</td><td>0.99</td><td>1</td></tr><tr><td>S7</td><td>0.99</td><td>0.97</td><td>0.92</td><td>0.99</td><td>0.99</td><td>0.99</td><td>1</td><td>0.99</td><td>0.99</td><td>0.96</td><td>0.85</td><td>0.88</td><td>0.99</td><td>0.99</td><td>0.99</td></tr><tr><td>S8</td><td>0.95</td><td>0.98</td><td>0.90</td><td>0.99</td><td>0.94</td><td>0.99</td><td>1</td><td>0.99</td><td>0.99</td><td>0.48</td><td>0.26</td><td>0.99</td><td>0.99</td><td>0.99</td><td>0.93</td></tr></tbody></table>

In conclusion, the AD test results enable us to reduce the three trials in each riding instruction by taking their average. Thus, for a given riding instruction, we can characterize each subject by an average behavior during bend taking, as shown in Figure 6. For the eight subjects and three riding instructions, we obtain a dataset in a matrix of size 24 × 50. The number of rows 24 is obtained by multiplying the average behavior of the 8 subjects by the number of riding instructions, and the number of columns 50 is obtained by multiplying $L = 10$ the window size chosen in the PAA algorithm by the number of features which is equal to 5.
<figure><img src="image_10.png" /><figcaption>**Figure 6.** Average of three trials of subject S1 in the FR instruction.</figcaption></figure>

### 3.3. Clustering of Subjects Using Hierarchical Agglomerative Clustering (HAC)

The final step of the methodology consists of using hierarchical agglomerative clustering (HAC) [33] on the data matrix $Y_{ij}$ composed by the average behaviors of the subjects during the three riding instructions. To use this algorithm, two methods must be implemented to calculate the distance between two observations and the distance between two clusters (linkage criterion). To calculate the distance between two observations $Y_{ij}$ and $Y_{mj}$ with fixed $i, m$, where $i \neq m$, we use Pearson correlation coefficient [34] and Dynamic Time Warping (DTW) [35] as metrics. Following are the main steps of HAC applied in our study (see Algorithm 1):

**Algorithm 1** The hierarchical agglomerative clustering (HAC) algorithm applied on the $Y_{ij}$ data matrix.

**Require:** $Y_{ij}$ data matrix.

1. Calculate the similarity : using the Pearson correlation and DTW metrics between each row pair of the data matrix $Y_{ij}$.

2. Compute the proximity matrix.

3. Merge the closest clusters according to the proximity matrix.

4. Update the proximity between new and original clusters.

5. Repeat steps 3 and 4 until a single cluster is obtained.

The distance between two clusters is measured by the single linkage criterion, which is the minimum distance among cluster data points.

## 4. Results and Discussion

In this section, we present the results obtained from the HAC algorithm. Choosing the number of clusters is the main difficulty for any clustering algorithm and has always been a subject of discussion, particularly with respect to unsupervised automatic clustering methods. Several criteria have been proposed in the literature, such as the Bayesian information criterion (BIC) [36], elbow rule, and gap statistic [37]. However, in our case, the number of clusters is chosen in relation to the experiment. If each subject respects the driving instructions while being consistent with their statements of riding preference (HR or BR), the ideal number of clusters is 2.

Figure 7 shows the dendrogram of the HAC results using the DTW and Pearson correlation metrics. The first analysis of the dendrograms shows that a clustering with 2 clusters is not a good idea. Indeed, the analysis will compare subject S3 (FR, HR, and BR) against the remaining subjects. Thus, we cut the dendrograms to have 3 clusters, as listed in Table 4. This table shows that we obtain the same results for both the DTW and Pearson correlation metrics.
<figure><img src="image_11.png" /><figcaption>**Figure 7.** Clustering results of the HAC algorithm with the Dynamic Time Warping (DTW) and Pearson correlation metrics.</figcaption></figure>

<table><caption>**Table 4.** Clustering results using the Dynamic Time Warping (DTW) and Pearson correlation metrics.</caption><thead><tr><th rowspan="2">Clusters</th><th colspan="2">Metrics</th></tr><tr><th>DTW</th><th>Pearson Correlation</th></tr></thead><tbody><tr><td>Cluster 1 (Red)</td><td>S3 (FR, HR, BR)</td><td>S3 (FR, HR, BR)</td></tr><tr><td rowspan="4">Cluster 2 (Green)</td><td>S1 (FR, HR, BR)</td><td>S1 (FR, HR, BR)</td></tr><tr><td>S2 (FR, HR, BR)</td><td>S2 (FR, HR, BR)</td></tr><tr><td><strong>S4 (FR, BR)</strong></td><td><strong>S4 (FR, BR)</strong></td></tr><tr><td>S8 (FR, HR, BR)</td><td>S8 (FR, HR, BR)</td></tr><tr><td rowspan="4">Cluster 3 (Blue)</td><td><strong>S4 (HR)</strong></td><td><strong>S4 (HR)</strong></td></tr><tr><td>S5 (FR, HR, BR)</td><td>S5 (FR, HR, BR)</td></tr><tr><td>S6 (FR, HR, BR)</td><td>S6 (FR, HR, BR)</td></tr><tr><td>S7 (FR, HR, BR)</td><td>S7 (FR, HR, BR)</td></tr></tbody></table>

## 4.1. Clusters Analysis

In this paragraph, we conduct cross-referencing for the analysis results of riding behavior extracted from the data mining (cf. Table 4) and stated riding preference (handlebar steering or body movements)
from Table 2. This cross-analysis enables us to determine the extent to which the subjects comply with the riding instructions.

*   **Cluster 1 analysis:** This cluster is essentially composed of three riding instructions of subject S3. An in-depth analysis of the dendrograms shows that in this cluster, the similarity is more observed between HR and BR instructions. Therefore, subject S3 has difficulty in differentiating between HR and BR instructions. Thus, for this subject, contradictions can be observed in terms of his riding preference (BR) and compliance with the HR and BR instructions.

*   **Cluster 2 analysis:** Subjects S2, S4, and S8 have their riding instructions in this cluster, whereas subject S4 has only FR and BR instructions. Therefore, subjects S2, S4, and S8 did not respect the HR and BR riding instructions. Subject S4, who declared a preference for BR instruction, is consistent with his declaration of preference while respecting the riding instructions.

*   **Cluster 3 analysis:** For the same reasons mentioned in the analysis of Cluster 2, subjects S5, S6, and S7 did not respect the HR and BR riding instructions. Subject S4, who has his HR instruction in this cluster, appears to respect this riding instruction.

In this cross-referencing, we detect an atypical rider, namely, subject S3, who behaves very differently from the other riders. This behavior is justified by the fact that this subject is obtaining a motorcycle license (LA is 0). Subject S4 appears to be the only subject to comply with the HR and BR riding instructions. Based on this hypothesis, we say that cluster 2 is composed of subjects using the body more than the handlebars, while cluster 3 is composed of subjects using the handlebars more than the body.

## 4.2. Interpretation of the Clustering Results

The riders who likely have a preference for handlebar riding exhibit significant variability of the handlebar steering angle. Similarly, riders with less variability in the handlebar steering angle are potential candidates to prefer riding with body control. To measure the variability of the handlebar steering angle, we use the standard deviation, which is a good and simple indicator. We use the average signal of three trials of our raw data (cf. Figure 2) to calculate the standard deviation. Figure 8a shows that the variability of the handlebar steering angle contributes the most to distinguishing the two main clusters (2 and 3). To verify this observation, we perform the AD-test between the standard deviation distributions of clusters 2 and 3. For the handlebar steering angle and curvature variables, we reject the null hypothesis at the threshold $\alpha = 5\%$. Thus, the existing difference between clusters 2 and 3 is significant for the variability of the handlebar steering angle and curvature variables, as depicted in Table 5.

The elements in cluster 2 have low variability of the handlebar steering angle compared those in cluster 3. Thus, cluster 2 gathers the elements riding with the body, while those of cluster 3 ride with the handlebars. Therefore, subjects S1, S2, and S8 are in contradiction with their riding preference, while subjects S5, S6, and S7 comply with their statements. This hypothesis does not contradict the behavior of subject S4, who appeared to be the only one to comply with the HR and BR instructions.

<table><caption>**Table 5.** *p*-value of the AD test between clusters 2 and 3.</caption><thead><tr><td>Riding Variables</td><td>p-Value</td></tr></thead><tbody><tr><td>Handlebar steering angle</td><td><strong>0.001</strong></td></tr><tr><td>Velocity</td><td>0.25</td></tr><tr><td>Normal acceleration</td><td>0.25</td></tr><tr><td>Jerk</td><td>0.25</td></tr><tr><td>Curvature</td><td><strong>0.001</strong></td></tr></tbody></table>
<figure><img src="image_12.png" /><figcaption>**Figure 8.** Distribution of the standard deviation by cluster for the five riding variables: (a) Handlebar steering angle, (b) Velocity, (c) Normal acceleration, (d) Jerk, and (e) Curvature.</figcaption></figure>

Based on the results in this section, we conclude that our initial hypothesis has been verified: the driving profile resulting from the interaction between the rider-vehicles-infrastructure can be categorized into classes. These classes represent different riding strategies, namely, handlebar riding and body riding. These different strategies have been identified using discriminative behavioral markers, which were calculated based on data obtained from physical sensors such as the handlebar steering angle and velocity as well as calculated features such as the normal acceleration, jerk, and curvature. Our explanation is that through this characterization of the driving profile, we can capture the sensorimotor processes used by the central nervous system of a rider to control the vehicle: visual anticipation, vestibular perception of the gravito-inertial forces, management of the handlebars, and postural adaptations required to control the vehicle tilt. These sensorimotor processes may differ from one rider to another. This difference is mainly related to her/his riding experience. Nevertheless, at a certain level, there is a common behavior. Another issue can arise from this data analysis: Is there a difference between methods of controlling the motorcycle (steering handlebar vs. rider's lean) in the bend-taking maneuver? In the following section, we will highlight these differences by performing an in-depth analysis of subject S4 by using rider action measurements made by pressure sensors.

## 5. In-Depth Analysis of the Behavior of Subject S4

In the case of PTWs, the driver is heavily involved in the driving task. The driver is continuously trying to ensure the stability of his/her vehicle to counter the gyroscopic effect, which plays an essential role in motorcycle driving, particularly during a turn [38]. For cornering, we note a significant variability in the real practices of drivers. Taking a curve requires the PTW driver to tilt the vehicle to compensate for the inertia caused by a change in direction. To perform this action, it is necessary to produce a tilt on the roll axis towards the side that one wants to approach during the change of direction. This change in direction may be initiated depending on the practices of each driver. The methods to control the motorcycle during a turn can be categorized into two main categories. In the first category, the rider applies a force on the handlebar, which results in a moment around
the steering axis, i.e., steering torque. In the second category, the rider applies a moment directly in the roll direction by basically changing his/her center of gravity in the lateral direction. This shift in body weight results in action and reaction forces in the vertical direction on the foot-pegs and saddle. The present work aims to analyze and compare different motorcycle riders in operating the motorcycle and especially how they initiate a curve in terms of actions measured by pressure sensors. The previous data analysis was performed based on the dynamical behavior of the motorcycle. To perform this in-depth analysis, each maneuver is manually labeled using proprietary software, which is composed of two synchronized applications: one can browse the dynamical data, and the other browses the video recording. The cornering maneuver is divided into four phases (see Figure 9):

* Phase 0 corresponds to when the rider travels in a straight line before entering the curve.

* Phase 1 corresponds to the time when the rider initiates the curve and reaches the middle of the curve. In the following figures, the time interval corresponding to phase 1 is represented by a shaded area.

* Phase 2 corresponds to the elapsed time from the moment that the rider reaches the middle of the curve and that when he/she exits the turn.

* Phase 3 corresponds to when the rider is traveling in a straight line after exiting the curve.

<figure><img src="image_13.png" /><figcaption>**Figure 9.** Average of three trials of subject S4 in the FR instruction. The curve maneuver is divided into four phases: Phase 0, Phase 1, Phase 2, and Phase 3.</figcaption></figure>

Before beginning our in-depth study, we will present the conventions adopted in this work regarding the orientation of roll angles, steering angles, and forces applied on the handlebar steering, see Figure 10a,b.

Rider “S4” is used to carry out an in-depth analysis of the turn maneuver with three instructions, namely, free riding (FR), riding principally using body movements (BR), and riding principally using the handlebar (HR). The previous data analysis study shows that this rider complied with the given instructions. We also find that the FR instruction complies with the BR instruction, which is coherent with the rider’s declared preference. This in-depth analysis of how rider S4 operates a motorcycle based on the average behavior made in the three trials of subject S4 with three instructions: HR, BR, and FR. Therefore, in this in-depth analysis, we expect to find that this rider will initiate the curve through a torque on the handlebars in the case of HR instruction. We do not have a sensor that measures the force-torque applied by the rider on the handlebars; instead, we have sensors that measure the forces applied by the rider on the right handlebar and left handlebar halves, knowing that the torque is the result of the application of the latter. In free riding and body riding cases, the rider will initiate the turn by the movement of his/her body by either moving her/his center of gravity or pressing on the footrest. The next figure shows the handlebar steering and roll angles of the average behavior of subject S4 during the riding instructions, see Figure 11a-c.
<figure><img src="image_14.png" /><figcaption>Figure 10. (a) Roll ($\phi$) and (b) steering ($\delta$) angles and forces ($fl_x, fr_x$) applied on the handlebar steering conventions.</figcaption></figure>

<figure><img src="image_15.png" /><figcaption>Figure 11. Handlebar steering angle $\delta$ (-) and roll angle $\phi$ (-) expressed in degrees (°) of subject S4 for the FR, BR, and HR instructions.</figcaption></figure>
In phase 0, the roll angle is approximately equal to 0, which implies that the motorcycle is upward. When phase 1 begins, the motorcycle starts to roll in the direction of the corner. The roll angle slowly increases until it reaches its maximum, which corresponds to the middle of the curve; this corresponds to when phase 1 ends and phase 2 begins. From this moment, the centrifugal force becomes too large and forces the motorcycle to roll back until the motorcycle returns to the upward position. In the handlebar riding instruction case (see Figure 12c), we notice that the steering angle starts decreasing before the beginning of phase 1, which corresponds to the curve initiation time, not the other instructions; the steering angle remains constant. In that case, we can hypothesize that by applying a steering torque, the rider initiates the curve. The rider creates this torque when he/she applies forces on the handlebar along the longitudinal axis, namely, $Fr_x$ and $Fl_x$. As depicted in Figure 12a-c, only in the HR case is $Fr_x$ significantly decreasing at the beginning of phase 1, which corresponds to the curve initiation.

<figure><img src="image_16.png" /><figcaption>**Figure 12.** Applied forces on the left (Flx) and right (Frx) handlebars along the longitudinal axis expressed in newtons (N), and the handlebar steering angle (δ) expressed in newtons (°) in the FR, BR, and HR instructions.</figcaption></figure>

Note that this turn is a right turn. Therefore, if the driver tries to initiate the corner with his body and not through the handlebars, he/she will either move his center of gravity or press on the right footrest. We will focus on the actions of the driver immediately before the start of phase 1. As shown in Figure 13a-c, there is a change in value of the force applied by the driver on the right footrest
during the FR case. The change is approximately 274 N, which corresponds to t = 2 s. From this instant, this force continuously increases until the end of phase 1 and the beginning of phase 2, where it reaches the value of 286 N. Meanwhile, there is no significant variation in the force applied by the driver on the left footrest. We notice the same behavior on the force applied by the right buttock, which means that simultaneously, the rider is trying to move his center of gravity and pressing the right footrest. This observation can be explained by the fact that the turn in this study is a right turn.

<figure><img src="image_17.png" /><figcaption>**Figure 13.** Left and right foot-pegs ($Fl_{f-p}$, $Fr_{f-p}$), buttocks forces ($FlB$, $FrB$) expressed in newtons (N) with the roll angle ($\phi$) expressed in degrees (°) in the FR instruction.</figcaption></figure>

In the BR case, as shown in Figure 14a-c, we notice a change in value of the force applied by the driver on the right footrest at t = 3 s before starting phase 1, which corresponds to the instant of the start of the turn. From that moment, this force continuously increases until the end of phase 1 and the beginning of phase 2. In addition, there is no large variation in the force applied by the driver to the left footrest. We also notice the same behavior on the force applied by the right buttock, which implies that simultaneously, the rider is trying to move his/her center of gravity and pressing the right footrest.
<figure><img src="image_18.png" /><figcaption>**Figure 14.** Left and right foot-pegs ($Fl_{f-p}$, $Fr_{f-p}$), buttocks forces ($FlB$, $FrB$) expressed in newtons (N) with the roll angle ($\phi$) expressed in degrees (°) in the BR instruction.</figcaption></figure>

In the HR case, as shown in Figure 15a-c, we observe no change in value of the force applied by the driver on the right footrest, and it is even kept constant, and during phase 1, no tendency emerges in its evolution. A similar trend is noted for the force applied by the right buttock evolution. Thus, there is no voluntary action on the driver's part in this condition.

Thus, we highlighted the relationship between the behavioral representation of the driving practices characterized by sensors and the actions of the rider.
<figure><img src="image_19.png" /><figcaption>**Figure 15.** Left and right foot-pegs ($Fl_{f-p}$, $Fr_{f-p}$), buttocks forces ($FIB$, $FrB$) expressed in newtons (N) with the roll angle ($\phi$) expressed in degrees (°) in the HR instruction.</figcaption></figure>

**6. Conclusions**

In this article, the objectives were to use data collected in a real-life experiment to analyze the behavior of several riders during bend taking, identify riders who respect handlebar riding and body riding instructions, and verify those who are consistent with the preferences declared in a questionnaire. We have developed a data mining approach based on feature extraction (using a dimensionality reduction technique) and automatic clustering. Through the standard deviation, the variability of the handlebar is a discriminating feature for rider clustering. The partitioning of provided riders essentially results in two groups according to whether high or low demands are placed on the handlebar during bend riding. Subjects with atypical behaviors were also highlighted. With three riding instructions provided to the riders, we identified the subjects who respect the riding instruction and those whose riding is consistent with the questionnaire. Finally, we initiated an in-depth analysis of bend-taking practices of one subject based on measurements from his inputs (pressures on the handlebar, tank, and footrests). The preliminary results presented in this study are promising. As this work is a part of ongoing research work, future trials should assess this framework methodology's effectiveness. In particular, additional experiments, including a higher number of subjects, should be considered.

This work calls for several perspectives. The bend-taking analysis in this work considers the curve in its entirety and depends on a parameterization of the signals regarding physical knowledge. A possible alternative is to formalize the problem as automatic segmentation of multivariate time series [39,40]. This segmentation will enable a differentiated and in-depth analysis for the entrance,
middle, and end of the bend. Finally, measurements from sensors installed on the rider are included to study the visual anticipation, postural aspects, and their relationship with trajectory control.

**Author Contributions:** A.B. and S.E. designed the main idea of this study; A.B., L.O., and S.E. supervised this study; M.D., A.B., L.O., and S.E. analyzed and discussed the experimental results; M.D. wrote the original draft preparation; A.B., L.O., and S.E. wrote and reviewed the manuscript. All authors have read and agreed to the published version of the manuscript.

**Funding:** This work was supported by the VIROLO++ project funded by Agence Nationale de la Recherche (ANR-15-CE22-0008-02).

**Acknowledgments:** The authors would like to thank Le Centre national de formation à la sécurité routière-Fontainebleau for their contribution to the experimentation. We would like to thank SATIE/MOSS team for their contribution to the instrumentation. The authors would also like to express their deepest gratitude to the reviewers for their constructive comments and time spent to review this paper.

**Conflicts of Interest:** The authors declare no conflicts of interest.

## References

1. 2018 Road Safety Annual Report. 2019. Available online: https://www.onisr.securite-routiere.gouv.fr/en/road-safety-performance/annual-road-safety-reports/2018-road-safety-annual-report (accessed on 22 November 2020).

2. Popov, A.; Rowell, S.; Meijaard, J. A review on motorcycle and rider modelling for steering control. *Veh. Syst. Dyn.* **2010**, *48*, 775–792.

3. Kooijman, J.; Schwab, A. A review on bicycle and motorcycle rider control with a perspective on handling qualities. *Veh. Syst. Dyn.* **2013**, *51*, 1722–1764.

4. Rice, R.S. Rider Skill Influences on Motorcycle Maneuvering. In *SAE Technical Paper*; SAE International: Detroit, MI, USA, 1978.

5. Bocciolone, M.; Cheli, F.; Leo, E.; Pezzola, M. Experimental identification of kinematic coupled effects between driver-motorcycle. In Proceedings of the IMAC-XXV: A Conference & Exposition on Structural Dynamics, Orlando, FL, USA, 19-22 February 2007; pp. 19-22.

6. Evertse, M.V.C. Rider Analysis Using a Fully Instrumented Motorcycle. Master's Thesis, Faculty of Aerospace Engineering-Delft University of Technology, Delft, The Netherlands, 2010.

7. Loiseau, P.; Boultifat, C.N.E.; Chevrel, P.; Claveau, F.; Espie, S.; Mars, F. Rider model identification using dynamic neural networks. In Proceedings of the IFAC World Congress, Berlin, Germany, 11-17 July 2020.

8. Loiseau, P.; Boultifat, C.N.E.; Chevrel, P.; Claveau, F.; Espie, S.; Mars, F. Rider model identification: Neural networks and quasi-LPV models. *IET Intell. Transp. Syst.* **2020**, *14*, 1259–1264.

9. Spiegel, B. *The Upper Half of the Motorcycle: On the Unity of Rider and Machine*; Meredith H., Translator; Whitehorse Press: Center Conway, NH, USA, 2010.

10. Kalsoom, R.; Halim, Z. Clustering the driving features based on data streams. In Proceedings of the INMIC, Lahore, Pakistan, 19-20 December 2013; pp. 89-94. [CrossRef]

11. Higgs, B.; Abbas, M. A two-step segmentation algorithm for behavioral clustering of naturalistic driving styles. In Proceedings of the 16th International IEEE Conference on Intelligent Transportation Systems (ITSC 2013), The Hague, The Netherlands, 6-9 October 2013; pp. 857-862.

12. Wu, C.; Sun, C.; Chu, D.; Huang, Z.; Ma, J.; Li, H. Clustering of Several Typical Behavioral Characteristics of Commercial Vehicle Drivers Based on GPS Data Mining: Case Study of Highways in China. *Transp. Res. Rec.* **2016**, *2581*, 154–163.

13. Xu, J.; Luo, X.; Shao, Y.M. Vehicle trajectory at curved sections of two-lane mountain roads: A field study under natural driving conditions. *ETRR* **2018**, *10*, 12.

14. Zardosht, M.; Beauchemin, S.; Bauer, M. Identifying Driver Behavior in Preturning Maneuvers Using In-Vehicle CANbus Signals. *J. Adv. Transp.* **2018**, *2018*, 1-10. [CrossRef]

15. Crundall, E.; Crundall, D.; Stedmon, A. Negotiating Left-Hand and Right-Hand Bends: A Motorcycle Simulator Study to Investigate Experiential and Behaviour Differences Across Rider Groups. *PLoS ONE* **2012**, *7*, e29978. [CrossRef]

16. Lobjois, R.; Mars, F. Effects of motorcycle simulator configurations on steering control and gaze behavior in bends. *J. Exp. Psychol. Appl.* **2019**. [CrossRef]
17. Smaiah, S.; Sadoun, R.; Elouardi, A.; Larnaudie, B.; Bouaziz, S.; Boubezoul, A.; Vincke, B.; Espié, S. A Practical Approach for High Precision Reconstruction of a Motorcycle Trajectory Using a Low-Cost Multi-Sensor System. *Sensors* **2018**, *18*, 2282. [CrossRef]

18. Available online: https://www.sick.com/fr/fr/capteurs-de-distance/capteurs-de-distance-a-moyenne-portee/dx35/dt35-b15251/p/p295353 (accessed on 22 November 2020).

19. Available online: https://ams.com/AS5047P (accessed on 22 November 2020).

20. Available online: https://www.xsens.com/products/mti-100-series (accessed on 22 November 2020).

21. Available online: https://www.mesurex.fr/en/ (accessed on 22 November 2020).

22. Available online: https://www.xsensor.com/body-pressure-sensors (accessed on 22 November 2020).

23. Available online: https://www.teaergo.com/products/tea-captiv-t-sens-motion-imu/ (accessed on 22 November 2020).

24. Available online: https://www.septentrio.com/en/products/gnss-receivers/rover-base-receivers/smart-antennas/aps3g (accessed on 22 November 2020).

25. Lemonakis, P.V.; Eliou, N.E.; Karakasidis, T.; Botzoris, G. A new methodology for approaching motorcycle riders' behavior at curved road sections. *Eur. Transp. Res. Rev.* **2014**, *6*, 303–314.

26. Nayak, A.; Murigendrappa, S.; Suresh, A. Stability enhancement of a powered two wheeler vehicle under curve negotiation. *Veh. Syst. Dyn.* **2016**, *7*, 85–95.

27. Keogh, E.; Chakrabarti, K.; Pazzani, M.; Mehrotra, S. Dimensionality Reduction for Fast Similarity Search in Large Time Series Databases. *Knowl. Inf. Syst.* **2002**, *3*. [CrossRef]

28. Senin, P.; Lin, J.; Wang, X.; Oates, T.; Gandhi, S.; Boedihardjo, A.P.; Chen, C.; Frankenstein, S. GrammarViz 3.0: Interactive Discovery of Variable-Length Time Series Patterns. *ACM Trans. Knowl. Discov. Data* **2018**, *12*, 10:1–10:28. [CrossRef]

29. Massey, F.J.M., Jr. The Kolmogorov-Smirnov Test for Goodness of Fit. *J. Am. Stat. Assoc.* **1951**, *46*, 68–78.

30. Scholz, F.W.; Stephens, M.A. K-Sample Anderson-Darling Tests. *J. Am. Stat. Assoc.* **1987**, *82*, 918–924.

31. Mohd Razali, N.; Yap, B. Power Comparisons of Shapiro-Wilk, Kolmogorov-Smirnov, Lilliefors and Anderson-Darling Tests. *J. Stat. Model. Analytics* **2011**, *2*, 21–33.

32. Engmann, S.; Cousineau, D. Comparing distributions: The two-sample Anderson-Darling test as an alternative to the Kolmogorov-Smirnov test. *J. Appl. Quant. Methods* **2011**, *6*, 1–17.

33. Zhang, Z.; Murtagh, F.; Van Poucke, S.; Lin, S.; Lan, P. Hierarchical cluster analysis in clinical research with heterogeneous study population: Highlighting its visualization with R. *Ann. Transl. Med.* **2017**, *5*, 75. [CrossRef]

34. Freedman, D.; Pisani, R.; Purves, R. *Purves*, 4th ed.; Statistics (international student edition); Pisani, R., Ed.; WW Norton & Company: New York, NY, USA, 2007.

35. Senin, P. Dynamic time warping algorithm review. *Inf. Comput. Sci. Dep. Univ. Hawaii Manoa Honolulu USA* **2008**, *855*, 40.

36. Weakliem, D.L. A Critique of the Bayesian Information Criterion for Model Selection. *Sociol. Methods Res.* **1999**, *27*, 359–397.

37. Zambelli, A. A Data-Driven Approach to Estimating the Number of Clusters in Hierarchical Clustering. *F1000Research* **2016**, *5*. [CrossRef]

38. Cossalter, V.; Doria, A.; Lot, R. Steady turning of two-wheeled vehicles. *Veh. Syst. Dyn.* **1999**, *31*, 157–181.

39. Chamroukhi, F.; Samé, A.; Govaert, G.; Aknin, P. Time series modeling by a regression approach based on a latent process. *Neural Netw.* **2009**, *22*, 593–602. [PubMed]

40. Chamroukhi, F.; Trabelsi, D.; Mohammed, S.; Oukhellou, L.; Amirat, Y. Joint segmentation of multivariate time series with hidden process regression for human activity recognition. *Neurocomputing* **2013**, *120*, 633–644.

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_20.png" /></figure>

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).