I Introduction
Graphbased approaches are widely considered the most popular solution to the simultaneous localization and mapping (SLAM) problem [1] due to their reliability, flexibility, and favorable computational cost [2, Sec. II]. In lifelong SLAM [3], we assume that an area is not just mapped once but revisited continuously over a long period of time. This leads to a number of additional challenges compared to a standard SLAM problem where an area is only covered once or twice during a short period of time [2, Sec. III & IV]. This paper focuses on one of these challenges: scalability. The graph grows indefinitely over time, which makes the SLAM algorithm prohibitively slow and implies virtually unlimited memory requirements.
To illustrate the problem, imagine an environment of a limited size. A robot keeps moving around in this environment for a long time while running a naive graphbased SLAM algorithm that continuously adds new measurements to its map. The number of vertices in the graph as well as the memory requirements will typically increase roughly linearly w.r.t. the time spent driving. As parts of the environment are traversed again and again, more and more loop closure edges will get added between nearby vertices. In the worst case, this leads to a quadratic number of edges w.r.t. time and to a graph whose structure is not sparse anymore. This impacts the runtime severely, because commonly used optimizers such as g2o [4] or gtsam [5] are only efficient if the graph is sparse and become extremely slow for dense graphs.
The goal of graph pruning is to carefully remove redundant vertices and edges, while minimizing the loss in localization and mapping accuracy. The algorithm ensures that the size of the graph only scales with the size of the environment and becomes independent of the time spent there or the length of the robot’s trajectory. An example of the effects of pruning is shown in Fig. 1.
In this paper, we present a novel approach that performs vertex pruning based on geometric criteria. Specifically, we make the following contributions:

An efficient approach for vertex pruning based on a scaleinvariant density that does not depend on the data associated to a given vertex and leads to an even distribution of vertices in the pruned graph.

A marginalization algorithm that is robust to incorrect loop closures.

A thorough evaluation on two realworld datasets using different pruning parameters.

An evaluation scheme that does not require ground truth and is independent of the particular SLAM algorithm by comparing the graph with and without pruning.
Ii Related work
Several other authors have considered the problem of graph pruning before. Kretschmar et al. [6] proposed a pruning scheme where the information gain of each lidar scan is computed using a probabilistic occupancy grid. Then, the vertices whose lidar scans provide the least additional information are removed. Marginalization of vertices is performed by replacing all edges of the pruned vertex with edges connected to one of its neighbors. An improved version of the paper uses a Chow–Liu tree [7] to approximate the clique obtained by marginalization [8].
Eade et al. proposed a similar graph pruning approach [9, Sec. VIII]
in the context of visual SLAM. Nodes in the graph are divided into view and pose nodes, the latter of which can be pruned after a while because they only encode the robot’s trajectory and constraints in the graph but are not relevant for the map itself. Marginalization is performed by inserting binary edges between all vertices adjacent to a pruned vertex, which leads to the creation of numerous new edges. Hence, the authors proposed an edge pruning algorithm that uses a simple heuristic to remove superfluous edges.
Lazaro et al. suggested a somewhat different approach where multiple scans are aggregated into a local map, which corresponds to a single vertex in the graph [10]. Later, several local maps can be merged into a single vertex by building on the idea of condensed measurements previously proposed by Grisetti et al. [11]. Their approach effectively combines pruning and change detection into a single operation. A disadvantage is that newly created vertices are not located on the robot’s trajectory anymore.
CarlevarisBianco et al. have further investigated the issue of marginalizing vertices from the graph [12, 13]. They propose an exact marginalization using nary factors and a sparse approximation, which is based on a Chow–Liu tree. This approach is different from [8], where the treebased approximation is computed from binary factors. The benefit of considering nary factors is that correlations between binary factors are not ignored and that the approach can also be applied to factors that do not constrain the full state (e.g., bearing only measurements, range only measurements). Later works by Mazuran et al. [14] and Ta et al. [15] generalize and improve this approach further.
A pruning approach for featurebased SLAM was proposed by [16]
. They use the Kullback–Leibler divergence between the graph before and after pruning to decide which vertex to prune. For marginalization, they consider different ways to approximate the clique using fewer edges. Unfortunately, their work is not directly applicable to nonfeature based SLAM, e.g., lidar SLAM based on scan matching.
Iii Vertex Pruning






In order to keep the number of vertices inside the graph limited, it is necessary to carefully remove vertices over time. In general, there are two strategies, selecting one or more vertices and 1) marginalizing them [6, 8, 9], or 2) fusing them into one [10]. In lidar SLAM, the second option has the disadvantage that the lidar scans have to be fused and that pose of the new vertex does not match the pose where the lidar scans were originally observed. This complicates or invalidates potential successive processing steps such as raytracing of free space using scan poses as, e.g., proposed for change detection in [17]. In addition, all further information that may be associated with the involved vertices has to be fused as well. Finally, the fused vertex is not part of the original trajectory and may be located on impassable terrain. For these reasons, we only consider the first option in the following.
Hence, the key question is: How do we choose which vertices to remove? In general, one could consider the simultaneous removal of multiple vertices and analyze the combined effect of this operation. However, this is difficult to analyze and for the sake of simplicity, we only consider removing a single vertex at a time. When choosing which vertex to remove, we would like to pick a vertex that is of low importance to the overall map quality and whose removal degrades the expected future performance of the SLAM algorithm as little as possible. In order to quantify the importance of a vertex, geometric or informationtheoretic measures can be used.
Informationtheoretic measures try to quantify the information contained in each vertex or, as proposed by Kretschmar et al. [6, Sec. 4.1], [8, Sec. IV]
, its associated laser scan. The idea is then to remove vertices in such a way that the loss of information is minimized. The main problem with this approach is that the actual information within a scan is hard to quantify. Kretschmar et al. compute the information based on the occupation probabilities of cells in a probabilistic grid map. However, this approach tends to consider scans with high noise, false measurements, or poor alignment to be high in information because they affect the probabilities inside the grid map more strongly than well aligned scans with little noise that agree well with other scans of the environment. Also, this information measure is fairly expensive to compute because all points in scans within range of the current vertex contributed to the result. This can be problematic especially for lidars with long range and high resolutions.
For these reasons, we propose to use geometric methods for vertex pruning. The basic idea is to keep the density of vertices below a certain threshold across the entire map, i.e., we seek to remove vertices in places where a lot of vertices cluster within a small area. For this purpose, we propose a novel density measure, which we call scaleinvariant density.
Iiia Scaleinvariant density
Consider vertices with unique positions , where for .
Definition 1.
For a radius , define the density of a vertex as
where is the Euclidean norm. This corresponds to the number of vertices (excluding ) inside the circle of radius around divided by the circle’s area.
However, this measure is strongly dependent on the choice of the radius . In particular, after pruning vertices inside the radius , there may be a lot vertices whose distance from is just slightly larger than , which leads to underestimating the actual density of vertices around .
In a pruning algorithm, this can leave behind clusters of vertices that are just slightly more than apart. This issue is illustrated in Fig. 2, where we consider the effect of vertex pruning on two synthetic examples.
To address the deficiencies of the radiusbased density , we propose the scaleinvariant density, which is independent of any fixed scale . It is obtained by integrating over all .
Definition 2.
The scaleinvariant density of is given by
As a result, this measure is not susceptible to the clustering problem mentioned above. In order to efficiently compute the density measure, we use the following theorem.
Theorem 1.
It holds that
Proof.
Wlog., we consider the parameter and assume that all other for are sorted by distance according to . Then, we can simplify
For practical use, this sum can be truncated to the closest neighbors, which speeds up computation and leads to a more localized version of the density. Based on Theorem 1, it is easy to show that scaling the positions of the vertices corresponds to inversely scaling the density.
IiiB Vertex Pruning Algorithm
Based on the density function introduced above, we propose a vertex pruning algorithm. This algorithm tries to find the vertex with highest density, marginalizes it and repeats until the graph size is below a predefined threshold or no more vertices have a density below the density threshold . Furthermore, it can be desirable to exclude certain vertices from pruning. For example, we always keep the most recent vertices to prevent pruning vertices prematurely, which could lead to poor robustness. We define the set of prunable vertices as and only consider theses vertices as candidates for pruning. Note that the density is still calculated using all vertices including nonprunable ones. Pseudocode of the proposed algorithm is given in Algorithm 1.
It should be noted that the proposed algorithm does not consider the direction the robot was facing at each vertex. This is fine for robots whose field of view is 360°or close to that. However, for robots with a small field of view it may be necessary to first partition the vertices by viewing angle and prune each partition separately.
IiiC Marginalization
Once we have chosen a vertex to prune from the graph, we need to consider how to perform the removal of a given vertex. The graph encodes the joint probability distribution of the random variables represented by all vertices. Thus, removing a vertex corresponds to marginalizing a set of random variables.
Marginalization is typically performed by creating an nary constraint between all neighbors of the pruned vertex [13, Sec. II]. Some authors approximate this nary constraint by connecting all neighbors of the pruned vertex pairwise, i.e., vertices adjacent to the pruned vertex will form a clique. This leads to a huge number of edges if the pruned vertex has a lot of neighbors or if several vertices are pruned successively [8, Sec. VA]. As a result, this method for marginalization is infeasible in practice, at least unless an aggressive edge pruning scheme is employed to remove most of these edges later (see Sec. IV).
To obtain a more sparse approximation of the graph in the neighborhood of the pruned vertex, several authors have proposed using a Chow–Liu tree [7] either directly from an nary constraint [13] or by selecting a subset of the fully connected graph of binary constraints among all neighbors [8, Sec. VC]. Either way, the Chow–Liu tree approximates the neighborhood of the marginalized vertex using a spanning tree weighted by the mutual information between the vertices involved. As the local tree structure might be too sparse in practice, other subgraphs have been considered, e.g., the Cliquey approach [14].
IiiC1 Impact of Wrong Loop Closures on Marginalization
An important issue that has not been sufficiently considered in related work on graph pruning is the influence of incorrect loop closures during marginalization. Typical SLAM algorithms can be made resilient to a small number of incorrect loop closures by using robust kernels in the error function or by specifically determining which loop closures are outliers using an expectation maximization algorithm
[18]. However, marginalization can amplify wrong loop closures. After multiple marginalization steps, even a single wrong loop closure might be amplified enough to lead to divergence. This is particularly problematic because wrong loop closures often do not have higher uncertainty according to their information matrix than correct loop closures because the scan matching algorithm only provides a local uncertainty estimate.
Consider the following example. We would like to prune a vertex with edges to . Let us assume that a single edge is incorrect, i.e., the ratio of incorrect edges is . After marginalization of , the combined nary constraint would be incorrect. We now approximate this constraint by creating pairwise edges for . This yields newly created edges, of which are wrong edges. Therefore, the number of incorrect edges increased from to and the ratio of incorrect edges increased from to . In addition, if there already was an edge between one of ’s neighbors for before marginalization, the existing and new edges would get fused. This leads to an incorrect edge if either of the original edges was incorrect, i.e., the ratio of incorrect edges would be even larger.
In practice however, we make the observation that odometry edges are typically much more reliable than loop closure edges even if they possess the same uncertainty. While odometry can drift over time and may be somewhat inaccurate, it is usually not affected by outliers with huge errors like loop closures. As a result, we propose a marginalization algorithm that is founded on the reliability of the odometry chain and is fairly robust to wrong loop closures.
IiiC2 Marginalization Algorithm
In [6, Sec. 4.2], an approximate marginalization was proposed where all edges of the pruned vertex are collapsed into a single neighboring vertex. This vertex is chosen such that the sum of lengths of all resulting edges is minimized. We propose a similar method that also collapses edges into adjacent vertices. However, we do not collapse all edges into a single vertex but instead move each loop closure forward or backward along the odometry chain depending on the distance between the vertices of the newly created edge. Pseudocode of this procedure is given in Algorithm 2.
To concatenate edges, we need three elementary operations, composition of two edges, combination of two edges and inversion of an edge, see [9, Sec. VIIB]. The equations for these operations depend on the exact error function used in the optimizer (in particular the frame covariances are represented in) and can be computed in closedform.
In the proposed algorithm, the number of wrong loop closures does not change during marginalization. Considering the same example as above, our algorithm reduces the number of edges from to , and thus the ratio of incorrect loop closures rises only slightly from to . When combining edges, we can detect contradictory edges by considering the differences in their relative poses or the Mahalanobis distance. If an odometry edge and a loop closure contradict, we trust the odometry edge and delete the loop closure. If two loop closures contradict, we delete them both as losing a correct loop closure is usually preferable to having an incorrect loop closure in the graph.
IiiD Evaluation
To investigate the behavior of the proposed algorithm in the presence of wrong loop closures, we created a synthetic example based on a grid similar to Fig. 2 and randomly replaced a given percentage of the loop closures with incorrect loop closures. Then, we applied pruning using different marginalization methods, optimized the graph again and compared the vertex positions to the true positions. In Fig. 3, we show the results from 50 Monte Carlo runs with different random broken loop closures. It can be seen that the proposed algorithm even slightly outperforms the original SLAM algorithm without pruning, which is mainly due to the detection and removal of some of the false loop closures when they contradict other edges. The ChowLiu approach performs significantly worse and diverges completely in some cases.
Iv Edge Pruning
In addition to vertex pruning and the removal of edges during marginalization, we also propose an independent edge pruning step, which removes superfluous edges without touching vertices at all. This is mostly beneficial in areas where vertices are not particularly dense but a lot of loop closures have been found. Also, depending on the particular algorithm, a lot of edges may be created during marginalization. In general, removing an edge is equivalent to removing one constraint from the graph and means losing the information encoded in that edge.
To choose which edges to prune, we consider the vertex with the largest number of edges and remove one of its edges via an edge selection criterion. This step is repeated until all vertices have less than a predefined number of edges . The key question is how to choose the edge selection criterion. Eade et al. employed pruning based on the edge residuals. Edges with the smallest residual are pruned first [9, Sec. VIIIC]. However, using the residual (or the error) means that edges with large errors (e.g., noisy edges, poor scan matches, wrong loop closures) will be kept and edges with small errors (good matches) will be discarded.
We therefore propose to use the trace or determinant of the edges’ covariance/information matrices as an alternative criterion. As the differences between these measures seem to be small in practice, we decided to use the trace of the information matrix, which is efficient to compute. We favor keeping edges with more certain estimates over edges with less certain estimates, i.e., edges with little information will be pruned.
Removing edges comes with the risk of removing important global loop closures or even disconnecting the graph. To avoid pruning edges that keep the graph together, we consider the factor by which the shortest path gets longer when removing an edge. Assume we want to prune the edge with length . If we removed this edge, the length of the shortest path from to would become , which can be computed efficiently using the A* algorithm [19]. We therefore compute the ratio and only allow pruning if it does not exceed a predefined threshold . Furthermore, we only prune loop closure but no odometry edges.
V Evaluation
Config Name  
10  50  50  5  5.0  
10  50  50  5  5.0  
No pruning 
To evaluate our graph pruning method, we apply it as part of a graphbased SLAM algorithm with 360 LIDAR data and evaluate the SLAM performance depending on different graph pruning parameters (see Table I) on two public datasets (see Table II).
Our SLAM algorithm uses g2o [4]
for graph optimization. Vertices in our graph have a single laser scan associated to them. New vertices are added in regular distances as keyframes; they are connected by odometry edges that encode incremental localization information. To obtain an odometry edge, we perform early fusion of scantoscan matches, wheel odometry, and IMU, if available. This makes our odometry edges quite reliable w.r.t. outliers and motivates our focus on the odometry chain for this work. Additionally, we add loop closure edges by matching scans of nonsubsequent vertices. For scan matching, we use the Normal Distributions Transform (NDT)
[20]. Place recognition techniques are employed to merge graphs of previous SLAM runs into the current graph. This allows us to build a consistent map using recordings from multiple days without knowing the robot’s starting pose in a shared reference frame beforehand. Dynamic objects are handled with a technique inspired by [17].Va Datasets
Dataset  Run  Travelled  Duration 
(m)  (hr)  
LCAS Strands Care Home  20161125  2634  7.68 
20161202  2092  9.32  
20161216  2081  8.00  
Sum  6807  25.00  
MIT Stata Center  20120118090907  683  0.60 
20120125121425  348  0.33  
20120125123329  239  0.23  
20120128111201  635  0.60  
20120202104408  1003  0.87  
Sum  2908  2.63 
As our first evaluation dataset, we use the Care Home^{1}^{1}1https://lcas.lincoln.ac.uk/nextcloud/shared/datasets/aaf.html dataset from the LCASSTRANDS longterm dataset collection, initially recorded for [3]. We picked the first three days with long trajectories () that had no recording issues (see Table II). Note that this dataset contains a lot of dynamic objects. Since it does not come with ground truth information, we run our SLAM algorithm without pruning () and use it as a baseline for comparison of trajectories and vertex positions of the pruned graph. This way, we can ascertain how much pruning changes the results of the SLAM algorithm and thus degrades the performance. This experiment design measures how pruning impacts relative SLAM performance rather than the SLAM algorithm’s performance as a whole.
As our second evaluation dataset, we use five mapping sessions of the second floor of the MIT Stata Center^{2}^{2}2http://projects.csail.mit.edu/stata/index.php. We have chosen the same datasets as [10] (see Table II). We use the laser data from the base of the PR2 robot (a Hokuyo UTM30LX Laser), the Microstrain 3DMGX2 IMU and the robot’s raw wheel odometry. This dataset comes with ground truth information, so in addition to the pruning/nopruning comparison, we can also evaluate our SLAM algorithm’s performance as a whole.
VB Effects on Runtime and Memory


We run our SLAM algorithm on both datasets (see Figure 4) and measure the time it takes to process a single scan (see Figure 5). This includes the duration of scan matching, loop closure search, and the graph optimization. Note that, especially for bigger graphs, the graph optimization makes up the largest portion of the measured time per scan. Considering time per scan and the LIDAR’s frequency allows us to judge whether the algorithm can run in realtime on a robot. All presented results were obtained on a single core of an Intel® i78650U @ 1.90GHz CPU.
In order to assess how pruning improves the memory requirements of our SLAM algorithm, we keep track of the number of vertices in its graph. Only laserscans that are associated to a vertex are needed longterm. So, whenever a vertex gets pruned, the associated laserscan can also be deleted. As visualized on the right side of Figure 5, the number of vertices grows linearly with time when no pruning is performed (). However, if pruning is performed, the number of vertices begins to converge to a level dependent on the and the size of the environment. When using the configuration for the MIT Stata Center dataset, this barely starts to happen. The effect is more evident for the longer LCASStrands Care Home dataset, or when using the configuration. Note that we start with a new graph for each individual bag file of a dataset, so the number of vertices drops down to zero at that point in time. After place recognition succeeds, we merge the new and the old graphs and the number of vertices jumps back up.
VC Effects on Trajectory Error and Map Quality
Dataset  Config  TE Mean SD  ME Mean SD  RME Mean SD  Runtime  Speedup  
name  (m)  (deg)  (m)  (deg)  (m)  (deg)  (s)  
MIT Stata Center  Comparison with ground truth  
—  
Comparison with SLAM ()  
—  —  —  —  —  —  —  
LCAS Strands Care Home  Comparison with SLAM ()  
—  —  —  —  —  —  —  
SLAM results vs. the ground truth information (top) and vs. reference SLAM configuration (center and bottom). Used metrics: Mean trajectory error (TE), map error (ME) and relative map error (RME) and their respective standard deviation (SD).
In the following, we consider several different error measures to assess the accuracy of the SLAM algorithm. The trajectory error (TE) measures the absolute position errors of the robot’s trajectory during execution. We calculate this metric by capturing our SLAM algorithm’s pose estimate after each scan. This pose estimate gets compared to the reference pose estimate or to the ground truth. This metric does not include corrections to the estimated trajectory that become available later, e.g., through loop closures. The map error (ME) shows the absolute position errors of the final map, after all data was processed. We calculate this metric by comparing the final poses of vertices in the SLAM graph to a reference graph or the ground truth. Both the trajectory error and the map error assess global consistency. For example, these metrics will punish small errors in orientation severely, especially if the map is big. We argue that for most usecases of lifelong SLAM, local consistency is sufficient. Therefore, we suggest to focus on the relative map error (RME) as proposed by [21, Section IIIA]. We calculate this metric by comparing the relative poses between subsequent vertices in the final SLAM graph with the respective relative pose from a reference graph or the ground truth. Note that this procedure may degenerate to assess global consistency, instead. This can happen if so many vertices get pruned that subsequent vertices are not close to each other anymore. In practice, this was not an issue throughout our evaluation, because our pruning thresholds are not that strict.
In Table III, we compare SLAM results obtained with different configurations (see Table I) to the ground truth of the MIT Stata Center dataset. Our experiments show that graph pruning results in a speedup of at the cost of slightly increased errors on this dataset. If no pruning is performed, the resource requirements of the SLAM algorithm do not stop growing. This is naturally more evident for the longer LCAS Strands Care Home dataset, where the speedup is . Note that the algorithm without pruning is not able to maintain realtime performance towards the end of the run. When using our pruning method, however, this is not the case. Also, we compare the results when performing graph pruning with the results of the configuration, which does not perform pruning. This evaluation is helpful when tuning pruning parameters for specific robots or environments that have different resource requirements, as it does not require ground truth information. Especially when considering the RME, the errors w.r.t. roughly fit the increase in errors w.r.t. the ground truth for MIT Stata Center.
Vi Conclusion
We presented a novel graph pruning algorithm for use in lifelong SLAM. Our key idea is to compute the vertex density using a special density function and to prune vertices in areas where their density is too high. Furthermore, we proposed an approach to marginalize vertices that is robust to wrong loop closures, which was experimentally validated on synthetic examples where it performs significantly better than the Chow–Liu tree method used by [8]. Our pruning algorithm was evaluated on two public realworld datasets. The evaluation shows a significant speedup (2.69 for the smaller and 40.68 for the larger dataset) at a moderate cost in terms of accuracy compared to standard graphbased SLAM without pruning ( and map error, respectively). Compared to the ground truth, the trajectory error increased from to (or ) when using pruning.
While the proposed algorithm is limited to 2D SLAM, generalization to 3D is straightforward. Future work may include the consideration of the combined effect of pruning multiple vertices/edges as well as the development of more sophisticated marginalization algorithms that are still robust to incorrect loop closures.
References
 [1] G. Grisetti, R. Kümmerle, C. Stachniss, and W. Burgard, “A tutorial on graphbased SLAM,” IEEE Intelligent Transportation Systems Magazine, vol. 2, no. 4, pp. 31–43, 2010.
 [2] C. Cadena, L. Carlone, H. Carrillo, Y. Latif, D. Scaramuzza, J. Neira, I. Reid, and J. J. Leonard, “Past, present, and future of simultaneous localization and mapping: Toward the robustperception age,” IEEE Transactions on Robotics, vol. 32, no. 6, pp. 1309–1332, dec 2016.
 [3] T. Krajnik, J. P. Fentanes, M. Hanheide, and T. Duckett, “Persistent localization and lifelong mapping in changing environments using the frequency map enhancement,” in 2016 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS), oct 2016.
 [4] R. Kümmerle, G. Grisetti, H. Strasdat, K. Konolige, and W. Burgard, “g^{2}o: A general framework for graph optimization,” in Proc. IEEE Int. Conf. Robotics and Automation, May 2011, pp. 3607–3613.
 [5] F. Dellaert, “Factor graphs and gtsam: A handson introduction,” Georgia Institute of Technology, Tech. Rep. GTRIMCP&R2012002, Sept. 2012.
 [6] H. Kretzschmar, G. Grisetti, and C. Stachniss, “Lifelong map learning for graphbased slam in static environments,” KIKünstliche Intelligenz, vol. 24, no. 3, pp. 199–206, 2010.
 [7] C. Chow and C. Liu, “Approximating discrete probability distributions with dependence trees,” IEEE Transactions on Information Theory, vol. 14, no. 3, pp. 462–467, may 1968.
 [8] H. Kretzschmar, C. Stachniss, and G. Grisetti, “Efficient informationtheoretic graph pruning for graphbased SLAM with laser range finders,” in IEEE/RSJ Int. Conf. Intelligent Robots and Systems, 2011, pp. 865–871.
 [9] E. Eade, P. Fong, and M. E. Munich, “Monocular graph SLAM with complexity reduction,” in Proc. IEEE/RSJ Int. Conf. Intelligent Robots and Systems, 2010, pp. 3017–3024.
 [10] M. T. Lázaro, R. Capobianco, and G. Grisetti, “Efficient longterm mapping in dynamic environments,” in Proc. IEEE/RSJ Int. Conf. Intelligent Robots and Systems (IROS), 2018, pp. 153–160.
 [11] G. Grisetti, R. Kummerle, and K. Ni, “Robust optimization of factor graphs by using condensed measurements,” in 2012 IEEE/RSJ International Conference on Intelligent Robots and Systems, oct 2012.
 [12] N. CarlevarisBianco and R. M. Eustice, “Generic factorbased node marginalization and edge sparsification for posegraph SLAM,” in 2013 IEEE International Conference on Robotics and Automation, may 2013.
 [13] N. CarlevarisBianco, M. Kaess, and R. M. Eustice, “Generic node removal for factorgraph SLAM,” IEEE Transactions on Robotics, vol. 30, no. 6, pp. 1371–1385, dec 2014.
 [14] M. Mazuran, W. Burgard, and G. D. Tipaldi, “Nonlinear factor recovery for longterm SLAM,” The International Journal of Robotics Research, vol. 35, no. 13, pp. 50–72, jun 2015.
 [15] D.N. Ta, N. Banerjee, S. Eick, S. Lenser, and M. E. Munich, “Fast nonlinear approximation of pose graph node marginalization,” in 2018 IEEE International Conference on Robotics and Automation (ICRA). IEEE, 2018, pp. 2494–2501.
 [16] Y. Wang, R. Xiong, and S. Huang, “A pose pruning driven solution to pose feature GraphSLAM,” Advanced Robotics, vol. 29, no. 10, pp. 683–698, jan 2015.
 [17] J. P. Underwood, D. Gillsjö, T. Bailey, and V. Vlaskine, “Explicit 3D change detection using raytracing in spherical coordinates,” in Proc. IEEE Int. Conf. Robotics and Automation, 2013, pp. 4735–4741.
 [18] G. H. Lee, F. Fraundorfer, and M. Pollefeys, “Robust posegraph loopclosures with expectationmaximization,” in 2013 IEEE/RSJ International Conference on Intelligent Robots and Systems, nov 2013.
 [19] P. Hart, N. Nilsson, and B. Raphael, “A formal basis for the heuristic determination of minimum cost paths,” IEEE Transactions on Systems Science and Cybernetics, vol. 4, no. 2, pp. 100–107, 1968.
 [20] P. Biber and W. Strasser, “The normal distributions transform: a new approach to laser scan matching,” in 2003 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS 2003), 2003.
 [21] W. Burgard, C. Stachniss, G. Grisetti, B. Steder, R. Kümmerle, C. Dornhege, M. Ruhnke, A. Kleiner, and J. D. Tardös, “A comparison of slam algorithms based on a graph of relations,” in 2009 IEEE/RSJ International Conference on Intelligent Robots and Systems, 2009, pp. 2089–2095.