<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.0 20120330//EN" "JATS-archivearticle1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink">
  <front>
    <journal-meta />
    <article-meta>
      <title-group>
        <article-title>Deep Learning for Terrain Surface Classification: Vibration-based Approach</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <string-name>Marcos Concon</string-name>
          <email>marcos@respiree.com</email>
          <xref ref-type="aff" rid="aff0">0</xref>
          <xref ref-type="aff" rid="aff2">2</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>W. K. Wong</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
          <xref ref-type="aff" rid="aff2">2</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Filbert H. Juwono</string-name>
          <email>filbert@ieee.org</email>
          <xref ref-type="aff" rid="aff0">0</xref>
          <xref ref-type="aff" rid="aff2">2</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Catur Apriono</string-name>
          <email>catur@eng.ui.ac.id</email>
          <xref ref-type="aff" rid="aff1">1</xref>
          <xref ref-type="aff" rid="aff2">2</xref>
        </contrib>
        <aff id="aff0">
          <label>0</label>
          <institution>Curtin University Malaysia</institution>
          ,
          <addr-line>CDT 250, Miri 98009, Sarawak</addr-line>
          ,
          <country country="MY">Malaysia</country>
        </aff>
        <aff id="aff1">
          <label>1</label>
          <institution>University of Indonesia</institution>
          ,
          <addr-line>Kampus Baru UI, Depok, West Java 16424</addr-line>
          ,
          <country country="ID">Indonesia</country>
        </aff>
        <aff id="aff2">
          <label>2</label>
          <institution>using Inertial Measurement Unit (IMU) sensors. In or-</institution>
        </aff>
      </contrib-group>
      <pub-date>
        <year>2021</year>
      </pub-date>
      <fpage>237</fpage>
      <lpage>243</lpage>
      <abstract>
        <p>As robots become more pervasive in the service sector, control in dynamic environment has become an important element in optimising the deployment of mobile robots. A mobile robot should be knowledgeable not only of the barriers, but also of the surface on which the robot navigates to estimate slippage and adaptive control. We note that various terrains/surfaces have diferent characteristics, which can directly influence the handling, driving, eficiency, and stability of the robot vehicle. Knowledge of the terrain can provide valuable information for establishing efective and secure navigation strategies. We built a mobile robot prototype equipped by Inertial Measurement Unit (IMU) to obtain the terrain data and applied deep learning models to classify the terrain using the data. Three deep learning configurations have been proposed in this paper, i.e. long short-term memory (LSTM), 1D convolutional network (1D CNN), and convolutional neural network-long short-term memory network (CNN-LSTM). The deep learning architectures were trained and evaluated based on the data collected from ifve diferent surfaces. It is shown that the CNN-LSTM performs the best with an F1 score of 98.49%. The other two networks also generalize relatively well with the unseen vibration sequences with F1 scores of 97.47% and 95.98% for the 1D CNN and LSTM, respectively. Finally, we investigate the efect of varying input sequence to find the optimal length, so that we are able to obtain the highest accuracy and generalization of the deep learning networks.</p>
      </abstract>
      <kwd-group>
        <kwd>Classification</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>1. Introduction</title>
      <p>
        Intelligent robotics have seen rapid advancement in their
scope of operations such as in military reconnaissance
in hostile environments [
        <xref ref-type="bibr" rid="ref6">1</xref>
        ], unmanned surveillance for
disaster management [
        <xref ref-type="bibr" rid="ref8">2</xref>
        ], and telemedicine robot used
for examining remote patients [
        <xref ref-type="bibr" rid="ref10">3</xref>
        ], and in factories. It
is necessary for a robot to acquire a clear
understanding of its current environment in order to successfully
manoeuvre and accomplish its planned operation, while
preventing any damage to itself and creating hazards to
others. As service robots have achieved broad adoption
in the above-mentioned industries, precise navigation
and surrounding awareness have become crucial issues
to improve the capacity of the device to deploy. An
significant consideration for the robot’s eficient navigation
is the motion control algorithm, based on the type of
terrain being travelled. Thus, a detailed classification of
the type of terrain is required for the robot to adapt its
nEvelop-O
      </p>
    </sec>
    <sec id="sec-2">
      <title>2. Related</title>
    </sec>
    <sec id="sec-3">
      <title>Works</title>
      <p>
        The problems of adaptive control in mobile robots have
been constantly researched. The challenges presents
various opportunities for researchers to develop methods in
predicting the dynamic changes in the environment. In
which depend on the characteristics of the terrain.
speed of navigation and the parameters of route planning, results to recognize the surrounding environment of the
[
        <xref ref-type="bibr" rid="ref12">4</xref>
        ], the authors investigated the use of kinematics-based to some factors such as environmental noise and the
analytic for wheel slippage calculation. The results were robot’s internal motor noise as described in [13]. A deep
validated using collected data on a mobile platform. Sim- learning approach was applied in [14] where a CNN was
ilar work was found in [
        <xref ref-type="bibr" rid="ref15">5</xref>
        ] where the authors applied developed and trained using the short time Fourier
transrolling resistance torque without using any additional form (STFT) spectograms extracted from the raw terrain
sensors. Rolling resistance torque in multiple terrains audio signals. It was demonstrated that the network was
can be acquired by reaction torque observer. Proposed robust even when the terrain audio signal was corrupted
concept was verified by using a diferential drive mobile with the white Gaussian noise.
robot. In [
        <xref ref-type="bibr" rid="ref18">6</xref>
        ], wheel slips were estimated based on the Haptic-based classification uses ground contact forces
odometric data. The collected data were analyzed using between a legged robot and terrain to describe diferent
two diferent approaches, which were instantaneous esti- terrain properties. Typically, features such as the robot’s
mator and temporal window approach. Results showed stride frequency, peak and average motor torque in a
that temporal window approach yielded a better result. single stride are used to train an SVM classifier [ 15]. In
      </p>
      <p>
        In [
        <xref ref-type="bibr" rid="ref20">7</xref>
        ], researchers presented a solution of using laser- [16], a 1-dimensional CNN and an RNN architecture were
based point cloud generation to detect robot traversal implemented and evaluated when raw force/torque
sigsurface. The researchers explored several terrains includ- nals from a hexapod robot were passed to them. There
ing carpet, coated asphalt, and asphalt. The solution was was a significant improvement of about 15% in
classificahighly precise with high computation cost as it gener- tion accuracy when compared to the SVM method with
ated point clouds which needed to be further processed a Gaussian kernel.
digitally. The authors stated that there were opportuni- The last reaction-based technique is based on the
vibraties to further investigate how a mobile robotic platform tion characteristics of the terrain. It was first suggested
could provide reliable and accurate surface prediction in [17] where the vibration signal was measured using
of the terrain for improving the navigation with prior an accelerometer during the robot’s traversal. In terms
knowledge of the surface. These research works have of performance, SVM has proven to be the best when
shown that there is strong motivation in investigating trained on hand crafted time domain features such as
methods to enable service robots to have perception on skewness, impulse factor, and root mean square (RMS),
the terrain and traversal surface. along with frequency-domain features from the discrete
      </p>
      <p>
        Several sensing methodologies have been developed to Fourier transform (DFT). Experiments using a CNN for
tackle the problem of terrain classification. The method- vibrational wheel slip estimation in ground robotics was
ology is typically categorized into two main groups: vision- carried out in [
        <xref ref-type="bibr" rid="ref4">11</xref>
        ]. The wheel torque, vertical
accelerabased and reaction-based techniques. Traditional vi- tion and degree of pitch were used to train the classifier.
sual feature engineering approaches include the scale- The diference of 10% for classification was obtained
beinvariant feature transform (SIFT) [
        <xref ref-type="bibr" rid="ref21">8</xref>
        ], speeded-up robust fore and after filtering the input data for the CNN which
features (SURF) [
        <xref ref-type="bibr" rid="ref2">9</xref>
        ], and the bag of visual words (BOVW) reinforces the generality of deep learning frameworks
[
        <xref ref-type="bibr" rid="ref3">10</xref>
        ], among many others. These algorithms pass the in extracting meaningful information directly from raw
useful features of the images obtained from light detec- input vibration data.
tion and ranging (lidar) or stereo camera to a classifier to
be trained and classified. In [
        <xref ref-type="bibr" rid="ref4">11</xref>
        ], raw grayscale terrain
images were trained using deep convolutional network 3. Methodology
and the accuracy was 6% less than the support vector
machine (SVM) classifier used jointly with the histogram The mobile robot used in this research is a two wheel
difof gradients (HOG) feature extractor. ferential drive with an attached 6-axis accelerometer that
      </p>
      <p>
        While vision-based approaches are useful because of measures six vibrational terrain signatures. The setup
their high accuracy, they are vulnerable to distortion is shown in Fig. 2. The form is very similar to the
concaused by lighting changes and other factors such as ventional indoor service robots such as robotic vaccuum
the realization of the surface’s physical properties (e.g. cleaner. The vibration characteristics are all dependent
material type and degree of hardness) [
        <xref ref-type="bibr" rid="ref5">12</xref>
        ]. Reaction- on the terrain’s texture/material and robot movement.
based techniques, on the other hand, utilize sensor mea- This study primarily aims to address terrain
classificasurements to obtain either the acoustics, haptics, or the tion by utilising raw time-series vibration data as input
vibration profiles for the classification. Acoustic-based to three implemented deep learning frameworks: LSTM,
classification relies on the use of microphone to record CNN, and a CNN-LSTM architectures. An overview of
the sound of signal generated between the robot and the experiment workflow is given in Fig. 1.
terrain during traversal. Noise removal and smoothing The data set used in this study contains a total of 24000
techniques are necessary in traditional acoustic-based samples distributed evenly from five diferent terrain
classification to achieve satisfactory results. This is due sources. The six features includes the lateral,
longitu(  ,   ,   ,   ,   ,   ) of the traversing robot. The setup
is shown in Fig. 2. Fig. 3 illustrates the five diferent
ifxed windows of 1.5 seconds (75 samples). An overlap
rate of 20% was applied between two consecutive 1.5
vibration signals corresponding to the surface type. The
second segments to conserve the temporal dependencies
vibration samples were collected via I2C using IMU unit
between the time steps in the vibration sequence.
OneMPU-6050 containing both an accelerometer and
gyrohot encoding was then performed to map the diferent
scope integrated in a single chip. The controlled condi- labelled surfaces numerically. Lastly, the vibration data
tions for the wheeled robot are: 50 Hz sampling rate, 1.6
set was split into training, validation, and testing sets to
minutes traversal time per surface, and circular motion
allow the neural networks to generalize with the unseen
of the robot.
sequence,  is the average, and  is the standard
devia(1)
4.
vibration characteristics. These data set partition were
set to be 70%, 15% and 15% for training, validation, and
testing, respectively.
3.1. Implementation
      </p>
      <sec id="sec-3-1">
        <title>LSTM is a type of recurrent neural network (RNN) that</title>
        <p>is typically used for sequence prediction. In
particular, LSTM solves the issues of the disappearing gradient
present in the RNNs while allowing the long-term
temporal dynamics of the series to be exploited. In contrast,</p>
      </sec>
      <sec id="sec-3-2">
        <title>CNNs have been commonly used for 2D problems (e.g.</title>
        <p>image classification task); however, it can be modified
to classify the 1D vibrational problem. The
dimensionality of the convolutional layers is reduced to match the
model’s 1D input.</p>
      </sec>
      <sec id="sec-3-3">
        <title>The CNN-LSTM model leverages the robustness of</title>
      </sec>
      <sec id="sec-3-4">
        <title>CNN in extracting spatial features and LSTM in exploit</title>
        <p>ing the temporal dependencies of the vibration sequence.</p>
      </sec>
      <sec id="sec-3-5">
        <title>In this paper, the time-series vibration is downsampled</title>
        <p>by the 1D CNN to extract the higher level features. This
can be considered as the pre-processing step which
allows the LSTM to interpret the features extracted at each
block of the sequence. The concept is illustrated in Fig.</p>
      </sec>
      <sec id="sec-3-6">
        <title>The three models were built and trained using a Ten</title>
        <p>Layer
LSTM (20 units)
Dropout (25%)
LSTM (70 units)
Dropout (40%)
Dense
Dense
Conv1D (80@6×1)
Dropout (50%)
Conv1D
(128@6×1)
Dropout(50%)
Max pooling
Flatten
Dense
Dense
Conv1D (96@6×1)
Conv1D (48@6×1)
Dropout (30%)
Max pooling
Flatten
LSTM (20 units)
Dropout (20%)
Dense
Dense
sorflow backend with the Keras API. A detailed overview
of the three models was summarized in Table 1. The
hyperband algorithm was used to select the
hyperparameters allowing for the best balance between training time
and accuracy. The learning rate, batch size, and number
of epochs were set at 0.001, 64, and 30, respectively.
Additionally, early stopping regularization was implemented
to avoid overfitting during model training. Further, the CNN-LSTM
Adam optimization algorithm based on the Stochastic
gradient descent was used as the optimizer.</p>
        <p>The implemented CNN-LSTM architecture is shown
in Fig. 5. For both the LSTM and 1D CNN networks,
the data length of a vibration training sample was a flat
vector of 75 time steps. In a stacked LSTM network, the
input sequence to the first LSTM layer returns a shape to retain its temporal representation during the
convoof (timestep, unit) to be passed on to the next layer. The lution process. The time distributed layer expects a 3D
output from the last LSTM layer returns only the unit. input and so the input sequence was reshaped from 75u
For the 1D CNN, the input shape to the network is repre- time steps into 3 subsequences of 25 time steps. The
consented as (timestep, features). In the case for the CNN- volutional layer used the ReLU activation and consisted
LSTM network, a time distributed wrapper is first used of a 6 × 1 kernel that moves across in one dimension
before the LSTM layers to allow the input vibration signal during the convolution operation.</p>
        <p>The dropout layers were then added to tackle overfit- Fig. 7 illustrates the performance of the models on the
ting issues by arbitrarily setting a fraction rate of input 5-class vibration test data set. In the worst case, about
units to zero. The pooling layer was added to reduce the 7% (average) of the wood class was mistakenly classified
spatial size of the output representation into half. Note as tiles across the three models. Overall, it can be seen
that both the dropout and pooling layer allows for faster that the three models exhibited good performance and
training time due to the reduced parameter size. The generalized well with the unseen data. Further, the
CNNlfatten layer was used to transform the input from the LSTM architecture has the best performance with F1
previous layers as input to the LSTM layer where the score of 98.49% (average). The 1D CNN follows at the
temporal characteristics of the vibration sequence were second place with F1 score of 97.49% (average). We note
extracted. Lastly, the fully connected layers with the soft- that the slight improvement of the CNN-LSTM model
max activation function was used to structure the outputs compared to the 1D CNN may suggest that the temporal
of the previous layer for the final classification task. In characteristics of the LSTM is less important than the
this experiment, the categorical cross-entropy loss func- feature generation capability of the CNN-LSTM. Table 4
tion was used to address the 5-class terrain classification summarizes the overall performance of the three models.
problem.</p>
      </sec>
    </sec>
    <sec id="sec-4">
      <title>4. Results</title>
      <sec id="sec-4-1">
        <title>The confusion matrix of the three models are depicted in</title>
        <p>Fig. 6. From the confusion matrix, we can calculate the
F1 score, the precision   , and the recall,   . The F1 score,
which is calculated using   and   , has been commonly
used to analyze the performance of the models. We used
macro-averaging technique to expand these benchmarks
towards multi-class terrain classification. The equations
for the precision, recall, and F1 score, respectively are
given by
  =
  =
 1 =
,
,
  +  
 
 
  +  
2    ,
  +  
where   shows the outcome where the model correctly
classifies the positive class,   is the outcome where the
model incorrectly classifies the positive class,   is the
outcome where the model correctly classifies the negative
class, and   is the outcome where the model incorrectly
classifies the negative class.
(2)
(3)
(4)</p>
        <p>One factor influencing the performance of the models
is the sequence length of the input vibration. To
further validate the performance of the three models, we
analyzed the F1 score with varying segment lengths as
shown in Fig. 8. It can be seen that a longer sequence
length results in a better accuracy with the cost of
performance saturation at a certain length. It can be shown
that the average F1 score rises as the duration of the
vibration series increases from 30 to 60 samples but
decreases afterwards. This may be caused by the lack of the
training data after the segmentation process of the given
length. Therefore, we can consider an optimal sequence
length of 75 (1.5 seconds) for obtaining high accuracy and
generalization of the models. Furthermore, the proposed
CNN-LSTM architecture slightly outperformed the CNN
and LSTM models across the varying segment lengths.</p>
      </sec>
    </sec>
    <sec id="sec-5">
      <title>5. Conclusion and Future Work</title>
      <sec id="sec-5-1">
        <title>In this paper, we have demonstrated the application of</title>
        <p>IMU-based surface classification task. We have compared
three candidates for classifying the IMU data, i.e. LSTM,
1D CNN, and a combination of CNN and LSTM. By
comparing the results, CNN-LSTM provided the best results
(F1 score of 98.49%). However, we can further observe
that the 1D CNN presented favorable results although
slightly lower than the CNN-LSTM. The results suggest
that 1D CNN is able to map the classification better when
compared to the LSTM on standalone basis. CNN and
LSTM works on diferent principle in which the latter is
based on the temporal dynamics of the data. On the other
hand, 1D CNN is based on static convolution, similar to
the 2D counterparts. This implies that there is a clear
static pattern when the IMU data enabling well defined
mapping to their respective classes.</p>
        <p>The results, despite counter intuitive, may prompt
furFigure 6: Confusion matrices of (a) LSTM, (b) CNN, and (c) ther research in the this direction. With the growing
CNN-LSTM on the vibration test dataset. of edge computing and capacity of embedded system,
enabling robots to recognize surface would enable
further applications for indoor or industrial applications.</p>
      </sec>
      <sec id="sec-5-2">
        <title>Reducing the complexity of the machine learning models to further benefit in terms of computation reduction is required. This is possible given that there is a clear static pattern demonstrated from the results using 1D CNN.</title>
      </sec>
    </sec>
    <sec id="sec-6">
      <title>Acknowledgments References</title>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          <article-title>This research was supported by Ministry of Research and Technology/National Agency for Research and Innovation, Republic of Indonesia through Penelitian Dasar Unggulan Perguruan Tinggi (PDUPT) Grant, contract number: NKB-2838/UN2</article-title>
          .RST/HKP.05.00/
          <year>2020</year>
          , year
          <year>2020</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          [9]
          <string-name>
            <surname>Seung-Youn Lee</surname>
            ,
            <given-names>D.</given-names>
          </string-name>
          <string-name>
            <surname>Kwak</surname>
          </string-name>
          ,
          <article-title>A terrain classification method for ugv autonomous navigation based on surf</article-title>
          ,
          <source>in: 2011 8th International Conference on Ubiquitous Robots and Ambient Intelligence (URAI)</source>
          ,
          <year>2011</year>
          , pp.
          <fpage>303</fpage>
          -
          <lpage>306</lpage>
          .
          <source>doi:1 0 . 1 1 0 9 / U R A I . 2 0</source>
          <volume>1 1 . 6 1 4 5 9 8 1 .</volume>
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          [10]
          <string-name>
            <given-names>H.</given-names>
            <surname>Wu</surname>
          </string-name>
          ,
          <string-name>
            <given-names>B.</given-names>
            <surname>Liu</surname>
          </string-name>
          ,
          <string-name>
            <given-names>W.</given-names>
            <surname>Su</surname>
          </string-name>
          ,
          <string-name>
            <given-names>Z.</given-names>
            <surname>Chen</surname>
          </string-name>
          ,
          <string-name>
            <given-names>W.</given-names>
            <surname>Zhang</surname>
          </string-name>
          ,
          <string-name>
            <given-names>X.</given-names>
            <surname>Ren</surname>
          </string-name>
          ,
          <string-name>
            <given-names>J.</given-names>
            <surname>Sun</surname>
          </string-name>
          ,
          <article-title>Optimum pipeline for visual terrain classiifcation using improved bag of visual words and fusion methods</article-title>
          ,
          <source>Journal of Sensors</source>
          <year>2017</year>
          (
          <year>2017</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          [11]
          <string-name>
            <given-names>R.</given-names>
            <surname>González</surname>
          </string-name>
          ,
          <string-name>
            <given-names>K.</given-names>
            <surname>Iagnemma</surname>
          </string-name>
          ,
          <article-title>Deepterramechanics: Terrain classification and slip estimation for ground robots via deep learning</article-title>
          , CoRR abs/
          <year>1806</year>
          .07379 (
          <year>2018</year>
          ). URL: http://arxiv.org/abs/
          <year>1806</year>
          .07379.
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          [12]
          <string-name>
            <given-names>P.</given-names>
            <surname>Roy</surname>
          </string-name>
          ,
          <string-name>
            <given-names>S.</given-names>
            <surname>Ghosh</surname>
          </string-name>
          ,
          <string-name>
            <given-names>S.</given-names>
            <surname>Bhattacharya</surname>
          </string-name>
          , U. Pal,
          <article-title>Efects of degradations on deep neural network architectures</article-title>
          , CoRR abs/
          <year>1807</year>
          .10108 (
          <year>2018</year>
          ). URL: http://arxiv.org/
        </mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          [1]
          <string-name>
            <given-names>J. G.</given-names>
            <surname>Bellingham</surname>
          </string-name>
          ,
          <string-name>
            <given-names>K.</given-names>
            <surname>Rajan</surname>
          </string-name>
          , Robotics in remote abs/
          <year>1807</year>
          .10108.
          <article-title>a r X i v : 1 8 0 7 . 1 0 1 0 8</article-title>
          .
        </mixed-citation>
      </ref>
      <ref id="ref7">
        <mixed-citation>
          <article-title>and hostile environments</article-title>
          ,
          <source>Science</source>
          <volume>318</volume>
          (
          <year>2007</year>
          ) [13]
          <string-name>
            <given-names>J.</given-names>
            <surname>Libby</surname>
          </string-name>
          ,
          <string-name>
            <given-names>A. J.</given-names>
            <surname>Stentz</surname>
          </string-name>
          ,
          <article-title>Using sound to classify vehicle1098-1102. terrain interactions in outdoor environments</article-title>
          , in:
        </mixed-citation>
      </ref>
      <ref id="ref8">
        <mixed-citation>
          [2]
          <string-name>
            <given-names>V.</given-names>
            <surname>Jorge</surname>
          </string-name>
          ,
          <string-name>
            <given-names>R.</given-names>
            <surname>Granada</surname>
          </string-name>
          ,
          <string-name>
            <given-names>R.</given-names>
            <surname>Maidana</surname>
          </string-name>
          ,
          <string-name>
            <given-names>D.</given-names>
            <surname>Jurak</surname>
          </string-name>
          , G. Heck, 2012 IEEE International Conference on Robotics A.
          <string-name>
            <surname>Negreiros</surname>
          </string-name>
          , D. dos
          <string-name>
            <surname>Santos</surname>
            ,
            <given-names>L.</given-names>
          </string-name>
          <string-name>
            <surname>Gonçalves</surname>
            ,
            <given-names>A.</given-names>
          </string-name>
          <string-name>
            <surname>Amory</surname>
          </string-name>
          , and
          <string-name>
            <surname>Automation</surname>
          </string-name>
          ,
          <year>2012</year>
          , pp.
          <fpage>3559</fpage>
          -
          <lpage>3566</lpage>
          . doi:
          <article-title>1 0 . 1 1 0 9 / A survey on unmanned surface vehicles for</article-title>
          disas
          <string-name>
            <surname>- I C R</surname>
          </string-name>
          <article-title>A</article-title>
          .
          <volume>2 0 1 2 . 6 2 2 5 3 5 7 .</volume>
        </mixed-citation>
      </ref>
      <ref id="ref9">
        <mixed-citation>
          <article-title>ter robotics: Main challenges and directions</article-title>
          , Sen- [14]
          <string-name>
            <given-names>A.</given-names>
            <surname>Valada</surname>
          </string-name>
          ,
          <string-name>
            <given-names>L.</given-names>
            <surname>Spinello</surname>
          </string-name>
          , W. Burgard, Deep feature sors
          <volume>19</volume>
          (
          <year>2019</year>
          )
          <article-title>702</article-title>
          . URL: http://dx.doi.org/10.3390/
          <article-title>learning for acoustics-based terrain classification, s19030702</article-title>
          .
          <source>doi:1 0 . 3 3 9 0 / s 1 9</source>
          <volume>0 3 0 7 0 2</volume>
          . in: International Symposium on Robotics Research
        </mixed-citation>
      </ref>
      <ref id="ref10">
        <mixed-citation>
          [3]
          <string-name>
            <given-names>K. K.</given-names>
            <surname>Chung</surname>
          </string-name>
          ,
          <string-name>
            <given-names>K. W.</given-names>
            <surname>Grathwohl</surname>
          </string-name>
          ,
          <string-name>
            <given-names>R. K.</given-names>
            <surname>Poropatich</surname>
          </string-name>
          ,
          <source>(ISRR)</source>
          ,
          <year>2015</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref11">
        <mixed-citation>
          <string-name>
            <given-names>S. E.</given-names>
            <surname>Wolf</surname>
          </string-name>
          ,
          <string-name>
            <given-names>J. B.</given-names>
            <surname>Holcomb</surname>
          </string-name>
          , Robotic telepresence: [15]
          <string-name>
            <given-names>X. A.</given-names>
            <surname>Wu</surname>
          </string-name>
          ,
          <string-name>
            <given-names>T. M.</given-names>
            <surname>Huh</surname>
          </string-name>
          ,
          <string-name>
            <given-names>R.</given-names>
            <surname>Mukherjee</surname>
          </string-name>
          ,
          <string-name>
            <given-names>M.</given-names>
            <surname>Cutkosky</surname>
          </string-name>
          , Inpast, present, and
          <article-title>future, Journal of cardiothoracic tegrated ground reaction force sensing and terrain and vascular anesthesia 21 (</article-title>
          <year>2007</year>
          )
          <fpage>593</fpage>
          -
          <lpage>596</lpage>
          .
          <article-title>classification for small legged robots</article-title>
          , IEEE Robotics
        </mixed-citation>
      </ref>
      <ref id="ref12">
        <mixed-citation>
          [4]
          <string-name>
            <given-names>R.</given-names>
            <surname>Chaichaowarat</surname>
          </string-name>
          , W. Wannasuphoprasit,
          <source>Wheel and Automation Letters</source>
          <volume>1</volume>
          (
          <year>2016</year>
          )
          <fpage>1125</fpage>
          -
          <lpage>1132</lpage>
          . doi:
          <volume>1</volume>
          <fpage>0</fpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref13">
        <mixed-citation>
          <article-title>slip angle estimation of a planar mobile platform, 1 1 0 9</article-title>
          / L R A .
          <volume>2 0 1 6 . 2 5 2 4 0 7 3 .</volume>
        </mixed-citation>
      </ref>
      <ref id="ref14">
        <mixed-citation>
          in: 2019 First International Symposium on In- [16]
          <string-name>
            <given-names>J.</given-names>
            <surname>Bednarek</surname>
          </string-name>
          ,
          <string-name>
            <given-names>M.</given-names>
            <surname>Bednarek</surname>
          </string-name>
          ,
          <string-name>
            <given-names>L.</given-names>
            <surname>Wellhausen</surname>
          </string-name>
          ,
          <string-name>
            <given-names>M.</given-names>
            <surname>Hutstrumentation</surname>
          </string-name>
          , Control,
          <string-name>
            <given-names>Artificial</given-names>
            <surname>Intelligence</surname>
          </string-name>
          , and ter, K. Walas,
          <article-title>What am i touching? learning to Robotics (ICA-SYMP</article-title>
          ),
          <year>2019</year>
          , pp.
          <fpage>163</fpage>
          -
          <lpage>166</lpage>
          . doi:
          <article-title>1 0 . classify terrain via haptic sensing</article-title>
          ,
          <source>in: 2019 In1 1</source>
          <volume>0</volume>
          <fpage>9</fpage>
          <string-name>
            <surname>/ I C A - S Y M P .</surname>
          </string-name>
          <article-title>2 0 1 9 . 8 6 4 6 1 9 8</article-title>
          . ternational Conference on Robotics and Automa-
        </mixed-citation>
      </ref>
      <ref id="ref15">
        <mixed-citation>
          [5]
          <string-name>
            <given-names>S. D. A. P.</given-names>
            <surname>Senadheera</surname>
          </string-name>
          ,
          <string-name>
            <given-names>A. M. H. S.</given-names>
            <surname>Abeykoon</surname>
          </string-name>
          , Sen- tion
          <source>(ICRA)</source>
          ,
          <year>2019</year>
          , pp.
          <fpage>7187</fpage>
          -
          <lpage>7193</lpage>
          .
          <source>doi:1 0 . 1 1</source>
          <volume>0</volume>
          <fpage>9</fpage>
          <string-name>
            <surname>/ I C R</surname>
          </string-name>
          <article-title>A</article-title>
          .
        </mixed-citation>
      </ref>
      <ref id="ref16">
        <mixed-citation>
          <source>sorless terrain estimation for a wheeled mobile 2 0</source>
          <volume>1 9 . 8 7 9 4 4 7 8 .</volume>
        </mixed-citation>
      </ref>
      <ref id="ref17">
        <mixed-citation>
          robot, in: 2017 IEEE International Conference on [17]
          <string-name>
            <given-names>C. A.</given-names>
            <surname>Brooks</surname>
          </string-name>
          ,
          <string-name>
            <given-names>K.</given-names>
            <surname>Iagnemma</surname>
          </string-name>
          ,
          <source>Vibration-based terIndustrial and Information Systems (ICIIS)</source>
          ,
          <year>2017</year>
          ,
          <article-title>rain classification for planetary exploration rovers</article-title>
          , pp.
          <fpage>1</fpage>
          -
          <lpage>6</lpage>
          .
          <source>doi:1 0 . 1 1</source>
          <volume>0</volume>
          <fpage>9</fpage>
          <string-name>
            <surname>/ I C I I N F S</surname>
          </string-name>
          .
          <volume>2 0 1 7 . 8 3 0 0 4 2</volume>
          <fpage>2</fpage>
          .
          <source>IEEE Transactions on Robotics</source>
          <volume>21</volume>
          (
          <year>2005</year>
          )
          <fpage>1185</fpage>
          -
          <lpage>1191</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref18">
        <mixed-citation>
          [6]
          <string-name>
            <given-names>D.</given-names>
            <surname>Masha</surname>
          </string-name>
          , M. Burke, b. Twala, Slip estimation meth-
          <source>doi:1 0 . 1 1 0 9 / T R O . 2 0</source>
          <volume>0 5 . 8 5 5 9 9 4 .</volume>
        </mixed-citation>
      </ref>
      <ref id="ref19">
        <mixed-citation>
          <article-title>ods for proprioceptive terrain classification using tracked mobile robots</article-title>
          , in: International Conference (
          <string-name>
            <surname>PRASA-RobMech</surname>
            <given-names>)</given-names>
          </string-name>
          ,
          <year>2017</year>
          , pp.
          <fpage>150</fpage>
          -
          <lpage>152</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref20">
        <mixed-citation>
          [7]
          <string-name>
            <given-names>S.</given-names>
            <surname>Wilson</surname>
          </string-name>
          ,
          <string-name>
            <given-names>J.</given-names>
            <surname>Potgieter</surname>
          </string-name>
          ,
          <string-name>
            <given-names>K.</given-names>
            <surname>Arif</surname>
          </string-name>
          ,
          <article-title>Floor surface mapping using mobile robot and 2d laser scanner</article-title>
          ,
          <source>in: 2017 24th International Conference on Mechatronics and Machine Vision in Practice (M2VIP)</source>
          ,
          <year>2017</year>
          , pp.
          <fpage>1</fpage>
          -
          <lpage>6</lpage>
          . doi:
          <article-title>1 0 . 1 1 0 9 / M 2 V I P</article-title>
          .
          <volume>2 0 1 7 . 8 2 1 1 5 0 8 .</volume>
        </mixed-citation>
      </ref>
      <ref id="ref21">
        <mixed-citation>
          [8]
          <string-name>
            <given-names>S.</given-names>
            <surname>Zenker</surname>
          </string-name>
          ,
          <string-name>
            <given-names>E. E.</given-names>
            <surname>Aksoy</surname>
          </string-name>
          ,
          <string-name>
            <given-names>D.</given-names>
            <surname>Goldschmidt</surname>
          </string-name>
          ,
          <string-name>
            <given-names>F.</given-names>
            <surname>Wörgötter</surname>
          </string-name>
          ,
          <string-name>
            <given-names>P.</given-names>
            <surname>Manoonpong</surname>
          </string-name>
          ,
          <article-title>Visual terrain classification for selecting energy eficient gaits of a hexapod robot</article-title>
          ,
          <source>in: 2013 IEEE/ASME International Conference on Advanced Intelligent Mechatronics</source>
          ,
          <year>2013</year>
          , pp.
          <fpage>577</fpage>
          -
          <lpage>584</lpage>
          .
          <source>doi:1 0 . 1 1 0 9 / A I M . 2 0</source>
          <volume>1 3 . 6 5 8 4 1 5 4 .</volume>
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>