CN103278159B - 机载2d激光测距机获取3d点云的方法 - Google Patents
机载2d激光测距机获取3d点云的方法 Download PDFInfo
- Publication number
- CN103278159B CN103278159B CN201310195723.7A CN201310195723A CN103278159B CN 103278159 B CN103278159 B CN 103278159B CN 201310195723 A CN201310195723 A CN 201310195723A CN 103278159 B CN103278159 B CN 103278159B
- Authority
- CN
- China
- Prior art keywords
- range finder
- laser
- angle
- data
- airborne
- 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
Links
Abstract
本发明提出一种机载2D激光测距机获取3D点云的方法,包括:通过EKF算法对惯性测量设备所测量的角速度进行处理,以得到横滚角和俯仰角;通过机载激光测距机获取飞行器的高度信息,并根据横滚角和俯仰角以及高度信息得到飞行器的高度坐标;根据横滚角和俯仰角对激光测距机所采集的扫描数据进行处理以得到第一扫描数据和第二扫描数据;根据第一扫描数据的数据进行数据间的匹配以得到飞行器的水平位置和偏航角;以及根据第二扫描数据、高度坐标、横滚角和俯仰角、水平位置和偏航角得到飞行器的世界坐标集,并根据世界坐标集构成飞行器所在环境的3D点云。根据本发明实施例的方法,可大幅提升无人机的自主性能。
Description
技术领域
本发明涉及飞行器自主飞行技术领域,特别涉及一种机载2D激光测距机获取3D点云的方法。
背景技术
随着微电子技术和机器人技术的发展,使飞行器具备自主飞行的能力成为可能,同时在民用和军用领域对飞行器的自主性能也提出了更高的需求,特别是在勘察,救援,监控等方面的需求也日益增加。其中,对环境的感知能力是飞行器自主飞行的重要基础,如果能获取飞行器所在环境的信息并实时估计自身位置,将会大幅度提升无人机的自主性能。不同于普通的地面机器人,飞行器在飞行过程中具有6自由度,因此需要获取所在环境的3D信息,但是现有技术中大部分的传感器都难以获取有效的3D环境信息,而成为自主飞行器发展的瓶颈。
发明内容
本发明的目的旨在至少解决上述的技术缺陷之一。
为此,本发明的目的在于提出一种机载2D激光测距机获取3D点云的方法。
为达到上述目的,本发明的实施例提出一种机载2D激光测距机获取3D点云的方法,包括以下步骤:通过EKF算法对惯性测量设备所测量的角速度进行处理,以得到横滚角和俯仰角;通过机载激光测距机获取飞行器的高度信息,并根据所述横滚角和俯仰角以及所述高度信息得到飞行器的高度坐标;根据所述横滚角和俯仰角对激光测距机所采集的扫描数据进行处理以得到第一扫描数据和第二扫描数据;根据所述第一扫描数据的数据进行数据间的匹配以得到所述飞行器的水平位置和偏航角;以及根据所述第二扫描数据、所述高度坐标、所述横滚角和俯仰角、所述水平位置和偏航角得到所述飞行器的世界坐标集,并根据所述世界坐标集构成所述飞行器所在环境的3D点云。
本发明的一个实施例中,通过如下方式获得所述飞行器的世界坐标:将所述第二扫描数据转化为测量点机载坐标系下的坐标;根据所述机载坐标系下的坐标、所述高度坐标、所述横滚角和俯仰角、所述水平位置和偏航角得到所机载坐标系与世界坐标系的转换关系;以及根据所述机载坐标系与世界坐标系的转换关系将所述机载坐标系下的坐标转化为所述世界坐标。
本发明的一个实施例中,所述激光测距机的光束与所述激光测距机反射装置的平面镜成45°角。
本发明的一个实施例中,所述第一扫描数据是通过对所述第二扫描数据进行投影处理得到的。
本发明的一个实施例中,所述激光测距机为2D激光测距机。
本发明的一个实施例中,所述机载坐标系转换为世界坐标系的公式如下,所述公式为:其中,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)为激光测距机中心的世界坐标。
本发明的一个实施例中,所述第一扫描数据通过如下公式表示,所述公式为:
其中,为投影后的测量距离,为投影后的测量角度,ri为投影前的测量距离,αi为投影前的测量角度,ψ为飞行器俯仰角,为飞行器横滚角。
根据本发明实施例的方法具有如下优点:
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。另外,由于飞行器存在水平方向的倾斜角度横滚角俯仰角ψ,因此当前得到的hL与飞行器在世界坐标系下的坐标z之间存在如下转换关系,
步骤S103,根据横滚角和俯仰角对激光测距机所采集的扫描数据进行处理以得到第一扫描数据和第二扫描数据。第一扫描数据是通过对第二扫描数据进行投影处理得到的。
具体地,通过去噪处理可以去除扫描数据中的杂散点,保留扫描数据中的有效信息,扫描数据S0表示为{(r1,α1),(r2,α2),...,(rn,αn)},其中,(ri,αi)为第i个扫描线的测量数据,ri表示角度为αi方向到扫描点的距离,扫描线的扫描方向是已确定的方向,扫描方向αi可以表示为αi=αi-1+Δ,其中,Δ为固定的扫描间隔,设扫描的起始角度为α0。去噪处理的过程为:计算当前扫描(ri,αi)的相邻的前m次扫描的距离的均值以及相邻的后m次扫描的距离的均值若|ri-rl|>rthr,且|ri-rr|>rthr,则视当前扫描数据(ri,αi)为噪点,将其ri置为0,从而去除噪点。
由于飞行器飞行过程中会有倾斜角,即横滚角仰角ψ,影响水平位置(x,y)和偏航角θ的估计,因此需要使用投影处理去除横滚角俯仰角ψ,将倾斜的测量数据投影为与OW-XW,YW平行的数据,在横滚角俯仰角ψ,(ri,αi)投影后的新的扫描数据为,
其中,为投影后的测量距离,为投影后的测量角度,ri为投影前的测量距离,αi为投影前的测量角度,ψ为飞行器俯仰角,为飞行器横滚角。
在本发明的一个实施例中,激光测距机的光束与激光测距机反射装置的平面镜成45°角。
步骤S104,根据第一扫描数据的数据进行数据间的匹配以得到飞行器的水平位置和偏航角。
具体地,ICP算法对连续两次扫描的数据,即第一扫描数据S1S2的数据进行数据间的匹配,从而计算得到相邻两次扫描间飞行器水平位置和偏航角的变化(Δx,Δy)和Δθ。使用从初始时刻到的(Δx,Δy)和Δθ进行累加即可得到当前的水平位置(x,y)和偏航角θ。
步骤S105,根据第二扫描数据、高度坐标、横滚角和俯仰角、水平位置和偏航角得到飞行器的世界坐标集,并根据世界坐标集构成飞行器所在环境的3D点云。
具体而言,将激光测距机的第二扫描数据S2转化为测量点机载坐标系下的坐标,并根据高度坐标z、横滚角和俯仰角ψ、水平位置(x,y)和偏航角θ得到所机载坐标系与世界坐标系的转换关系,并将测量点在机载坐标系下的坐标转化为世界坐标。将扫描数据S1中的测量点(ri,αi)转化为机载坐标系下的坐标的换算公式为,xi L=ricosαi,yi L=risinαi,其中,xi L为第i个测量点的x轴局部坐标,yi L为第i个测量点的y轴局部坐标,ri为第i个测量点的测量距离,αi为第i个测量点的测量角度。
由机载坐标系转换为世界坐标系的公式为,
其中,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)为激光测距机中心的世界坐标。
根据多个将多次扫描数据所生成的世界坐标集即可得到飞行器周围环境的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.如权利要求2所述的机载2D激光测距机获取3D点云的方法,其特征在于,所述机载坐标系转换为世界坐标系的公式如下,所述公式为:
其中,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.如权利要求4所述的机载2D激光测距机获取3D点云的方法,其特征在于,所述第一扫描数据通过如下公式表示,所述公式为:
其中,为投影后的测量距离,为投影后的测量角度,ri为投影前的测量距离,αi为投影前的测量角度,ψ为飞行器俯仰角,为当前激光测距机的横滚角。
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 CN103278159A (zh) | 2013-09-04 |
CN103278159B true CN103278159B (zh) | 2016-01-20 |
Family
ID=49060747
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310195723.7A Active CN103278159B (zh) | 2013-05-23 | 2013-05-23 | 机载2d激光测距机获取3d点云的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103278159B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105425809B (zh) * | 2015-12-02 | 2018-01-23 | 深圳市易飞行科技有限公司 | 一种无人机避障方法及系统 |
CN105549616B (zh) * | 2016-01-05 | 2018-02-16 | 深圳市易飞行科技有限公司 | 一种基于激光阵列的多轴无人机避障系统及其避障方法 |
CN106969763B (zh) * | 2017-04-07 | 2021-01-01 | 百度在线网络技术(北京)有限公司 | 用于确定无人驾驶车辆的偏航角的方法和装置 |
CN108415459A (zh) * | 2018-05-23 | 2018-08-17 | 宜昌快马仕网络科技有限公司 | 一种无人机绕目标点环绕飞行的控制方法及装置 |
CN110866927B (zh) * | 2019-11-21 | 2021-07-20 | 哈尔滨工业大学 | 一种基于垂足点线特征结合的ekf-slam算法的机器人定位与构图方法 |
CN111366155B (zh) * | 2020-03-27 | 2022-10-14 | 西安应用光学研究所 | 基于机载光电系统的局域扫描方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101161151A (zh) * | 2007-11-08 | 2008-04-16 | 浙江理工大学 | 基于线结构光传感器自动生成鞋底喷胶轨迹的方法及系统 |
CN102607457A (zh) * | 2012-03-05 | 2012-07-25 | 西安交通大学 | 基于惯性导航技术的大尺寸三维形貌测量装置及方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101688774A (zh) * | 2006-07-13 | 2010-03-31 | 威力登音响公司 | 高精确度激光雷达系统 |
-
2013
- 2013-05-23 CN CN201310195723.7A patent/CN103278159B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101161151A (zh) * | 2007-11-08 | 2008-04-16 | 浙江理工大学 | 基于线结构光传感器自动生成鞋底喷胶轨迹的方法及系统 |
CN102607457A (zh) * | 2012-03-05 | 2012-07-25 | 西安交通大学 | 基于惯性导航技术的大尺寸三维形貌测量装置及方法 |
Non-Patent Citations (2)
Title |
---|
姿态角扰动对机载激光雷达点云数据的影响;王建军等;《仪器仪表学报》;20110815;第32卷(第8期);第1810-1817页 * |
机载激光3D探测成像系统的关键技术;孟庆季等;《中国光学》;20110615;第4卷(第3期);第327-339页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103278159A (zh) | 2013-09-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103278159B (zh) | 机载2d激光测距机获取3d点云的方法 | |
CN107655473B (zh) | 基于slam技术的航天器相对自主导航系统 | |
JP7274674B1 (ja) | 無人航空機による3次元再構成の実行 | |
WO2018205119A1 (zh) | 基于激光雷达扫描的路沿检测方法和系统 | |
CN110609570A (zh) | 一种基于无人机的自主避障巡检方法 | |
CN103424112B (zh) | 一种基于激光平面辅助的运动载体视觉导航方法 | |
CN103697855B (zh) | 一种基于海天线检测的船体水平姿态测量方法 | |
US20160363451A1 (en) | Multi-sensor merging based super-close distance autonomous navigation apparatus and method | |
CN110889808A (zh) | 一种定位的方法、装置、设备及存储介质 | |
CN109724586B (zh) | 一种融合深度图和点云的航天器相对位姿测量方法 | |
CN103994755A (zh) | 一种基于模型的空间非合作目标位姿测量方法 | |
CN103697883B (zh) | 一种基于天际线成像的飞行器水平姿态确定方法 | |
CN105023254A (zh) | 一种合成孔径雷达图像的高度重建方法 | |
CN112068152A (zh) | 使用3d扫描仪同时进行2d定位和2d地图创建的方法和系统 | |
JP2019504418A (ja) | 移動物体の位置を判定するための方法およびシステム | |
CN112154303B (zh) | 高精度地图定位方法、系统、平台及计算机可读存储介质 | |
JP5375249B2 (ja) | 移動経路計画装置、移動体制御装置及び移動体 | |
JP5451457B2 (ja) | 3次元モデル生成装置およびコンピュータプログラム | |
Wang et al. | Micro aerial vehicle navigation with visual-inertial integration aided by structured light | |
CN104075710A (zh) | 一种基于航迹预测的机动扩展目标轴向姿态实时估计方法 | |
CN113610910B (zh) | 一种移动机器人避障方法 | |
Ma et al. | Robust visual-inertial odometry with point and line features for blade inspection UAV | |
CN114689046A (zh) | 一种无人机巡检隧道的方法与系统 | |
CN114089376A (zh) | 一种基于单激光雷达的负障碍物检测方法 | |
Liu et al. | 6-DOF motion estimation using optical flow based on dual cameras |
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 |