CN111538029A - 一种视觉与雷达融合的测量方法及终端 - Google Patents

一种视觉与雷达融合的测量方法及终端 Download PDF

Info

Publication number
CN111538029A
CN111538029A CN202010329732.0A CN202010329732A CN111538029A CN 111538029 A CN111538029 A CN 111538029A CN 202010329732 A CN202010329732 A CN 202010329732A CN 111538029 A CN111538029 A CN 111538029A
Authority
CN
China
Prior art keywords
carrier
feature points
radar
binocular camera
position error
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.)
Pending
Application number
CN202010329732.0A
Other languages
English (en)
Inventor
余平
项崴
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jiangsu Shenghai Intelligent Technology Co ltd
Original Assignee
Jiangsu Shenghai Intelligent Technology Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Jiangsu Shenghai Intelligent Technology Co ltd filed Critical Jiangsu Shenghai Intelligent Technology Co ltd
Priority to CN202010329732.0A priority Critical patent/CN111538029A/zh
Publication of CN111538029A publication Critical patent/CN111538029A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S17/00Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
    • G01S17/86Combinations of lidar systems with systems other than lidar, radar or sonar, e.g. with direction finders
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S17/00Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
    • G01S17/88Lidar systems specially adapted for specific applications
    • G01S17/93Lidar systems specially adapted for specific applications for anti-collision purposes
    • G01S17/931Lidar systems specially adapted for specific applications for anti-collision purposes of land vehicles
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/02Control of position or course in two dimensions
    • G05D1/021Control of position or course in two dimensions specially adapted to land vehicles
    • G05D1/0231Control of position or course in two dimensions specially adapted to land vehicles using optical position detecting means
    • G05D1/0246Control of position or course in two dimensions specially adapted to land vehicles using optical position detecting means using a video camera in combination with image processing means
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/02Control of position or course in two dimensions
    • G05D1/021Control of position or course in two dimensions specially adapted to land vehicles
    • G05D1/0257Control of position or course in two dimensions specially adapted to land vehicles using a radar
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/70Determining position or orientation of objects or cameras
    • G06T7/73Determining position or orientation of objects or cameras using feature-based methods

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Automation & Control Theory (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Multimedia (AREA)
  • Theoretical Computer Science (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种视觉与雷达融合的测量方法及终端,将双目摄像机与毫米波雷达放置于载体上;提取双目摄像机所拍摄图像上的特征点,对所述特征点进行追踪并匹配,确定所述特征点的第一位置误差;通过所述毫米波雷达跟踪载体运动时周边的静态物体,确定所述跟踪的静态物体的第二位置误差;联合所述第一位置误差和第二位置误差,估计所述载体的位姿;通过毫米波雷达和双目摄像机共同测量目标的深度信息,提升了测量系统的鲁棒性,计算双目摄像机所拍摄图像上特征点的匹配误差,及毫米波雷达跟踪载体周边静态物体时的匹配误差,并将其进行联合优化,做出载体的最优位姿估计,确保了载体位姿估计结果的准确性,实现在降低成本的同时提高定位精度。

Description

一种视觉与雷达融合的测量方法及终端
技术领域
本发明涉及无人驾驶测量领域,尤其涉及一种视觉与雷达融合的测量方法及终端。
背景技术
目前,在无人驾驶领域,主要通过应用激光雷达融合视觉进行局部定位,其中视觉摄像机主要是单目,现有的方法下单目能保证计算的效率,但是单目自身无法通过时差直接估计测量物的尺度信息,间接进行估计得到的数据又无法保证精确度,直接使用会影响无人驾驶的安全性,一般采用激光雷达对单目视觉进行辅助,但激光雷达的价格昂贵,需要在车辆外安装、影响车辆外观美观,并且容易受外界环境制约,如雨天、雾天和风沙大的天气,都会较大影响激光雷达的测量精度从而影响到驾驶安全。
发明内容
本发明所要解决的技术问题是:提供一种视觉与雷达融合的测量方法及终端,在降低成本的同时提高定位的精度。
为了解决上述技术问题,本发明采用的一种技术方案为:
一种视觉与雷达融合的测量方法,包括步骤:
S1、将双目摄像机与毫米波雷达放置于载体上;
S2、提取双目摄像机所拍摄图像上的特征点,对所述特征点进行追踪并匹配,确定所述特征点的第一位置误差;
S3、通过所述毫米波雷达跟踪载体运动时周边的静态物体,确定所述跟踪的静态物体的第二位置误差;
S4、联合所述第一位置误差和第二位置误差,估计所述载体的位姿。
为了解决上述技术问题,本发明采用的另一种技术方案为:
一种视觉与雷达融合的测量终端,包括存储器、处理器及存储在存储器上并可在所述处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现以下步骤:
S1、将双目摄像机与毫米波雷达放置于载体上;
S2、提取双目摄像机所拍摄图像上的特征点,对所述特征点进行追踪并匹配,确定所述特征点的第一位置误差;
S3、通过所述毫米波雷达跟踪载体运动时周边的静态物体,确定所述跟踪的静态物体的第二位置误差;
S4、联合所述第一位置误差和第二位置误差,估计所述载体的位姿。
本发明的有益效果在于:利用双目视觉和毫米波雷达共同提供信息对载体的位姿状态进行估计,测量数据能够互相验证,保证了测量的精度,并且在一方出现问题时,另一方可以单独完成基本数据的测量,同时毫米波雷达的造价相对较低,能够有效节约成本,并且毫米波雷达适用范围较广,可以在一些特殊环境下正常工作,保证了自动驾驶的安全性,提升了测量方法的鲁棒性;将双目摄像机测量得到的第一位置误差和雷达跟踪得到的第二位置误差进行联合后估计载体的位姿,实现在降低成本的同时提高定位精度。
附图说明
图1为本发明实施例的一种视觉与雷达融合的测量方法的步骤流程图;
图2为本发明实施例的一种视觉与雷达融合的测量终端的结构示意图;
图3为本发明实施例的一种视觉与雷达融合的测量方法的总流程图;
标号说明:
1、一种视觉与雷达融合的测量终端;2、处理器;3、存储器。
具体实施方式
为详细说明本发明的技术内容、所实现目的及效果,以下结合实施方式并配合附图予以说明。
请参照图1,一种视觉与雷达融合的测量方法,包括步骤:
S1、将双目摄像机与毫米波雷达放置于载体上;
S2、提取双目摄像机所拍摄图像上的特征点,对所述特征点进行追踪并匹配,确定所述特征点的第一位置误差;
S3、通过所述毫米波雷达跟踪载体运动时周边的静态物体,确定所述跟踪的静态物体的第二位置误差;
S4、联合所述第一位置误差和第二位置误差,估计所述载体的位姿。
本发明的有益效果在于:利用双目视觉和毫米波雷达共同提供信息对载体的位姿状态进行估计,测量数据能够互相验证,保证了测量的精度,并且在一方出现问题时,另一方可以单独完成基本数据的测量,同时毫米波雷达的造价相对较低,能够有效节约成本,并且毫米波雷达适用范围较广,可以在一些特殊环境下正常工作,保证了自动驾驶的安全性,提升了测量方法的鲁棒性;将双目摄像机测量得到的第一位置误差和雷达跟踪得到的第二位置误差进行联合后估计载体的位姿,实现在降低成本的同时提高定位精度。
进一步的,所述S2中所述提取双目摄像机所拍摄图像上的特征点包括:
提取双目摄像机拍摄的图像上的角点作为图像特征点;
通过KLT光流法进行不同时刻的摄像机拍摄的图像上同个角点的跟踪,对不同时刻图像上的特征点进行匹配;
计算匹配上的特征点在不同时刻相对于摄像机的深度信息;
根据所述匹配上的特征点的位置信息和深度信息确定其第一位置误差。
由上述描述可知,以角点作为图像特征点进行跟踪,角点的特征较明显,方便实现匹配,并且保证了匹配的精度,且匹配角点位置之后,能够较方便得到其他点的匹配的结果,利用KLT光流法进行不同时刻特征点的匹配,进一步加大了匹配的准确度。
进一步的,所述S2中确定所述特征点的第一位置误差具体为:
Figure BDA0002464506460000031
其中,ucj和vcj是j时刻双目摄像机测量到的归一化深度的特征点相对于双目摄像机的位置信息,xcj,ycj和zcj是由j时刻摄像头测量值计算得出,与i时刻双目摄像机测量到的归一化深度的特征点相对于双目摄像机的位置信息的关系为:
Figure BDA0002464506460000041
由上式估计出Rbjn
其中,λ为逆深度,其值为
Figure BDA0002464506460000042
z为j时刻双目摄像机测量出的所述特征点对应的深度信息,Rcb为摄像机相对于载体坐标系的位姿矩阵,Rbin和Rbjn分别为i时刻和j时刻所述载体相对于导航坐标系的位姿矩阵。
由上述描述可知,建立本时刻双目摄像机测量到的归一化深度的特征点与其估计值的差值,并将所述差值与载体的位姿矩阵相关联,实现对载体位姿矩阵的最优估计。
进一步的,所述S3中确定所述跟踪的静态物体的第二位置误差具体为:
Figure BDA0002464506460000043
上式中,xmrj和ymrj为毫米波雷达对于所述周边静态物体的水平面扫描测量得到的位置信息,
Figure BDA0002464506460000044
Figure BDA0002464506460000045
是由j时刻毫米波雷达测量值计算得出,为所述周边静态物体的水平面位置信息,其与i时刻周边静态物体的水平面位置信息的关系如下:
Figure BDA0002464506460000046
由上式估计出Rbjn,hori
上式中,Rrb,hori表示毫米波雷达相对于载体坐标系的水平方向的位姿矩阵,Rbin,hori和Rbjn,hori分别为i时刻和j时刻所述载体相对于导航坐标系的水平方向位姿矩阵。
由上述描述可知,建立本时刻毫米波雷达对于周边静态物体的水平面扫描测量得到的位置信息及其估计值的差值,并将所述差值与载体位姿相关联,实现在差值最小的情况下估计出最优的载体位姿。
进一步的,所述S4具体为:
设定目标函数:
Figure BDA0002464506460000051
上式中resprior=ρ(||rsw,c||)+ρ(||rsw,r||);
rescam=ρ(residualc);
resmrad=ρ(residualr);
resprior为先验误差,rsw,c为滑动窗口包含的起始时刻前的相机定位误差,rsw,r为滑动窗口包含的起始时刻前的雷达定位误差;
rescam为第一位置误差通过柯西鲁棒核函数后的值,resmard第二位置误差通过柯西鲁棒核函数后的值,ρ表示柯西鲁棒核函数;
计算所述目标函数的雅可比矩阵,通过LM优化器,得到载体相对于导航坐标系的位姿矩阵
Figure BDA0002464506460000052
由上述描述可知,将第一位置误差、第二位置误差、先验误差都考虑到位姿估计中,保证最终的估计结果误差最小,计算雅可比矩阵并通过LM优化器对所述目标函数进行优化,确保最终估计的位姿矩阵与实际值的差距最小,提高位姿估计的准确度。
请参照图2,一种视觉与雷达融合的测量终端,包括存储器、处理器及存储在存储器上并可在所述处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现以下步骤:
S1、将双目摄像机与毫米波雷达放置于载体上;
S2、提取双目摄像机所拍摄图像上的特征点,对所述特征点进行追踪并匹配,确定所述特征点的第一位置误差;
S3、通过所述毫米波雷达跟踪载体运动时周边的静态物体,确定所述跟踪的静态物体的第二位置误差;
S4、联合所述第一位置误差和第二位置误差,估计所述载体的位姿。
本发明的有益效果在于:利用双目视觉和毫米波雷达共同提供信息对载体的位姿状态进行估计,保证了测量的精度,并且在一方出现问题时,另一方可以单独完成基本数据的测量,可以在一些特殊环境下正常工作,保证了自动驾驶的安全性,提升了测量方法的鲁棒性,将双目摄像机测量得到的第一信息和雷达跟踪得到的第二信息进行匹配,融合了两种传感器的深度测量数据,大大提高了定位精度;双目摄像头可以直接进行深度的估计,误差较小,不用次次都引入雷达的测量数据进行矫正,提高了效率,毫米波雷达的成本大大小于激光雷达,节约了成本。
进一步的,所述S2中所述提取双目摄像机所拍摄图像上的特征点包括:
提取双目摄像机拍摄的图像上的角点作为图像特征点;
通过KLT光流法进行不同时刻的摄像机拍摄的图像上同个角点的跟踪,对不同时刻图像上的特征点进行匹配;
计算匹配上的特征点在不同时刻相对于摄像机的深度信息;
根据所述匹配上的特征点的位置信息和深度信息确定其第一位置误差。
由上述描述可知,以角点作为图像特征点进行跟踪,角点的特征较明显,方便实现匹配,并且保证了匹配的精度,且匹配角点位置之后,能够较方便得到其他点的匹配的结果,利用KLT光流法进行不同时刻特征点的匹配,进一步加大了匹配的准确度。
进一步的,所述S2中确定所述特征点的第一位置误差具体为:
Figure BDA0002464506460000061
其中,ucj和vcj是j时刻双目摄像机测量到的归一化深度的特征点相对于双目摄像机的位置信息,xcj,ycj和zcj是由j时刻摄像头测量值计算得出,与i时刻双目摄像机测量到的归一化深度的特征点相对于双目摄像机的位置信息的关系为:
Figure BDA0002464506460000071
由上式估计出Rbjn
其中,λ为逆深度,其值为
Figure BDA0002464506460000072
z为j时刻双目摄像机测量出的所述特征点对应的深度信息,Rcb为摄像机相对于载体坐标系的位姿矩阵,Rbin和Rbjn分别为i时刻和j时刻所述载体相对于导航坐标系的位姿矩阵。
由上述描述可知,建立本时刻双目摄像机测量到的归一化深度的特征点与其估计值的差值,并将所述差值与载体的位姿矩阵相关联,实现对载体位姿矩阵的最优估计。
进一步的,所述S3中确定所述跟踪的静态物体的第二位置误差具体为:
Figure BDA0002464506460000073
上式中,xmrj和ymrj为毫米波雷达对于所述周边静态物体的水平面扫描测量得到的位置信息,
Figure BDA0002464506460000074
Figure BDA0002464506460000075
是由j时刻毫米波雷达测量值计算得出,为所述周边静态物体的水平面位置信息,其与i时刻周边静态物体的水平面位置信息的关系如下:
Figure BDA0002464506460000076
由上式估计出Rbjn,hori
上式中,Rrb,hori表示毫米波雷达相对于载体坐标系的水平方向的位姿矩阵,Rbin,hori和Rbjn,hori分别为i时刻和j时刻所述载体相对于导航坐标系的水平方向位姿矩阵。
由上述描述可知,建立本时刻毫米波雷达对于周边静态物体的水平面扫描测量得到的位置信息及其估计值的差值,并将所述差值与载体位姿相关联,实现在差值最小的情况下估计出最优的载体位姿。
进一步的,所述S4具体为:
设定目标函数:
Figure BDA0002464506460000081
上式中resprior=ρ(||rsw,c||)+ρ(||rsw,r||);
rescam=ρ(residualc);
resmrad=ρ(residualr);
resprior为先验误差,rsw,c为滑动窗口包含的起始时刻前的相机定位误差,rsw,r为滑动窗口包含的起始时刻前的雷达定位误差;
rescam为第一位置误差通过柯西鲁棒核函数后的值,resmard第二位置误差通过柯西鲁棒核函数后的值,ρ表示柯西鲁棒核函数;
计算所述目标函数的雅可比矩阵,通过LM优化器,得到载体相对于导航坐标系的位姿矩阵
Figure BDA0002464506460000082
由上述描述可知,将第一位置误差、第二位置误差、先验误差都考虑到位姿估计中,保证最终的估计结果误差最小,计算雅可比矩阵并通过LM优化器对所述目标函数进行优化,确保最终估计的位姿矩阵与实际值的差距最小
请参照图1,本发明的实施例一为:
下文所述导航坐标系:以载体的初始坐标为原点,过所述原点且指向所述载体在所述原点时的东方向的射线为xn轴,北方向的射线为yn轴,过所述原点且垂直于xnyn平面并指向所述载体顶部的射线为zn轴;
载体坐标系:以载体的重心为原点,过所述原点且指向所述载体左侧的射线为xb轴,过所述原点且指向所述载体前方的射线为yb轴,过所述原点且垂直于xbyb平面并指向所述载体顶部的射线为zb轴;
雷达坐标系:以雷达测量零点为原点,过所述原点且指向所述雷达左侧的射线为xr轴,过所述原点且指向所述雷达前方的射线为yr轴,过所述原点且垂直于xryr平面并指向所述雷达上方的射线为zr轴;
相机坐标系:以相机重心为原点,过所述原点且指向所述相机右侧的射线为xc轴,过所述原点且指向所述相机下方的射线为yc轴,过所述原点且垂直于xcyc平面并指向所述相机前方的射线为zc轴;
一种视觉与雷达融合的测量方法,具体包括步骤:
S1、将双目摄像机与毫米波雷达放置于载体上;
进行毫米波雷达与双目摄像头外参数的标定,得到摄像机坐标系与雷达坐标系之间的转化关系、摄像机坐标系与载体坐标系之间的转化关系及雷达坐标系与载体坐标系之间的转化关系;
S2、提取双目摄像机所拍摄图像上的特征点,对所述特征点进行追踪并匹配,确定所述特征点的第一位置误差;
所述提取双目摄像机所拍摄图像上的特征点包括:
提取双目摄像机拍摄的图像上的角点作为图像特征点;
通过KLT光流法进行不同时刻的摄像机拍摄的图像上同个角点的跟踪,对不同时刻图像上的特征点进行匹配;
计算匹配上的特征点在不同时刻相对于摄像机的深度信息;
根据所述匹配上的特征点的位置信息和深度信息确定其第一位置误差;
S3、通过所述毫米波雷达跟踪载体运动时周边的静态物体,确定所述跟踪的静态物体的第二位置误差;
S4、联合所述第一位置误差和第二位置误差,估计所述载体的位姿。
本发明的实施例二为:
一种视觉与雷达融合的测量方法,其与实施例一的不同之处在于,所述S2中,确定所述特征点的第一位置误差具体为:
Figure BDA0002464506460000091
其中,ucj和vcj是j时刻双目摄像机测量到的归一化深度的特征点相对于双目摄像机的位置信息,xcj,ycj和zcj是由j时刻摄像头测量值计算得出,与i时刻双目摄像机测量到的归一化深度的特征点相对于双目摄像机的位置信息的关系为:
Figure BDA0002464506460000101
由上式估计出Rbjn
其中,λ为逆深度,其值为
Figure BDA0002464506460000102
z为i时刻双目摄像机测量出的所述特征点对应的深度信息,Rcb为摄像机相对于载体坐标系的位姿矩阵,可由标定后的摄像机坐标系与载体坐标系之间的转换关系确定,Rbin和Rbjn分别为i时刻和j时刻所述载体相对于导航坐标系的位姿矩阵,Rbin为已估计出的i时刻载体相对于导航坐标系的位姿矩阵,若i为初始时刻,则Rbin已知;
所述S3中,确定所述跟踪的静态物体的第二位置误差具体为:
Figure BDA0002464506460000103
上式中,xmrj和ymrj为毫米波雷达对于所述周边静态物体的水平面扫描测量得到的位置信息,
Figure BDA0002464506460000104
Figure BDA0002464506460000105
是由j时刻毫米波雷达测量值计算得出,为所述周边静态物体的水平面位置信息,其与i时刻周边静态物体的水平面位置信息的关系如下:
Figure BDA0002464506460000106
由上式估计出Rbjn,hori
上式中,Rrb,hori表示毫米波雷达相对于载体坐标系的水平方向的位姿矩阵,可以由标定后的雷达坐标系与载体坐标系之间的转换关系确定,Rbin,hori和Rbjn,hori分别为i时刻和j时刻所述载体相对于导航坐标系的水平方向位姿矩阵,Rbin,hori为已估计出的i时刻所述载体相对于导航坐标系的水平方向位姿矩阵,若i为初始时刻,则Rbin,hori已知;
所述S4具体为:
视觉估计的是载体的三维位姿Rbn,毫米波雷达估计的是载体的二维位姿Rbn,hori,Rbn,hori被Rbn所包含,因此,所需估计的载体位姿为:
Figure BDA0002464506460000111
上式中,X表示滑动窗口,N为滑动窗口内待估计系数矩阵起始的时刻,M为滑动窗口内待估计系数矩阵由起始到最后时刻的时间间隔,也可以表示所述滑动窗口的数目;
设定目标函数:
Figure BDA0002464506460000112
上式中resprior=ρ(||rsw,c||)+ρ(||rsw,r||);
rescam=ρ(residualc);
resmrad=ρ(residualr);
resprior为先验误差,rsw,c为滑动窗口包含的起始时刻前的相机定位误差,rsw,r为滑动窗口包含的起始时刻前的雷达定位误差;
rescam为第一位置误差通过柯西鲁棒核函数后的值,resmard第二位置误差通过柯西鲁棒核函数后的值,ρ表示柯西鲁棒核函数,参照下式:
Figure BDA0002464506460000113
其中,c为系数,一般设定为2.3849,s为输入值;
使用柯西鲁棒核函数,在误差过大时(可能特征提取错误,匹配错误,或者测量值误差过大等造成的误差过大),可以减小其误差值,并且不会对误差值小的数据有影响,这样可以通过这个柯西鲁棒核函数减低误差过大的值对于优化过程中的正确收敛点产生影响,从而影响定位精度;
计算所述目标函数的雅可比矩阵(为代价函数对于带估参数的求导),通过LM优化器,得到载体相对于导航坐标系的位姿矩阵
Figure BDA0002464506460000121
LM优化器:
(JTJ+μI)ΔXlm=-JTf withμ≥0
上式中,J为其代价函数
Figure BDA0002464506460000122
的雅可比矩阵,f为resprior+rescam+resmrad;的值,μ为阻尼系数,设定初值后,一直迭代计算,μ会根据每次迭代后的ΔXlm的变化量进行调整,当ΔXlm变化值小于额定阈值,并且其值小于特定阈值后,收敛,得到估计出f中带估参数
Figure BDA0002464506460000123
的值;
然后滑动窗口进行滑动,剔除头一个点,把新的点的观测值加入窗口内,然后把剔除点的误差加入到先验误差内,再次计算加入新点的滑动窗口内代价函数的误差值和雅可比矩阵的值,然后再次优化;
重复上述步骤直至结束对载体位置的确定。
请参照图2,本发明的实施例三为:
一种视觉与雷达融合的测量终端1,所述终端1包括处理器2、存储器3及存储在存储器3上并可在所述处理器2上运行的计算机程序,所述处理器2执行所述计算机程序时实现实施例一或实施例二中的各个步骤。
综上所述,本发明提供了一种视觉与雷达融合的测量方法及终端,通过毫米波雷达和双目摄像机共同对目标的深度信息进行测量,提升了系统的鲁棒性,通过计算双目摄像机所拍摄图像上的特征点的匹配误差,及毫米波雷达跟踪载体周边静态物体时的匹配误差,并将所述两项误差与载体的位姿相关联,从而实现在进行误差最小的优化时,做出载体的最优位姿估计,同时将雷达扫描获得的水平面位姿与摄像机获得的三维位姿进行联合,确保了载体位姿估计结果的准确性,同时在最后采用LM方法进一步进行优化,保证最后获得的载体位姿与实际值的误差最小。
以上所述仅为本发明的实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书及附图内容所作的等同变换,或直接或间接运用在相关的技术领域,均同理包括在本发明的专利保护范围内。

Claims (10)

1.一种视觉与雷达融合的测量方法,其特征在于,包括步骤:
S1、将双目摄像机与毫米波雷达放置于载体上;
S2、提取双目摄像机所拍摄图像上的特征点,对所述特征点进行追踪并匹配,确定所述特征点的第一位置误差;
S3、通过所述毫米波雷达跟踪载体运动时周边的静态物体,确定所述跟踪的静态物体的第二位置误差;
S4、联合所述第一位置误差和第二位置误差,估计所述载体的位姿。
2.根据权利要求1所述的一种视觉与雷达融合的测量方法,其特征在于,所述S2中所述提取双目摄像机所拍摄图像上的特征点包括:
提取双目摄像机拍摄的图像上的角点作为图像特征点;
通过KLT光流法进行不同时刻的摄像机拍摄的图像上同个角点的跟踪,对不同时刻图像上的特征点进行匹配;
计算匹配上的特征点在不同时刻相对于摄像机的深度信息;
根据所述匹配上的特征点的位置信息和深度信息确定其第一位置误差。
3.根据权利要求2所述的一种视觉与雷达融合的测量方法,其特征在于,所述S2中确定所述特征点的第一位置误差具体为:
Figure FDA0002464506450000011
其中,ucj和vcj是j时刻双目摄像机测量到的归一化深度的特征点相对于双目摄像机的位置信息,xcj,ycj和zcj是由j时刻摄像头测量值计算得出,与i时刻双目摄像机测量到的归一化深度的特征点相对于双目摄像机的位置信息的关系为:
Figure FDA0002464506450000012
由上式估计出Rbjn
其中,λ为逆深度,其值为
Figure FDA0002464506450000021
z为j时刻双目摄像机测量出的所述特征点对应的深度信息,Rcb为摄像机相对于载体坐标系的位姿矩阵,Rbin和Rbjn分别为i时刻和j时刻所述载体相对于导航坐标系的位姿矩阵。
4.根据权利要求3所述的一种视觉与雷达融合的测量方法,其特征在于,所述S3中确定所述跟踪的静态物体的第二位置误差具体为:
Figure FDA0002464506450000022
上式中,xmrj和ymrj为毫米波雷达对于所述周边静态物体的水平面扫描测量得到的位置信息,
Figure FDA0002464506450000023
Figure FDA0002464506450000024
是由j时刻毫米波雷达测量值计算得出,为所述周边静态物体的水平面位置信息,其与i时刻周边静态物体的水平面位置信息的关系如下:
Figure FDA0002464506450000025
由上式估计出Rbjn,hori
上式中,Rrb,hori表示毫米波雷达相对于载体坐标系的水平方向的位姿矩阵,Rbin,hori和Rbjn,hori分别为i时刻和j时刻所述载体相对于导航坐标系的水平方向位姿矩阵。
5.根据权利要求4所述的一种视觉与雷达融合的测量方法,其特征在于,所述S4具体为:
设定目标函数:
Figure FDA0002464506450000026
上式中resprior=ρ(||rsw,c||)+ρ(||rsw,r||);
rescam=ρ(residualc);
resmrad=ρ(residualr);
resprior为先验误差,rsw,c为滑动窗口包含的起始时刻前的相机定位误差,rsw,r为滑动窗口包含的起始时刻前的雷达定位误差;
rescam为第一位置误差通过柯西鲁棒核函数后的值,resmard第二位置误差通过柯西鲁棒核函数后的值,ρ表示柯西鲁棒核函数;
计算所述目标函数的雅可比矩阵,通过LM优化器,得到载体相对于导航坐标系的位姿矩阵
Figure FDA0002464506450000031
6.一种视觉与雷达融合的测量终端,包括存储器、处理器及存储在存储器上并可在所述处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现以下步骤:
S1、将双目摄像机与毫米波雷达放置于载体上;
S2、提取双目摄像机所拍摄图像上的特征点,对所述特征点进行追踪并匹配,确定所述特征点的第一位置误差;
S3、通过所述毫米波雷达跟踪载体运动时周边的静态物体,确定所述跟踪的静态物体的第二位置误差。
7.根据权利要求6所述的一种视觉与雷达融合的测量终端,其特征在于,所述S2中所述提取双目摄像机所拍摄图像上的特征点包括:
提取双目摄像机拍摄的图像上的角点作为图像特征点;
通过KLT光流法进行不同时刻的摄像机拍摄的图像上同个角点的跟踪,对不同时刻图像上的特征点进行匹配;
计算匹配上的特征点在不同时刻相对于摄像机的深度信息;
根据所述匹配上的特征点的位置信息和深度信息确定其第一位置误差。
8.根据权利要求7所述的一种视觉与雷达融合的测量终端,其特征在于,所述S2中确定所述特征点的第一位置误差具体为:
Figure FDA0002464506450000041
其中,ucj和vcj是j时刻双目摄像机测量到的归一化深度的特征点相对于双目摄像机的位置信息,xcj,ycj和zcj是由j时刻摄像头测量值计算得出,与i时刻双目摄像机测量到的归一化深度的特征点相对于双目摄像机的位置信息的关系为:
Figure FDA0002464506450000042
由上式估计出Rbjn
其中,λ为逆深度,其值为
Figure FDA0002464506450000043
z为j时刻双目摄像机测量出的所述特征点对应的深度信息,Rcb为摄像机相对于载体坐标系的位姿矩阵,Rbin和Rbjn分别为i时刻和j时刻所述载体相对于导航坐标系的位姿矩阵。
9.根据权利要求8所述的一种视觉与雷达融合的测量终端,其特征在于,所述S3中确定所述跟踪的静态物体的第二位置误差具体为:
Figure FDA0002464506450000044
上式中,xmrj和ymrj为毫米波雷达对于所述周边静态物体的水平面扫描测量得到的位置信息,
Figure FDA0002464506450000045
Figure FDA0002464506450000046
是由j时刻毫米波雷达测量值计算得出,为所述周边静态物体的水平面位置信息,其与i时刻周边静态物体的水平面位置信息的关系如下:
Figure FDA0002464506450000047
由上式估计出Rbjn,hori
上式中,Rrb,hori表示毫米波雷达相对于载体坐标系的水平方向的位姿矩阵,Rbin,hori和Rbjn,hori分别为i时刻和j时刻所述载体相对于导航坐标系的水平方向位姿矩阵。
10.根据权利要求9所述的一种视觉与雷达融合的测量终端,其特征在于,所述S4具体为:
设定目标函数:
Figure FDA0002464506450000051
上式中resprior=ρ(||rsw,c||)+ρ(||rsw,r||);
rescam=ρ(residualc);
resmrad=ρ(residualr);
resprior为先验误差,rsw,c为滑动窗口包含的起始时刻前的相机定位误差,rsw,r为滑动窗口包含的起始时刻前的雷达定位误差;
rescam为第一位置误差通过柯西鲁棒核函数后的值,resmard第二位置误差通过柯西鲁棒核函数后的值,ρ表示柯西鲁棒核函数;
计算所述目标函数的雅可比矩阵,通过LM优化器,得到载体相对于导航坐标系的位姿矩阵
Figure FDA0002464506450000052
CN202010329732.0A 2020-04-24 2020-04-24 一种视觉与雷达融合的测量方法及终端 Pending CN111538029A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010329732.0A CN111538029A (zh) 2020-04-24 2020-04-24 一种视觉与雷达融合的测量方法及终端

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010329732.0A CN111538029A (zh) 2020-04-24 2020-04-24 一种视觉与雷达融合的测量方法及终端

Publications (1)

Publication Number Publication Date
CN111538029A true CN111538029A (zh) 2020-08-14

Family

ID=71977222

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010329732.0A Pending CN111538029A (zh) 2020-04-24 2020-04-24 一种视觉与雷达融合的测量方法及终端

Country Status (1)

Country Link
CN (1) CN111538029A (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112230211A (zh) * 2020-10-15 2021-01-15 长城汽车股份有限公司 车辆定位方法、装置、存储介质及车辆
CN113030944A (zh) * 2021-04-16 2021-06-25 深圳市众云信息科技有限公司 一种雷达目标追踪方法
CN113625271A (zh) * 2021-07-29 2021-11-09 中汽创智科技有限公司 基于毫米波雷达和双目相机的同时定位与建图方法
WO2023142814A1 (zh) * 2022-01-30 2023-08-03 中国第一汽车股份有限公司 目标识别方法、装置、设备及存储介质

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107869989A (zh) * 2017-11-06 2018-04-03 东北大学 一种基于视觉惯导信息融合的定位方法及系统
CN109360228A (zh) * 2018-10-18 2019-02-19 清华大学苏州汽车研究院(吴江) 单目摄像头与毫米波雷达间的位姿配准方法
CN110126824A (zh) * 2019-05-22 2019-08-16 河南工业大学 一种集成双目摄像头和毫米波雷达的商用车aebs系统
CN110345944A (zh) * 2019-05-27 2019-10-18 浙江工业大学 融合视觉特征和imu信息的机器人定位方法
CN110390697A (zh) * 2019-07-11 2019-10-29 浙江大学 一种基于lm算法的毫米波雷达与摄像头联合标定方法
CN110517303A (zh) * 2019-08-30 2019-11-29 的卢技术有限公司 一种基于双目相机和毫米波雷达的融合slam方法及系统

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107869989A (zh) * 2017-11-06 2018-04-03 东北大学 一种基于视觉惯导信息融合的定位方法及系统
CN109360228A (zh) * 2018-10-18 2019-02-19 清华大学苏州汽车研究院(吴江) 单目摄像头与毫米波雷达间的位姿配准方法
CN110126824A (zh) * 2019-05-22 2019-08-16 河南工业大学 一种集成双目摄像头和毫米波雷达的商用车aebs系统
CN110345944A (zh) * 2019-05-27 2019-10-18 浙江工业大学 融合视觉特征和imu信息的机器人定位方法
CN110390697A (zh) * 2019-07-11 2019-10-29 浙江大学 一种基于lm算法的毫米波雷达与摄像头联合标定方法
CN110517303A (zh) * 2019-08-30 2019-11-29 的卢技术有限公司 一种基于双目相机和毫米波雷达的融合slam方法及系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
刘振彬 等: "基于单目及惯导的同步定位与建图方案", 《测绘科学》 *
王录涛 等: "基于图优化的视觉SLAM研究进展与应用分析", 《计算机应用研究》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112230211A (zh) * 2020-10-15 2021-01-15 长城汽车股份有限公司 车辆定位方法、装置、存储介质及车辆
CN113030944A (zh) * 2021-04-16 2021-06-25 深圳市众云信息科技有限公司 一种雷达目标追踪方法
CN113030944B (zh) * 2021-04-16 2024-02-02 深圳市众云信息科技有限公司 一种雷达目标追踪方法
CN113625271A (zh) * 2021-07-29 2021-11-09 中汽创智科技有限公司 基于毫米波雷达和双目相机的同时定位与建图方法
CN113625271B (zh) * 2021-07-29 2023-10-27 中汽创智科技有限公司 基于毫米波雷达和双目相机的同时定位与建图方法
WO2023142814A1 (zh) * 2022-01-30 2023-08-03 中国第一汽车股份有限公司 目标识别方法、装置、设备及存储介质

Similar Documents

Publication Publication Date Title
CN109345588B (zh) 一种基于Tag的六自由度姿态估计方法
CN111538029A (zh) 一种视觉与雷达融合的测量方法及终端
US11830216B2 (en) Information processing apparatus, information processing method, and storage medium
CN111795686B (zh) 一种移动机器人定位与建图的方法
CN111524194B (zh) 一种激光雷达和双目视觉相互融合的定位方法及终端
WO2018142900A1 (ja) 情報処理装置、データ管理装置、データ管理システム、方法、及びプログラム
CN112102458A (zh) 基于激光雷达点云数据辅助的单镜头三维图像重构方法
US20140300736A1 (en) Multi-sensor camera recalibration
CN109727278B (zh) 一种机载LiDAR点云数据与航空影像的自动配准方法
CN107330927B (zh) 机载可见光图像定位方法
CN114474056B (zh) 一种面向抓取操作的单目视觉高精度目标定位方法
CN112465877B (zh) 一种基于运动状态估计的卡尔曼滤波视觉追踪稳定方法
CN115187658B (zh) 多相机视觉的大目标定位方法、系统及设备
CN114714356A (zh) 基于双目视觉的工业机器人手眼标定误差精确检测方法
CN114494629A (zh) 一种三维地图的构建方法、装置、设备及存储介质
CN113052897A (zh) 定位初始化方法和相关装置、设备、存储介质
JP6922348B2 (ja) 情報処理装置、方法、及びプログラム
CN113450334B (zh) 一种水上目标检测方法、电子设备及存储介质
CN112113564B (zh) 一种基于图像传感器和惯性传感器的定位方法及系统
CN114078163A (zh) 激光雷达与可见光相机的精确标定方法
CN114485648B (zh) 一种基于仿生复眼惯性系统的导航定位方法
CN115018922A (zh) 畸变参数标定方法、电子设备和计算机可读存储介质
CN113504385B (zh) 复数相机测速方法及测速装置
CN112767482B (zh) 一种多传感器融合的室内外定位方法及系统
CN113379840B (zh) 一种基于共面目标的单目视觉位姿估计方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20200814