CN106443641A - 一种激光雷达扫描均一性测量方法 - Google Patents

一种激光雷达扫描均一性测量方法 Download PDF

Info

Publication number
CN106443641A
CN106443641A CN201610857173.4A CN201610857173A CN106443641A CN 106443641 A CN106443641 A CN 106443641A CN 201610857173 A CN201610857173 A CN 201610857173A CN 106443641 A CN106443641 A CN 106443641A
Authority
CN
China
Prior art keywords
scan line
laser
point
scanning
scan
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
CN201610857173.4A
Other languages
English (en)
Other versions
CN106443641B (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.)
INSTITUTE OF SOURCE INFORMATION CHINESE ACADEMY OF FORESTRY
Original Assignee
INSTITUTE OF SOURCE INFORMATION CHINESE ACADEMY OF FORESTRY
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 INSTITUTE OF SOURCE INFORMATION CHINESE ACADEMY OF FORESTRY filed Critical INSTITUTE OF SOURCE INFORMATION CHINESE ACADEMY OF FORESTRY
Priority to CN201610857173.4A priority Critical patent/CN106443641B/zh
Publication of CN106443641A publication Critical patent/CN106443641A/zh
Application granted granted Critical
Publication of CN106443641B publication Critical patent/CN106443641B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO 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/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/48—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
    • G01S7/497—Means for monitoring or calibrating

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Radar Systems Or Details Thereof (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

本发明公开一种激光雷达扫描均一性测量方法，包括以下步骤：(1)判断末次回波；(2)判断扫描线；(3)计算扫描线内激光回波点之间的水平距离及其统计量；(4)计算扫描线0度扫描角激光回波点的中心点；(5)计算扫描线间的水平距离及其统计量；(6)计算激光扫描均一指数；该激光雷达扫描均一性测量方法能够判别激光扫描线、判别激光扫描线中线点以及判别激光扫描均一性。

Description

一种激光雷达扫描均一性测量方法
技术领域
本发明涉及一种激光雷达扫描均一性测量方法。
背景技术
对于机载激光雷达扫描系统来说，激光雷达按照一定的频率发射激光脉冲，并接收来自地物的回波。激光雷达扫描装置按照一定的扫描频率改变激光脉冲的发射方向，扫描方向一般垂直于飞行平台的运动方向。通过飞行平台的运动和扫描装置的扫描可以得到覆盖了一定宽度的条带状激光雷达点云数据。
飞行平台的运动速度受大气湍流、平台稳定性等因素影响，出现一定的波动和非均一性，导致激光扫描线之间的间隔出现变化。另外，不同的激光扫描模式将产生不同的激光回波点间隔，对于匀速的旋转棱镜扫描模式，随着扫描角的增加，激光回波点之间间隔而变大；对于变速的震荡镜扫描模式，随着扫描角的增加，激光回波点之间间隔呈现变小趋势。
激光扫描模式一般分为单向扫描和双向扫描，单向扫描是指激光沿着单个扫描方向扫描，扫描线通常为平行线，双向扫描是指激光沿着往返扫描方向扫描，扫描线可以为之字形、椭圆形等。通过激光扫描角的变化可以判别激光扫描模式，对于单向扫描模式来说，前后扫描线之间存在扫描角突变，即前一扫描线的结束点扫描角和后一扫描线的开始点扫描角之间存在突变，在扫描线内前后扫描角表现为递增或递减规律；对于双向扫描模式来说，前后扫描线的扫描角存在重叠现象，即前一扫描线的结束点扫描角等于后一扫描线的开始点扫描角，在扫描线内前后扫描角表现为递增与递减交互变换规律。
激光雷达点云数据中激光扫描角存储为整型，单位为度，由于激光脉冲发射频率较高，在激光扫描线上，1度扫描角范围内可能存在多个激光回波点，即多个激光回波点具有相同的扫描角整型值，因此，判断激光扫描线时需要考虑扫描角取整产生的影响。
扫描线之间的距离采用扫描线中心点之间的距离，由于存在多个0度角激光回波点，将0度角激光回波点的中心点作为扫描线中心点。另外，为了减小激光损失点产生的影响，采用GPS时间判别规则计算0度角激光回波点的中心点。
为了评价激光扫描间隔的均一性，这里提出了一种激光雷达扫描均一指数，并给出了测量方法。
发明内容
本发明要解决的技术问题：(1)根据前后激光回波点的扫描角判断激光扫描线，计算激光扫描线内前后激光回波点之间的水平距离；(2)根据激光回波点的0度扫描角判断激光扫描线中心点，计算激光扫描线之间的水平距离；(3)根据扫描线内激光回波点之间的水平距离和扫描线之间的水平距离计算激光扫描均一指数，用于判别激光扫描的均一性。
为解决上述问题，本发明采用如下技术方案：
一种激光雷达扫描均一性测量方法，包括以下步骤：
(1)判断末次回波：根据每个激光脉冲的回波个数NR和回波编号RN判断末次回波，若NR等于RN，则表示该回波点为末次回波点；
(2)判断扫描线：根据前后激光回波点扫描角计算扫描角差值dSA，以及扫描角差值的累计和sum(dSA)，根据dSA和sum(dSA)判断扫描线是否结束，具体判别规则如下：若dSA＝0，表示扫描线内扫描角相同的激光回波点；若dSA＞0并且sum(dSA)＜0，表示扫描线内扫描角的单调递减方向改变，当前扫描线结束，新扫描线开始，并将扫描角差值的累计和清零；若dSA＜0并且sum(dSA)＞0，表示扫描线内扫描角的单调递增方向改变，当前扫描线结束，新扫描线开始，并将扫描角差值的累计和清零；
(3)计算扫描线内激光回波点之间的水平距离及其统计量：根据前后激光回波点位置计算回波点之间的水平距离Da，并统计水平距离Da的平均值Mean(Da)和标准差StdDev(Da)：
(4)计算扫描线0度扫描角激光回波点的中心点：对于扫描线内多个0度扫描角激光回波点，根据第1个回波点的GPS时间和最后一个回波点的GPS时间，计算GPS时间中值，依次计算每个回波点GPS时间与GPS时间中值之间差值的绝对值dT，将dT最小值对应的回波点作为0度扫描角激光回波点的中心点；
(5)计算扫描线间的水平距离及其统计量：根据前后激光扫描线中心点位置计算扫描线间的水平距离Db，并统计水平距离Db的平均值Mean(Db)和标准差StdDev(Db)：
(6)计算激光扫描均一指数：根据激光扫描线内激光回波点之间水平距离Da的统计量计算扫描线内激光扫描的均一性Ua：根据激光扫描线间的水平距离Db的统计量计算扫描线间激光扫描的均一性Ub：根据扫描线内激光扫描的均一性和扫描线间激光扫描的均一性计算激光扫描的均一性U：
基本原理：激光雷达航带数据存储为点云形式，点云数据包含了每个回波点的空间坐标、扫描角等信息，仅使用末次回波点，首先根据前后两个激光回波点扫描角之间的关系判断扫描线，计算扫描线内激光回波点之间的水平距离，并统计扫描线内水平距离的平均值和标准差；根据GPS时间判断扫描线内0度扫描角激光回波点的中心点，计算激光扫描线中心点之间的水平距离，并统计扫描线间水平距离的平均值和标准差；根据扫描线内和扫描线间的统计量计算激光扫描均一指数，用于判断激光扫描的均一性。
本发明的有益效果为：通过该技术方案，能够对激光雷达扫描数据进行处理，计算激光扫描线内激光回波点之间的水平距离统计量及激光扫描线间的水平距离统计量，进而计算激光扫描均一指数，用于评价激光扫描数据的均一性，激光扫描均一指数越接近于1，表示激光扫描数据越均一。
具体实施方式
一种激光雷达扫描均一性测量方法，包括以下步骤：
(1)判断末次回波：根据每个激光脉冲的回波个数NR和回波编号RN判断末次回波，若NR等于RN，则表示该回波点为末次回波点；
(2)判断扫描线：根据前后激光回波点扫描角计算扫描角差值dSA，以及扫描角差值的累计和sum(dSA)，根据dSA和sum(dSA)判断扫描线是否结束，具体判别规则如下：若dSA＝0，表示扫描线内扫描角相同的激光回波点；若dSA＞0并且sum(dSA)＜0，表示扫描线内扫描角的单调递减方向改变，当前扫描线结束，新扫描线开始，并将扫描角差值的累计和清零；若dSA＜0并且sum(dSA)＞0，表示扫描线内扫描角的单调递增方向改变，当前扫描线结束，新扫描线开始，并将扫描角差值的累计和清零；
(3)计算扫描线内激光回波点之间的水平距离及其统计量：根据前后激光回波点位置计算回波点之间的水平距离Da，并统计水平距离Da的平均值Mean(Da)和标准差StdDev(Da)：
(4)计算扫描线0度扫描角激光回波点的中心点：对于扫描线内多个0度扫描角激光回波点，根据第1个回波点的GPS时间和最后一个回波点的GPS时间，计算GPS时间中值，依次计算每个回波点GPS时间与GPS时间中值之间差值的绝对值dT，将dT最小值对应的回波点作为0度扫描角激光回波点的中心点；
(5)计算扫描线间的水平距离及其统计量：根据前后激光扫描线中心点位置计算扫描线间的水平距离Db，并统计水平距离Db的平均值Mean(Db)和标准差StdDev(Db)：
(6)计算激光扫描均一指数：根据激光扫描线内激光回波点之间水平距离Da的统计量计算扫描线内激光扫描的均一性Ua：根据激光扫描线间的水平距离Db的统计量计算扫描线间激光扫描的均一性Ub：根据扫描线内激光扫描的均一性和扫描线间激光扫描的均一性计算激光扫描的均一性U：
本发明的有益效果为：通过该技术方案，能够对激光雷达扫描数据进行处理，计算激光扫描线内激光回波点之间的水平距离统计量及激光扫描线间的水平距离统计量，进而计算激光扫描均一指数，用于评价激光扫描数据的均一性，激光扫描均一指数越接近于1，表示激光扫描数据越均一。
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何不经过创造性劳动想到的变化或替换，都应涵盖在本发明的保护范围之内。

Claims (1)

1.一种激光雷达扫描均一性测量方法，其特征在于：包括以下步骤：
(1)判断末次回波：根据每个激光脉冲的回波个数NR和回波编号RN判断末次回波，若NR等于RN，则表示该回波点为末次回波点；
(2)判断扫描线：根据前后激光回波点扫描角计算扫描角差值dSA，以及扫描角差值的累计和sum(dSA)，根据dSA和sum(dSA)判断扫描线是否结束，具体判别规则如下：若dSA＝0，表示扫描线内扫描角相同的激光回波点；若dSA＞0并且sum(dSA)＜0，表示扫描线内扫描角的单调递减方向改变，当前扫描线结束，新扫描线开始，并将扫描角差值的累计和清零；若dSA＜0并且sum(dSA)＞0，表示扫描线内扫描角的单调递增方向改变，当前扫描线结束，新扫描线开始，并将扫描角差值的累计和清零；
(3)计算扫描线内激光回波点之间的水平距离及其统计量：根据前后激光回波点位置计算回波点之间的水平距离Da，并统计水平距离Da的平均值Mean(Da)和标准差StdDev(Da)：
(4)计算扫描线0度扫描角激光回波点的中心点：对于扫描线内多个0度扫描角激光回波点，根据第1个回波点的GPS时间和最后一个回波点的GPS时间，计算GPS时间中值，依次计算每个回波点GPS时间与GPS时间中值之间差值的绝对值dT，将dT最小值对应的回波点作为0度扫描角激光回波点的中心点；
(5)计算扫描线间的水平距离及其统计量：根据前后激光扫描线中心点位置计算扫描线间的水平距离Db，并统计水平距离Db的平均值Mean(Db)和标准差StdDev(Db)：
(6)计算激光扫描均一指数：根据激光扫描线内激光回波点之间水平距离Da的统计量计算扫描线内激光扫描的均一性Ua：根据激光扫描线间的水平距离Db的统计量计算扫描线间激光扫描的均一性Ub：根据扫描线内激光扫描的均一性和扫描线间激光扫描的均一性计算激光扫描的均一性U：
CN201610857173.4A 2016-09-28 2016-09-28 一种激光雷达扫描均一性测量方法 Active CN106443641B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610857173.4A CN106443641B (zh) 2016-09-28 2016-09-28 一种激光雷达扫描均一性测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610857173.4A CN106443641B (zh) 2016-09-28 2016-09-28 一种激光雷达扫描均一性测量方法

Publications (2)

Publication Number Publication Date
CN106443641A true CN106443641A (zh) 2017-02-22
CN106443641B CN106443641B (zh) 2019-03-08

Family

ID=58169541

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610857173.4A Active CN106443641B (zh) 2016-09-28 2016-09-28 一种激光雷达扫描均一性测量方法

Country Status (1)

Country Link
CN (1) CN106443641B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106772318A (zh) * 2016-11-24 2017-05-31 东台银信钢结构工程有限公司 一种钢结构内部激光雷达扫描密度的测量方法
CN109540142A (zh) * 2018-11-27 2019-03-29 达闼科技（北京）有限公司 一种机器人定位导航的方法、装置、计算设备

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011196916A (ja) * 2010-03-23 2011-10-06 Mitsubishi Electric Corp 計測車両、道路地物計測システム
CN104091321A (zh) * 2014-04-14 2014-10-08 北京师范大学 适用于地面激光雷达点云分类的多层次点集特征的提取方法
CN104360339A (zh) * 2014-10-31 2015-02-18 中国林业科学研究院资源信息研究所 激光雷达航带扫描边界识别方法及装置
CN105738909A (zh) * 2016-02-02 2016-07-06 江苏大学 一种适用于低空植保的无人直升机作业边界提取方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011196916A (ja) * 2010-03-23 2011-10-06 Mitsubishi Electric Corp 計測車両、道路地物計測システム
CN104091321A (zh) * 2014-04-14 2014-10-08 北京师范大学 适用于地面激光雷达点云分类的多层次点集特征的提取方法
CN104360339A (zh) * 2014-10-31 2015-02-18 中国林业科学研究院资源信息研究所 激光雷达航带扫描边界识别方法及装置
CN105738909A (zh) * 2016-02-02 2016-07-06 江苏大学 一种适用于低空植保的无人直升机作业边界提取方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
ODED RAZ等: "Submicrosecond Scan一Angle Switching Photonic Beamformer With Flat RF Response in the C and X B ands", 《JOURNAL OF LIGHTWAVE TECHNOLOGY》 *
ZHI-WEI XIA等: "Quality Metrics of Simulated Intensity Images of Cohereht Ladar", 《2012 INTERNATIONAL CONFERENCE ON OPTOELECTRONICS AND MICROELECTRONICS》 *
刘清旺: "机载激光雷达森林参数估测方法研究", 《中国博士学位论文全文数据库 信息科技辑》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106772318A (zh) * 2016-11-24 2017-05-31 东台银信钢结构工程有限公司 一种钢结构内部激光雷达扫描密度的测量方法
CN109540142A (zh) * 2018-11-27 2019-03-29 达闼科技（北京）有限公司 一种机器人定位导航的方法、装置、计算设备
CN109540142B (zh) * 2018-11-27 2021-04-06 达闼科技（北京）有限公司 一种机器人定位导航的方法、装置、计算设备
US11638997B2 (en) 2018-11-27 2023-05-02 Cloudminds (Beijing) Technologies Co., Ltd. Positioning and navigation method for a robot, and computing device thereof

Also Published As

Publication number Publication date
CN106443641B (zh) 2019-03-08

Similar Documents

Publication Publication Date Title
CN104597430B (zh) 机载气象雷达低空风切变三维回波仿真方法
CN107153193A (zh) 多线激光雷达和多线激光雷达控制方法
CN108182660B (zh) 一种区域气象雷达网数据融合方法及装置
EP3309513A1 (en) Three-dimensional topographic mapping system and mapping method
CN109032157A (zh) 无人机仿地作业方法、装置、设备及存储介质
JP6037762B2 (ja) 風計測装置
CN111624623B (zh) 基于激光雷达非均匀扫描的风场反演方法
CN106019286B (zh) 一种机载气象雷达多扫描多重频全空域气象目标探测方法
CN106066154A (zh) 一种适用于快速扫描场景的靶标及其控制点的提取方法
KR102555054B1 (ko) 멀티포인트 레이더 유속계로부터 획득된 정보와 다양한 기상 관측 정보를 빅데이터화 및 데이터 마이닝하여 하천 유량을 측정하는 브레이니 rtu 및 상기 멀티포인트 레이더 유속계와 상기 브레이니 rtu를 포함하는 하천 유량 측정 시스템
CN106443641A (zh) 一种激光雷达扫描均一性测量方法
CN106291584B (zh) 一种基于激光雷达扫描特征的点云抽稀方法
US11513197B2 (en) Multiple-pulses-in-air laser scanning system with ambiguity resolution based on range probing and 3D point analysis
CN106772318A (zh) 一种钢结构内部激光雷达扫描密度的测量方法
CN101923158B (zh) 用于高斯分解气象雷达数据以供传送的系统和方法
CN104880699B (zh) 一种合成孔径雷达的脉冲重复频率下限值的确定方法
US20220075038A1 (en) Apparatus and methods for long range, high resolution lidar
JP2009002921A (ja) レーダ装置
EP3517998B1 (en) Airborne lidar pulse rate modulation
CN102621557B (zh) 一种动态目标二维捕获的方法
CN112731355A (zh) 计算激光雷达安装角度偏差的方法、装置、终端和介质
CN113376625A (zh) 目标物体的偏离角度获得方法、装置、电子设备及存储介质
CN111190178A (zh) 一种相频扫三坐标雷达的多波束顺序比幅测高方法
CN106405527A (zh) 一种可自适应补偿被测地形高程变化的机载LiDAR装置
CN116381731A (zh) 一种基于相干多普勒激光雷达的气溶胶参数反演方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant