一步一步学习S-MSCKF(五)观测更新

1 单个特征对单个双目的观测模型

考虑一个特征点 f j f_j fj被位置在 ( G C i q , G p C i ) (^{C_i}_Gq,{^{G}p}_{C_i}) (GCiq,GpCi)的双目相机观测到。双目的左右相机位姿可以表示为 ( G C i , 1 q , G p C i , 1 ) (^{C_{i,1}}_Gq,{^{G}p_{C_{i,1}}}) (GCi,1q,GpCi,1) ( G C i , 2 q , G p C i , 2 ) (^{C_{i,2}}_Gq,{^{G}p_{C_{i,2}}}) (GCi,2q,GpCi,2)。虽然状态向量仅仅包含左相机的位姿,但是相机左目与右目的位姿可以通过双目相机的外参产生联系。

双目的观测值(归一化坐标)可以表达为:
z i j = [ u i , 1 j v i , 1 j u i , 2 j v i , 2 j ] = [ 1 C i , 1 Z j I 2 × 2 0 2 × 2 0 2 × 2 1 C i , 2 Z j I 2 × 2 ] [ C i , 1 X j C i , 1 Y j C i , 2 X j C i , 2 Y j ] + n i j z^{j}_i= \left[\begin{matrix} u_{i,1}^{j} \\ v_{i,1}^{j} \\ u_{i,2}^{j} \\ v_{i,2}^{j} \end{matrix}\right]= \left[\begin{matrix} \frac{1}{^{C_{i,1}}Z_j}I_{2\times 2} & 0_{2\times 2} \\ 0_{2\times 2} & \frac{1}{^{C_{i,2}}Z_j}I_{2\times 2} \end{matrix}\right] \left[\begin{matrix} {^{C_{i,1}}X_j} \\ {^{C_{i,1}}Y_j} \\ {^{C_{i,2}}X_j} \\ {^{C_{i,2}}Y_j} \end{matrix}\right]+n_i^j zij=ui,1jvi,1jui,2jvi,2j=[Ci,1Zj1I2×202×202×2Ci,2Zj1I2×2]Ci,1XjCi,1YjCi,2XjCi,2Yj+nij
左右目的观测值与相机的位姿的关系为:
C i , 1 p j = [ C i , 1 X j C i , 1 Y j C i , 1 Z j ] = C ( G C i , 1 q ) ( G p j − G p C i , 1 ) ^{C_{i,1}}p_j=\left[\begin{matrix}{^{C_{i,1}}X_j} \\ {^{C_{i,1}}Y_j} \\ {^{C_{i,1}}Z_j}\end{matrix}\right]=C({^{C_{i,1}}_Gq})({^Gp_j}-{^Gp_{C_{i,1}}}) Ci,1pj=Ci,1XjCi,1YjCi,1Zj=C(GCi,1q)(GpjGpCi,1)
C i , 2 p j = [ C i , 2 X j C i , 2 Y j C i , 2 Z j ] = C ( G C i , 2 q ) ( G p j − G p C i , 2 ) = C ( C i , 1 C i , 2 q ) ( C i , 1 p j − C i , 1 p C i , 2 ) ^{C_{i,2}}p_j=\left[\begin{matrix}{^{C_{i,2}}X_j} \\ {^{C_{i,2}}Y_j} \\ {^{C_{i,2}}Z_j}\end{matrix}\right]=C({^{C_{i,2}}_Gq})({^Gp_j}-{^Gp_{C_{i,2}}})= C(^{C_{i,2}}_{C_{i,1}}q)(^{C_{i,1}}p_j-{^{C_{i,1}}p_{C_{i,2}}}) Ci,2pj=Ci,2XjCi,2YjCi,2Zj=C(GCi,2q)(GpjGpCi,2)=C(Ci,1Ci,2q)(Ci,1pjCi,1pCi,2)
从上面可以发现右目观测量可以用左目位姿来表示。

z i j z_i^j zij在估计值 ( x ^ C i , G p ^ j ) ({\hat x_{C_i}},{^G\hat p_j}) (x^Ci,Gp^j)处泰勒展开得到:
z i j = [ u ^ i , 1 j v ^ i , 1 j u ^ i , 2 j v ^ i , 2 j ] + ∂ z i j ∂ x C i , 1 x ~ C i + ∂ z i j ∂ G p j G p ~ j + n i j z_i^j= \left[\begin{matrix} \hat u_{i,1}^{j} \\ \hat v_{i,1}^{j} \\ \hat u_{i,2}^{j} \\ \hat v_{i,2}^{j} \end{matrix}\right]+ \frac{\partial z^j_i}{\partial x_{C_{i,1}}}\tilde x_{C_i}+ \frac{\partial z^j_i}{\partial {^Gp_j}}{^G \tilde p_j}+n_i^j zij=u^i,1jv^i,1ju^i,2jv^i,2j+xCi,1zijx~Ci+GpjzijGp~j+nij
所以测量残差约为:
r i j = z i j − z ^ i j = ∂ z i j ∂ x C i , 1 x ~ C i + ∂ z i j ∂ G p j G p ~ j + n i j = H C i j x ~ C i + H f i j G p ~ j + n i j r^j_i= z_i^j-{\hat z^j_i}=\frac{\partial z^j_i}{\partial x_{C_{i,1}}}\tilde x_{C_i}+ \frac{\partial z^j_i}{\partial {^Gp_j}}{^G \tilde p_j}+n_i^j= H^j_{C_i}\tilde x_{C_i}+H^j_{f_i}{^G\tilde p_j}+n^j_i rij=zijz^ij=xCi,1zijx~Ci+GpjzijGp~j+nij=HCijx~

  • 4
    点赞
  • 9
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值