Creating a robot localization monitor using particle filter and machine learning approaches

Robot localization is a fundamental capability of all mobile robots. Because of uncertainties in acting and sensing, and environmental factors such as people flocking around robots, there is always the risk that a robot loses its localization. Very often behaviors of robots rely on a reliable position estimation. Thus, for dependability of robot systems it is of great interest for the system to know the state of its localization component. In this paper we present an approach that allows a robot to asses if the localization is still correct. The approach assumes that the underlying localization approach is based on a particle filter. We use deep learning to identify temporal patterns in the particles in the case of losing/lost localization. These patterns are then combined with weak classifiers from the particle set and sensor perception for boosted learning of a localization estimator. Through the extraction of features generated by neural networks and its usage for training strong classifiers, the robots localization accuracy can be estimated. The approach is evaluated in a simulated transport robot environment where a degraded localization is provoked by disturbances cased by dynamic obstacles. Results show that it is possible to monitor the robots localization accuracy using convolutional as well as recurrent neural networks. The additional boosting using Adaboost also yields an increase in training accuracy. Thus, this paper directly contributes to the verification of localization performance.


Introduction
Robot localization capability of all mobile intelligent robots, aiming to determine the position and orientation of an object in space, which is also called an objects pose.There is a huge corpus of research that has been 1 Institute of Software Technology, Graz University of Technology, Graz, Austria 2 IncubedIT GmbH, Hart bei Graz, Austria conducted to allow reliable estimation of the current pose using uncertain sensor data and action execution.Typically, these approaches perform some kind of state estimation to determine the current robot pose, using sensors like 2D lasers [1].A typical approach used for localization are particle filters [2].Particle filters are a popular method for representing arbitrary probability distributions and solving state estimation problems.The technique behind particle filters is the Monte Carlo method [3][4][5] which has existed for over five decades.The key idea of particle filters is to spread particles in space which represent the posterior distribution [1].Instead of using a parametric form for representing the distribution, particle filters generate samples based on its own distribution.Figure 1 shows a general overview on the localization with a particle filter.First, an environment must be mapped.This is often done using a simultaneous localization and mapping (SLAM) approach [6], which needs takes the robot movement (twist) and scanning data on the environment into consideration.The result is a map of the environment that can be used in addition with scanner data and the robot's twist to estimate the pose of the robot.The main problem with particle filters is that they are computationally expensive.Particle filters allow the analysis of complex systems which are non-linear and non-Gaussian.The goal is to deal with arbitrary probability distributions and model them correctly [7,8].Since the computational power has increased in the last years and particle filters are a non-parametric approach for solving complex models, they were applied in many different fields like neuroscience [9], biochemical networks [10], signal processing [11], economics [12] and robotics [13].
Because of uncertainties in sensor data and environmental factors such as people flocking around robots there is always the risk that the error of the pose estimation increases significantly or that it diverges completely.Very often navigation and other behaviors of robots rely on a valid pose estimation.Thus, for performance as well as dependability of robot systems it is of great interest for the system to know the state of its localization component.Moreover, if a robot gets delocalized near an object this can result in undesirable behaviors, such as a crash.To avoid this, it is aimed to detect the delocalization before it occurs to be able to perform a predefined error-correcting behavior.Thus undesired behaviors can be prevented.
In this paper we present an approach that allows a robot to asses if the localization is still valid and thus contribute to the topic of localization performance verification.The approach assumes that the underlying localization approach is based on a particle filter.Following the assumption that the particle set representing a probability distribution of the robot's location bears information about the localization state, we propose to use deep learning to identify temporal patterns in the particle set to be able to act in pre-caution before losing localization.Thus, we investigated different deep neural networks (non-recurrent and recurrent) for their ability to learn identifying such patterns.Moreover, we propose to combine these networks with weak classifiers obtained from statistical information in the particle set and the actual robot perception for boosted learning of a more reliable localization monitor.
For the training and evaluation of the approach we use a realistic simulation of an industrial transport robot environment where due to access to ground truth a large amount of labeled data can be generated.Moreover, a simulation allows to provoke a degraded localization easily by disturbing the robot's localization system by randomly adding dynamic obstacles.By disturbing the sensor readings, the localization algorithm cannot fit the sensor data to the prerecorded map and thus might delocalize the robot.
The remainder of the paper is organized as follows.In the next section we discuss briefly related research.In Section 3 we introduce the proposed approach for a localization monitor.An experimental evaluation of the proposed approach is presented in the succeeding section.Finally, in Section 5 we draw some conclusions and provide some ideas for future work.

Related research
In this section we briefly discuss related research concerning robot localization and machine learning.First, problems of current particle filter approaches are presented before discussing recent advances in the topics of robot localization and machine learning.

Dilemma of particle filters and its limitations in robot localization
Localization approach based on particle filters are a good method to overcome common issues in the field of robot localization [1,2].For example they are able to solve the global localization problem as well as the kidnapping problem.Another advantage of particle filters is that they are non-parametric [1] which means that they get an arbitrary probability density function (PDF) without adjusting parameters by hand.However, there also occur some minor problems.The first is that the accuracy and the performance depend on the number of particles used to represent the PDF.To keep this algorithm efficient the number of particles should be small while to keep it accurate the number of particles should be high [14,15].
Another issue is the inflation of the particle set.When sensor data can not be clearly assigned to a pose within the environment, the particle set starts to inflate because the uncertainty about the current pose increases.It is hard to determine if the robot is still localized or if it is already delocalized.The best example for such an inflation are a long straight corridor.When navigating through the corridor the particle filters pose estimation of the robot becomes imprecise.This is due to the issue that for the particle filter every pose within this corridor looks similar.The current laser scan fits fine to the current robots pose but also fits to various poses ahead and behind the robot due to ambiguous measurements.This paper aims to solve this inflation problem by using the proposed method and detect uncertainty increases.This way, an impending delocalization can be detected in advance and the robot can start a predefined error-correcting behaviour.

Recent approaches in robot localization and its accuracy
A publication which was presented in 2012 by Röwekämper et al. evaluates the position accuracy of a mobile robot localization method that is based on particle filtering and laser scan matching [16].For evaluation they used a motion capture system that tracks the position of a robot within its environment with high accuracy [17].The results show, that a high accuracy of +/-0.5cm can be achieved in a static environment.The localization system which was used in their evaluation was a combination of basic state-of-theart approaches.Therefore Monte Carlo Localization, a scan matching procedure [18] and the Kullback-Leibler distance sampling (KLD-sampling) [15] were used.Although the proposed methods yields a high localization accuracy, the combination of particle filter and laser scan matching have the potential the delocalize the robot in a dynamic environment.To overcome this, a similar approach is used in this paper in combination with the presented approach.
A novel approach is proposed by Kallasi et al. who introduced a new method for detecting features to localize and navigate robots [19].They propose two new feature detectors called Fast Adaptive Laser Keypoint Orientation-Invariant (FALKO) and Orthogonal Corner (OC).Those two detectors are an improvement of the Fast Laser Interest Region Transform (FLIRT) [20] approach that can be used to detect high curvature points in laser scan images.While the FLIRT method searches for general features which depend on the viewpoint of the robot, FALKO and OC are designed to detect stable features like corner walls.The difference between the two proposed methods is that FALKO detects features by selecting meaningful neighbours and scoring the cornerness of a feature and OC uses orthogonal alignments to detect important features.
Since particle filters are among the most popular methods for state estimation, a lot of research is done to improve their accuracy and performance.In general the state of a complex system can be estimated correctly by using (infinite) many particles.But with an increasing number of particles the efficiency of the particle filter decreases.Thus, researcher aim at decreasing the number of particles while offering the same performance.One attempt is called adaptive particle filtering as described in [13].
Oliveira et al. present a visual localization approach to estimate a robots pose [21].In contrast to Lidar systems, visual pose estimation is more cost-effective but the accuracy is usually inferior to Lidar systems.To overcome this issue, Oliveira et al. present a visionbased localization approach that uses visual odometry and topological localization which is used to train deep neural networks.Using their approach, they manage to get up to 10 times smaller localization errors than with traditional vision-based localization approaches.
Havangi recently proposed a localization method, using particle swarm optimization (PSO) estimators to overcome the localization weaknesses which occur when using particle filters [22].In the proposed method, dynamic optimization is converted to find a pose estimate that suits best to the robots location.The advantage compared to PF is, that no noise distribution and resampling step is required.
A different localization approach is the use of ultrawideband (UWB) sensors, such as proposed in [23].In this paper, Güleret al. propose an on-board multi-robot localization approach which does not require infrastructure.therefore, three UWB sensors are mounted on an anchor robot which estimates the others robot location without inter-robot communication.However, the usage of UWB for localization has the major disadvantage that is is usually not very accurate, especially then obstacles are within the way.
Other recent localization monitoring approaches can be found in [24,25].In there, intregity risk metrics and methodologies are presented to verify a robots localization performance.

Review on machine learning approaches
The main application of convolutional neural networks (CNN) is the detection of features within images.Therefore a given image is used as input of a neural network and sent into convolutional layers which store single features in feature maps.Those features can then be used to recognize patterns for classification.A simple and modern example is presented by Shamov and Shelest [26].They present the main features of convolutional neural networks and show how they are applied for feature detection.For this purpose they created the task of detecting tower lighthouses from a video stream.
A new application of convolutional neural networks is the detection of dynamic obstacles in grid maps.Piewak et al. use a grid map and a deep neural network to detect whether grid cells within a map are moving [27].The difference to a normal tracking approach, like particle filters, is to use the complete map grid with a top view as input image.They also proposed an approach which is optimized for real time applications.
A Long-Short Term Memory (LSTM) network is good at classifying, predicting and processing time series.Its special power is that it is also good in learning and evaluating time series which have a time lag of unknown size between two events.It is well suited for training data that really has to remember longer time series and also focuses on the current input.This is a reason why in many fields LSTM performs better than other recurrent architectures [28].Gensler et al. investigted how good LSTM networks can predict a future event [29].Therefore they compared different types of deep learning networks on a defined task.The task was to forecast the power production of solar panels.With this task they want to provide a reliable power forecasting method that allows to efficiently operate a solar power station.The results of the experiment for solar power forecasting showed that in nearly all error measurements an auto-encoded LSTM performed the best.This also shows that recurrent networks and especially long-short term memory networks excel at their task of predicting future events.
Pattern recognition is an important task in visual computing.It is about detecting specific patterns within an given input.In visual computing these inputs are images which are used to detect certain content like humans.Modern pattern recognition methods are based on convolutional neural networks which aim to identify patterns within the input data by training the network on a big set of data.This method can be used to detect licence plates on a car [30] and is also applied in vineyards to detect birds which needs to be scared away [31].
Morales et al. presented a method for object tracking by using a 3D occupancy grid as environment representation and a particle filter based approach for detecting and tracking obstacles [32].
In [33], a classification approach is presented to detect the invasiveness of pulmonary subsolid nodules in CT images.Altough the topic is not in direct relation with robot localization, the training approach using adaptiveboost deep learning shows how to train a strong classifier using weak CNN-based classifiers.The results show, that it is possible to achieve an acceptable accuracy when training features extracted from a CNN for a binary classifier.

Localization accuracy estimation approach
In order to judge if the estimation of the robot's position is accurate we follow the approach depicted in Fig. 2. In our application we need to solve a binary classification problem: either the robot is well localized or delocalized.We declare the robot delocalized if the estimated position and orientation diverged form the true position and orientation more than a given threshold.
In order to have access to a ground truth of the position and orientation and to be able to provoke delocalization we simulate the robot's motion and sensing (i.e.laser scans) as well as the environment.The data from the perception as well as the data from the particle filter-based pose estimation (i.e.particles) are recorded and converted to labeled training data for the following classification steps.For the classification we follow two approaches which we also use in combination.The first approach follows the assumption that the probability distribution of the robot's pose modelled by the particle set and its typical growth of uncertainty bears information about the localization accuracy.Thus, we train non-recurrent and recurrent neural networks for the classification based on the particle set.The second approach follows the assumption that statistical features (weak classifiers) derived from the perception and the particle set can be used to learn a strong classifier using boosting, since the particle set represents a distribution Fig. 2 The overall approach for estimating the accuracy of robot localization in the environment.By analyzing this distribution with statistical tools and using the received scan data to match it to the environment, information on the localization state might be revealed.To see, how promising a feature is, the KL-divergence between a localized and delocalized set is calculated, as described in Section 3.3.Finally, the classification from the trained neural network can be used in addition to the weak classifiers in a boosting step, to see if the neural network output can further improve the classifier.
In the remainder of the section we introduce the following steps of the proposed process in more detail: (1) generation of training data, (2) selection and validation of neural network architectures, (3) validation and selection of features, and (4) learning classifiers.

Generation of training data
This section describes how the training data for the neural networks is generated.
Figure 3 shows the overall procedure of generating a training sample and finding the correct label.To create a sample one needs the particle set from the particle filter which is used to create a training image x.A particle represents a potential pose (position and orientation) along with a weight representing the importance of the particle x, y, θ, w .As the performance of the particle filter depends on the size of the particle set, a sufficient number of particles has to be selected.In this work, a particle set size of M = 1000 was selected.We then convert the particle set into a binary image because it is able to represent the distribution of the particle set and is suitable for state-of-the-art machine learning tools.The image itself then contains black pixels which represent a particle of the particle cloud.The image is then labelled by comparing the exact pose x GT , y GT , θ GT (retrieved from the simulation) and the estimated pose x P F , y P F , θ P F (retrieved by the distribution of the particles).Using a defined distance α, β as threshold between them one can classify the robots localization state and thus also create a label y for the training sample.the values for the thresholds are defined in Section 4.2.
Usually the particles are distributed over a large state space (i.e. the entire environment map).Since this map does differ in size and the particles accumulate at one area, it is not necessary to store the complete map as data sample.In order to focus on relevant areas of the particle distribution, the image is centered around the mean of the positions represented in the particle set.Moreover, the image is oriented along the mean of the orientations represented in the particles.Then the area of the distribution is cut and cropped into a default size.By analysing the distribution of the particle, the optimal area size was estimated.This was done by cutting the the area to a size s × s meters that when plotting the particle set, in 95% of the cases all particles are in the specified area.To determine the optimal side, sample data was recorded which collected the information of the particle cloud.The result was a cutting area with s = 1.5m.Using only a focused representation has the advantage that the trained network can be used for environments of different size and shape.The cut area is then cropped into an image of size p × p, where p is the number of pixels on the axis.The size of p depends on the details which are represented by the particle distribution and on the desired performance of the neural network training.After analyzing different values for p by pre training small neural networks, a size of p = 36 pixel was selected as this showed the best results with acceptable training performance.This values had been selected based on the results of run-time and accuracy.
Figure 4 shows generated images which are created during the movement of the robot, stored at different time steps.The first three images illustrate a particle set of a Fig. 3 The procedure of generating a data sample and labelling for training well localized robot while the fourth image is an example of a delocalized robot.To provoke the delocalization of the robot, dynamic obstacles where included, so the particle filter algorithm got disturbed.
The transformation and inflation of the particle cluster are indicators for the localization state of the robot.Thus a neural network can be trained to detect relevant patterns within this transformation.
To generate enough distinctive data samples the robot has to drive randomly within a predefined environment.This leads to different robot locations within the environment and samples representing different localization situation in the environment.Figure 6a shows an example environment.Due the facts that methods based on particle filters are quite robust and we need positive as well as negative examples we need to provoke delocalization.One of the simplest and most realistic attempts is to disturb the lidar sensor using dynamic obstacles that are not represented in the map.For this purpose we randomly generate additional obstacles around the robot in the simulation.Such obstacles cause unexpected measurements and may lead to delocalization.

Selecting and validating the neural networks
Having created a training set for neural networks one has to define an appropriate network structure (type, layout) which can be used for localization classification.Thus, several possible network structures need to be trained and validated.
When using data from a particle filters as input a priori it is hard to specify the best network type and structure since it is rather unclear which information is extracted by a neural network for estimating the robot's localization accuracy.A possible network structure is a convolutional neural network (CNN) which is used for pattern recognition.Another possibility is a recurrent network that may learn the transformation of the particle set over a certain time period.
To evaluate a wide range of possibilities, both network types are considered.A CNN for feature extraction and a recurrent network structure for learning the transformation over time.For learning a recurrent network a long-short term memory (LSTM) structure is used since a better performance is expected [28].Additionally a combination of both, CNN and LSTM, is trained.This is called a Long-Term Recurrent Convolutional Network (LRCN).By using this network structure it is evaluated if the advantages of both former types can be combined.
For each different network types (CNN, LSTM and LRCN) three different network layouts are created and trained.We generate a simple, a moderate, and a complex layout.Those three layouts are chosen to evaluate whether a simple network structure leads to underfitting or a complex structure leads to overfitting.This results in a total number of nine networks to be trained, which are summarized in Table 1.
To validate the usefulness of a network structure one needs to define a metric that allows to compare different network structures in respect to localization monitoring.A suitable metric is accuracy [34].Accuracy is a value that describes how well a given data set is classified.If a classification task with two classes is considered it can be used to calculate the accuracy by dividing the predictions of a data sample into a binary classification.
Having two classes, one can assign to a class either the positive or negative.To validate a data sample it is fed into, the network structure and the output is recorded.The output is then assigned to one of the two classes, depending on the result.By comparing the assigned class with the expected class label one can determine whether a class was correctly assigned.If both, the result and the label match, then the outcome is true.When a data sample is classified as positive and the expected class is also positive the data sample is true positive.Vice versa, if the data sample is classified as negative and the result is negative, it is categorized as true negative.If the two classifications do not match the outcome is false.If the data sample is classified as negative and the expected result is positive the outcome is false negative and false positive if the data sample is labelled as positive and the prediction is negative.By counting the results for the examples in the data set one can determine the accuracy as where t x indicates the number of correct classification and f x the number of the incorrect ones with x either p (positive) or n (negative).

Feature selection
Once an appropriate network structure has been selected useful features for the boosting step needs to be identified.By using features which can be extracted from the particle set and use them in a boosting step, it is aimed to improve the accuracy for estimating the localization state of the robot.In the context of the company in which the project was carried out, a localization state estimation called localization scoring was used prior this implementation.This localization scoring approach uses features that are extracted by comparing the scan data to the robots map.For example a Hough line transform [35] was performed on the map and scan points and compared to each other.Using features form the existing localization scoring mechanism and expanding it with relevant data on the particle set, such as the center of gravity, the features were investigated to see whether they reveal information on the robots localization estimate.
In general it is hard to specify which features may be important for localization scoring.Using a simple evaluation one can search for various possible features that might be useful as weak classifiers in scoring the accuracy of a robots localization.By using the thresholds α and β as above the robots localization status can be easily determined and every feature can be inspected on its own.We identified 33 possible features that are based on information from the robot's perception and data from the particle filter.Table 2 summarizes the selection of these features.For a complete description of the list we refer to [36].During the recording of the training data feature estimations are recorded and labelled with their localization state, as described in (1).
To determine which features are important the robot was randomly driven around in a simulation environment and the features were recorded.After generating roughly 100.000 feature samples every feature is analysed individually.These recordings are split into a localized X Q and delocalized X P data set.These sets are then represented as a discrete distribution over k bins to determine the divergence between the two classes.However, as continuous variables reveal more information about the robots location estimate, the original continuous values of the were used again for boosting after selecting the relevant features.The size k of the bin is selected by observing the prerecorded data set and evaluating its maximum and minimum values.Then, a preselected bin size was chosen considering the unit.for meters the bin size was 1cm, for percentage the bin size was 1% and for degree the bin size was 1 • .Each bin contains a number of samples which lie within that bin.A discrete probability distribution P is then represented as the number of samples within that bin divided by the total number of samples n The same holds for the distribution Q.Having two discrete distributions P , Q for the localized and delocalized set one can calculate the Kullback-Leibler divergence [37] as where k is the number of bins.The KL-divergence is only defined if ∀i : Q(i) = 0 → P (i) = 0 applies.If P (i) = 0 the contribution of the i-th bin is also 0. It is defined that D(P ||Q) ≥ 0 for all distributions and D(P ||Q) = 0 if Using the KL-divergence between the probability distribution of the feature values for the positive and negative cases allows to select well discriminating features.The higher the KL-divergence the more different the distributions are and the more informative is the underlying feature.On this basis, promising features can be selected for the boosting step.As criterion for the feature selection a threshold of D(P ||Q) ≥ 0.1 was selected.
Table 2 shows the KL-divergence of some example features calculated using the same simulated disturbed navigation as discussed above.The selected features which were used for boosting are highlighted in grey.Figure 5 depicts the probability distribution of feature values for a discriminating and a non-discriminating feature.

Adaptive boosting
Adaptive Boosting (AdaBoost) is a machine learning approach which uses supervised data for classification.Since a single feature by itself is a rather bad classifier for estimating the localization state, the idea is to combine multiple weak classifier which do not provide enough information about a certain class itself and combine them into a strong classifier that can be used for identifying classes [38,39].In general an adaptive boost algorithm takes N training samples (x i , y i ), 1 ≤ i ≤ N with x i ∈ R K and y i ∈ {−1, +1}.x i is the input vector which contains different components (features) that are used for training.y i is the desired label which is either −1 or +1.There exist several variants of boosting algorithm which all have a similar structure [40].In this work the standard discrete Ada Boosting algorithm which uses two classes is used.It uses its N-sized input set and initializes the weights for each input sample with w i = 1/N .Then a weak classifier f m (x), the weighted training error m , and the scaling factor c m is computed.For a detailed formula on how to compute these values, we refer to [41].Then the weights are increased for input samples that have been wrongly classified.After this step the weights are normalized and the steps for finding a new weak classifier are repeated M times.At the end a final classifier F (x) is found which uses the sign of the weighted sum of the input set.This can then be used to estimate the state of the input data, in our case the localization state.
In the proposed approach we use an appropriate selection of features from the previous section and the same examples from training trajectory used in training the neural network.Moreover, we treat the binary classification of the trained neural networks as additional weak classifier.

Support vector machines
A Support Vector Machine (SVM) is a popular tool for classification problems [42,43].It searches for an optimal hyperplane which can be used to separate two classes.It takes training samples as input which are classified and outputs an optimal hyperplane which can be used to categorize new samples.To find the optimal plane which separates the two classes best is not a trivial task.Consider m training samples (x (i) , y (i) ), 1 ≤ i ≤ m where x (i) is a sample which is labelled with y (i) ∈ {−1, 1}.If those samples are linearly separable, multiple possible solutions exist which could be applied.The task is now to find the best hyperplane, determined with w o and b o , that maximizes the separation space between the two classes w T o x + b o = 0.If data samples are not linearly separable, the so-called kernel methods can be applied.Those methods take the training samples which are not linearly separable and map them into a higher dimensional space using some non-linear transformation ϕ(x) [44].Kernels allow to operate in highdimensional feature spaces and allow a linear separation of the classes within a different state space.This is also called the kernel trick and has better performance than explicit computation methods.
As for the approach with AdaBoost, we use the features selected above in the previous section and the binary classification result from the neural networks as classifiers for the SVM.At the end, the two methods, AdaBoost and SVM, are then to be compared to each.First, the extracted features from Section 3.3 are used to estimate the localization state of the robot.Then, the binary classification of the trained neural network is used as additional feature to see, whether this further improves the performance of the SVM or Adaptive Boost.

Experimental evaluation
In the experimental evaluation we are interested to investigate how well the proposed approaches are able to identify a lost localization of a robot.Moreover, we are interested which network structure minimizes the uncertainty and if the approach is dependent on the environmental structure.In the remainder of the section we introduce the experimental setup, give details on the preparation of the training data, and discuss the results with respect to the used network structure and the used features as well as the boosting step.

Experimental setup
The training and evaluation is based on a simulation of a team of transport robots in industrial environments based on Stage [45] and provided by an industrial partner.The used navigation software is the same as embedded into the real transport robots and uses laser scans, odometry, a gyro, and a 2D gridmap for the localization based on a particle filter.The simulation is based on ROS Indigo [46] and runs on The training and classification is based on the libraries Caffe 1.0.0,OpenCV 3.3, and CUDA 7.5.17 and runs on a standard PC equipped with an Intel i7 (2.1 GHz), 8 GB RAM, Nvida GPU, and Ubuntu 16.04.
Figure 6 depicts the environment map used where the left map was used to collect data for training and the basic evaluation and the right map was used to evaluate if the trained classifiers generalize to other environments.

Training data
As described in Section 3 the robot is randomly driven trough the training environment which is simulated.As stated earlier, the simulated environment allows to add random obstacles within the environment that disturb the scan data and thus provoke delocalization.Also the precision of the scan data can be defined.In our simulation environment, the precision of +/-3cm are used, as this is the standard error of currently existing sensor systems.During this run all relevant output from the particle filter and the perception were recorded.Using the thresholds α = 1.2 m and β = 1.2 rad all data points were labeled.Localized samples got the label 1 and delocalized samples the label 0. The thresholds were selected in the way that they maximize the KL-distance of the hand-crafted localization scoring approach.
This leads to a sequence of 150.000 data points where 50% are labeled delocalized.For the training of the nonrecurrent networks this data are equally split into a training and evaluation set and shuffled randomly.For the training of the recurrent networks individual data points are useless.Thus, we split the data points in positive and negative sequences of the length of 2 × 7 data samples.The negative sequences are organized around a time point t where the robot become delocalized and 7 data samples from before and after t are added.Using a frequency of 5Hz this results in a sequence of about 3 seconds.This lead to 10.714 sequences.Moreover, the same number of positive sequences were generated where the robot stayed localized.Table 3 gives an overview on how the data was separated for training, validation and test purposes.The same data sets were also generated for the second environment which are used to evaluate the generalization of the approach and not used in training.

Validation of the neural networks
In order to evaluate the most appropriate network structure for the estimation of localization accuracy we trained each network structure with the proper (individual data points or sequence of data points) training data and evaluated the results using the evaluation set of the first environment as well as the data set of the second environment.As criterion we use the metric defined in Section 3.2.
Table 4 summarizes the results of this evaluation.The number of iterations in training was 100.000 for all networks.In general it can be seen that all networks perform quite well on the training data where we have equally many positive and negative examples.The mid-complex CNN performs best.For the validation set in the same environment the performance drops by about 2-13%.The mid-complex CNN is still in the lead but shows the largest performance loss showing a tendency for overfitting.For the data set from novel environment the performance drops for The number of data samples is equally distributed over localized (label 1) and delocalized (label 0) samples The accuracy is obviously a compressed metric.In order to analyze the performance of the mid-complex CNN, Table 5 shows the detailed results of its evaluation.It can be seen that the true and false results for the positive and negative cases are almost balanced.Thus, in contract to other networks this network is not biased in a sense of being too optimistic (many false negatives -missing delocalizations) or too pessimistic (many false positives -raising a lot of false alarms).Figure 7 shows the loss function during the training of the mid complex network.It can be seen, that the loss is decreased within the first 60.000 iterations and then does not decrease further.Thus, a total of 60.000 iterations would have been sufficient to train this neural network.We refer to [36] for more detailed results.

Result of SVM
The first method for classifying the localization state of a robot with various features are Support Vector Machines.SVMs are used in this paper to find an optimal separation hyperplane that separates the features introduced in Section 3.3 into localized and delocalized states.Therefore, two different feature sets are used for training.At first the selected features from Section 3.3 are used without the binary classification of the selected neural network output.This is then compared with a SVM that trains all the features including the neural network output.For finding an optimal separation hyperplane the library offered by OpenCV is used.This library offers an automatic trainer that automatically searches for the optimal separation hyperplane and also tries to detect the perfect Kernel for it.Since this SVM algorithm is well tested by OpenCV and its automatic mode decreases the computational effort this method was used for training.
Table 6a shows the outcome of the boosting with a SVM.It appears that this method is not applicable for this task since the testing accuracy in the same environment is with 61.94% lower than the results of the neural networks themselves.The accuracy in a different environment decreases further to 58.11% which is shockingly close to a random distribution of 50%.However, since this is the base line for comparison the neural network feature is still added and evaluated with a Support Vector Machine to see if it can be improved too.By adding the neural network feature one can improve the accuracy of a SVM in the same environment by about 2% as shown in Table 6b.Although this is an improvement compared to the SVM without the network feature it is still an insufficient result.Also, it can be seen, that the SVM fitted the network to detect more false negatives than false positives within the training environment.Due to the low accuracy it can be assumed that a SVM is not suitable for solving the task of estimating a robots localization state.

Result of Ada boosting
First we evaluate how well Ada Boosting using the features introduced in Section 3.3 performs on the classification task.The features represent information extracted from the particle filter and the actual perception.Using the discrimination metric 26 promising features out of 33 were selected for boosting.In this evaluation no output of the neural networks had been used.
In order to find the optimal number of weak classifiers learned in boosting we varied the maximum number between 26 and 300. Figure 8a shows the performance of the learned classification on the training and evaluation data set.We selected 276 weak classifiers as an optimal number.Table 7a shows the performance of the so learned classifier on the training and the evaluation sets.The performance for the unknown environment is with 81.83% significantly better than the trained neural networks.But this comes at a price as the false negative rate is significantly higher than the false positive rate.This can be a problem for robots acting in an open environment.
Figure 8b shows the same learning and evaluation process where additionally the output of the mid-complex CNN was used as a feature.Here the optimal number of weak classifiers is 229 as the accuracy on the validation set does not significantly increase with further classifiers.Table 7b shows the performance of the so learned classifier on the training and the evaluation sets.It is clearly visible that the additional feature allows only a marginal performance improvement of about 1% on the data from the novel environment.However, the distribution of false negatives and false positives is now more balanced.

Conclusion and future work
Estimating its own position in an environment is a crucial capability of intelligent mobile robots.Knowing reliably the state of the localization process is important for the dependability of such robot systems.
Following the observation that most robots use distance sensors, environment maps, and particle filter based approaches we investigate if information generated from these components (patterns and statistical features of the particle set) can be used to train a reliable localization monitor.Using training data from test runs in a simulation environment where delocalization was provoked we trained several deep networks (non-recurrent and recurrent) and a boosting approach.While the trained deep networ ks showed moderate classification rates in particular in an environment not used for training and the SVM delivers insufficient results, the adaptive boosting approach using the network output as an additional feature showed detection rates of more than 80%.In conclusio, the evaluation of different training settings showed that it is possible to use information obtained from the particle set for scoring the localization accuracy.
For future work it needs to be investigated if delocalization cased by other uncertainties such as slipping can be detected as well.Moreover, a more detailed investigation of networks structures needs to be done.As the results show a significant decrease in accuracy when using the approach in a different environment, an investigation is to be conducted to see whether the approach can be generalized.Finally, so far only the position of particles is used for training.
It would be interesting if their orientation bears additional useful information.

Fig. 1
Fig. 1 General overivew on localization using particle filters

Fig. 4
Fig. 4 Example images generated for training the neural network.a Localized particle set at time step t. b Localized set at time t + 40 s.c Localized set at time step t + 90 s. d Delocalized image at time step t + 110 s

Fig. 5
Fig.5 Probability distributions of feature values for localized and delocalized cases for feature 18 (left, non-discriminating) and feature31 (right, discriminating)

Fig. 6
Fig. 6 Environment maps used for training and evaluation.Based on an real industrial environment enriched with additional clutter

Table 1
Overview of the tested network types and their complexity

Table 2
List of features which were used as weak classifiers

Table 3
Separation of training data, validation data and test data

Table 5
Result of the mid-complex CNN

Table 6
Separation results using a SVM (a) without the neural network feature (b) with the neural network feature

Table 7
Boosting result using AdaBoost (a) without the neural network feature and 276 weak classifiers (b) with the neural network feature and 229 weak classifiers