CN113093155A - 一种激光雷达联合标定方法及系统 - Google Patents

一种激光雷达联合标定方法及系统 Download PDF

Info

Publication number
CN113093155A
CN113093155A CN202110228154.6A CN202110228154A CN113093155A CN 113093155 A CN113093155 A CN 113093155A CN 202110228154 A CN202110228154 A CN 202110228154A CN 113093155 A CN113093155 A CN 113093155A
Authority
CN
China
Prior art keywords
marker
ground
point
laser radar
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.)
Granted
Application number
CN202110228154.6A
Other languages
English (en)
Other versions
CN113093155B (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.)
Shanghai New Era Robot Co ltd
Original Assignee
Shanghai New Era Robot 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 Shanghai New Era Robot Co ltd filed Critical Shanghai New Era Robot Co ltd
Priority to CN202110228154.6A priority Critical patent/CN113093155B/zh
Publication of CN113093155A publication Critical patent/CN113093155A/zh
Application granted granted Critical
Publication of CN113093155B publication Critical patent/CN113093155B/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
    • 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
    • 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/02Systems using the reflection of electromagnetic waves other than radio waves
    • G01S17/06Systems determining position data of a target
    • G01S17/46Indirect determination of position data

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Electromagnetism (AREA)
  • Optical Radar Systems And Details Thereof (AREA)

Abstract

本发明提供了一种激光雷达联合标定方法及系统,包括步骤:分别通过第一激光雷达和第二激光雷达采集棍棒状的标志物附近的点云数据;分别对两组点云数据进行地面与标志物分割,将地面点拟合成平面方程,将标志物点拟合成直线方程;分别计算标志物的朝向及标志物与地面的交点;分别根据标志物的朝向及对应的标志物与地面的交点计算标志物上预设一点的坐标;分别将第一激光雷达和第二激光雷达观测到标志物与地面的交点,以及标志物上预设一点的坐标组合成一组三维点集;根据两组三维点集计算联合标定结果。该方案的自动化程度更高,标定精度更高,且在标定时不需要看到标志物全部,适用于各种视场的激光雷达的标定。

Description

一种激光雷达联合标定方法及系统
技术领域
本发明涉及雷达标定技术领域,尤指一种激光雷达联合标定方法及系统。
背景技术
激光雷达,是激光探测及测距系统的简称。激光雷达是激光技术与雷达技术相结合的产物,由发射机、天线、接收机、跟踪架及信息处理等部分组成。发射机是各种形式的激光器,如二氧化碳激光器、掺钕钇铝石榴石激光器、半导体激光器及波长可调谐的固体激光器等;天线是光学望远镜;接收机采用各种形式的光电探测器,如光电倍增管、半导体光电二极管、雪崩光电二极管、红外和可见光多元探测器件等。
为了实现机器人的障碍物感知、定位和建图的功能,需要对无人驾驶车辆的双激光雷达传感器进行联合标定,即估计激光雷达2到激光雷达1的位置和旋转的变化,以对准两个激光雷达。现有的激光雷达联合标定方法大多自动化程度低,需要手动标定,且在标定时,需要激光雷达看到标志物全部,不利于小视场的激光雷达的标定。因此,需要一种自动化程度更高,且不需要看到标志物全部,适用范围更广,精度更高的联合标定方法。
发明内容
本发明的目的是提供一种激光雷达联合标定方法及系统,该方案的自动化程度更高,标定精度更高,且在标定时不需要看到标志物全部,适用于各种视场的激光雷达的标定。
本发明提供的技术方案如下:
本发明提供一种激光雷达联合标定方法,包括步骤:
分别通过第一激光雷达和第二激光雷达采集棍棒状的标志物附近的点云数据;
分别对两组所述点云数据进行地面与标志物分割,并分别将地面点拟合成平面方程,将标志物点拟合成直线方程;
分别根据所述平面方程和其对应的所述直线方程计算所述标志物的朝向及所述标志物与地面的交点;
分别根据所述标志物的朝向及对应的所述标志物与地面的交点计算所述标志物上预设一点的坐标;
将所述第一激光雷达观测到所述标志物与地面的交点,以及所述标志物上预设一点的坐标组合成一组三维点集,将所述第二激光雷达观测到的所述标志物与地面的交点,以及所述标志物上预设一点的坐标组合成另一组三维点集;
根据两组所述三维点集计算联合标定结果。
分别通过第一激光雷达和第二激光雷达采集棍棒状的标志物附近的点云数据,并对两组点云数据进行地面与标志物分割,能够将地面点拟合成平面方程,将标志物点拟合成直线方程,通过平面方程和其对应的直线方程能够计算获得标志物的朝向及标志物与地面的交点,进而获得标志物上预设一点的坐标,将第一激光雷达观测到标志物与地面的交点,以及标志物上预设一点的坐标组合成一组三维点集,将第二激光雷达观测到的标志物与地面的交点,以及标志物上预设一点的坐标组合成另一组三维点集,便能够根据两组三维点集计算联合标定结果。由于该方案选用棍棒状的标志物,在采集数据时,对数据的要求较低,仅需要看到部分标志物,便能进行数据计算,适用于各种视场的激光雷达的标定,且数据的分割、拟合和计算都能够自动进行,使得标定的自动化程度更高,有利于提高标定精度。
另外,本方案是用于两个激光雷达的联合标定,在激光雷达的数量超过三个时,可以按照上述方法两两进行标定,得到全部的联合标定结果,从而使得标定更为灵活通用。
具体的,在本方案中,由于对标志物是进行直线拟合,因此,对标志物的设置没有具体要求,可以竖直设置,也可以倾斜设置。
进一步地,所述的分别对两组所述点云数据进行地面与标志物分割,并分别将地面点拟合成平面方程,将标志物点拟合成直线方程,具体包括:
分别对所述第一激光雷达和所述第二激光雷达采集的所述点云数据进行ROI过滤;
采用随机采样一致算法致(RANSAC)分割所述点云数据的地面点和标志物点;
将地面点拟合成平面方程:nT(r-p)=0,
将标志物点拟合成直线方程:r=λd+q,
其中,p=(px,py,pz)为地面所在的平面上的一点,n=(nx,ny,nz)为平面法向量,q=(qx,qy,qz)为标志物所在的直线上的一点,d=(dx,dy,dz)为标志物所在的直线的朝向,r=(rx,ry,rz)为地面所在的平面与标志物所在的直线的交点。
通过对点云数据进行ROI过滤,能够便于选择确定两个激光雷达的感兴趣区域(ROI),保证区域内有较多的地面点和标志物点,从而方便后续的数据分割和计算。
另外,若分割失败,则抛弃本次点云数据,待车辆位置改变后,再次尝试分割。
进一步地,根据所述平面方程和其对应的所述直线方程计算所述第一激光雷达观测到的所述标志物的朝向d1,所述标志物与地面的交点r1,以及所述第二激光雷达观测到的所述标志物的朝向d2,所述标志物与地面的交点r2,其中:
Figure BDA0002957655850000031
进一步地,根据所述第一激光雷达观测到的所述标志物与地面的交点r1,计算所述标志物上的预设三维点r1+kd1,根据所述第二激光雷达观测到的所述标志物与地面的交点r2,计算所述标志物上的预设三维点r2+kd2,其中k为预设参数,用于控制预设点的位置,k的具体数值可以根据实际测量需求进行调整,且为了保证标定的精准度,通过第一激光雷达和第二激光雷达计算的标志物上的预设点应为同一点,即两个预设三维点的k值相等。
进一步地,将当前标定获得的所述第一激光雷达观测到的所述标志物与地面的交点r1,以及所述标志物上的预设三维点r1+kd1,记为x1,将当前标定获得的所述第二激光雷达观测到的所述标志物与地面的交点r2,以及所述标志物上的预设三维点r2+kd2,记为y1
重复标定过程,获得若干组所述第一激光雷达和所述第二激光雷达分别观测到的所述标志物与地面的交点,以及对应的所述标志物上的三维点;
将所有所述第一激光雷达观测到的所述标志物与地面的交点,以及对应的所述标志物上的三维点组成三维点集X,将所有所述第二激光雷达观测到的所述标志物与地面的交点,以及对应的所述标志物上的三维点组成三维点集Y,其中,X={x1,x2,…,xN},Y={y1,y2,...,yN},N≥4,如果N小于4,那么本次标定失败,可以亮起红色标定指示灯并启动蜂鸣器发声3秒,以提醒工作人员。
分别计算X和Y的平均值和
Figure BDA0002957655850000044
Figure BDA0002957655850000045
其中,
Figure BDA0002957655850000041
根据
Figure BDA0002957655850000046
Figure BDA0002957655850000047
获得协方差矩阵H,其中,
Figure BDA0002957655850000042
对所述协方差矩阵进行SVD分解,获得H=UΛVT
计算联合标定的旋转矩阵R,其中,R=VUT∈O(3);
计算平移向量t,其中,
Figure BDA0002957655850000043
存储所述旋转矩阵和所述平移向量,将它们和所标定的两个激光雷达的ID号保存起来,并亮起绿色标定指示灯。
由于旋转矩阵R在正交群中,即R∈O(3),需要进一步检验R的行列式,如果det(R)=+1,则表示标定成功,R在特殊正交群中R∈SO(3),即不存在翻转变换;如果det(R)=-1,则表示本次标定失败,亮起红色标定指示灯并启动蜂鸣器发声3秒。
进一步地,所述的重复标定过程,具体包括:
移动车辆,从而调整所述第一激光雷达和所述第二激光雷达的相对位置,获得预设数量的标定数据,或在预设时间段内获得若干组标定数据,如在本实施例中,可采集400组标定数据,或在预设时间段300秒内,采集若干标定数据,在其它实施例中,可以根据实际需求调整数据数量或预设时间段。
另外,本发明还提供一种激光雷达联合标定系统,包括:
第一激光雷达;
第二激光雷达;
处理器,所述处理器与所述第一激光雷达、所述第二激光雷达通信连接;
其中,所述处理器包括:
点云采集模块,用于分别通过第一激光雷达和第二激光雷达采集棍棒状的标志物附近的点云数据;
特征分割及拟合模块,用于分别对两组所述点云数据进行地面与标志物分割,并分别将地面点拟合成平面方程,将标志物点拟合成直线方程;
交点及朝向计算模块,用于分别根据所述平面方程和其对应的所述直线方程计算所述标志物的朝向及所述标志物与地面的交点;
三维点集存储模块,分别根据所述标志物的朝向及对应的所述标志物与地面的交点计算所述标志物上预设一点的坐标,并将所述第一激光雷达观测到所述标志物与地面的交点,以及所述标志物上预设一点的坐标组合成一组三维点集,将所述第二激光雷达观测到的所述标志物与地面的交点,以及所述标志物上预设一点的坐标组合成另一组三维点集;
联合标定计算模块,用于根据两组所述三维点集计算联合标定结果。
分别通过第一激光雷达和第二激光雷达采集棍棒状的标志物附近的点云数据,并对两组点云数据进行地面与标志物分割,能够将地面点拟合成平面方程,将标志物点拟合成直线方程,通过平面方程和其对应的直线方程能够计算获得标志物的朝向及标志物与地面的交点,进而获得标志物上预设一点的坐标,将第一激光雷达观测到标志物与地面的交点,以及标志物上预设一点的坐标组合成一组三维点集,将第二激光雷达观测到的标志物与地面的交点,以及标志物上预设一点的坐标组合成另一组三维点集,便能够根据两组三维点集计算联合标定结果。由于该方案选用棍棒状的标志物,在采集数据时,对数据的要求较低,仅需要看到部分标志物,便能进行数据计算,适用于各种视场的激光雷达的标定,且数据的分割、拟合和计算都能够自动进行,使得标定的自动化程度更高,有利于提高标定精度。
进一步地,所述处理器还包括ROI过滤模块,用于对所述点云数据进行ROI过滤。
通过对点云数据进行ROI过滤,能够便于选择确定两个激光雷达的感兴趣区域(ROI),保证区域内有较多的地面点和标志物点,从而方便后续的数据分割和计算
进一步地,所述标志物竖直设置或倾斜设置。
进一步地,还包括:
标定指示灯,用于在标定成功时点亮绿灯,或在标定失败时点亮红灯;
标定蜂鸣器,用于在标定失败时进行发声。
具体的,第一激光雷达和第二激光雷达采集的数据可以通过路由器与处理器连接,在本方案中,处理器可以为工控机处理器,标定指示灯和标定蜂鸣器与路由器连接,用于在工控机处理器的控制下工作。
根据本发明提供的一种激光雷达联合标定方法及系统,分别通过第一激光雷达和第二激光雷达采集棍棒状的标志物附近的点云数据,并对两组点云数据进行地面与标志物分割,能够将地面点拟合成平面方程,将标志物点拟合成直线方程,通过平面方程和其对应的直线方程能够计算获得标志物的朝向及标志物与地面的交点,进而获得标志物上预设一点的坐标,将第一激光雷达观测到标志物与地面的交点,以及标志物上预设一点的坐标组合成一组三维点集,将第二激光雷达观测到的标志物与地面的交点,以及标志物上预设一点的坐标组合成另一组三维点集,便能够根据两组三维点集计算联合标定结果。由于该方案选用棍棒状的标志物,在采集数据时,对数据的要求较低,仅需要看到部分标志物,便能进行数据计算,适用于各种视场的激光雷达的标定,且数据的分割、拟合和计算都能够自动进行,使得标定的自动化程度更高,有利于提高标定精度。
附图说明
下面将以明确易懂的方式,结合附图说明优选实施方式,对本方案的上述特性、技术特征、优点及其实现方式予以进一步说明。
图1是本发明实施例的整体流程示意图;
图2是本发明实施例的标志物设置示意图;
图3是本发明另一个实施例的流程示意图;
图4是本发明实施例的标志物朝向示意图;
图5是本发明实施例的系统结构示意图;
图6是本发明实施例的处理器结构示意图。
图中标号:1-第一激光雷达;2-第二激光雷达;3-处理器;31-点云采集模块;32-特征分割及拟合模块;33-交点及朝向计算模块;34-三维点集存储模块;35-联合标定计算模块;36-ROI过滤模块;4-标定指示灯;5-标定蜂鸣器。
具体实施方式
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对照附图说明本发明的具体实施方式。显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图,并获得其他的实施方式。
为使图面简洁,各图中只示意性地表示出了与本发明相关的部分,它们并不代表其作为产品的实际结构。另外,以使图面简洁便于理解,在有些图中具有相同结构或功能的部件,仅示意性地绘示了其中的一个,或仅标出了其中的一个。在本文中,“一个”不仅表示“仅此一个”,也可以表示“多于一个”的情形。
实施例1
本发明的一个实施例,如图1所示,本发明提供一种激光雷达联合标定方法,包括步骤:
S1、分别通过第一激光雷达和第二激光雷达采集棍棒状的标志物附近的点云数据。
具体的,在本方案中,由于对标志物是进行直线拟合,因此,对标志物的设置没有具体要求,可以竖直设置,也可以倾斜设置。第一激光雷达、第二激光雷达以及标志物的设置方式可以如图2所示。
S2、分别对两组点云数据进行地面与标志物分割,并分别将地面点拟合成平面方程,将标志物点拟合成直线方程。
S3、分别根据平面方程和其对应的直线方程计算标志物的朝向及标志物与地面的交点。
S4、分别根据标志物的朝向及对应的标志物与地面的交点计算标志物上预设一点的坐标。
S5、将第一激光雷达观测到标志物与地面的交点,以及标志物上预设一点的坐标组合成一组三维点集,将第二激光雷达观测到的标志物与地面的交点,以及标志物上预设一点的坐标组合成另一组三维点集。
S6、根据两组三维点集计算联合标定结果。
分别通过第一激光雷达和第二激光雷达采集棍棒状的标志物附近的点云数据,并对两组点云数据进行地面与标志物分割,能够将地面点拟合成平面方程,将标志物点拟合成直线方程,通过平面方程和其对应的直线方程能够计算获得标志物的朝向及标志物与地面的交点,进而获得标志物上预设一点的坐标,将第一激光雷达观测到标志物与地面的交点,以及标志物上预设一点的坐标组合成一组三维点集,将第二激光雷达观测到的标志物与地面的交点,以及标志物上预设一点的坐标组合成另一组三维点集,便能够根据两组三维点集计算联合标定结果。由于该方案选用棍棒状的标志物,在采集数据时,对数据的要求较低,仅需要看到部分标志物,便能进行数据计算,适用于各种视场的激光雷达的标定,且数据的分割、拟合和计算都能够自动进行,使得标定的自动化程度更高,有利于提高标定精度。
另外,本方案是用于两个激光雷达的联合标定,在激光雷达的数量超过三个时,可以按照上述方法两两进行标定,得到全部的联合标定结果,从而使得标定更为灵活通用。
实施例2
本发明的一个实施例,如图3所示,在实施例1的基础上,分别对两组点云数据进行地面与标志物分割,并将地面点拟合成平面方程,将分别标志物点拟合成直线方程,具体包括:
S21、分别对第一激光雷达和第二激光雷达采集的点云数据进行ROI过滤。
通过对点云数据进行ROI过滤,能够便于选择确定两个激光雷达的感兴趣区域(ROI),保证区域内有较多的地面点和标志物点,从而方便后续的数据分割和计算。
S22、采用随机采样一致算法致(RANSAC)分割点云数据的地面点和标志物点。
S23、将地面点拟合成平面方程:nT(r-p)=0,将标志物点拟合成直线方程:r=λd+q,其中,p=(px,py,pz)为地面所在的平面上的一点,n=(nx,ny,nz)为平面法向量,q=(qx,qy,qz)为标志物所在的直线上的一点,d=(dx,dy,dz)为标志物所在的直线的朝向,标志物朝向示意如图4所示,r=(rx,ry,rz)为地面所在的平面与标志物所在的直线的交点。
另外,若分割失败,则抛弃本次点云数据,待车辆位置改变后,再次尝试分割。
优选的,根据平面方程和其对应的直线方程计算第一激光雷达观测到的标志物的朝向d1,标志物与地面的交点r1,以及第二激光雷达观测到的标志物的朝向d2,标志物与地面的交点r2,其中:
Figure BDA0002957655850000101
优选的,根据第一激光雷达观测到的标志物与地面的交点r1,计算标志物上的预设三维点r1+kd1,根据第二激光雷达观测到的标志物与地面的交点r2,计算标志物上的预设三维点r2+kd2,其中k为预设参数,用于控制预设点的位置,k的具体数值可以根据实际测量需求调整,且为了保证标定的精准度,通过第一激光雷达和第二激光雷达计算的标志物上的预设点应为同一点,即两个预设三维点的k值相等。
进一步优选的,将当前标定获得的第一激光雷达观测到的标志物与地面的交点r1,以及标志物上的预设三维点r1+kd1,记为x1,将当前标定获得的第二激光雷达观测到的标志物与地面的交点r2,以及标志物上的预设三维点r2+kd2,记为y1
重复标定过程,获得若干组第一激光雷达和第二激光雷达分别观测到的标志物与地面的交点,以及对应的标志物上的三维点。
优选的,重复标定过程,具体包括:移动车辆,从而调整所述第一激光雷达和所述第二激光雷达的相对位置,获得预设数量的标定数据,或在预设时间段内获得若干组标定数据,如在本实施例中,可采集400组标定数据,或在预设时间段300秒内,采集若干标定数据,在其它实施例中,可以根据实际需求调整数据数量或预设时间段。
具体的,将所有所述第一激光雷达观测到的所述标志物与地面的交点,以及对应的所述标志物上的三维点组成三维点集X,将所有所述第二激光雷达观测到的所述标志物与地面的交点,以及对应的所述标志物上的三维点组成三维点集Y,其中,X={x1,x2,…,xN},Y={y1,y2,…,yN},N≥4,如果N小于4,那么本次标定失败,可以亮起红色标定指示灯并启动蜂鸣器发声3秒,以提醒工作人员。
分别计算X和Y的平均值和
Figure BDA0002957655850000111
Figure BDA0002957655850000112
其中,
Figure BDA0002957655850000113
根据
Figure BDA0002957655850000116
Figure BDA0002957655850000117
获得协方差矩阵H,其中,
Figure BDA0002957655850000114
对协方差矩阵进行SVD分解,获得H=UΛVT
计算联合标定的旋转矩阵R,其中,R=VUT∈O(3)。
计算平移向量t,其中,
Figure BDA0002957655850000115
存储旋转矩阵和平移向量,将它们和所标定的两个激光雷达的ID号保存起来,并亮起绿色标定指示灯。
由于旋转矩阵R在正交群中,即R∈O(3),需要进一步检验R的行列式,如果det(R)=+1,则表示标定成功,R在特殊正交群中R∈SO(3),即不存在翻转变换;如果det(R)=-1,则表示本次标定失败,亮起红色标定指示灯并启动蜂鸣器发声3秒。
实施例3
本发明的一个实施例,如图5和图6所示,本发明还提供一种激光雷达联合标定系统,包括:第一激光雷达(即图中的激光雷达1)1、第二激光雷达(即图中的激光雷达2)2和处理器3,处理器3与第一激光雷达1、第二激光雷达2通信连接。
具体的,第一激光雷达1和第二激光雷达2采集的数据可以通过路由器与处理器3连接,在本方案中,处理器可以为工控机处理器。
优选的,该系统还包括标定指示灯4和标定蜂鸣器5。
标定指示灯4用于在标定成功时点亮绿灯,或在标定失败时点亮红灯;标定蜂鸣器5用于在标定失败时进行发声。标定指示灯4和标定蜂鸣器5与路由器连接,用于在工控机处理器的控制下工作。
其中,处理器3包括点云采集模块31、特征分割及拟合模块32、交点及朝向计算模块33、三维点集存储模块34和联合标定计算模块35。
点云采集模块31用于分别通过第一激光雷达1和第二激光雷达2采集棍棒状的标志物附近的点云数据。
具体的,在本方案中,由于对标志物是进行直线拟合,因此,对标志物的设置没有具体要求,可以竖直设置,也可以倾斜设置。第一激光雷达、第二激光雷达以及标志物的设置方式可以如图2所示。
特征分割及拟合模块32用于分别对两组点云数据进行地面与标志物分割,并分别将地面点拟合成平面方程,将标志物点拟合成直线方程。
交点及朝向计算模块33用于分别根据平面方程和其对应的直线方程计算标志物的朝向及标志物与地面的交点。
三维点集存储模块(即图中的3维点集存储模块)34分别根据标志物的朝向及对应的标志物与地面的交点计算标志物上预设一点的坐标,并将第一激光雷达观测到标志物与地面的交点,以及标志物上预设一点的坐标组合成一组三维点集,将第二激光雷达观测到的标志物与地面的交点,以及标志物上预设一点的坐标组合成另一组三维点集。
联合标定计算模块(即图中的联合标定SVD计算模块)35用于根据两组三维点集计算联合标定结果。
优选的地,处理器3还包括ROI过滤模块36,用于对点云数据进行ROI过滤。
通过对点云数据进行ROI过滤,能够便于选择确定两个激光雷达的感兴趣区域(ROI),保证区域内有较多的地面点和标志物点,从而方便后续的数据分割和计算。
分别通过第一激光雷达1和第二激光雷达2采集棍棒状的标志物附近的点云数据,并对两组点云数据进行地面与标志物分割,能够将地面点拟合成平面方程,将标志物点拟合成直线方程,通过平面方程和其对应的直线方程能够计算获得标志物的朝向及标志物与地面的交点,进而获得标志物上预设一点的坐标,将第一激光雷达1观测到标志物与地面的交点,以及标志物上预设一点的坐标组合成一组三维点集,将第二激光雷达2观测到的标志物与地面的交点,以及标志物上预设一点的坐标组合成另一组三维点集,便能够根据两组三维点集计算联合标定结果。由于该方案选用棍棒状的标志物,在采集数据时,对数据的要求较低,仅需要看到部分标志物,便能进行数据计算,适用于各种视场的激光雷达的标定,且数据的分割、拟合和计算都能够自动进行,使得标定的自动化程度更高,有利于提高标定精度。
具体的,联合标定的具体计算方式如实施例2所述,在此不累赘叙述。
另外,本方案是用于两个激光雷达的联合标定,在激光雷达的数量超过三个时,可以按照上述方法两两进行标定,得到全部的联合标定结果,从而使得标定更为灵活通用。
应当说明的是,上述实施例均可根据需要自由组合。以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。

Claims (10)

1.一种激光雷达联合标定方法,其特征在于,包括步骤:
分别通过第一激光雷达和第二激光雷达采集棍棒状的标志物附近的点云数据;
分别对两组所述点云数据进行地面与标志物分割,并分别将地面点拟合成平面方程,将标志物点拟合成直线方程;
分别根据所述平面方程和其对应的所述直线方程计算所述标志物的朝向及所述标志物与地面的交点;
分别根据所述标志物的朝向及对应的所述标志物与地面的交点计算所述标志物上预设一点的坐标;
将所述第一激光雷达观测到所述标志物与地面的交点,以及所述标志物上预设一点的坐标组合成一组三维点集,将所述第二激光雷达观测到的所述标志物与地面的交点,以及所述标志物上预设一点的坐标组合成另一组三维点集;
根据两组所述三维点集计算联合标定结果。
2.根据权利要求1所述的一种激光雷达联合标定方法,其特征在于,所述的分别对两组所述点云数据进行地面与标志物分割,并分别将地面点拟合成平面方程,将标志物点拟合成直线方程,具体包括:
分别对所述第一激光雷达和所述第二激光雷达采集的所述点云数据进行ROI过滤;
采用随机采样一致算法分割所述点云数据的地面点和标志物点;
将地面点拟合成平面方程:nT(r-p)=0,
将标志物点拟合成直线方程:r=λd+q,
其中,p=(px,py,pz)为地面所在的平面上的一点,n=(nx,ny,nz)为平面法向量,q=(qx,qy,qz)为标志物所在的直线上的一点,d=(dx,dy,dz)为标志物所在的直线的朝向,r=(rx,ry,rz)为地面所在的平面与标志物所在的直线的交点。
3.根据权利要求2所述的一种激光雷达联合标定方法,其特征在于:
根据所述平面方程和其对应的所述直线方程计算所述第一激光雷达观测到的所述标志物的朝向d1,所述标志物与地面的交点r1,以及所述第二激光雷达观测到的所述标志物的朝向d2,所述标志物与地面的交点r2
4.根据权利要求3所述的一种激光雷达联合标定方法,其特征在于:根据所述第一激光雷达观测到的所述标志物与地面的交点r1,计算所述标志物上的预设三维点r1+kd1,根据所述第二激光雷达观测到的所述标志物与地面的交点r2,计算所述标志物上的预设三维点r2+kd2,其中k为预设参数。
5.根据权利要求4所述的一种激光雷达联合标定方法,其特征在于:
将当前标定获得的所述第一激光雷达观测到的所述标志物与地面的交点r1,以及所述标志物上的预设三维点r1+kd1,记为x1,将当前标定获得的所述第二激光雷达观测到的所述标志物与地面的交点r2,以及所述标志物上的预设三维点r2+kd2,记为y1
重复标定过程,获得若干组所述第一激光雷达和所述第二激光雷达分别观测到的所述标志物与地面的交点,以及对应的所述标志物上的三维点;
将所有所述第一激光雷达观测到的所述标志物与地面的交点,以及对应的所述标志物上的三维点组成三维点集X,将所有所述第二激光雷达观测到的所述标志物与地面的交点,以及对应的所述标志物上的三维点组成三维点集Y,其中,X={x1,x2,…,xN},Y={y1,y2,…,yN},N≥4,
分别计算X和Y的平均值和
Figure FDA0002957655840000021
Figure FDA0002957655840000022
其中,
Figure FDA0002957655840000031
根据
Figure FDA0002957655840000034
Figure FDA0002957655840000035
获得协方差矩阵H,其中,
Figure FDA0002957655840000032
对所述协方差矩阵进行SVD分解,获得H=UΛVT
计算联合标定的旋转矩阵R,其中,R=VUT∈O(3);
计算平移向量t,其中,
Figure FDA0002957655840000033
存储所述旋转矩阵和所述平移向量。
6.根据权利要求5所述的一种激光雷达联合标定方法,其特征在于,所述的重复标定过程,具体包括:
移动车辆,从而调整所述第一激光雷达和所述第二激光雷达的相对位置,获得预设数量的标定数据,或在预设时间段内获得若干组标定数据。
7.一种激光雷达联合标定系统,其特征在于,包括:
第一激光雷达;
第二激光雷达;
处理器,所述处理器与所述第一激光雷达、所述第二激光雷达通信连接;
其中,所述处理器包括:
点云采集模块,用于分别通过第一激光雷达和第二激光雷达采集棍棒状的标志物附近的点云数据;
特征分割及拟合模块,用于分别对两组所述点云数据进行地面与标志物分割,并分别将地面点拟合成平面方程,将标志物点拟合成直线方程;
交点及朝向计算模块,用于分别根据所述平面方程和其对应的所述直线方程计算所述标志物的朝向及所述标志物与地面的交点;
三维点集存储模块,用于分别根据所述标志物的朝向及对应的所述标志物与地面的交点计算所述标志物上预设一点的坐标,并将所述第一激光雷达观测到所述标志物与地面的交点,以及所述标志物上预设一点的坐标组合成一组三维点集,将所述第二激光雷达观测到的所述标志物与地面的交点,以及所述标志物上预设一点的坐标组合成另一组三维点集;
联合标定计算模块,用于根据两组所述三维点集计算联合标定结果。
8.根据权利要求7所述的一种激光雷达联合标定系统,其特征在于,所述处理器还包括ROI过滤模块,用于对所述点云数据进行ROI过滤。
9.根据权利要求7所述的一种激光雷达联合标定系统,其特征在于:所述标志物竖直设置或倾斜设置。
10.根据权利要求7所述的一种激光雷达联合标定系统,其特征在于,还包括:
标定指示灯,用于在标定成功时点亮绿灯,或在标定失败时点亮红灯;
标定蜂鸣器,用于在标定失败时进行发声。
CN202110228154.6A 2021-03-02 2021-03-02 一种激光雷达联合标定方法及系统 Active CN113093155B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110228154.6A CN113093155B (zh) 2021-03-02 2021-03-02 一种激光雷达联合标定方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110228154.6A CN113093155B (zh) 2021-03-02 2021-03-02 一种激光雷达联合标定方法及系统

Publications (2)

Publication Number Publication Date
CN113093155A true CN113093155A (zh) 2021-07-09
CN113093155B CN113093155B (zh) 2022-12-23

Family

ID=76667866

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110228154.6A Active CN113093155B (zh) 2021-03-02 2021-03-02 一种激光雷达联合标定方法及系统

Country Status (1)

Country Link
CN (1) CN113093155B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113671454A (zh) * 2021-08-16 2021-11-19 中汽创智科技有限公司 一种车载雷达的位置参数标定方法、装置及存储介质

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103837869A (zh) * 2014-02-26 2014-06-04 北京工业大学 基于向量关系的单线激光雷达和ccd相机标定方法
CN105300379A (zh) * 2015-10-13 2016-02-03 上海新纪元机器人有限公司 一种基于加速度的卡尔曼滤波姿态估计方法及系统
CN206523630U (zh) * 2017-02-13 2017-09-26 杨华军 一种激光雷达标定装置
CN107229043A (zh) * 2017-05-22 2017-10-03 中国农业科学院农业资源与农业区划研究所 一种距离传感器外参数标定方法和系统
CN206740970U (zh) * 2017-05-09 2017-12-12 吉林大学 一种多线激光雷达的标定装置
US10282720B1 (en) * 2018-07-16 2019-05-07 Accel Robotics Corporation Camera-based authorization extension system
US10282852B1 (en) * 2018-07-16 2019-05-07 Accel Robotics Corporation Autonomous store tracking system
CN110031824A (zh) * 2019-04-12 2019-07-19 杭州飞步科技有限公司 激光雷达联合标定方法及装置
CN110221276A (zh) * 2019-05-31 2019-09-10 文远知行有限公司 激光雷达的标定方法、装置、计算机设备和存储介质
CN111458697A (zh) * 2020-05-27 2020-07-28 北京百度网讯科技有限公司 激光雷达标定方法、装置和系统
CN112180348A (zh) * 2020-11-27 2021-01-05 深兰人工智能(深圳)有限公司 车载多线激光雷达的姿态标定方法和装置

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103837869A (zh) * 2014-02-26 2014-06-04 北京工业大学 基于向量关系的单线激光雷达和ccd相机标定方法
CN105300379A (zh) * 2015-10-13 2016-02-03 上海新纪元机器人有限公司 一种基于加速度的卡尔曼滤波姿态估计方法及系统
CN206523630U (zh) * 2017-02-13 2017-09-26 杨华军 一种激光雷达标定装置
CN206740970U (zh) * 2017-05-09 2017-12-12 吉林大学 一种多线激光雷达的标定装置
CN107229043A (zh) * 2017-05-22 2017-10-03 中国农业科学院农业资源与农业区划研究所 一种距离传感器外参数标定方法和系统
US10282720B1 (en) * 2018-07-16 2019-05-07 Accel Robotics Corporation Camera-based authorization extension system
US10282852B1 (en) * 2018-07-16 2019-05-07 Accel Robotics Corporation Autonomous store tracking system
CN110031824A (zh) * 2019-04-12 2019-07-19 杭州飞步科技有限公司 激光雷达联合标定方法及装置
CN110221276A (zh) * 2019-05-31 2019-09-10 文远知行有限公司 激光雷达的标定方法、装置、计算机设备和存储介质
CN111458697A (zh) * 2020-05-27 2020-07-28 北京百度网讯科技有限公司 激光雷达标定方法、装置和系统
CN112180348A (zh) * 2020-11-27 2021-01-05 深兰人工智能(深圳)有限公司 车载多线激光雷达的姿态标定方法和装置

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
FUMIO ITAMI: "A Simple Calibration Procedure for a 2D LiDAR With Respect to a Camera", 《IEEE SENSORS JOURNAL》 *
敖平等: "基于线激光车体旋转扫描系统标定方法", 《新技术新工艺》 *
沈小军等: "基于地面激光雷达的输电线路铁塔倾斜度测量", 《电子测量与仪器学报》 *
陆承超: "基于并行差分进化算法的相量测量单元优化配置", 《系统仿真技术》 *
陈贵宾等: "车载三维激光雷达外参数的分步自动标定算法", 《中国激光》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113671454A (zh) * 2021-08-16 2021-11-19 中汽创智科技有限公司 一种车载雷达的位置参数标定方法、装置及存储介质
CN113671454B (zh) * 2021-08-16 2024-04-26 中汽创智科技有限公司 一种车载雷达的位置参数标定方法、装置及存储介质

Also Published As

Publication number Publication date
CN113093155B (zh) 2022-12-23

Similar Documents

Publication Publication Date Title
US8712144B2 (en) System and method for detecting crop rows in an agricultural field
US8737720B2 (en) System and method for detecting and analyzing features in an agricultural field
Barawid Jr et al. Development of an autonomous navigation system using a two-dimensional laser scanner in an orchard application
US8855405B2 (en) System and method for detecting and analyzing features in an agricultural field for vehicle guidance
Turner et al. Development of an Unmanned Aerial Vehicle (UAV) for hyper resolution vineyard mapping based on visible, multispectral, and thermal imagery
Teng et al. Development of uncut crop edge detection system based on laser rangefinder for combine harvesters
CN113093155B (zh) 一种激光雷达联合标定方法及系统
Vrochidou et al. Computer vision in self-steering tractors
CN111583342A (zh) 一种基于双目视觉的目标快速定位方法及装置
Fasiolo et al. Towards autonomous mapping in agriculture: A review of supportive technologies for ground robotics
CN107255446A (zh) 一种矮化密植果树冠层三维地图构建系统和方法
Thamrin et al. Simultaneous localization and mapping based real-time inter-row tree tracking technique for unmanned aerial vehicle
Rovira-Más et al. Design parameters for adjusting the visual field of binocular stereo cameras
Bhanu et al. Inertial navigation sensor integrated motion analysis for obstacle detection
Zhou et al. Point cloud registration for agriculture and forestry crops based on calibration balls using Kinect V2
Paturkar et al. Overview of image-based 3D vision systems for agricultural applications
Escolà et al. Mobile terrestrial laser scanner vs. UAV photogrammetry to estimate woody crop canopy parameters–Part 1: Methodology and comparison in vineyards
CN115112115A (zh) 一种机器人果园高精度实时定位方法
CN112489118B (zh) 一种无人机机载传感器组外参快速标定方法
JP2019185213A (ja) 自律移動ロボット及びその制御方法
Reiser et al. Clustering of laser scanner perception points of maize plants
Ramirez et al. 3D maps of vegetation indices generated onboard a precision agriculture UAV
Izere Plant Height Estimation Using RTK-GNSS Enabled Unmanned Aerial Vehicle (UAV) Photogrammetry
Bayat et al. The Architecture of a stereo image based system to measure tree geometric parameters
Tang et al. A Fruit-Tree Mapping System for Semi-Structured Orchards based on Multi-Sensor-Fusion SLAM

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