CN109270534B - 一种智能车激光传感器与相机在线标定方法 - Google Patents

一种智能车激光传感器与相机在线标定方法 Download PDF

Info

Publication number
CN109270534B
CN109270534B CN201810426821.XA CN201810426821A CN109270534B CN 109270534 B CN109270534 B CN 109270534B CN 201810426821 A CN201810426821 A CN 201810426821A CN 109270534 B CN109270534 B CN 109270534B
Authority
CN
China
Prior art keywords
laser
camera
image
data
calibration
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.)
Active
Application number
CN201810426821.XA
Other languages
English (en)
Other versions
CN109270534A (zh
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.)
Xian Jiaotong University
Original Assignee
Xian Jiaotong University
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 Xian Jiaotong University filed Critical Xian Jiaotong University
Priority to CN201810426821.XA priority Critical patent/CN109270534B/zh
Publication of CN109270534A publication Critical patent/CN109270534A/zh
Application granted granted Critical
Publication of CN109270534B publication Critical patent/CN109270534B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88Radar or analogous systems specially adapted for specific applications
    • G01S13/93Radar or analogous systems specially adapted for specific applications for anti-collision purposes
    • G01S13/931Radar or analogous systems specially adapted for specific applications for anti-collision purposes of land vehicles
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/10Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration
    • G01C21/12Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning
    • G01C21/16Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation
    • G01C21/165Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation combined with non-inertial navigation instruments
    • 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
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/86Combinations of radar systems with non-radar systems, e.g. sonar, direction finder
    • G01S13/867Combination of radar systems with cameras
    • 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
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/48Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
    • G01S7/497Means for monitoring or calibrating

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Automation & Control Theory (AREA)
  • Electromagnetism (AREA)
  • Length Measuring Devices By Optical Means (AREA)
  • Optical Radar Systems And Details Thereof (AREA)
  • Traffic Control Systems (AREA)

Abstract

本发明公开了一种智能车激光传感器与相机在线标定方法,通过相机标定、三维激光传感器与图像传感器离线标定、激光数据与图像数据时间对齐、激光传感器与图像传感器的在线对齐等方法步骤将激光数据与图像数据进行精确的标定。该标定系统可适应于多种不同路况和场景,实现激光和图像的在线的高精度标定。将激光和图像标定起来以后,可利用两种传感器的信息对某一障碍物进行综合分析并作出准确决策,对智能车的感知技术有着重要的意义,因此这一技术可广泛应用于无人驾驶汽车视觉导航、智能车视觉辅助驾驶等领域。

Description

一种智能车激光传感器与相机在线标定方法
技术领域
本发明属于计算机视觉与智能交通领域,特别涉及一种智能车激光传感器与相机在线标定方法。
背景技术
智能车辆(Intelligent Vehicle,IV)是集环境感知、规划决策、多等级辅助驾驶等功能于一体的综合系统,它集中运用了计算机、现代传感、信息融合、通讯、人工智能及自动控制等技术,是典型的高新技术综合体,是一个国家科研实力和产业水平的重要标志。环境感知技术是智能车三大关键技术之一,主动传感激光系统,即三维激光雷达是其中的重要传感器,可以感知障碍物的深度信息,并取得了部分成功的应用,但该传感器存在分辨率低、无法感知RGB信息等问题,制约了其在智能车技术研发和应用上的推广。而被动可见光传感,即相机,在分辨率、感知RGB信息等方面具有显著优势,但无法感知深度信息。近年来,大量研究团队和机构在利用可见光传感完成交通场景感知方面做出了很多卓有成效的研究,基于视觉的交通要素和障碍物感知成为该领域的研究热点。
将三维激光雷达和相机标定起来以后,就可以充分利用两种传感器的优势,综合RGB信息和深度信息,对所感知的物体类别和位置做出准确判断,使无人车做出相关决策。但传统的离线标定方法难以获得精确的标定参数,在车辆运动过程中也存在传感器的偏移问题使得标定参数发生改变。因此,如何设计一套可靠的标定系统使得无人车在长时间的运行中能够高精度的自动在线修正标定参数在无人车的感知技术中非常关键。
发明内容
本发明的目的在于提供一种智能车激光传感器与相机在线标定方法,以解决上述问题。
为实现上述目的,本发明采用以下技术方案:
一种智能车激光传感器与相机在线标定方法,包括以下步骤:
步骤1,将GPS/INS惯性导航系统、工业相机和三维激光雷达安装在智能车上,并通过GPS/INS惯性导航系统、工业相机和三维激光雷达采集依次采集位姿数据、图像数据和激光数据;
步骤2,步骤1采集的数据通过工控机存储,并连接千兆网交换机;手动选取激光数据和图像数据中的标定板,然后对相机内参进行标定,得到精确相机内参;对激光雷达与相机外参进行标定,得到激光与相机初始外参;
步骤3,步骤1的位姿数据,利用坐标变换对自然场景激光数据和图像数据进行时间对齐,并校正激光数据的畸变;
步骤4,对步骤3时间对齐后的激光数据和图像数据分别提取激光数据和图像数据的边缘信息,并对提取的激光数据和图像数据进行预处理;并对图像边界图进行反向距离变换,通过梯度下降的方法优化外部标定参数使得激光与图像的边缘得到精确对齐,优化得到的外部标定参数就是所需的最终标定参数。
进一步的,对相机内参进行标定,根据标定板的实际大小和标定板中黑白格在图像中的位置计算出标定板在相机坐标系中的位置,再通过最小化重投影误差的方式计算出所需标定相机的内部参数;相机的内部参数包括相机的焦距、主点、偏态系数以及畸变系数;标定板为黑白相间的棋盘格。
进一步的,对激光雷达与相机外参进行标定,计算标定板在相机坐标系中的位置和激光坐标系中的位置,通过计算两个坐标系的外参使得激光中的标定板和图像中的标定板对齐;激光与相机的外部参数包括相机坐标系与激光坐标系的旋转关系和平移关系;具体方法如下:
利用标定相机内参是使用的标定板完成相机与激光坐标系的联合标定;将标定板摆放不同的姿态和位置,保证其在两个坐标系均可观测,然后采集对于那个的图像和激光数据组成数据对;对于图像数据,利用张正友标定法标定时定义标定板黑白格图形的两条边缘为X轴和Y轴,垂直于标定板的平面方向为Z轴;对标定板图像数据进行处理得到标定板坐标系到相机坐标系的旋转平移变换关系,并以标定板平面法向量θc,i和其到摄像机坐标系距离αc,i表示第i帧图像中标定板相对于摄像机坐标系的位姿;对于激光数据,手动选取出位于标定板平面的激光数据点后,利用RANSAC方法和最小二乘法鲁棒的估计第i帧标定板平面相对于激光坐标系的位姿,其中标定板平面法向量用θl,i表示,平面到激光坐标系原点的距离为αl,i;由此,标定激光坐标系与摄像机坐标系之间的位姿变换即求解使标定板通过两种不同的数据观测方式在同一坐标系中的表示差别最小的旋转矩阵R与平移向量t。
进一步的,步骤3中,获得位姿数据,并获得图像时间戳和激光时间戳,使用线性插值的方法得到获取每个激光点时刻车体的位姿和获取图像时刻的车体位姿;通过坐标转换对三维激光点的值进行重新计算,从而校正激光的畸变并将其和图像实现时间对齐。
进一步的,激光数据的预处理包括使用车辆位姿校正激光数据和提取激光点云中的边界点;通过时间戳得到采集到图像的时刻t0和获取到某个激光点的时刻t1;通过线性插值可以获得t0时刻对应的位姿和t1时刻对应的位姿;这两个时刻位姿对应的旋转矩阵和车辆坐标分别为R0,p0和R1,p1;由于车辆的运动,车辆上的激光坐标系也在不断运动,激光校正就是将每个激光点在t1时刻的激光坐标系下的坐标(x1,y1,z1)T转换到t0时刻下激光坐标系下的坐标(x0,y0,z0)T,转换公式如下:
Figure BDA0001652344260000031
激光点云数据共有64线,每个点都有一个标号表明其属于第几条线,提取激光边界点需要单独在每束线上分别进行提取;在对激光点进行校正后,首先使用高斯过程去除地面上的激光点;分割掉路面激光点后,对每束线上的激光点按角度排列,再用聚类算法在该束线上进行聚类,将距离较近的激光点聚为一类,然后就可以得到若干的类;某个类的左/右类间距离定义为该类的最左/右侧点与其他类中的点的最小距离;然后选择所含激光点数较多的类且左或右类间距离大于阈值θ的类,如果激光的左/右间距大于θ那么就令左/右侧激光点即为激光边界点;令激光边界点的值为1,非边界点的值为0;得到激光点云的二值边界图像X;用数学语言描述即为Xx,y,z∈(0,1);其中(x,y,z)为激光点的三维坐标。
进一步的,图像数据预处理包括提取图像边界和对图像边界进行反向距离变换;图像边界提取使用结构化随机森林提取图像边界点,得到边界图像,它提取图像中轮廓信息,并给每个图像中每个像素点赋予不同的权重;然后对这些点按阈值分割,得到一个二值边界图E,其中边界点为1,非边界点为0;随后对图像进行反向距离变换,得到预处理所需的图像D;D中坐标为(i,j)的像素点的灰度值的计算公式如下式所示:
Figure BDA0001652344260000041
其中x,y为图像中任意一点的坐标值,γ为一常数。
进一步的,步骤4中优化外部参数采用迭代的方法来优化目标函数;目标函数如下式所示:
Figure BDA0001652344260000042
其中f遍历选择的所有帧,p遍历Xf中所有的三维激光点,(i,j)是第f帧激光点云第p个激光点通过离线方法得到的相机内参所映射在图像上的坐标;假设外参对应的3*3的旋转矩阵和3*1的平移向量分别为R和t,设相机内参中的焦距为f,主点为u0,v0,像素点在x方向和y方向的长度分别为dx和dy,激光点坐标为(x,y,z);采用齐次坐标的形式的话,(i,j)可用下式计算得到:
Figure BDA0001652344260000043
其中k为尺度因子;
外参包括6个参数,有三个旋转角度,即偏航角
Figure BDA0001652344260000044
俯仰角
Figure BDA0001652344260000045
以及横滚角
Figure BDA0001652344260000046
并且这三个角度唯一的对应着一个旋转矩阵R;有三个方向的平移量,在坐标轴x方向的平移量tx,y方向的平移量ty和z方向的平移tz,这三个平移量唯一对应着一个平移向量t;对当前外参进行扰动,即对每个角度加减0.1度或不变,对每个平移量加减0.1m或不变,即可得到729个新的邻域参数;然后在这些参数里找一个对应目标函数最大的参数作为新的外部参数,如此循环,直到外参不再改变,停止迭代,所得外参即为所要求的参数。
与现有技术相比,本发明有以下技术效果:
本发明实现对激光雷达与相机的高精度标定,利用离线标定模块提供的初始值,并进一步进行在线优化,从而完成激光雷达与相机的高精度标定。
本发明采用低成本的标定物和计算平台,对所使用的自然场景没有过高的要求,可广泛应用于使用了各种配置方案无人车的激光雷达和相机的高精度动态标定,对无人车的感知技术有着十分重要的意义。
附图说明
图1为本发明的流程框图。
图2为本发明的系统硬件关系图。
图3为本发明激光坐标系和摄像机坐标系的定义示意图,其中:(a)侧视图,(b)俯视图。
图4为本发明的反向距离变换示意图。其中上部为原图,下部为变换后的图片。
具体实施方式
以下结合附图对本发明进一步说明:
一种智能车激光传感器与相机在线标定方法,包括以下步骤:
步骤1,将GPS/INS惯性导航系统、工业相机和三维激光雷达安装在智能车上,并通过GPS/INS惯性导航系统、工业相机和三维激光雷达采集依次采集位姿数据、图像数据和激光数据;
步骤2,步骤1采集的数据通过工控机存储,并连接千兆网交换机;手动选取激光数据和图像数据中的标定板,然后对相机内参进行标定,得到精确相机内参;对激光雷达与相机外参进行标定,得到激光与相机初始外参;
步骤3,步骤1的位姿数据,利用坐标变换对自然场景激光数据和图像数据进行时间对齐,并校正激光数据的畸变;
步骤4,对步骤3时间对齐后的激光数据和图像数据分别提取激光数据和图像数据的边缘信息,并对提取的激光数据和图像数据进行预处理;并对图像边界图进行反向距离变换,通过梯度下降的方法优化外部标定参数使得激光与图像的边缘得到精确对齐,优化得到的外部标定参数就是所需的最终标定参数。
该激光雷达与相机高精度标定系统工作原理:主要包括离线标定和在线标定两部分,离线标定法利用张正友标定法得到精确的相机内参和CMU激光相机联合标定工具箱完成激光雷达和相机的外部参数;在线标定法首先获得车辆的全局位姿数据,并对三维激光数据进行坐标变换,完成畸变校正和时间对齐。然后使用聚类的算法提取激光中的边界点,使用结构化随机森林的方法提取对应于激光的图像中的轮廓,并对其进行二值化得到边界图像,再对边界图像进行反向距离变换从而为距离边界不同的像素赋予不同的灰度值。最后,利用先前得到的相机内参和投射投影原理将激光边界点映射在图像中,并迭代优化激光雷达与相机的外部标定参数使得映射在图像中的激光点的对应的像素灰度值最大。最后优化得到的外部标定参数和预先标定好的相机内部参数即为所需的高精度标定参数。
目前该系统已应用于西安交通大学人工智能与机器人研究所视觉认知计算与智能车实验室的“夸父一号”无人驾驶平台上激光雷达和相机的标定结果,并顺利在各种不同场景下完成测试,每次均能得到可靠有效的标定参数,表明该算法是稳定鲁棒的。
(一)实验平台硬件介绍
1)相机和镜头:
采集道路场景使用的相机为加拿大PointGrey公司的GrassHopper2高分辨率CCD摄像机,其具体参数如表所示。摄像机配套的镜头为Kowa手动光圈百万像素镜头。本发明中选择的Kowa镜头焦距为25mm,其详细参数如表所示。
表1 GrassHopper2摄像机规格指标
Figure BDA0001652344260000061
Figure BDA0001652344260000071
表2 Kowa镜头规格指标
Figure BDA0001652344260000072
2)GPS/INS惯导系统:
实验平台安装的惯导系统为加拿大NovAtel公司的SPAN-LCI。由惯性测量单元(IMU)、GPS接收机和GPS天线等组成,同时,配备了数据传输单元(DTU)传输差分改正数,组合多种定位设备实现高精度的全局定位。
使用的惯导系统的各项性能指标如表所示:
表3 SPAN-LCI惯导系统性能指标
Figure BDA0001652344260000073
Figure BDA0001652344260000081
3)工控机:
实验平台使用的工控机为凌华公司的车载刀片服务器ADLINK嵌入式电脑,体积小、可挂式安装、运算速度快,同时配置了固态硬盘,其指标如表4所示:
表4 ADLINK工控机性能指标
Figure BDA0001652344260000082
4)三维激光雷达:
实验平台上安装的三维激光雷达是美国Velodyne公司生产的64线激光雷达HDL-64E,其扫描线数多、覆盖范围广、数据量丰富,可以较逼真的还原场景。
各硬件的连接关系如图1所示。系统计算负载在车载刀片服务器ADLINK嵌入式电脑(分为位姿工控机、交通灯工控机以及激光工控机)上进行;所有模块间通过千兆以太网,以UDP数据包方式进行通信。系统成本低廉、功耗低,且具有较高的整体移植性,适合批量推广应用。
激光雷达与相机在线标定时,需要对齐激光数据和图像数据。激光数据从三维激光雷达获取,原始数据通过千兆以太网传输到激光工控机,经过数据解析从极坐标系转换到激光坐标系下,同时也存储了经过时间戳同步后的激光数据对应时间戳。图像数据通过摄像机抓取,然后通过千兆网接口传输到交通灯工控机上,同时保存对准后的图像时间戳。获取两种传感器的原始数据后,开始进行对齐操作,完成标定后保存相机与激光标定的结果到本地硬盘中供其他模块使用。
(二)本发明的软件系统介绍如下:
1)开发平台
本发明采用windows 10操作系统作为软件的运行平台,Windows操作系统具有支持图形化显示和多任务机制,操作简单,实时性高,运行稳定等优点。程序主要使用Matlab语言开发,同时部分开发在Visual Studio 2015编译环境下使用C++语言完成,并使用了OpenCV库辅助图像处理过程。
2)系统框架设计
整个系统如图2所示,输入为GPS数据、图像序列数据和激光序列数据,设计工作包含两大部分:(1)相机内参以及激光雷达和相机之间外参的离线标定;(2)在线标定激光雷达与相机。
离线过程包括相机内参标定和相机外参标定。
(1)相机内参标定
本发明采用广泛使用的张正友标定法标定相机内部参数(即相机内参)。标定的思路为,首先只考虑相机的4个线性参数进行初步标定,计算得到内参初始值,在此基础上标定畸变参数,然后以线性参数为初始值,对标定参数进行重新优化,直至得到收敛的线性和非线性参数。内参标定主要用于坐标系之间的旋转平移参数的求解。
(2)激光雷达与相机外参标定
对于激光雷达与相机之间外参的离线标定,本发明利用标定相机内参是使用的标定板完成相机与激光坐标系的联合标定。将标定板摆放不同的姿态和位置,保证其在两个坐标系均可观测,然后采集对于那个的图像和激光数据组成数据对。对于图像数据,利用张正友标定法标定时定义标定板黑白格图形的两条边缘为X轴和Y轴,垂直于标定板的平面方向为Z轴。对标定板图像数据进行处理,不仅可以返回相机的内部参数,还可以得到标定板坐标系到相机坐标系的旋转平移变换关系,并以标定板平面法向量θc,i和其到摄像机坐标系距离αc,i表示第i帧图像中标定板相对于摄像机坐标系的位姿。对于激光数据,手动选取出位于标定板平面的激光数据点后,利用RANSAC方法和最小二乘法鲁棒的估计第i帧标定板平面相对于激光坐标系的位姿,其中标定板平面法向量用θl,i表示,平面到激光坐标系原点的距离为αl,i。由此,标定激光坐标系与摄像机坐标系之间的位姿变换即求解使标定板通过两种不同的数据观测方式在同一坐标系中的表示差别最小的旋转矩阵R与平移向量t。
在线过程中主要包括激光预处理、图像预处理和优化标定参数三个重要过程。
(1)激光预处理
激光预处理包括使用车辆位姿校正激光数据和提取激光点云中的边界点。通过时间戳可知采集到图像的时刻t0和获取到某个激光点的时刻t1。通过线性插值可以获得t0时刻对应的位姿和t1时刻对应的位姿。这两个时刻位姿对应的旋转矩阵和车辆坐标分别为R0,p0和R1,p1。由于车辆的运动,车辆上的激光坐标系也在不断运动,激光校正就是将每个激光点在t1时刻的激光坐标系下的坐标(x1,y1,z1)T转换到t0时刻下激光坐标系下的坐标(x0,y0,z0)T,转换公式如下:
Figure BDA0001652344260000101
激光点云数据共有64线,每个点都有一个标号表明其属于第几条线,提取激光边界点需要单独在每束线上分别进行提取。在对激光点进行校正后,首先使用高斯过程去除地面上的激光点,这一方法由谌彤童在2012年提出。分割掉路面激光点后,对每束线上的激光点按角度排列,再用聚类算法在该束线上进行聚类,将距离较近的激光点聚为一类,然后就可以得到若干的类。某个类的左/右类间距离定义为该类的最左/右侧点与其他类中的点的最小距离。然后选择所含激光点数较多的类且左或右类间距离大于阈值θ的类,如果激光的左/右间距大于θ那么就令左/右侧激光点即为激光边界点。令激光边界点的值为1,非边界点的值为0。如此我们即可得到激光点云的二值边界图像X。用数学语言描述即为Xx,y,z∈(0,1)。其中(x,y,z)为激光点的三维坐标。
(2)图像预处理
图像预处理包括提取图像边界和对图像边界进行反向距离变换。图像边界提取使用结构化随机森林提取图像边界点,得到边界图像,这一方法由Piotr在2013年提出,它可以提取图像中轮廓信息,并给每个图像中每个像素点赋予不同的权重。然后对这些点按阈值分割,得到一个二值边界图E,其中边界点为1,非边界点为0。随后对图像进行反向距离变换,得到预处理所需的图像D。D中坐标为(i,j)的像素点的灰度值的计算公式如下式所示:
Figure BDA0001652344260000102
其中x,y为图像中任意一点的坐标值,γ为一常数,本发明设其为0.98。
(3)优化外部参数
优化外部参数采用迭代的方法来优化目标函数。目标函数如下式所示:
Figure BDA0001652344260000111
其中f遍历选择的所有帧,p遍历Xf中所有的三维激光点,(i,j)是第f帧激光点云第p个激光点通过离线方法得到的相机内参所映射在图像上的坐标。假设外参对应的3*3的旋转矩阵和3*1的平移向量分别为R和t,设相机内参中的焦距为f,主点为u0,v0,像素点在x方向和y方向的长度分别为dx和dy,激光点坐标为(x,y,z)。采用齐次坐标的形式的话,(i,j)可用下式计算得到:
Figure BDA0001652344260000112
其中k为尺度因子。
外参包括6个参数,有三个旋转角度,即偏航角
Figure BDA0001652344260000113
俯仰角
Figure BDA0001652344260000114
以及横滚角
Figure BDA0001652344260000115
并且这三个角度唯一的对应着一个旋转矩阵R;有三个方向的平移量,在坐标轴x方向的平移量tx,y方向的平移量ty和z方向的平移tz,这三个平移量唯一对应着一个平移向量t。对当前外参进行扰动,即对每个角度加减0.1度或不变,对每个平移量加减0.1m或不变,即可得到729个新的邻域参数。然后在这些参数里找一个对应目标函数最大的参数作为新的外部参数,如此循环,直到外参不再改变,停止迭代,所得外参即为所要求的参数。
与传统的激光雷达相机标定方法比,本发明充分利用了自然场景中的信息,使得标定精度具有显著的提高,并且可以自动地更新标定参数以校准传感器的偏移。利用位姿数据对激光点云进行畸变校正,有效改善了激光点云在车辆高速运动场景下的畸变问题。可广泛用于配备了三维激光雷达和摄像机的无人驾驶汽车中。

Claims (5)

1.一种智能车激光传感器与相机在线标定方法,其特征在于,包括以下步骤:
步骤1,将GPS/INS惯性导航系统、工业相机和三维激光雷达安装在智能车上,并通过GPS/INS惯性导航系统、工业相机和三维激光雷达采集依次采集位姿数据、图像数据和激光数据;
步骤2,步骤1采集的数据通过工控机存储,并连接千兆网交换机;手动选取激光数据和图像数据中的标定板,然后对相机内参进行标定,得到精确相机内参;对激光雷达与相机外参进行标定,得到激光与相机初始外参;
步骤3,步骤1的位姿数据,利用坐标变换对自然场景激光数据和图像数据进行时间对齐,并校正激光数据的畸变;
步骤4,对步骤3时间对齐后的激光数据和图像数据分别提取激光数据和图像数据的边缘信息,并对提取的激光数据和图像数据进行预处理;并对图像边界图进行反向距离变换,通过梯度下降的方法优化外部标定参数使得激光与图像的边缘得到精确对齐,优化得到的外部标定参数就是所需的最终标定参数;
激光数据的预处理包括使用车辆位姿校正激光数据和提取激光点云中的边界点;通过时间戳得到采集到图像的时刻t0和获取到某个激光点的时刻t1;通过线性插值可以获得t0时刻对应的位姿和t1时刻对应的位姿;这两个时刻位姿对应的旋转矩阵和车辆坐标分别为R0,p0和R1,p1;由于车辆的运动,车辆上的激光坐标系也在不断运动,激光校正就是将每个激光点在t1时刻的激光坐标系下的坐标(x1,y1,z1)T转换到t0时刻下激光坐标系下的坐标(x0,y0,z0)T,转换公式如下:
Figure FDA0002599384980000011
激光点云数据共有64线,每个点都有一个标号表明其属于第几条线,提取激光边界点需要单独在每束线上分别进行提取;在对激光点进行校正后,首先使用高斯过程去除地面上的激光点;分割掉路面激光点后,对每束线上的激光点按角度排列,再用聚类算法在该束线上进行聚类,将距离较近的激光点聚为一类,然后就可以得到若干的类;某个类的左/右类间距离定义为该类的最左/右侧点与其他类中的点的最小距离;然后选择所含激光点数较多的类且左或右类间距离大于阈值θ的类,如果激光的左/右间距大于θ那么就令左/右侧激光点即为激光边界点;令激光边界点的值为1,非边界点的值为0;得到激光点云的二值边界图像X;用数学语言描述即为Xx,y,z∈(0,1);其中(x,y,z)为激光点的三维坐标;
图像数据预处理包括提取图像边界和对图像边界进行反向距离变换;图像边界提取使用结构化随机森林提取图像边界点,得到边界图像,它提取图像中轮廓信息,并给每个图像中每个像素点赋予不同的权重;然后对这些点按阈值分割,得到一个二值边界图E,其中边界点为1,非边界点为0;随后对图像进行反向距离变换,得到预处理所需的图像D;D中坐标为(i,j)的像素点的灰度值的计算公式如下式所示:
Figure FDA0002599384980000021
其中x,y为图像中任意一点的坐标值,γ为一常数。
2.根据权利要求1所述的一种智能车激光传感器与相机在线标定方法,其特征在于,对相机内参进行标定,根据标定板的实际大小和标定板中黑白格在图像中的位置计算出标定板在相机坐标系中的位置,再通过最小化重投影误差的方式计算出所需标定相机的内部参数;相机的内部参数包括相机的焦距、主点、偏态系数以及畸变系数;标定板为黑白相间的棋盘格。
3.根据权利要求1所述的一种智能车激光传感器与相机在线标定方法,其特征在于,对激光雷达与相机外参进行标定,计算标定板在相机坐标系中的位置和激光坐标系中的位置,通过计算两个坐标系的外参使得激光中的标定板和图像中的标定板对齐;激光与相机的外部参数包括相机坐标系与激光坐标系的旋转关系和平移关系;具体方法如下:
利用标定相机内参是使用的标定板完成相机与激光坐标系的联合标定;将标定板摆放不同的姿态和位置,保证其在两个坐标系均可观测,然后采集对应的图像和激光数据组成数据对;对于图像数据,利用张正友标定法标定时定义标定板黑白格图形的两条边缘为X轴和Y轴,垂直于标定板的平面方向为Z轴;对标定板图像数据进行处理得到标定板坐标系到相机坐标系的旋转平移变换关系,并以标定板平面法向量θc,i和其到摄像机坐标系距离αc,i表示第i帧图像中标定板相对于摄像机坐标系的位姿;对于激光数据,手动选取出位于标定板平面的激光数据点后,利用RANSAC方法和最小二乘法鲁棒的估计第i帧标定板平面相对于激光坐标系的位姿,其中标定板平面法向量用θl,i表示,平面到激光坐标系原点的距离为αl,i;由此,标定激光坐标系与摄像机坐标系之间的位姿变换即求解使标定板通过两种不同的数据观测方式在同一坐标系中的表示差别最小的旋转矩阵R与平移向量t。
4.根据权利要求1所述的一种智能车激光传感器与相机在线标定方法,其特征在于,步骤3中,获得位姿数据,并获得图像时间戳和激光时间戳,使用线性插值的方法得到获取每个激光点时刻车体的位姿和获取图像时刻的车体位姿;通过坐标转换对三维激光点的值进行重新计算,从而校正激光的畸变并将其和图像实现时间对齐。
5.根据权利要求1所述的一种智能车激光传感器与相机在线标定方法,其特征在于,步骤4中优化外部参数采用迭代的方法来优化目标函数;目标函数如下式所示:
Figure FDA0002599384980000031
其中f遍历选择的所有帧,p遍历Xf中所有的三维激光点,(i,j)是第f帧激光点云第p个激光点通过离线方法得到的相机内参所映射在图像上的坐标;假设外参对应的3*3的旋转矩阵和3*1的平移向量分别为R和t,设相机内参中的焦距为f,主点为u0,v0,像素点在x方向和y方向的长度分别为dx和dy,激光点坐标为(x,y,z);采用齐次坐标的形式的话,(i,j)可用下式计算得到:
Figure FDA0002599384980000032
其中k为尺度因子;
外参包括6个参数,有三个旋转角度,即偏航角
Figure FDA0002599384980000033
俯仰角
Figure FDA0002599384980000034
以及横滚角
Figure FDA0002599384980000035
并且这三个角度唯一的对应着一个旋转矩阵R;有三个方向的平移量,在坐标轴x方向的平移量tx,y方向的平移量ty和z方向的平移tz,这三个平移量唯一对应着一个平移向量t;对当前外参进行扰动,即对每个角度加减0.1度或不变,对每个平移量加减0.1m或不变,即可得到729个新的邻域参数;然后在这些参数里找一个对应目标函数最大的参数作为新的外部参数,如此循环,直到外参不再改变,停止迭代,所得外参即为所要求的参数。
CN201810426821.XA 2018-05-07 2018-05-07 一种智能车激光传感器与相机在线标定方法 Active CN109270534B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810426821.XA CN109270534B (zh) 2018-05-07 2018-05-07 一种智能车激光传感器与相机在线标定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810426821.XA CN109270534B (zh) 2018-05-07 2018-05-07 一种智能车激光传感器与相机在线标定方法

Publications (2)

Publication Number Publication Date
CN109270534A CN109270534A (zh) 2019-01-25
CN109270534B true CN109270534B (zh) 2020-10-27

Family

ID=65152624

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810426821.XA Active CN109270534B (zh) 2018-05-07 2018-05-07 一种智能车激光传感器与相机在线标定方法

Country Status (1)

Country Link
CN (1) CN109270534B (zh)

Families Citing this family (48)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109765901A (zh) * 2019-02-18 2019-05-17 华南理工大学 基于线激光与双目视觉的动态代价地图导航方法
CN111650598A (zh) * 2019-02-19 2020-09-11 北京京东尚科信息技术有限公司 一种车载激光扫描系统外参标定方法和装置
CN111665470A (zh) * 2019-03-07 2020-09-15 阿里巴巴集团控股有限公司 一种定位方法及装置和机器人
CN111678533A (zh) * 2019-03-11 2020-09-18 武汉小狮科技有限公司 一种惯性导航设备与激光雷达的标定方法
CN111678534A (zh) * 2019-03-11 2020-09-18 武汉小狮科技有限公司 一种结合rgbd双目深度相机、imu和多线激光雷达的联合标定平台和方法
CN109919893B (zh) * 2019-03-20 2021-04-23 湖北亿咖通科技有限公司 点云修正方法、装置和可读存储介质
CN110068814B (zh) * 2019-03-27 2021-08-24 东软睿驰汽车技术(沈阳)有限公司 一种测量障碍物距离的方法及装置
CN109828250B (zh) * 2019-03-28 2020-07-21 森思泰克河北科技有限公司 一种雷达标定方法、标定装置及终端设备
CN110097620A (zh) * 2019-04-15 2019-08-06 西安交通大学 基于图像和三维激光的高精度地图创建系统
CN111862224B (zh) * 2019-04-17 2023-09-19 杭州海康威视数字技术股份有限公司 确定相机与激光雷达之间外参的方法和装置
CN110044379B (zh) * 2019-04-19 2020-10-16 武汉地震计量检定与测量工程研究院有限公司 一种移动测量系统标定方法
CN112214009B (zh) * 2019-06-25 2022-07-26 上海商汤临港智能科技有限公司 传感器数据处理方法、装置、电子设备及系统
CN112147598A (zh) * 2019-06-27 2020-12-29 武汉爱速达机器人科技有限公司 一种基于直角墙面的激光标定方法
CN110310339B (zh) * 2019-07-04 2021-05-04 厦门金龙联合汽车工业有限公司 一种标定装置及一种雷达与摄像头联合标定方法及系统
CN110244284B (zh) * 2019-07-29 2023-05-05 南通润邦重机有限公司 一种用于多线激光雷达和gps\ins标定的标定板及其方法
US11846731B2 (en) 2019-07-31 2023-12-19 Canon Kabushiki Kaisha Distance detection device and imaging apparatus
CN112396630A (zh) * 2019-08-15 2021-02-23 纳恩博(北京)科技有限公司 目标对象状态的确定方法、装置、存储介质及电子装置
CN110501036A (zh) * 2019-08-16 2019-11-26 北京致行慕远科技有限公司 传感器参数的标定检查方法及装置
CN110470333B (zh) * 2019-08-16 2022-05-24 纳恩博(常州)科技有限公司 传感器参数的标定方法及装置、存储介质和电子装置
CN110766758B (zh) * 2019-09-12 2022-08-09 浙江大华技术股份有限公司 标定方法、装置、系统及存储装置
CN110780285B (zh) * 2019-10-24 2022-10-18 深圳市镭神智能系统有限公司 激光雷达与组合惯导的位姿标定方法、系统及介质
CN110823252B (zh) * 2019-11-06 2022-11-18 大连理工大学 一种多线激光雷达和单目视觉的自动标定方法
CN111123912B (zh) * 2019-11-29 2023-01-31 苏州智加科技有限公司 行车定位坐标的标定方法和装置
CN111103578B (zh) * 2020-01-10 2022-02-15 清华大学 一种基于深度卷积神经网络的激光雷达在线标定方法
CN113494927A (zh) * 2020-03-20 2021-10-12 郑州宇通客车股份有限公司 一种车辆多传感器标定方法、装置及车辆
CN111443337B (zh) * 2020-03-27 2022-03-08 北京航空航天大学 一种基于手眼标定的雷达-imu标定方法
CN111640158B (zh) * 2020-06-11 2023-11-10 武汉斌果科技有限公司 基于对应掩模的端到端的摄像头与激光雷达外参标定方法
CN112017205B (zh) * 2020-07-27 2021-06-25 清华大学 一种激光雷达和相机传感器空间位置自动标定方法及系统
CN111928776A (zh) * 2020-07-31 2020-11-13 中国航空工业集团公司济南特种结构研究所 基于多传感器的数控机床非接触式在线测量系统与方法
CN112184828B (zh) * 2020-08-21 2023-12-05 阿波罗智联(北京)科技有限公司 激光雷达与摄像头的外参标定方法、装置及自动驾驶车辆
CN112146682B (zh) * 2020-09-22 2022-07-19 福建牧月科技有限公司 智能汽车的传感器标定方法、装置、电子设备及介质
CN114322751B (zh) * 2020-09-30 2024-01-23 广东博智林机器人有限公司 目标测量方法、装置、计算机设备和存储介质
CN112230241B (zh) * 2020-10-23 2021-07-20 湖北亿咖通科技有限公司 基于随机扫描型雷达的标定方法
CN112270311B (zh) * 2020-11-25 2023-12-19 武汉理工大学 一种基于车载环视逆投影的近目标快速检测方法及系统
CN112711263A (zh) * 2021-01-19 2021-04-27 未来机器人(深圳)有限公司 仓储自动导引车避障方法、装置、计算机设备和存储介质
CN113147738A (zh) * 2021-02-26 2021-07-23 重庆智行者信息科技有限公司 一种自动泊车定位方法和装置
CN113050108B (zh) * 2021-03-23 2024-01-09 湖南盛鼎科技发展有限责任公司 电子界址视觉测量系统及测量方法
CN112964291B (zh) * 2021-04-02 2023-07-14 清华大学 一种传感器标定的方法、装置、计算机存储介质及终端
CN113091771B (zh) * 2021-04-13 2022-09-23 清华大学 一种激光雷达-相机-惯导联合标定方法及系统
CN113284197B (zh) * 2021-07-22 2021-11-23 浙江华睿科技股份有限公司 Agv的tof相机外参标定方法及装置、电子设备
CN113340334B (zh) * 2021-07-29 2021-11-30 新石器慧通(北京)科技有限公司 用于无人车的传感器标定方法、装置及电子设备
CN113848541B (zh) * 2021-09-22 2022-08-26 深圳市镭神智能系统有限公司 一种标定的方法、装置、无人机及计算机可读存储介质
CN113884099B (zh) * 2021-12-07 2022-04-12 智道网联科技(北京)有限公司 一种路端移动物体位置测量方法及装置
CN114111811A (zh) * 2021-12-17 2022-03-01 奇瑞万达贵州客车股份有限公司 一种自动驾驶公交客车的导航控制系统和方法
CN115092254B (zh) * 2022-07-18 2023-08-18 苏州智加科技有限公司 车轮转向角的获取方法、装置、设备及存储介质
CN115100299B (zh) * 2022-08-29 2023-02-10 广州镭晨智能装备科技有限公司 一种标定方法、装置、设备和存储介质
CN117173257A (zh) * 2023-11-02 2023-12-05 安徽蔚来智驾科技有限公司 3d目标检测及其标定参数增强方法、电子设备、介质
CN117826129A (zh) * 2024-03-04 2024-04-05 南京航空航天大学 一种单目相机和激光雷达在轨外参标定方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106651964A (zh) * 2017-02-28 2017-05-10 嘉兴锐视智能科技有限公司 激光相机的标定系统及标定装置
CN107976668A (zh) * 2016-10-21 2018-05-01 法乐第(北京)网络科技有限公司 一种确定相机与激光雷达之间的外参数的方法

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101698303B (zh) * 2009-09-11 2011-02-02 大连理工大学 一种三维激光和单目视觉间的自动标定方法
CN103049912B (zh) * 2012-12-21 2015-03-11 浙江大学 一种基于任意三面体的雷达-相机系统外部参数标定方法
JP2016048172A (ja) * 2014-08-27 2016-04-07 株式会社トプコン 画像処理装置、画像処理方法、およびプログラム
JP6121063B1 (ja) * 2014-11-04 2017-04-26 エスゼット ディージェイアイ テクノロジー カンパニー リミテッドSz Dji Technology Co.,Ltd カメラ較正方法、デバイス及びシステム
CN104484887B (zh) * 2015-01-19 2017-07-07 河北工业大学 摄像机与二维激光测距仪联合使用时的外参数标定方法
CN105737849A (zh) * 2015-12-24 2016-07-06 武汉武大卓越科技有限责任公司 隧道车用激光扫描仪与相机相对位置的检校方法
CN105930819B (zh) * 2016-05-06 2019-04-12 西安交通大学 基于单目视觉和gps组合导航系统的实时城区交通灯识别系统

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107976668A (zh) * 2016-10-21 2018-05-01 法乐第(北京)网络科技有限公司 一种确定相机与激光雷达之间的外参数的方法
CN106651964A (zh) * 2017-02-28 2017-05-10 嘉兴锐视智能科技有限公司 激光相机的标定系统及标定装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
车载激光雷达外参数的标定方法;程金龙;《光电工程》;20131231;第40卷(第12期);第89-94页 *

Also Published As

Publication number Publication date
CN109270534A (zh) 2019-01-25

Similar Documents

Publication Publication Date Title
CN109270534B (zh) 一种智能车激光传感器与相机在线标定方法
EP4283515A1 (en) Detection method, system, and device based on fusion of image and point cloud information, and storage medium
CN108868268A (zh) 基于点到面距离和互相关熵配准的无人车位姿估计方法
KR20210042274A (ko) 장애물을 탐지하는 방법, 장치, 전자 기기, 저장 매체 및 프로그램
CN111340797A (zh) 一种激光雷达与双目相机数据融合检测方法及系统
CN108594245A (zh) 一种目标物运动监测系统及方法
CN110246175A (zh) 全景相机与云台相机结合的变电站巡检机器人图像检测系统及方法
CN110853037A (zh) 一种基于球面投影的轻量化彩色点云分割方法
CN107300382B (zh) 一种用于水下机器人的单目视觉定位方法
CN111191625A (zh) 一种激光-单目视觉融合的物体识别与定位方法
CN112967344B (zh) 相机外参标定的方法、设备、存储介质及程序产品
CN112288825B (zh) 相机标定方法、装置、电子设备、存储介质和路侧设备
CN106908064B (zh) 一种基于Kinect2传感器的室内夜间视觉导航方法
CN111998862B (zh) 一种基于bnn的稠密双目slam方法
WO2019144289A1 (en) Systems and methods for calibrating an optical system of a movable object
CN113409459A (zh) 高精地图的生产方法、装置、设备和计算机存储介质
CN113095154A (zh) 基于毫米波雷达与单目相机的三维目标检测系统及方法
Ma et al. Crlf: Automatic calibration and refinement based on line feature for lidar and camera in road scenes
CN113985405A (zh) 障碍物检测方法、应用于车辆的障碍物检测设备
CN115984766A (zh) 一种煤矿井下快速单目视觉三维目标检测方法
Mi et al. A vision-based displacement measurement system for foundation pit
CN113327296B (zh) 基于深度加权的激光雷达与相机在线联合标定方法
KR102490521B1 (ko) 라이다 좌표계와 카메라 좌표계의 벡터 정합을 통한 자동 캘리브레이션 방법
CN112580683B (zh) 一种基于互相关的多传感器数据时间对齐系统及其方法
CN103260008A (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
GR01 Patent grant
GR01 Patent grant