CN103278159A - 机载2d激光测距机获取3d点云的方法 - Google Patents

机载2d激光测距机获取3d点云的方法 Download PDF

Info

Publication number
CN103278159A
CN103278159A CN2013101957237A CN201310195723A CN103278159A CN 103278159 A CN103278159 A CN 103278159A CN 2013101957237 A CN2013101957237 A CN 2013101957237A CN 201310195723 A CN201310195723 A CN 201310195723A CN 103278159 A CN103278159 A CN 103278159A
Authority
CN
China
Prior art keywords
range finder
laser range
angle
aircraft
acquiring
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
CN2013101957237A
Other languages
English (en)
Other versions
CN103278159B (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.)
Tsinghua University
Original Assignee
Tsinghua 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 Tsinghua University filed Critical Tsinghua University
Priority to CN201310195723.7A priority Critical patent/CN103278159B/zh
Publication of CN103278159A publication Critical patent/CN103278159A/zh
Application granted granted Critical
Publication of CN103278159B publication Critical patent/CN103278159B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Length Measuring Devices By Optical Means (AREA)
  • Optical Radar Systems And Details Thereof (AREA)

Abstract

本发明提出一种机载2D激光测距机获取3D点云的方法,包括:通过EKF算法对惯性测量设备所测量的角速度进行处理,以得到横滚角和俯仰角;通过机载激光测距机获取飞行器的高度信息,并根据横滚角和俯仰角以及高度信息得到飞行器的高度坐标;根据横滚角和俯仰角对激光测距机所采集的扫描数据进行处理以得到第一扫描数据和第二扫描数据;根据第一扫描数据的数据进行数据间的匹配以得到飞行器的水平位置和偏航角;以及根据第二扫描数据、高度坐标、横滚角和俯仰角、水平位置和偏航角得到飞行器的世界坐标集,并根据世界坐标集构成飞行器所在环境的3D点云。根据本发明实施例的方法,可大幅提升无人机的自主性能。

Description

机载2D激光测距机获取3D点云的方法
技术领域
本发明涉及飞行器自主飞行技术领域,特别涉及一种机载2D激光测距机获取3D点云的方法。
背景技术
随着微电子技术和机器人技术的发展,使飞行器具备自主飞行的能力成为可能,同时在民用和军用领域对飞行器的自主性能也提出了更高的需求,特别是在勘察,救援,监控等方面的需求也日益增加。其中,对环境的感知能力是飞行器自主飞行的重要基础,如果能获取飞行器所在环境的信息并实时估计自身位置,将会大幅度提升无人机的自主性能。不同于普通的地面机器人,飞行器在飞行过程中具有6自由度,因此需要获取所在环境的3D信息,但是现有技术中大部分的传感器都难以获取有效的3D环境信息,而成为自主飞行器发展的瓶颈。
发明内容
本发明的目的旨在至少解决上述的技术缺陷之一。
为此,本发明的目的在于提出一种机载2D激光测距机获取3D点云的方法。
为达到上述目的,本发明的实施例提出一种机载2D激光测距机获取3D点云的方法,包括以下步骤:通过EKF算法对惯性测量设备所测量的角速度进行处理,以得到横滚角和俯仰角;通过机载激光测距机获取飞行器的高度信息,并根据所述横滚角和俯仰角以及所述高度信息得到飞行器的高度坐标;根据所述横滚角和俯仰角对激光测距机所采集的扫描数据进行处理以得到第一扫描数据和第二扫描数据;根据所述第一扫描数据的数据进行数据间的匹配以得到所述飞行器的水平位置和偏航角;以及根据所述第二扫描数据、所述高度坐标、所述横滚角和俯仰角、所述水平位置和偏航角得到所述飞行器的世界坐标集,并根据所述世界坐标集构成所述飞行器所在环境的3D点云。
本发明的一个实施例中,通过如下方式获得所述飞行器的世界坐标:将所述第二扫描数据转化为测量点机载坐标系下的坐标;根据所述机载坐标系下的坐标、所述高度坐标、所述横滚角和俯仰角、所述水平位置和偏航角得到所机载坐标系与世界坐标系的转换关系;以及根据所述机载坐标系与世界坐标系的转换关系将所述机载坐标系下的坐标转化为所述世界坐标。
本发明的一个实施例中,所述激光测距机的光束与所述激光测距机反射装置的平面镜成45°角。
本发明的一个实施例中,所述第一扫描数据是通过对所述第二扫描数据进行投影处理得到的。
本发明的一个实施例中,所述激光测距机为2D激光测距机。
本发明的一个实施例中,所述机载坐标系转换为世界坐标系的公式如下,所述公式为:
Figure BDA00003238024400021
其中,xi W为第i个观测点的x轴世界坐标,yi W为第i个观测点的y轴世界坐标,zi W为第i个观测点的z轴世界坐标,θ为当前激光测距机的偏航角,
Figure BDA00003238024400022
为当前激光测距机的横滚,ψ为当前激光测距机的俯仰角,xi L为第i个观测点的x轴局部坐标,yi L为第i个观测点的y轴局部坐标,x为激光测距机所在世界坐标系x坐标,y为激光测距机所在世界坐标系y坐标,z为激光测距机所在世界坐标系z坐标,(xi W,yi W,zi W)为第i个观测点的世界坐标,(xi L,yi L,0)为第i个观测点的局部坐标,(x,y,z)为激光测距机中心的世界坐标。
本发明的一个实施例中,所述第一扫描数据通过如下公式表示,所述公式为:
Figure BDA00003238024400023
Figure BDA00003238024400024
其中,为投影后的测量距离,
Figure BDA00003238024400026
为投影后的测量角度,ri为投影前的测量距离,αi为投影前的测量角度,ψ为飞行器俯仰角,
Figure BDA00003238024400027
为飞行器横滚角。
根据本发明实施例的方法具有如下优点:
1、可获取飞行器所在环境的3D点云数据,以弥补传统激光测距机只能获取2D数据的不足。
2、仅利用飞行器的惯导传感器和2D激光测距机,而不需要其余设备,因此有效地减轻了飞行器的载重负担。
3、通过激光测距机获取水平位置和偏航角,因此与使用传统的惯导传感器相比更具有鲁棒性。
4、通过激光测距的倾斜校正和噪声处理使得激光测距数据的利用更加有效。
5、可以根据需求获取感兴趣区域的3D点云信息,并将激光测距机扫描的数据累积到环境点云中,因此避免了不必要的测量以及由此带来的资源消耗。
本发明附加的方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
本发明上述的和/或附加的方面和优点从下面结合附图对实施例的描述中将变得明显和容易理解,其中:
图1为根据本发明一个实施例的机载2D激光测距机获取3D点云的方法的流程图;
图2为根据本发明一个实施例的机载2D激光测距机获取3D点云的方法的数据处理示意图;
图3为根据本发明一个实施例的2D激光测距机测量范围的示意图;以及
图4为根据本发明一个实施例的飞行器位姿示意图。
具体实施方式
下面详细描述本发明的实施例,实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,仅用于解释本发明,而不能解释为对本发明的限制。
图1为根据本发明一个实施例的机载2D激光测距机获取3D点云的方法的流程图。图2为根据本发明一个实施例的机载2D激光测距机获取3D点云的方法的数据处理示意图。如图1所示,根据本发明实施例的机载2D激光测距机获取3D点云的方法,包括以下步骤:
步骤S101,通过EKF算法对惯性测量设备所测量的角速度进行处理,以得到横滚角和俯仰角。
步骤S102,通过机载激光测距机获取飞行器的高度信息,并根据横滚角和俯仰角以及高度信息得到飞行器的高度坐标。激光测距机为2D激光测距机。
图3为根据本发明一个实施例的2D激光测距机测量范围的示意图。图4为根据本发明一个实施例的飞行器位姿示意图。使用平面镜将机载坐标系下和OL-XL,YL平面平行的激光光束中的一小部分反射成为和机载坐标系Z轴平行的方向,通过这一小部分激光光束所测量的距离得到机载坐标下飞行器的高度hL。其中,经过反射的激光光束测量得到的飞行器与障碍物的距离为l,并且测量得到的反射光线和OL-XL,YL平面平行段的长度为l0,因此测量到的飞行器的高度hL可以由下式计算得到,hL=l-l0。另外,由于飞行器存在水平方向的倾斜角度横滚角
Figure BDA00003238024400031
俯仰角ψ,因此当前得到的hL与飞行器在世界坐标系下的坐标z之间存在如下转换关系,
Figure BDA00003238024400032
步骤S103,根据横滚角和俯仰角对激光测距机所采集的扫描数据进行处理以得到第一扫描数据和第二扫描数据。第一扫描数据是通过对第二扫描数据进行投影处理得到的。
具体地,通过去噪处理可以去除扫描数据中的杂散点,保留扫描数据中的有效信息,扫描数据S0表示为{(r11),(r22),...,(rnn)},其中,(rii)为第i个扫描线的测量数据,ri表示角度为αi方向到扫描点的距离,扫描线的扫描方向是已确定的方向,扫描方向αi可以表示为αii-1+Δ,其中,Δ为固定的扫描间隔,设扫描的起始角度为α0。去噪处理的过程为:计算当前扫描(rii)的相邻的前m次扫描的距离的均值
Figure BDA00003238024400041
以及相邻的后m次扫描的距离的均值
Figure BDA00003238024400042
若|ri-rl|>rthr,且|ri-rr|>rthr,则视当前扫描数据(rii)为噪点,将其ri置为0,从而去除噪点。
由于飞行器飞行过程中会有倾斜角,即横滚角
Figure BDA00003238024400043
仰角ψ,影响水平位置(x,y)和偏航角θ的估计,因此需要使用投影处理去除横滚角
Figure BDA00003238024400044
俯仰角ψ,将倾斜的测量数据投影为与OW-XW,YW平行的数据,在横滚角俯仰角ψ,(rii)投影后的新的扫描数据为,
Figure BDA00003238024400046
Figure BDA00003238024400047
其中,为投影后的测量距离,
Figure BDA00003238024400049
为投影后的测量角度,ri为投影前的测量距离,αi为投影前的测量角度,ψ为飞行器俯仰角,
Figure BDA000032380244000410
为飞行器横滚角。
在本发明的一个实施例中,激光测距机的光束与激光测距机反射装置的平面镜成45°角。
步骤S104,根据第一扫描数据的数据进行数据间的匹配以得到飞行器的水平位置和偏航角。
具体地,ICP算法对连续两次扫描的数据,即第一扫描数据S1S2的数据进行数据间的匹配,从而计算得到相邻两次扫描间飞行器水平位置和偏航角的变化(Δx,Δy)和Δθ。使用从初始时刻到的(Δx,Δy)和Δθ进行累加即可得到当前的水平位置(x,y)和偏航角θ。
步骤S105,根据第二扫描数据、高度坐标、横滚角和俯仰角、水平位置和偏航角得到飞行器的世界坐标集,并根据世界坐标集构成飞行器所在环境的3D点云。
具体而言,将激光测距机的第二扫描数据S2转化为测量点机载坐标系下的坐标,并根据高度坐标z、横滚角
Figure BDA000032380244000411
和俯仰角ψ、水平位置(x,y)和偏航角θ得到所机载坐标系与世界坐标系的转换关系,并将测量点在机载坐标系下的坐标转化为世界坐标。将扫描数据S1中的测量点(rii)转化为机载坐标系下的坐标的换算公式为,xi L=ricosαi,yi L=risinαi,其中,xi L为第i个测量点的x轴局部坐标,yi L为第i个测量点的y轴局部坐标,ri为第i个测量点的测量距离,αi为第i个测量点的测量角度。
由机载坐标系转换为世界坐标系的公式为,
Figure BDA00003238024400051
其中,xi W为第i个观测点的x轴世界坐标,yi W为第i个观测点的y轴世界坐标,zi W为第i个观测点的z轴世界坐标,θ为当前激光测距机的偏航角,
Figure BDA00003238024400052
为当前激光测距机的横滚,ψ为当前激光测距机的俯仰角,xi L为第i个观测点的x轴局部坐标,yi L为第i个观测点的y轴局部坐标,x为激光测距机所在世界坐标系x坐标,y为激光测距机所在世界坐标系y坐标,z为激光测距机所在世界坐标系z坐标,(xi W,yi W,zi W)为第i个观测点的世界坐标,(xi L,yi L,0)为第i个观测点的局部坐标,(x,y,z)为激光测距机中心的世界坐标。
根据多个将多次扫描数据所生成的世界坐标集即可得到飞行器周围环境的3D点云数据。
根据本发明实施例的方法具有如下优点:
1、可获取飞行器所在环境的3D点云数据,以弥补传统激光测距机只能获取2D数据的不足。
2、仅利用飞行器的惯导传感器和2D激光测距机,而不需要其余设备,因此有效地减轻了飞行器的载重负担。
3、通过激光测距机获取水平位置和偏航角,因此与使用传统的惯导传感器相比更具有鲁棒性。
4、通过激光测距的倾斜校正和噪声处理使得激光测距数据的利用更加有效。
5、可以根据需求获取感兴趣区域的3D点云信息,并将激光测距机扫描的数据累积到环境点云中,因此避免了不必要的测量以及由此带来的资源消耗。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。

Claims (7)

1.一种机载2D激光测距机获取3D点云的方法,其特征在于,包括以下步骤:
通过EKF算法对惯性测量设备所测量的角速度进行处理,以得到横滚角和俯仰角;
通过机载激光测距机获取飞行器的高度信息,并根据所述横滚角和俯仰角以及所述高度信息得到飞行器的高度坐标;
根据所述横滚角和俯仰角对激光测距机所采集的扫描数据进行处理以得到第一扫描数据和第二扫描数据;
根据所述第一扫描数据的数据进行数据间的匹配以得到所述飞行器的水平位置和偏航角;以及
根据所述第二扫描数据、所述高度坐标、所述横滚角和俯仰角、所述水平位置和偏航角得到所述飞行器的世界坐标集,并根据所述世界坐标集构成所述飞行器所在环境的3D点云。
2.如权利要求1所述的机载2D激光测距机获取3D点云的方法,其特征在于,通过如下方式获得所述飞行器的世界坐标:
将所述第二扫描数据转化为测量点机载坐标系下的坐标;
根据所述机载坐标系下的坐标、所述高度坐标、所述横滚角和俯仰角、所述水平位置和偏航角得到所机载坐标系与世界坐标系的转换关系;以及
根据所述机载坐标系与世界坐标系的转换关系将所述机载坐标系下的坐标转化为所述世界坐标。
3.如权利要求1所述的机载2D激光测距机获取3D点云的方法,其特征在于,所述激光测距机的光束与所述激光测距机反射装置的平面镜成45°角。
4.如权利要求1所述的机载2D激光测距机获取3D点云的方法,其特征在于,所述第一扫描数据是通过对所述第二扫描数据进行投影处理得到的。
5.如权利要求1所述的机载2D激光测距机获取3D点云的方法,其特征在于,所述激光测距机为2D激光测距机。
6.如权利要求1所述的机载2D激光测距机获取3D点云的方法,其特征在于,所述机载坐标系转换为世界坐标系的公式如下,所述公式为:
Figure FDA00003238024300011
其中,xi W为第i个观测点的x轴世界坐标,yi W为第i个观测点的y轴世界坐标,zi W为第i个观测点的z轴世界坐标,θ为当前激光测距机的偏航角,为当前激光测距机的横滚,ψ为当前激光测距机的俯仰角,xi L为第i个观测点的x轴局部坐标,yi L为第i个观测点的y轴局部坐标,x为激光测距机所在世界坐标系x坐标,y为激光测距机所在世界坐标系y坐标,z为激光测距机所在世界坐标系z坐标,(xi W,yi W,zi W)为第i个观测点的世界坐标,(xi L,yi L,0)为第i个观测点的局部坐标,(x,y,z)为激光测距机中心的世界坐标。
7.如权利要求1或4所述的机载2D激光测距机获取3D点云的方法,其特征在于,所述第一扫描数据通过如下公式表示,所述公式为:
Figure FDA00003238024300022
其中,
Figure FDA00003238024300024
为投影后的测量距离,为投影后的测量角度,ri为投影前的测量距离,αi为投影前的测量角度,ψ为飞行器俯仰角,为飞行器横滚角。
CN201310195723.7A 2013-05-23 2013-05-23 机载2d激光测距机获取3d点云的方法 Expired - Fee Related CN103278159B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310195723.7A CN103278159B (zh) 2013-05-23 2013-05-23 机载2d激光测距机获取3d点云的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310195723.7A CN103278159B (zh) 2013-05-23 2013-05-23 机载2d激光测距机获取3d点云的方法

Publications (2)

Publication Number Publication Date
CN103278159A true CN103278159A (zh) 2013-09-04
CN103278159B CN103278159B (zh) 2016-01-20

Family

ID=49060747

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310195723.7A Expired - Fee Related CN103278159B (zh) 2013-05-23 2013-05-23 机载2d激光测距机获取3d点云的方法

Country Status (1)

Country Link
CN (1) CN103278159B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105425809A (zh) * 2015-12-02 2016-03-23 深圳市易飞行科技有限公司 一种无人机避障方法及系统
CN105549616A (zh) * 2016-01-05 2016-05-04 深圳市易飞行科技有限公司 一种基于激光阵列的多轴无人机避障系统及其避障方法
CN106969763A (zh) * 2017-04-07 2017-07-21 百度在线网络技术(北京)有限公司 用于确定无人驾驶车辆的偏航角的方法和装置
CN108415459A (zh) * 2018-05-23 2018-08-17 宜昌快马仕网络科技有限公司 一种无人机绕目标点环绕飞行的控制方法及装置
CN110866927A (zh) * 2019-11-21 2020-03-06 哈尔滨工业大学 一种基于垂足点线特征结合的ekf-slam算法的机器人定位与构图方法
CN111366155A (zh) * 2020-03-27 2020-07-03 西安应用光学研究所 基于机载光电系统的局域扫描方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101161151A (zh) * 2007-11-08 2008-04-16 浙江理工大学 基于线结构光传感器自动生成鞋底喷胶轨迹的方法及系统
US20100020306A1 (en) * 2006-07-13 2010-01-28 Velodyne Acoustics, Inc. High definition lidar system
CN102607457A (zh) * 2012-03-05 2012-07-25 西安交通大学 基于惯性导航技术的大尺寸三维形貌测量装置及方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100020306A1 (en) * 2006-07-13 2010-01-28 Velodyne Acoustics, Inc. High definition lidar system
CN101161151A (zh) * 2007-11-08 2008-04-16 浙江理工大学 基于线结构光传感器自动生成鞋底喷胶轨迹的方法及系统
CN102607457A (zh) * 2012-03-05 2012-07-25 西安交通大学 基于惯性导航技术的大尺寸三维形貌测量装置及方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
孟庆季等: "机载激光3D探测成像系统的关键技术", 《中国光学》, vol. 4, no. 3, 15 June 2011 (2011-06-15) *
王建军等: "姿态角扰动对机载激光雷达点云数据的影响", 《仪器仪表学报》, vol. 32, no. 8, 15 August 2011 (2011-08-15) *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105425809A (zh) * 2015-12-02 2016-03-23 深圳市易飞行科技有限公司 一种无人机避障方法及系统
CN105425809B (zh) * 2015-12-02 2018-01-23 深圳市易飞行科技有限公司 一种无人机避障方法及系统
CN105549616A (zh) * 2016-01-05 2016-05-04 深圳市易飞行科技有限公司 一种基于激光阵列的多轴无人机避障系统及其避障方法
CN105549616B (zh) * 2016-01-05 2018-02-16 深圳市易飞行科技有限公司 一种基于激光阵列的多轴无人机避障系统及其避障方法
CN106969763A (zh) * 2017-04-07 2017-07-21 百度在线网络技术(北京)有限公司 用于确定无人驾驶车辆的偏航角的方法和装置
CN108415459A (zh) * 2018-05-23 2018-08-17 宜昌快马仕网络科技有限公司 一种无人机绕目标点环绕飞行的控制方法及装置
CN110866927A (zh) * 2019-11-21 2020-03-06 哈尔滨工业大学 一种基于垂足点线特征结合的ekf-slam算法的机器人定位与构图方法
CN111366155A (zh) * 2020-03-27 2020-07-03 西安应用光学研究所 基于机载光电系统的局域扫描方法

Also Published As

Publication number Publication date
CN103278159B (zh) 2016-01-20

Similar Documents

Publication Publication Date Title
CN103278159B (zh) 机载2d激光测距机获取3d点云的方法
CN112347840B (zh) 视觉传感器激光雷达融合无人机定位与建图装置和方法
US10712747B2 (en) Data processing method, apparatus and terminal
CN107340522B (zh) 一种激光雷达定位的方法、装置及系统
AU2018282302B2 (en) Integrated sensor calibration in natural scenes
WO2017057041A1 (ja) 信号処理装置、信号処理方法、およびプログラム
CN103424112B (zh) 一种基于激光平面辅助的运动载体视觉导航方法
CN108692719B (zh) 物体检测装置
JP2018503913A (ja) 船舶補助ドッキング方法およびシステム
JP2012225806A (ja) 道路勾配推定装置及びプログラム
US10109074B2 (en) Method and system for inertial measurement having image processing unit for determining at least one parameter associated with at least one feature in consecutive images
JP6333396B2 (ja) モバイルプラットフォームの変位を計測する方法及び装置
WO2018056441A1 (ja) 軸ずれ推定装置
US11783507B2 (en) Camera calibration apparatus and operating method
CN113655453B (zh) 用于传感器标定的数据处理方法、装置及自动驾驶车辆
WO2017038659A1 (ja) 運動検出装置及びそれを用いた三次元形状測定装置
JP6736931B2 (ja) 立体物検出方法及び立体物検出装置
CN110533649B (zh) 一种无人机通用结构裂缝识别检测装置及方法
KR102132242B1 (ko) 스테레오 카메라 보정 시스템
CN112068152A (zh) 使用3d扫描仪同时进行2d定位和2d地图创建的方法和系统
JP2019102007A (ja) 勾配推定装置、勾配推定方法、プログラムおよび制御システム
JP5375249B2 (ja) 移動経路計画装置、移動体制御装置及び移動体
CN113340272B (zh) 一种基于无人机微群的地面目标实时定位方法
WO2021212297A1 (en) Systems and methods for distance measurement
JP7302966B2 (ja) 移動体

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160120

CF01 Termination of patent right due to non-payment of annual fee