标题:Bayesian Generalized Kernel Inference forTerrain Traversability Mapping
作者:Tixiao Shan,Jinkun Wang,Brendan Englot and Kevin Doherty
来源:CoRL 2018
代码:
https://github.com/TixiaoShan/traversability_mapping/tree/master/traversability_mapping
文章目录
1. 主要工作
本文提出贝叶斯广义核推理来解决可通行区域建图问题:
【1】将BGK高程推断应用于解决地形建图中遇到的稀疏数据问题;
【2】仅对选定位置的高程数据进行可通行性计算,从而减轻了计算负担;
【3】通过BGK可通行性推断来估计其他位置的可通行性。
这个框架可以使用稀疏激光雷达数据和与小型UGV兼容的硬件来实时建立可通行区域地图,是贝叶斯广义核推断在可通行性地图上的第一次应用。
2. 系统流程
3. 基于贝叶斯广义核推理的可通行区域建图
3.1 贝叶斯广义核推理
给定观测: D = { ( x i , y i ) i = 1 : N } \mathcal{D}=\left\{\left(\mathbf{x}_{i}, y_{i}\right)_{i=1: N}\right\} D={
(xi,yi)i=1:N},试图推断在潜空间上参数化的目标值的概率分布:
Θ : p ( y ∗ ∣ x ∗ , D ) ∝ ∫ p ( y ∗ ∣ θ ∗ ) p ( θ ∗ ∣ x ∗ , D ) d θ ∗ \Theta: p\left(y^{*} \mid \mathbf{x}^{*}, \mathcal{D}\right) \propto \int p\left(y^{*} \mid \boldsymbol{\theta}^{*}\right) p\left(\boldsymbol{\theta}^{*} \mid \mathbf{x}^{*}, \mathcal{D}\right) d \boldsymbol{\theta}^{*} Θ:p(y∗∣x∗,D)∝∫p(y∗∣θ∗)p(θ∗∣x∗,D)dθ∗
其中,
p ( θ ∗ ∣ x ∗ , D ) ∝ ∫ θ 1 ⋅ N ∏ i = 1 N p ( y i ∣ θ i ) p ( θ 1 : N , θ ∗ ∣ x 1 : N , x ∗ ) d θ 1 : N p\left(\boldsymbol{\theta}^{*} \mid \mathbf{x}^{*}, \mathcal{D}\right) \propto \int_{\boldsymbol{\theta}_{1 \cdot N}} \prod_{i=1}^{N} p\left(y_{i} \mid \boldsymbol{\theta}_{i}\right) p\left(\boldsymbol{\theta}_{1: N}, \boldsymbol{\theta}^{*} \mid \mathbf{x}_{1: N}, \mathbf{x}^{*}\right) d \boldsymbol{\theta}_{1: N} p(θ∗∣x∗,D)∝∫θ1⋅N∏i=1Np(yi∣θi)p(θ1:N,θ∗∣x1:N,x∗)dθ1:N,
为与目标输入相关的潜参数的后验分布。
在高斯过程中,假设所有参数 θ \theta θ 是相关的,而在贝叶斯广义核推理中,假设在给定目标输入的情况下,与观测输入相关的参数条件独立:
p ( θ 1 : N , θ ∗ ∣ x 1 : N , x ∗ ) = p\left(\boldsymbol{\theta}_{1: N}, \boldsymbol{\theta}^{*} \mid \mathbf{x}_{1: N}, \mathbf{x}^{*}\right)= p(θ1:N,θ∗∣x1:N,x∗)= ∏ i = 1 N p ( θ i ∣ x i , θ ∗ , x ∗ ) p ( θ ∗ ∣ x ∗ ) , \prod_{i=1}^{N} p\left(\boldsymbol{\theta}_{i} \mid \mathbf{x}_{i}, \boldsymbol{\theta}^{*}, \mathbf{x}^{*}\right) p\left(\boldsymbol{\theta}^{*} \mid \mathbf{x}^{*}\right), ∏i=1Np(θi∣xi,θ∗,x∗)p(θ∗∣x∗),
将其代入上式可得:
p ( θ ∗ ∣ x ∗ , D ) ∝ p\left(\boldsymbol{\theta}^{*} \mid \mathbf{x}^{*}, \mathcal{D}\right) \propto p(θ∗∣x∗,D)∝ ∏ i = 1 N p ( y i ∣ θ ∗ , x ∗ , x i ) p ( θ ∗ ∣ x ∗ ) \prod_{i=1}^{N} p\left(y_{i} \mid \boldsymbol{\theta}^{*}, \mathbf{x}^{*}, \mathbf{x}_{i}\right) p\left(\boldsymbol{\theta}^{*} \mid \mathbf{x}^{*}\right) ∏i=<