Course:CPSC522/Rao Blackwellized Particle Filtering
Rao Blackwellized Particle Filtering for Grid Mapping
Monte Carlo methods are a common approach for large dimensional problems such as grid mapping, RaoBlackwellized Particle Filters aim to do the needed sampling as efficiently as possible.
Principal Author: Jocelyn Minns
Collaborators:
Abstract
RaoBlackwellized Particle Filters (RBPF) incorporates the Rao–Blackwell theorem to improve the sampling done in a particle filter by marginalizing out some variables. Particle filters alone struggle in high dimensional spaces since a large number of particles must be used, but with improved RBPFs the amount of particles needed decreases greatly.
The two papers discussed are:
 RaoBlackwellised Particle Filtering for Dynamic Bayesian Networks^{[1]}
 Improved Techniques for Grid Mapping With RaoBlackwellized Particle Filters ^{[2]}
Builds on
Particle fillers are a Monte Carlo methods that are used to update nonlinear system that can be modelled as a Hidden Markov network.
Related Pages
RaoBlackwell theorem, Monte Carlo methods, Particle fillers, Hidden Markov networks, SLAM.
Background
Particle Filters
Main article: Particle fillers
Particle filtering is used to estimate the updated state of a system given noisy sensor data and noisy state transition information. Unlike other filters, such as a Kalman Filter, particle filtering works with nonlinear systems and nonGaussian noise distributions.
The system must first be modelled as a hidden Markov model where the observed value are the measurements and the hidden variables are the true values . It is also needed to define the transition probability density with an initial probability density and the conditional distribution .
Initially, a set of N particles are sampled from the initial probability density . The size of the set of particles will depend on the complexity of the problem and the accuracy that is desired. For the initial set, all particles will be given equal weight. Meaning each particle is equally likely to be the true status of the system.
The algorithm for a particle filter can be described in three steps:
 Sample a new set of particles given the weights of the previous particles. Particles with a heaver weight will be sampled more often.
 Sample the updated state given the previous state and the current measurement data
 Update the weights of each particle and normalize the weights.
Given a large enough set of particles, this method has been shown to be effective in estimating state transitions for a variety of nonlinear problems.
RaoBlackwell theorem
Main article: RaoBlackwell theorem
The RaoBlackwell theorem states that given any kind of estimator of a parameter then the conditional expectation of given a sufficient statistic is typically a better estimator of .
While this may seem to be a simple statement, in practice the improvement of the estimator can be very large.
Papers
Paper 1: RaoBlackwellised Particle Filtering for Dynamic Bayesian Networks
The aim of RaoBlackwellised Particle Filtering is to find an estimator of the conditional distribution such that less particles will be required to reach the same accuracy as a typical particle filter.
The objective remains the same as a particle filter: to estimate the distribution .
The posterior gives the recession:
However, this cannot be solved analytically hence why sampling is typically used. To decrease the state space into something that will be easier to sample from, they break the hidden variables into two groups: and such that . The original problem assumes no structure to the hidden variables, but by breaking them into two groups as stated above, they can determine part of the posterior analytically then marginalize out the variables .
The posterior can be decompose as such:
The dimensions of are smaller than the original so they can obtain better results.
RaoBlackwellisation
For sampling, It is needed to define an estimator for .
First they define an estimate of any function f of the hidden variables as
The paper goes on to use the Central Limit theorem to prove why this is an estimate, but that has been excluded from this page as it is unnecessary to understand the implementation of RBPF.
Since each particle may not have even weights, they replace the factor with the normalized weights of each particle giving:
However, since they have shown that they can marginalize out , they can write a RaoBlackwell estimator as
Using this estimator can be proven to reduce the number of particles needed to reach a desired level of accuracy.
Algorithm
The algorithm for RaoBlackwellised Particle Filtering can be described as three steps:
 Sequential importance sampling step where they sample N particles from our given estimator, then compute and normalize their importance weights
 Selection step where resample N particles based on the weights calculated above
 Apply the transition step based on the probability distribution
Application for Grid Mapping
This form of filtering has been shown to be useful in Simultaneous Localization and Mapping (SLAM) of a grid. Consider a robot that can move in a discrete, dimensional grid. Both the motors of the robot and the sensors are imperfect and have some noise to their readings. Say you want to colour the map grid as the robot moves around the room to indicate obstacle that it should avoid. Let be the location at time t, be the colour of grid cell i and time t. The hidden Markov network is show in the figure for a 3 cell grid. is the observation model.
Murphy describes using RBPF for this purpose further in his paper^{[3]}. For this problem, if is known then that they can marginalize out analytically and only sample the reduced state space .
This paper only tested the problem on a grid map which in comparison to real world problems, is unrealistically small. This experiment does show to give equivalent results as other methods while requiring few particles.
Paper 2: Improved Techniques for Grid Mapping With RaoBlackwellized Particle Filters
The original paper outlined the RBPF problem as a general solution for any problem that uses particle filters, even though they chose to use SLAM to evaluate the method. The second paper focuses only on the SLAM problem and improves the performance of the filter only for SLAM problems were the observation data is highly accurate as it would be if the robot was using a Laser Range Finder (LRF).
With this, they propose two approaches to improve RBPF:
 A proposal distribution that considers the accuracy of the observed readings
 An adaptive resampling technique that reduces the risk of particle depletion.
Improved Proposal Distribution
In the original RBPF algorithm, the particles are updated using the odometry motion model as the proposal distribution where is the odometry measurements at time t1. Then the weights of each particle are determined by both the sensor measurements and the odometry measurements. However, the sensor measurements are far more accurate than the odometry measurements, meaning particle weights can differ significantly from each other. The figure shows the difference in the models. To overcome this problem, they incorporate the sensor readings into the proposal such as where is the sensor readings at time t.
This model may then be hard to sample since is focused on such small area. To simplify the computation, they use the samples to estimate a Gaussian distribution with parameters
Where
This does limit the type of noise that can be expected from the sensor though, which was one of the main advantages of using a particle filter, but since they are only focusing on LRF sensors, this is not an unreasonable assumption.
Adaptive sampling
Since the samples are now highly focused on one area, they now must deal with the risk of particle depletion. This occurs when particles become overconfident in their positions and are incorrect due to poor observations or too small of an overlapping area between successive scans.
To mitigate this occurrence, they propose adding a resampling step based on a parameter such that
where are normalized weights. This means that as the weights grow and become more confident, will drop. Once drops below the threshold N/2 where N is the number of particles, then they will resample the particles. In the extensive testing done, this was shown to help decrease the risk of particle depletion.
Evaluation
To evaluate this improved RBFP, testing was done are three different locations: Intel Reseach Lab, Freibugrg Campus, and MIT Killian Court. These location varied in size, moving obstacles, and types of locations environments such as indoor vs outdoor maps. Each environment was tested multiple times with different sizes of particle sets to see the map accuracy that was able to be achieved.
Outdoor environments proved to be an issue since LRF have a limited range and open areas such as a parking lot caused the robot to get lost.
In all these experiments, the number of particles needed decrease by approximately one order of magnitude. As expected, smaller maps with fewer moving obstacles required the lease particles and larger maps required more. This was tested against the leading SLAM algorithm at the time known as FastSLAM ^{[4]}.
Conclusion
The difference between the two papers shows the transition from theory to application. The first paper proposed the RBPF providing a large amount of mathematical background and performing 'toy' experiments to show that the theory did have promise. The second paper took the theory and looked at what could be done to improve it in the real world environments. Doing this came with its own limitations such that this approach would not transfer easily to another similar application if it does not fit model given. This means that someone may choose to create their own improved RBPF for their own system and disregard the approaches proposed by the second paper. For further studies, it would be interesting to see if the RBPF could be improved for a general case without adding application specific limitation.
Annotated Bibliography
 ↑ Arnaud Doucet, Nando de Freitas, Kevin Murphy, and Stuart Russell. 2000. Raoblackwellised particle filtering for dynamic Bayesian networks. In Proceedings of the Sixteenth conference on Uncertainty in artificial intelligence (UAI'00), Craig Boutilier and Moisés Goldszmidt (Eds.). Morgan Kaufmann Publishers Inc., San Francisco, CA, USA, 176183.
 ↑ G. Grisetti, C. Stachniss and W. Burgard, "Improved Techniques for Grid Mapping With RaoBlackwellized Particle Filters," in IEEE Transactions on Robotics, vol. 23, no. 1, pp. 3446, Feb. 2007. doi: 10.1109/TRO.2006.889486
 ↑ Kevin P. Murphy. 1999. Bayesian map learning in dynamic environments. In Proceedings of the 12th International Conference on Neural Information Processing Systems (NIPS'99), S. A. Solla, T. K. Leen, and K. Müller (Eds.). MIT Press, Cambridge, MA, USA, 10151021.
 ↑ D. Hahnel, W. Burgard, D. Fox and S. Thrun, "An efficient fastSLAM algorithm for generating maps of largescale cyclic environments from raw laser range measurements," Proceedings 2003 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS 2003) (Cat. No.03CH37453), 2003, pp. 206211 vol.1.
