CN108362288A - 一种基于无迹卡尔曼滤波的偏振光slam方法 - Google Patents
一种基于无迹卡尔曼滤波的偏振光slam方法 Download PDFInfo
- Publication number
- CN108362288A CN108362288A CN201810128645.1A CN201810128645A CN108362288A CN 108362288 A CN108362288 A CN 108362288A CN 201810128645 A CN201810128645 A CN 201810128645A CN 108362288 A CN108362288 A CN 108362288A
- Authority
- CN
- China
- Prior art keywords
- coordinate system
- polarized light
- uav
- time
- landmark
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000005259 measurement Methods 0.000 claims abstract description 44
- 238000001914 filtration Methods 0.000 claims abstract description 13
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 6
- 239000011159 matrix material Substances 0.000 claims description 35
- 238000004364 calculation method Methods 0.000 claims description 9
- 230000010287 polarization Effects 0.000 claims description 8
- 230000001133 acceleration Effects 0.000 claims description 6
- 230000007704 transition Effects 0.000 claims description 6
- 230000009466 transformation Effects 0.000 claims description 4
- 238000013461 design Methods 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 238000010276 construction Methods 0.000 abstract description 4
- 230000008569 process Effects 0.000 description 4
- 230000007613 environmental effect Effects 0.000 description 3
- 238000013507 mapping Methods 0.000 description 3
- 238000009825 accumulation Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 102000008115 Signaling Lymphocytic Activation Molecule Family Member 1 Human genes 0.000 description 1
- 108010074687 Signaling Lymphocytic Activation Molecule Family Member 1 Proteins 0.000 description 1
- 238000000149 argon plasma sintering Methods 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/20—Instruments for performing navigational calculations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Mathematical Physics (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Algebra (AREA)
- Geometry (AREA)
- Computing Systems (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Automation & Control Theory (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
- Navigation (AREA)
Abstract
本发明公开了一种基于无迹卡尔曼滤波的偏振光SLAM方法,属于无人机自主导航的范畴。该方法结合无人机的状态模型和基于激光雷达传感器、偏振光传感器的量测模型,通过无迹卡尔曼滤波,即UKF算法,实现无人机位置的确定和周围环境地图的构建,利用偏振光信息和激光雷达信息匹配互补、不受其他外界干扰的特性,提高了无人机SLAM系统的稳定性和精度。
Description
技术领域
本发明涉及无人机同时定位与构图(SLAM)属于无人机自主导航的范畴,具体涉及一种基于无迹卡尔曼滤波的偏振光SLAM方法,对无人机如何确定自身位置及感知外界环境的问题,SLAM系统旨在通过无人机系统模型,结合相应的滤波方法完成无人机的定位与周围环境的绘制。
背景技术
SLAM是Simultaneous Localization and Mapping缩写,意为“同时定位与建图”。它是指运动体根据传感器的信息,一边计算自身位置,一边构建环境地图的过程。目前,SLAM技术已经被运用于无人机、无人驾驶、机器人、AR、智能家居等领域。
SLAM研究侧重于使用滤波器理论,最小化运动体位姿和地图的路标点的噪声,一般采用里程计的输入作为预测过程的输入,激光传感器的输出作为更新过程的输入。扩展卡尔曼滤波(Extended Kalman Filter,EKF)滤波算法是目前大多数SLAM采用的滤波算法,缺点是线性化处理时需要用雅克比(Jacobian)矩阵,其繁琐的计算过程导致该方法实现相对困难,而UKF使用无迹(UT)变换来处理均值和协方差的非线性传递,不需要用求解雅克比矩阵,精度较高。近年来无人机的发展非常迅速,而将SLAM导航用于无人机可以解决无人机自身位置不确定,航迹推算中位置误差累积的情况,利用自身携带的传感器,反复探测环境中的特征,从而完成自身位置及特征位置的校正,同时构建环境地图,无需预知地图信息或依靠外部辅助设备,即可完成无人机位置信息和周围环境地图的构建,但是这种SLAM方案的定位精度和构图精确度都不高,鲁棒性很差。
近年来对偏振光的研究越来越多,1871年英国著名物理学家瑞利提出了瑞利散射定律,揭示了光线散射特性,随后人们基于瑞利散射定律获得了全天域大气偏振分布模式。大气偏振分布模式相对稳定,其中蕴涵着丰富的导航信息,自然界中很多生物都能够利用天空偏振光进行导航或辅助导航。偏振光导航机制是一种非常有效的导航手段,具有无源、无辐射、隐蔽性好等特点,能够为复杂环境下的导航任务提供新的解决途径。
将偏振光用于无人机SLAM上可以提高无人确定自身位置的和构建周围环境地图的精确度,解决鲁棒性不高的问题。
发明内容
本发明主要解决的问题是:将自然界中广泛存在的偏振光信息应用到无人机SLAM中以解决无人机同时定位与构图中存在的自身位置确定困难、环境适应性差、构图不够精确的问题。
本发明采用的技术方案为:
一种基于无迹卡尔曼滤波的偏振光SLAM方法,包括以下步骤:
(1)选取无人机的姿态、速度、位置和路标点的位置为系统状态,建立无人机的动力学模型;
(2)建立激光雷达的量测模型;
(3)建立偏振光传感器的量测模型;
(4)系统初始化、地图初始化;
(5)路标点匹配;
(6)利用路标点的激光雷达数据和偏振传感器数据,设计UKF滤波器,估计无人机位置、速度、姿态和路标点的位置;
(7)地图更新;
其中,所述步骤(1)选取无人机的姿态、速度、位置和路标点为系统状态,建立无人机的动力学模型;以无人机起始位置为原点的世界坐标系,即w系,以正北方向为世界坐标系的x轴的正方向,以正西方向为世界坐标系的y轴的正方向,根据右手准则确定世界坐标系的z轴正方向;建立以无人机机体中心为原点的机体坐标系,即b系,以平行于机身纵轴指向机头方向为机体坐标系x轴的正方向,以平行于机身横轴指向左方为机体坐标系y轴的正方向,根据右手准则确定机体坐标系z轴正方向;选取无人机的位置、速度、姿态和路标点的位置作为系统状态;则系统状态转移方程如下:
其中
表示k-1时刻的系统状态,表示k-1时刻无人机在世界坐标系下的坐标,分别表示k-1时刻无人机在世界坐标系下x,y,z三个轴方向的坐标,表示k-1时刻无人机在世界坐标系下的速度,分别表示k-1时刻无人机在世界坐标系下指向x,y,z三个轴方向的速度,表示无人机的姿态,以欧拉角的形式表示,分别为滚转角、俯仰角和偏航角,表示k-1时刻世界坐标系下的m个路标点的坐标,表示第i,i=1,2,...,m个路标点在世界坐标系下的坐标,为k-1时刻机体坐标系下的比力加速度,为k-1时刻机体坐标系下的角速度,为k-1时刻世界坐标下的重加速度,△t为采样时间间隔,
为k-1时刻机体坐标系到世界坐标系的变换矩阵,nk-1为系统噪声,噪声设定为高斯白噪声,Qk-1为系统噪声协方差矩阵;
其中,所述步骤(2)建立激光雷达的量测模型;选取激光雷达测量的距离、角度作为测量值,针对激光雷达量测的数据给出路标点的量测模型:
其中
表示第k时刻对第i个路标点的量测方程中的非线性函数,yk,i,lidar=[dk,i θk,i λk,i]T表示k时刻系统的量测值,dk,i表示第k时刻激光雷达测量的无人机质心与第i个路标点的距离,θk,i表示第k时刻激光雷达测量的无人机质心与第i个路标点的距离俯仰角,λk,i表示第k时刻激光雷达测量的无人机质心与第i个路标点的方位角,uk,lidar表示激光雷达量测噪声,噪声设定为高斯白噪声,Rk,lidar为量测噪声协方差矩阵;
其中,所述步骤(3)建立偏振光传感器的量测模型:选取偏振光传感器量测的机体纵轴与太阳子午线的夹角量作为测值,针对偏振光给出的量测数据的观测模型如下:
其中为偏振光传感器获取的机体长轴与太阳子午线的夹角,为太阳方位角,为太阳高度角,为太阳赤纬,为观测点的纬度,t为太阳时角,uk,polar表示偏振光传感器量测噪声,噪声设定为高斯白噪声,Rk,polar为量测噪声协方差矩阵;
其中,所述步骤(4)系统初始化;将无人机坐标系转换到世界坐标系并将无人机,给定初始时刻的系统状态为初始协方差矩阵P0,并建立以无人机起始位置为原点的全局地图;
其中,所述步骤(5)路标点匹配;对k时刻激光雷达量测到的数据,和在全局地图中存储的路标点利用迭代最近点(ICP)算法进行匹配,对于匹配成功的路标点,则将该路标点送入相应的滤波器进行后续的运算;
其中,所述步骤(6)利用路标点的激光雷达数据和偏振光传感器数据,设计UKF滤波器,估计无人机的位置、速度、姿态和路标点的位置;首先将k-1时刻的估计状态和协方差矩阵,结合k-1时刻的IMU的输出数据,进行一步递推得到由IMU和状态转移矩阵推算出的机器人位置,然后根据预测值,加入激光传感器信息和偏振光信息进行量测更新的计算,输出无人机的位置、速度、姿态和路标点的位置,并同时输出相应的协方差矩阵,具体步骤如下:
时间更新
其中为状态向量维度,Pk-1|k-1为k-1时刻的误协方差矩阵,为尺度参数,计算方法为该式中ε为一小量取值范围为[10-4,1],常量可取0或者3-n;
量测更新
Kk=Pxy,k|k-1Pyy,k|k-1 -1
Pk|k=Pk|k-1-KkPyy,k|k-1Kk T
其中Pyy,k|k-1为自协方差矩阵,Pxy,k|k-1为互协方差矩阵;
其中,所述步骤(7)地图、状态更新;对于全局地图中已存在的路标点,将对应的滤波后的路标点作为最新路标点,对步骤(5)中未匹配成功的路标点,直接将其加入到全局地图中。
本发明与现有技术相比的优点在于:
本发明公开了一种基于无迹卡尔曼滤波的偏振光SLAM方法,属于无人机自主导航的范畴。该方法结合无人机的状态模型和基于激光雷达传感器、偏振光传感器的量测模型,通过无迹卡尔曼滤波算法实现无人机位置的确定和周围环境地图的构建,利用偏振光信息和激光雷达信息匹配互补、不受其他外界干扰的特性,提高了无人机SLAM系统的稳定性和精度。在保持激光雷达量测信息不变的情况下,利用偏振光量测误差不随时间累积的优点,降低了长时间情况下偏航角的误差的累积、提高了全局位姿的准确性。使无人机在位置环境中的定位与构图精度更高,对未知环境的适应能力增强。
附图说明
图1为一种基于无迹卡尔曼滤波的偏振光SLAM方法流程图;
图2为一种基于无迹卡尔曼滤波的偏振光SLAM方法仿真实验结果图。
具体实施方式
下面结合附图以及具体实施方式进一步说明本发明。
如图1所示,本发明一种基于无迹卡尔曼滤波的偏振光SLAM方法,包括以下步骤:
(1)选取无人机的姿态、速度、位置和路标点的位置为系统状态,建立无人机的动力学模型;
所述步骤(1)选取无人机的姿态、速度、位置和路标点为系统状态,建立无人机的动力学模型;以无人机起始位置为原点的世界坐标系,即w系,以正北方向为世界坐标系的x轴的正方向,以正西方向为世界坐标系的y轴的正方向,根据右手准则确定世界坐标系的z轴正方向;建立以无人机机体中心为原点的机体坐标系,即b系,以平行于机身纵轴指向机头方向为机体坐标系x轴的正方向,以平行于机身横轴指向左方为机体坐标系y轴的正方向,根据右手准则确定机体坐标系z轴正方向;选取无人机的位置、速度、姿态和路标点的位置作为系统状态;则系统状态转移方程如下:
其中
表示k-1时刻的系统状态,表示k-1时刻无人机在世界坐标系下的坐标,分别表示k-1时刻无人机在世界坐标系下x,y,z三个轴方向的坐标,表示k-1时刻无人机在世界坐标系下的速度,分别表示k-1时刻无人机在世界坐标系下指向x,y,z三个轴方向的速度,表示无人机的姿态,以欧拉角的形式表示,分别为滚转角、俯仰角和偏航角,表示k-1时刻世界坐标系下的m个路标点的坐标,表示第i(i=1,…,m)个路标点在世界坐标系下的坐标,为k-1时刻机体坐标系下的比力加速度,为k-1时刻机体坐标系下的角速度,为k-1时刻世界坐标下的重加速度,△t为采样时间间隔,
为k-1时刻机体坐标系到世界坐标系的变换矩阵,nk-1为系统噪声,噪声设定为高斯白噪声,Qk-1为系统噪声协方差矩阵;
(2)建立激光雷达的量测模型;
所述步骤(2)建立激光雷达的量测模型;选取激光雷达测量的距离、角度作为测量值,针对激光雷达量测的数据给出路标点的量测模型:
其中
表示第k时刻对第i个路标点的量测方程中的非线性函数,yk,i,lidar=[dk,i θk,i λk,i]T表示k时刻系统的量测值,dk,i表示第k时刻激光雷达测量的无人机质心与第i个路标点的距离,θk,i表示第k时刻激光雷达测量的无人机质心与第i个路标点的距离俯仰角,λk,i表示第k时刻激光雷达测量的无人机质心与第i个路标点的方位角,uk,lidar表示激光雷达量测噪声,噪声设定为高斯白噪声,Rk,lidar为量测噪声协方差矩阵;
(3)建立偏振光传感器的量测模型;
所述步骤(3)建立偏振光传感器的量测模型:选取偏振光传感器量测的机体纵轴与太阳子午线的夹角量作为测值,针对偏振光给出的量测数据的观测模型如下:
其中为偏振光传感器获取的机体长轴与太阳子午线的夹角,为太阳方位角,为太阳高度角,为太阳赤纬,为观测点的纬度,为太阳时角,uk,polar表示偏振光传感器量测噪声,噪声设定为高斯白噪声,Rk,polar为量测噪声协方差矩阵;
(4)系统初始化、地图初始化;
所述步骤(4)系统初始化;将无人机坐标系转换到世界坐标系并将无人机,给定初始时刻的系统状态为x0,初始协方差矩阵P0,并建立以无人机起始位置为原点的全局地图;
(5)路标点匹配;
所述步骤(5)路标点匹配;对k时刻激光雷达量测到的数据,和在全局地图中存储的路标点利用迭代最近点(ICP)算法进行匹配,对于匹配成功的路标点,则将该路标点送入相应的滤波器进行后续的运算;
(6)利用路标点的激光雷达数据和偏振传感器数据,设计UKF滤波器,估计无人机位置、速度、姿态和路标点的位置;
所述步骤(6)利用路标点的激光雷达数据和偏振光传感器数据,设计UKF滤波器,估计无人机的位置、速度、姿态和路标点的位置;首先将k-1时刻的估计状态和协方差矩阵,结合k-1时刻的IMU的输出数据,进行一步递推得到由IMU和状态转移矩阵推算出的机器人位置,然后根据预测值,加入激光传感器信息和偏振光信息进行量测更新的计算,输出无人机的位置、速度、姿态和路标点的位置,并同时输出相应的协方差矩阵,具体步骤如下:时间更新:
其中为状态向量维度,Pk-1|k-1为k-1时刻的误协方差矩阵,为尺度参数,计算方法为该式中ε为一小量取值范围为[10-4,1],常量可取0或者3-n;
量测更新:
Kk=Pxy,k|k-1Pyy,k|k-1 -1
Pk|k=Pk|k-1-KkPyy,k|k-1Kk T
其中Pyy,k|k-1为自协方差矩阵,Pxy,k|k-1为互协方差矩阵;
(7)地图更新;
所述步骤(7)地图、状态更新;对于全局地图中已存在的路标点,将对应的滤波后的路标点作为最新路标点,对步骤(5)中未匹配成功的路标点,直接将其加入到全局地图中。
如图2所示为仿真实验结果,*点为真实的无人机位姿,方框点为所提出方法估计的无人机位姿,五星点为设定的路标点。经过多次试验,无人机在x轴的平均误差为0.1602米,y轴平均误差为0.2554米,z轴平均误差为0.3532米,位置平均误差为0.5323米。
Claims (8)
1.一种基于无迹卡尔曼滤波的偏振光SLAM方法,其特征在于:包括以下步骤:
(1)选取无人机的姿态、速度、位置和路标点的位置为系统状态,建立无人机的动力学模型;
(2)建立激光雷达的量测模型;
(3)建立偏振光传感器的量测模型;
(4)系统初始化、地图初始化;
(5)路标点匹配;
(6)利用路标点的激光雷达数据和偏振传感器数据,设计UKF滤波器,估计无人机位置、速度、姿态和路标点的位置;
(7)地图更新。
2.根据权利要求1所述的一种基于无迹卡尔曼滤波的偏振光SLAM方法,其特征在于:所述步骤(1)选取无人机的姿态、速度、位置和路标点为系统状态,建立无人机的动力学模型;以无人机起始位置为原点的世界坐标系,即w系,以正北方向为世界坐标系的x轴的正方向,以正西方向为世界坐标系的y轴的正方向,根据右手准则确定世界坐标系的z轴正方向;建立以无人机机体中心为原点的机体坐标系,即b系,以平行于机身纵轴指向机头方向为机体坐标系x轴的正方向,以平行于机身横轴指向左方为机体坐标系y轴的正方向,根据右手准则确定机体坐标系z轴正方向;选取无人机的位置、速度、姿态和路标点的位置作为系统状态;则系统状态转移方程如下:
其中
表示k-1时刻的系统状态,表示k-1时刻无人机在世界坐标系下的坐标,分别表示k-1时刻无人机在世界坐标系下x,y,z三个轴方向的坐标,表示k-1时刻无人机在世界坐标系下的速度,分别表示k-1时刻无人机在世界坐标系下指向x,y,z三个轴方向的速度,表示无人机的姿态,以欧拉角的形式表示,分别为滚转角、俯仰角和偏航角,表示k-1时刻世界坐标系下的m个路标点的坐标,表示第i,i=1,2,...,m个路标点在世界坐标系下的坐标,为k-1时刻机体坐标系下的比力加速度,为k-1时刻机体坐标系下的角速度,为k-1时刻世界坐标下的重加速度,△t为采样时间间隔,
为k-1时刻机体坐标系到世界坐标系的变换矩阵,nk-1为系统噪声,噪声设定为高斯白噪声,Qk-1为系统噪声协方差矩阵。
3.根据权利要求1所述的一种基于无迹卡尔曼滤波的偏振光SLAM方法,其特征在于:所述步骤(2)建立激光雷达的量测模型;选取激光雷达测量的距离、角度作为测量值,针对激光雷达量测的数据给出路标点的量测模型:
其中
i=1,2,...,m表示第k时刻对第i个路标点的量测方程中的非线性函数,yk,i,lidar=[dk,i θk,i λk,i]T表示k时刻系统的量测值,dk,i表示第k时刻激光雷达测量的无人机质心与第i个路标点的距离,θk,i表示第k时刻激光雷达测量的无人机质心与第i个路标点的距离俯仰角,λk,i表示第k时刻激光雷达测量的无人机质心与第i个路标点的方位角,uk,lidar表示激光雷达量测噪声,噪声设定为高斯白噪声,Rk,lidar为量测噪声协方差矩阵。
4.根据权利要求1所述的一种基于无迹卡尔曼滤波的偏振光SLAM方法,其特征在于:所述步骤(3)建立偏振光传感器的量测模型:选取偏振光传感器量测的机体纵轴与太阳子午线的夹角量作为测值,针对偏振光给出的量测数据的观测模型如下:
其中为偏振光传感器获取的机体长轴与太阳子午线的夹角,为太阳方位角,为太阳高度角,为太阳赤纬,为观测点的纬度,为太阳时角,uk,polar表示偏振光传感器量测噪声,噪声设定为高斯白噪声,Rk,polar为量测噪声协方差矩阵。
5.根据权利要求1所述的一种基于无迹卡尔曼滤波的偏振光SLAM方法,其特征在于:所述步骤(4)系统初始化;将无人机坐标系转换到世界坐标系并将无人机,给定初始时刻的系统状态为初始协方差矩阵P0,并建立以无人机起始位置为原点的全局地图。
6.根据权利要求1所述的一种基于无迹卡尔曼滤波的偏振光SLAM方法,其特征在于:所述步骤(5)路标点匹配;对k时刻激光雷达量测到的数据,和在全局地图中存储的路标点利用迭代最近点(ICP)算法进行匹配,对于匹配成功的路标点,则将该路标点送入相应的滤波器进行后续的运算。
7.根据权利要求1所述的一种基于无迹卡尔曼滤波的偏振光SLAM方法,其特征在于:所述步骤(6)利用路标点的激光雷达数据和偏振光传感器数据,设计UKF滤波器,估计无人机的位置、速度、姿态和路标点的位置;首先将k-1时刻的估计状态和协方差矩阵,结合k-1时刻的IMU的输出数据,进行一步递推得到由IMU和状态转移矩阵推算出的机器人位置,然后根据预测值,加入激光传感器信息和偏振光信息进行量测更新的计算,输出无人机的位置、速度、姿态和路标点的位置,并同时输出相应的协方差矩阵,具体步骤如下:
时间更新:
其中为状态向量维度,Pk-1|k-1为k-1时刻的误协方差矩阵,为尺度参数,计算方法为该式中ε为一小量取值范围为[10-4,1],常量可取0或者3-n;
量测更新:
其中Pyy,k|k-1为自协方差矩阵,Pxy,k|k-1为互协方差矩阵。
8.根据权利要求1所述的一种基于无迹卡尔曼滤波的偏振光SLAM方法,其特征在于:所述步骤(7)地图、状态更新;对于全局地图中已存在的路标点,将对应的滤波后的路标点作为最新路标点,对步骤(5)中未匹配成功的路标点,直接将其加入到全局地图中。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810128645.1A CN108362288B (zh) | 2018-02-08 | 2018-02-08 | 一种基于无迹卡尔曼滤波的偏振光slam方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810128645.1A CN108362288B (zh) | 2018-02-08 | 2018-02-08 | 一种基于无迹卡尔曼滤波的偏振光slam方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108362288A true CN108362288A (zh) | 2018-08-03 |
CN108362288B CN108362288B (zh) | 2021-05-07 |
Family
ID=63004949
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810128645.1A Active CN108362288B (zh) | 2018-02-08 | 2018-02-08 | 一种基于无迹卡尔曼滤波的偏振光slam方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108362288B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109341705A (zh) * | 2018-10-16 | 2019-02-15 | 北京工业大学 | 智能探测机器人同时定位与地图构建系统 |
CN109613549A (zh) * | 2018-12-28 | 2019-04-12 | 芜湖哈特机器人产业技术研究院有限公司 | 一种基于卡尔曼过滤的激光雷达定位方法 |
CN109840517A (zh) * | 2019-03-08 | 2019-06-04 | 兰州交通大学 | 一种mems陀螺噪声估计和滤波方法 |
CN111136660A (zh) * | 2020-02-19 | 2020-05-12 | 清华大学深圳国际研究生院 | 机器人位姿定位方法及系统 |
CN112197741A (zh) * | 2020-12-04 | 2021-01-08 | 华南理工大学 | 基于扩展卡尔曼滤波的无人机slam技术测量倾斜角系统 |
CN112581616A (zh) * | 2020-12-18 | 2021-03-30 | 杭州电子科技大学 | 基于序贯分块滤波的最近邻ukf-slam方法 |
CN112697138A (zh) * | 2020-12-07 | 2021-04-23 | 北方工业大学 | 一种基于因子图优化的仿生偏振同步定位与构图的方法 |
CN113739795A (zh) * | 2021-06-03 | 2021-12-03 | 东北电力大学 | 一种基于偏振光/惯性/视觉组合导航的水下同步定位与建图方法 |
CN115574816A (zh) * | 2022-11-24 | 2023-01-06 | 东南大学 | 仿生视觉多源信息智能感知无人平台 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101201403A (zh) * | 2007-04-27 | 2008-06-18 | 北京航空航天大学 | 三维偏振成像激光雷达遥感器 |
EP2405281A1 (de) * | 2010-07-09 | 2012-01-11 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung | Verfahren und Vorrichtung zur Bestimmung der Position und Orientierung eines mobilen Senders |
CN103808316A (zh) * | 2012-11-12 | 2014-05-21 | 哈尔滨恒誉名翔科技有限公司 | 室内飞行智能体惯性系统与激光测距仪组合导航改进方法 |
CN106886030A (zh) * | 2017-03-24 | 2017-06-23 | 黑龙江硅智机器人有限公司 | 应用于服务机器人的同步式地图构建与定位系统及方法 |
CN107300697A (zh) * | 2017-06-07 | 2017-10-27 | 南京航空航天大学 | 基于无人机的运动目标ukf滤波方法 |
-
2018
- 2018-02-08 CN CN201810128645.1A patent/CN108362288B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101201403A (zh) * | 2007-04-27 | 2008-06-18 | 北京航空航天大学 | 三维偏振成像激光雷达遥感器 |
EP2405281A1 (de) * | 2010-07-09 | 2012-01-11 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung | Verfahren und Vorrichtung zur Bestimmung der Position und Orientierung eines mobilen Senders |
CN103808316A (zh) * | 2012-11-12 | 2014-05-21 | 哈尔滨恒誉名翔科技有限公司 | 室内飞行智能体惯性系统与激光测距仪组合导航改进方法 |
CN106886030A (zh) * | 2017-03-24 | 2017-06-23 | 黑龙江硅智机器人有限公司 | 应用于服务机器人的同步式地图构建与定位系统及方法 |
CN107300697A (zh) * | 2017-06-07 | 2017-10-27 | 南京航空航天大学 | 基于无人机的运动目标ukf滤波方法 |
Non-Patent Citations (2)
Title |
---|
LIZHONG SONGHAIMING QIXIANDE MENG: "Scheme of adaptive polarization filtering based on Kalman model", 《JOURNAL OF SYSTEMS ENGINEERING AND ELECTRONICS》 * |
王道斌: "基于天空偏振光的SLAM方法的研究", 《中国博士学位论文全文数据库 信息科技辑》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109341705A (zh) * | 2018-10-16 | 2019-02-15 | 北京工业大学 | 智能探测机器人同时定位与地图构建系统 |
CN109613549A (zh) * | 2018-12-28 | 2019-04-12 | 芜湖哈特机器人产业技术研究院有限公司 | 一种基于卡尔曼过滤的激光雷达定位方法 |
CN109613549B (zh) * | 2018-12-28 | 2023-04-07 | 芜湖哈特机器人产业技术研究院有限公司 | 一种基于卡尔曼过滤的激光雷达定位方法 |
CN109840517A (zh) * | 2019-03-08 | 2019-06-04 | 兰州交通大学 | 一种mems陀螺噪声估计和滤波方法 |
CN111136660B (zh) * | 2020-02-19 | 2021-08-03 | 清华大学深圳国际研究生院 | 机器人位姿定位方法及系统 |
CN111136660A (zh) * | 2020-02-19 | 2020-05-12 | 清华大学深圳国际研究生院 | 机器人位姿定位方法及系统 |
CN112197741A (zh) * | 2020-12-04 | 2021-01-08 | 华南理工大学 | 基于扩展卡尔曼滤波的无人机slam技术测量倾斜角系统 |
CN112697138A (zh) * | 2020-12-07 | 2021-04-23 | 北方工业大学 | 一种基于因子图优化的仿生偏振同步定位与构图的方法 |
CN112581616A (zh) * | 2020-12-18 | 2021-03-30 | 杭州电子科技大学 | 基于序贯分块滤波的最近邻ukf-slam方法 |
CN112581616B (zh) * | 2020-12-18 | 2024-06-07 | 杭州电子科技大学 | 基于序贯分块滤波的最近邻ukf-slam方法 |
CN113739795A (zh) * | 2021-06-03 | 2021-12-03 | 东北电力大学 | 一种基于偏振光/惯性/视觉组合导航的水下同步定位与建图方法 |
CN113739795B (zh) * | 2021-06-03 | 2023-10-20 | 东北电力大学 | 一种基于偏振光/惯性/视觉组合导航的水下同步定位与建图方法 |
CN115574816A (zh) * | 2022-11-24 | 2023-01-06 | 东南大学 | 仿生视觉多源信息智能感知无人平台 |
Also Published As
Publication number | Publication date |
---|---|
CN108362288B (zh) | 2021-05-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108362288A (zh) | 一种基于无迹卡尔曼滤波的偏振光slam方法 | |
You et al. | Data fusion of UWB and IMU based on unscented Kalman filter for indoor localization of quadrotor UAV | |
CN110243358B (zh) | 多源融合的无人车室内外定位方法及系统 | |
CN108387236B (zh) | 一种基于扩展卡尔曼滤波的偏振光slam方法 | |
Fan et al. | Data fusion for indoor mobile robot positioning based on tightly coupled INS/UWB | |
CN108871336B (zh) | 一种车辆位置估算系统及方法 | |
CN103697889B (zh) | 一种基于多模型分布式滤波的无人机自主导航与定位方法 | |
CN111426320B (zh) | 一种基于图像匹配/惯导/里程计的车辆自主导航方法 | |
CN112697138B (zh) | 一种基于因子图优化的仿生偏振同步定位与构图的方法 | |
Zheng et al. | An optimization-based UWB-IMU fusion framework for UGV | |
CN101762273A (zh) | 一种深空探测器软着陆自主光学导航方法 | |
CN107727101B (zh) | 基于双偏振光矢量的三维姿态信息快速解算方法 | |
CN110057356B (zh) | 一种隧道内车辆定位方法及装置 | |
CN113252038B (zh) | 基于粒子群算法的航迹规划地形辅助导航方法 | |
CN102788580A (zh) | 无人机视觉导航中的飞行路线合成方法 | |
CN108195376A (zh) | 小型无人机自主导航定位方法 | |
CN105333869A (zh) | 一种基于自适应ekf的无人侦察机同步定位与构图方法 | |
CN106885573A (zh) | 面向四旋翼飞行器的运动捕捉系统实时测姿方法 | |
CN106441372A (zh) | 一种基于偏振与重力信息的静基座粗对准方法 | |
CN104316058B (zh) | 一种采用交互多模型的移动机器人wsn/ins组合导航方法 | |
CN111238469A (zh) | 一种基于惯性/数据链的无人机编队相对导航方法 | |
CN110887472A (zh) | 一种偏振-地磁信息深度融合全自主姿态解算方法 | |
CN113129377A (zh) | 一种三维激光雷达快速鲁棒slam方法和装置 | |
Huang et al. | Design and performance analysis of landmark-based INS/Vision Navigation System for UAV | |
CN114994732B (zh) | 基于gnss载波相位的车载航向快速初始化装置及方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |