一种基于圆形扫描的机载激光雷达测距精度的检校方法技术

一种基于圆形扫描的机载激光雷达测距精度的检校方法，包括如下步骤：（1）特征点距离差值计算；（2）精度修正值计算。本发明专利技术所公开的检校方法，适用于基于圆形扫描的激光雷达测距精度检校，可有效的标定基于圆形扫描的激光雷达的测距误差，该误差值可在后续解算中修正系统误差，进而提高激光雷达的测距精度，便于激光雷达产品与高精度位置和姿态系统的整体集成，提高解算后的激光点云定位精度。此外该检校方法同样适用于对激光雷达产品测距精度的检验。

A calibration method for ranging accuracy of airborne lidar based on circular scanning

【技术实现步骤摘要】
一种基于圆形扫描的机载激光雷达测距精度的检校方法
本专利技术属于激光雷达测量系统的测距精度检校方法领域，特别涉及该领域中的一种基于圆形扫描的机载激光雷达测距精度的检校方法。
技术介绍
激光雷达的工作原理是向目标发射激光束，然后将接收到的目标回波信号与激光束发射时的时间进行解算，以获得目标的距离信息。基于圆形扫描的机载激光雷达采用单波束工作方式，激光器每次只发射一路光束，光束经楔形镜折射后按照固定的方向对外发射，并通过电机驱动楔形镜转动的方式在雷达的发射镜面形成圆形扫描。激光雷达测距原理为S＝v×t，v即光速为3×108m/s，t为激光雷达FPGA板卡发出激光束到接收到回波信号的时间间隔，影响激光雷达测距精度的主要原因是FPGA记录的时间t1和激光实际传输的时间t2不相等，从而导致距离的解算值与真值之间存在差异。引起时间t记录错误的原因主要有以下几个方面：1)FPGA发出激光发射指令到激光器发出激光束产生时间差Δt1；2)激光器发射激光束经过雷达楔形镜时速度下降，产生时间差Δt2；3)接收装置接收到激光回波信号到FPGA记录该回波时间产生时间差Δt3。另外，激光雷达产品在后期的系统集成中需要与组合导航(GNSS+IMU)协同工作，激光点云经与位姿数据联合解算后，才能输出带有大地坐标的三维空间点云。而激光雷达与GNSS、IMU联合标定时，无法测量雷达内部的光路原点的准确位置，只能对激光雷达外表面的特征点进行标定。
技术实现思路
本专利技术所要解决的技术问题就是提供一种基于圆形扫描的机载激光雷达测距精度的检校方法。本专利技术采用如下技术方案：一种基于圆形扫描的机载激光雷达测距精度的检校方法，其改进之处在于，包括如下步骤：(1)特征点距离差值计算：(11)特征点选取：机载激光雷达工作方式是圆形扫描式，即通过伺服电机带动楔形镜沿中心线旋转，使得激光束呈以中心线为轴心的圆锥形，在目标平面上体现为以中心线在目标平面上的投影为圆心的圆弧，选取圆弧上与激光雷达镜面中心在同一水平面和同一铅垂面的四个扫描点为距离差值计算的特征点，标记为A、B、C、D；(12)激光雷达定点观测特征点：将激光雷达水平安置在安置台上，设置楔形镜的旋转角度，使激光束水平打出到达A点，记录此时激光雷达数据解算软件解算的距离值SPA0，并标记A点，同理，设置楔形镜旋转角度，可得到B、C、D三点的距离解算值SPB0、SPC0、SPD0，并标记B、C、D三点，A、B、C、D四点的观测为一组，完成该组观测后标记激光雷达镜面中心点为P点；(13)全站仪观测记录：将全站仪架设在空旷无遮挡处，整平后标记全站仪中心点为O，开始测量激光雷达和特征点：(131)转动全站仪的照准部使其照准P点观测，记录OP距离SOP和竖角β1，并将全站仪水平度盘置零；(132)转动全站仪的照准部使其照准特征点A观测，记录OA距离SOA、竖角β2和水平角α1，水平角即∠POA；(133)将全站仪依次照准B、C、D三点，可依次得到B、C、D三点的观测值(SOB、β3、α2)，(SOC、β4、α3)，(SOD、β5、α4)；(14)计算激光雷达镜面中心到各特征点的距离观测值S：SPA12＝SOP2×cos2β1+SOA2×cos2β2-2×SOP×cosβ1×SOA×cosβ2×cosα1SPB12＝SOP2×cos2β1+SOB2×cos2β3-2×SOP×cosβ1×SOB×cosβ3×cosα2SPC12＝SOP2+SOC2-2×SOP×SOC×(cosβ1×cosβ4×cosα3+sinβ1×sinβ4)SPD12＝SOP2+SOD2-2×SOP×SOD×(cosβ1×cosβ5×cosα4-sinβ1×sinβ5)根据上述公式带入(13)中观测值，计算得到雷达镜面中心到四个特征点的一组距离值SPA1、SPB1、SPC1、SPD1；(15)计算距离差值ΔS：根据(14)中计算得到的一组距离值，结合(12)中记录的距离解算值，通过差值计算公式：ΔS1＝[(SPA0-SPA1)+(SPB0-SPB1)+(SPC0-SPC1)+(SPD0-SPD1)]/4计算得到第一组观测的距离差值ΔS1，重复步骤(12)—(14)，得到多组距离解算值和观测值，根据差值计算公式，得到多个距离差值ΔSi，对多个距离差值取平均作为最终距离差值ΔS；(2)精度修正值计算：(21)激光雷达定点观测：将激光雷达安置在任意位置，设置激光雷达转速为0，任意角度发射激光束，标记激光雷达镜面中心点为P，标记激光落点为E，记录激光雷达距离解算值SPE0；(22)全站仪观测记录：(221)将全站仪整平后，标记全站仪中心点为O；(222)转动全站仪的照准部使其照准P点观测，记录OP距离SOP和竖角β1，并将全站仪水平度盘置零；(223)转动全站仪的照准部使其照准E点观测，记录OE距离SOE、竖角β2和水平角α1，水平角即∠POE；(23)计算激光雷达镜面中心到激光落点的距离观测值S：根据公式SPE12＝SOP2+SOC2-2×SOP×SOC×(cosβ1×cosβ4×cosα3+sinβ1×sinβ4)利用(22)中观测值计算得到P、E两点的距离值SPE1；(24)计算距离真值：利用(23)中计算得到的距离值SPE1，结合(15)计算得到的距离差值ΔS，根据公式SPE＝SPE1+ΔS，计算得到激光雷达探测点的距离真值；(25)计算激光雷达测距误差：根据误差公式Δ1＝SPE0–SPE，计算得到激光雷达在E点的测距误差Δ1，旋转激光雷达使激光打在不同的位置，重复步骤(21)-(24)，得到多组距离解算值、距离真值，并计算得到多个测距误差Δi，对多个测距误差值取平均作为激光雷达测距误差Δ，Δ即为激光雷达测距精度的修正值和相对于镜面中心位置的距离值。本专利技术的有益效果是：本专利技术所公开的检校方法，适用于基于圆形扫描的激光雷达测距精度检校，可有效的标定基于圆形扫描的激光雷达的测距误差，该误差值可在后续解算中修正系统误差，进而提高激光雷达的测距精度，便于激光雷达产品与高精度位置和姿态系统的整体集成，提高解算后的激光点云定位精度。此外该检校方法同样适用于对激光雷达产品测距精度的检验。附图说明图1是激光雷达的测量原理图；图2是激光雷达的光路示意图；图3是全站仪的精度验证示意图；图4是距离计算原理图。具体实施方式为了使本专利技术的目的、技术方案及优点更加清楚明白，以下结合附图和实施例，对本专利技术进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本专利技术，并不用于限定本专利技术。实施例1，激光雷达发射的激光束经楔形镜折射后，沿着楔形镜固定的位置射出，本文档来自技高网...

【技术保护点】
1.一种基于圆形扫描的机载激光雷达测距精度的检校方法，其特征在于，包括如下步骤：/n(1)特征点距离差值计算：/n(11)特征点选取：/n机载激光雷达工作方式是圆形扫描式，即通过伺服电机带动楔形镜沿中心线旋转，使得激光束呈以中心线为轴心的圆锥形，在目标平面上体现为以中心线在目标平面上的投影为圆心的圆弧，选取圆弧上与激光雷达镜面中心在同一水平面和同一铅垂面的四个扫描点为距离差值计算的特征点，标记为A、B、C、D；/n(12)激光雷达定点观测特征点：/n将激光雷达水平安置在安置台上，设置楔形镜的旋转角度，使激光束水平打出到达A点，记录此时激光雷达数据解算软件解算的距离值S

【技术特征摘要】
1.一种基于圆形扫描的机载激光雷达测距精度的检校方法，其特征在于，包括如下步骤：
(1)特征点距离差值计算：
(11)特征点选取：
机载激光雷达工作方式是圆形扫描式，即通过伺服电机带动楔形镜沿中心线旋转，使得激光束呈以中心线为轴心的圆锥形，在目标平面上体现为以中心线在目标平面上的投影为圆心的圆弧，选取圆弧上与激光雷达镜面中心在同一水平面和同一铅垂面的四个扫描点为距离差值计算的特征点，标记为A、B、C、D；
(12)激光雷达定点观测特征点：
将激光雷达水平安置在安置台上，设置楔形镜的旋转角度，使激光束水平打出到达A点，记录此时激光雷达数据解算软件解算的距离值SPA0，并标记A点，同理，设置楔形镜旋转角度，可得到B、C、D三点的距离解算值SPB0、SPC0、SPD0，并标记B、C、D三点，A、B、C、D四点的观测为一组，完成该组观测后标记激光雷达镜面中心点为P点；
(13)全站仪观测记录：
将全站仪架设在空旷无遮挡处，整平后标记全站仪中心点为O，开始测量激光雷达和特征点：
(131)转动全站仪的照准部使其照准P点观测，记录OP距离SOP和竖角β1，并将全站仪水平度盘置零；
(132)转动全站仪的照准部使其照准特征点A观测，记录OA距离SOA、竖角β2和水平角α1，水平角即∠POA；
(133)将全站仪依次照准B、C、D三点，可依次得到B、C、D三点的观测值(SOB、β3、α2)，(SOC、β4、α3)，(SOD、β5、α4)；
(14)计算激光雷达镜面中心到各特征点的距离观测值S：
SPA12＝SOP2×cos2β1+SOA2×cos2β2-2×SOP×cosβ1×SOA×cosβ2×cosα1
SPB12＝SOP2×cos2β1+SOB2×cos2β3-2×SOP×cosβ1×SOB×cosβ3×cosα2
SPC12＝SOP2+SOC2-2×SOP×SOC×(cosβ1×cosβ4×cosα3+sinβ1×sinβ4)
SPD12＝SOP2+SOD2-2×SOP×SOD×(cosβ1×cosβ5×cosα4-sinβ1×sinβ5)
根...

【专利技术属性】
技术研发人员：张鹏，周亮，张迎春，郎锐，方淑燕，董建业，
申请(专利权)人：中国电波传播研究所中国电子科技集团公司第二十二研究所，
类型：发明
国别省市：山东;37