目录
Abstract
本文提出了一种称名为激光雷达-惯性导航仪-相机三合一(LIC-Fusion)的紧耦合多传感器融合算法,该算法可有效地融合IMU测量,稀疏的视觉特征和提取的激光雷达点。尤其是,提出的LIC-Fusion能够在三个异步传感器之间执行在线的空间和时间传感器标定,以补偿可能的标定偏差。
本方法的关键贡献在于,能够基于高效的MSCKF框架,利用从激光雷达检测和跟踪到的系数边缘/surf特征点,以及稀疏的视觉特征与IMU读数,完成理想的(根据线性误差)多模态传感器融合。
我们在室内和室外环境中均进行了广泛的实验,结果表明,所提出的LIC-Fusion在估计精度和对剧烈运动的鲁棒性方面均优于最新的视觉惯性里程计(VIO)和激光雷达里程计方法。
Ⅰ. INTRODUCTION AND RELATED WORK
- 惯性导航系统(INS)辅助单目摄像机对于跟踪自动驾驶车辆和移动感知系统的三维运动效果很好------>互补传感模式、低成本和小尺寸
- 相机受到光照条件的限制,在低光或夜间条件下无法提供高质量的信息
- 3D激光雷达传感器无论在何种光照条件下都能提供更稳健、更精确的距离测量---->存在点云稀疏、成本高、采集率低的问题
- 惯性测量单元(IMU)测量局部角速度和线性加速度,可以在动态轨迹中提供大量的信息,但如果不与其他信息融合,则会由于噪声而产生较大的漂移
- 所以,着重于激光-惯性-视觉里程计(LIC)。
本文在计算效率高的多状态约束卡尔曼滤波(MSCKF)框架下,提出了一种快速、紧耦合、单线程的激光-惯性-视觉(LIC)里程计算法,该算法具有在线时空多传感器标定功能。这项工作的主要贡献如下:
1.开发了一种紧密耦合的LIC里程计技术(称为LIC-Fusion),通过在线空间和时间校准实现了高效的6自由度姿态估计 。在MSCKF框架内,所提出的LIC-Fusion有效地结合了IMU测量值、稀疏视觉特征和两种不同的稀疏激光雷达特征(见图1)。明确地建立了标定的外部参数和对测量值的姿态估计两者之间的依赖关系模型,并进行了解析推导。
2. 对所提出的方法进行了大量的实验验证,包括室内和室外环境的真实实验,表明所提出的LIC-Fusion比最先进的方法更准确、更鲁棒。
![](https://i-blog.csdnimg.cn/blog_migrate/4ee90855322aaec406749081cd6acaad.png)
图1:LIC-Fusion方法融合了图像中所跟踪的稀疏视觉特征和点云中提取的激光雷达特征。红色和蓝色的激光雷达点分别是边缘和平面特征。估计的轨迹用绿色标记。
II. THE PROPOSED LIC-FUSION
A. State Vector
- 状态向量
- k时刻,IMU的状态: x I x_I xI
- IMU和相机之间的转换外参: x c a l i b _ C x_{calib\_C} xcalib_C
- IMU和激光雷达之间的转换外参: x c a l i b _ L x_{calib\_L} xcalib_L
- 克隆的滑动窗口,包括过去m个图像时间的局部IMU克隆状态 x C x_C xC和过去n个激光雷达扫描时间的局部IMU克隆状态 x L x_L xL
总状态向量为:
其中
G
I
k
q
^{I_k}_G q
GIkq表示在
t
k
t_k
tk时刻,IMU从世界坐标系{G}到IMU局部坐标系
I
k
{I_k}
Ik的旋转。
G
v
I
k
^Gv_{I_k}
GvIk和
G
p
I
k
^Gp_{I_k}
GpIk代表在世界坐标系中,IMU的速度和
t
k
t_k
tk时刻的位置。
b
g
b_g
bg和
b
a
b_a
ba是陀螺仪和加速度计的偏置。
I
C
q
ˉ
^C_I \bar{q}
ICqˉ和
C
p
I
^Cp_I
CpI表示相机坐标系{C}和IMU坐标系{I}之间的刚体变换。
I
L
q
ˉ
^L_I \bar{q}
ILqˉ和
L
p
I
^Lp_I
LpI就是激光雷达和IMU坐标系之间的三维刚性转换。
- 估计外感传感器和IMU之间的时间偏移量:
IMU通常存在于低成本设备中。以IMU时间为基准时间,将相机和激光雷达都建模为具有偏移 t d C t_{dC} tdC和 t d L t_{dL} tdL的基钟,可以修正如下的测量时间:
![](https://i-blog.csdnimg.cn/blog_migrate/daa182273cb04a61e13189107b685246.png)
t C t_C tC和 t L t_L tL分别是相机和激光雷达报告的时间。
- 状态误差的定义
我们定义状态的真实值为 x x x,状态估计值为 x ^ \hat{x} x^,对应的状态误差为 δ x δx δx,则其关系为:(都是些比较基础的知识了)
![](https://i-blog.csdnimg.cn/blog_migrate/b793cd7b1c544508a163fe20b5eecad6.png)
其操作是对于向量空间中的状态 v v v是简单的欧几里德加法,比如 v = v ^ + δ v v = \hat{v}+δv v=v^+δv,对于四元数有:
![](https://i-blog.csdnimg.cn/blog_migrate/8709d03978e7cbdd60d39d6d2217034f.png)
其操作为四元数乘法
B. IMU Propagation
IMU提供角速率和线性加速度测量,使用以下连续时间运动学模型进行建模:(以下基础知识应该比较熟悉)
![](https://i-blog.csdnimg.cn/blog_migrate/cffdda1a781a4e1e34b816161691a2db.png)
其中,
I k ω ^{I_k} ω Ikω 和 I k a ^{I_k} a Ika表示局部IMU坐标系下的真实角速度和线性加速度。 G g ^Gg Gg表示世界坐标系下的重力加速度。陀螺仪和加速度计的 b g b_g bg和 b a b_a ba分别受高斯白噪声 n w g n_{wg} nwg和 n w a n_{wa} nwa的影响。
然后对这个连续时间系统进行积分和线性化,使状态协方差矩阵在时间上向前传播。
C. State Augmentation
当系统接收到新的图像或激光雷达扫描时,IMU状态将向前传播到该时间瞬间,并将传播的惯性状态复制(克隆)到 x C x_C xC或 x L x_L xL状态向量中。为了校准不同传感器之间的时间偏移,我们将一直传播到IMU时间 t ^ I k \hat{t}_{I_k} t^Ik,这是目前对IMU时钟中测量采集时间的最佳估计。例如,如果以时间戳 t L k t_{L_k} tLk接收到一个新的激光雷达扫描,我们将一直传播到 t ^ I k = t L k + t ^ d L \hat{t}_{Ik} = t_{Lk}+\hat{t}_{dL} t^Ik=tLk+t^dL,并增加状态向量 x L x_L xL来包括这个新的克隆状态估计。
![](https://i-blog.csdnimg.cn/blog_migrate/44122371fbc598c01ff6bf2b4bcbc8ad.png)
![](https://i-blog.csdnimg.cn/blog_migrate/71919d65f416203a3e152e99747b6a52.png)
其中,J是新克隆的 x ^ L k \hat{x}_{Lk} x^Lk对当前状态的雅可比矩阵:
![](https://i-blog.csdnimg.cn/blog_migrate/af3797f3b9ca693677c27b2a6dfa8cbf.png)
在上面的表达式中, J I J_I JI是关于IMU状态 x I x_I xI的雅可比矩阵:
![](https://i-blog.csdnimg.cn/blog_migrate/6441c50449c3b79d4e5c0a60743eeccb.png)
J c a l i b L J_{calib_L} JcalibL为IMU与激光雷达关于外参(包括时间偏移)方面的雅可比矩阵:
![](https://i-blog.csdnimg.cn/blog_migrate/628f9414bf60e13ac5d77ad6b2e75705.png)
同理, J c a l i b C J_{calib_C} JcalibC、 J C J_C JC、 J L J_L JL分别是IMU与camera、camera时刻克隆、LiDAR时刻克隆在外参方面的雅可比,在本例中为0。
要注意的是,与激光雷达测量相对应的新克隆IMU状态对 t d L t_{dL} tdL的依赖性是通过IMU运动学建模的,因此我们的测量模型(见第II-D.1节)可以直接作为IMU时钟帧中真实测量时间的克隆函数。这个激光雷达克隆过程类似于一个新的相机测量过程。
D. Measurement Models
1. 激光雷达特征测量:
为了限制所需的计算成本,我们希望从原始激光雷达扫描中选择一组稀疏的高质量特征进行状态估计。
因此,提取激光雷达扫描环的高曲率部分和低曲率部分,分别对应于边缘和平面surf特征(见图1)。我们使用KD-树快速索引,将当前激光雷达扫描中提取的边缘和surf特征投影到之前的扫描中,并找到最近的对应特征。
例如,我们在激光雷达扫描{ L l + 1 L_{l+1} Ll+1}到{ L l L_l Ll}中投影一个特征点 L l + 1 p f i ^{L_{l+1}}p_{f i} Ll+1pfi,这个投影点记作 L l p f ^{L_l}p_f Llpf:
![](https://i-blog.csdnimg.cn/blog_migrate/936e486764185e174644b49407aef2b4.png)
其中 L l + 1 L l R ^{L_l}_{L_l+1}R Ll+1LlR和 L l p L l + 1 ^{L_l}p_{L_{l+1}} LlpLl+1是两个激光雷达帧间的相对旋转和平移,可由状态向量中的状态计算得到:
![](https://i-blog.csdnimg.cn/blog_migrate/8e3c3d6458e6edb739f227d83cfebdff.png)
跟踪之后,我们会在原先的扫描中找到两个边缘特征 L l p f j ^{L_l}p_{fj} Llpfj, L l p f k ^{L_l}p_{fk} Llpfk,对应于投影的边缘特征 L l p f i ^{L_l}p_{fi} Llpfi。我们假设它们是从与 L l p f i ^{L_l}p_{fi} Llpfi相同的物理边缘采样的。如果最近的边缘特征 L l p f j ^{L_l}p_{fj} Llpfj在第 r r r个扫描环上,则第二个最近的边缘特征 L l p f k ^{L_l}p_{fk} Llpfk应该在其邻近的环 r − 1 r-1 r−1或 r + 1 r+1 r+1上。因此,边缘特征 L l + 1 p f i ^{L_{l+1}}p_{fi} Ll+1pfi的测量残差是其投影特征点 L l p f i ^{L_l}p_{fi} Llpfi与 L l p f j ^{L_l}p_{fj} Llpfj和 L l p f k ^{L_l}p_{fk} Llpfk形成的直线之间的距离。
![](https://i-blog.csdnimg.cn/blog_migrate/5442151d571b6764326dd60569735cf0.png)
其中 ∣ ∣ . ∣ ∣ 2 ||.||_2 ∣∣.∣∣2为欧几里德范数,×为两个向量的外积。
在当前状态估计时,我们将上述边缘特征的距离测量线性化:
![](https://i-blog.csdnimg.cn/blog_migrate/63f592abdcdfa35ea439417ad9519933.png)
其中
H
x
H_x
Hx为距离相对于状态向量中的状态的雅可比矩阵,
n
r
n_r
nr为方差是
C
r
C_r
Cr的高斯白噪声。
H
x
H_x
Hx中的非零元素只与克隆的poses:
以及IMU和激光雷达之间的严格校准:
相关,所以得到:
为了进行EKF更新,我们需要知道距离测量的显式协方差
C
r
C_r
Cr。由于这种测量不是直接从激光雷达传感器中获得的,我们在激光雷达扫描
C
r
C_r
Cr中传播原始测量(点)的协方差。假设
L
l
+
1
p
f
i
^{L_{l+1}}p_{fi}
Ll+1pfi,
L
l
p
f
j
^{L_l}p_{fj}
Llpfj,
L
l
p
f
k
^{L_l}p_{fk}
Llpfk点对应的协方差分别为
C
i
,
C
j
,
C
k
C_i, C_j , C_k
Ci,Cj,Ck;则
C
r
C_r
Cr可如下计算:
![](https://i-blog.csdnimg.cn/blog_migrate/eaa4584f5ffec2043fb5d7f5f32f6a8e.png)
我们执行简单的基于Mahalanobis距离的概率离群拒绝:
![](https://i-blog.csdnimg.cn/blog_migrate/73ecc25abc720470e6cf8c1ce07bca1a.png)
其中, P x P_x Px表示相关状态的协方差矩阵, r m r_m rm应该服从于 χ 2 χ^2 χ2卡方分布,因此,如果它通过了这个测试, r ( L l + 1 p f i ) r(^{L_{l+1}} p_{fi}) r(Ll+1pfi)将在EKF更新中使用。
同样地,为了投影surf特征 L l p f i ^{L_l}p_{fi} Llpfi,我们将发现三个相应的surf特征 L l p f j ^{L_l}p_{fj} Llpfj、 L l p f k ^{L_l}p_{fk} Llpfk、 L l p f l ^{L_l}p_{fl} Llpfl,假设这些特征都是在与 L l p f i ^{L_l}p_{fi} Llpfi相同的物理平面上进行采样的。surf特征 L l + 1 p f i ^{L_{l+1}}p_{fi} Ll+1pfi的测量残差为其投影特征点 L l p f i ^{L_l}p_{fi} Llpfi与 L l p f j ^{L_l}p_{fj} Llpfj、 L l p f k ^{L_l}p_{fk} Llpfk、 L l p f l ^{L_l}p_{fl} Llpfl构成的平面之间的距离。距离测量的协方差传播特性,线性化和Mahalanobis距离测试与边缘特征相似。
2. 视觉特征测量:
给定一个新图像,我们同样地传播和增益状态。使用KLT光流从图像中提取FAST特征并跟踪到未来帧中。一旦某个视觉特征丢失或在整个滑动窗口中被跟踪,我们使用相机克隆的当前估计值在3D空间中对该特征进行三角化测量。更新中使用了标准的视觉特征重投影误差。给定某个三维视觉特征点 G p f i ^Gp_{fi} Gpfi 的一组特征方位测量值 z i z_i zi,一般的线性化残差为:
![](https://i-blog.csdnimg.cn/blog_migrate/c8b5e8c43640d031ce3773eb2bb0dd7e.png)
其中 H f H_f Hf是视觉特征测量相对于3D特征 G p f i ^Gp_{fi} Gpfi的雅可比矩阵,并且两个雅可比矩阵都是按当前的最佳估计值计算的。因为我们的测量值是关于 G p f i ^Gp_{fi} Gpfi的函数(公式29),我们利用MSCKF零空间投影去除这种依赖性。在零空间投影之后,我们有:
![](https://i-blog.csdnimg.cn/blog_migrate/cdb2bd31c9aab7478431170d0d4051f7.png)
需要注意的是,IMU和camera{ I C q ˉ , C p I ^C_I\bar{q}, ^Cp_I ICqˉ,CpI}之间的刚性变换的雅可比矩阵是非0的,也就是说IMU和camera之间的变换可以在线校准。
E. Measurement Compression
在当前状态估计下对激光雷达特征和视觉特征测量进行线性化之后,我们可以执行EKF更新,但是由于大量的视觉和激光雷达特征测量,这会带来很大的计算成本。考虑所有测量残差和雅可比(来自激光雷达或视觉特征)的叠加:
![](https://i-blog.csdnimg.cn/blog_migrate/250c5a0b6cd862505afa693a9c441546.png)
其中 r r r和 n n n是带有残差和噪声块元素的向量(公式25/30),通常假设所有的测量数据在统计上是独立的,那么噪声向量 n n n就不相关。为了降低计算复杂度,我们使用Givens旋转来执行thin QR来压缩测量值,即
![](https://i-blog.csdnimg.cn/blog_migrate/6db3532853eee2a6fcae0eebefa74c05.png)
其中 Q H 1 Q_{H1} QH1和 Q H 2 Q_{H2} QH2是酉矩阵。经过测量压缩得到:
![](https://i-blog.csdnimg.cn/blog_migrate/e8c441bc17d6b38a301269df6ded61b8.png)
其中压缩后的雅可比矩阵 T H T_H TH应该是状态向量x的维数的平方,压缩后的噪声由 n c = Q H 1 T n n_c = Q^T_{H_1}n nc=QH1Tn给出。然后利用这个压缩的线性测量残差来有效地更新状态估计和与标准EKF的一致性。
III. EXPERIMENTAL RESULTS
-
实验设备
自组装的激光-惯性-相机,由激光雷达,Xsens IMU,和单目相机组装而成 -
部分实验结果
IV. CONCLUSIONS AND FUTURE WORK
在本文中,我们开发了一种紧密耦合的高效多模态传感器融合算法,用于激光-惯性-相机里程计。在MSCKF框架内。在三个传感器之间进行在线时空校准,以补偿校准灵敏度,并简化传感器部署。该方法通过激光雷达扫描检测和跟踪稀疏边缘和平面特征点,并融合这些测量值和从单目图像中提取的视觉特征。因此,通过利用不同的传感模式,所提出的LIC-Fusion能够在不同的环境和主动运动下提供精确和鲁棒的6自由度运动跟踪。在未来,我们将研究如何有效地集成从激光雷达和摄像机获得的环路闭合约束,以约束导航误差。