低能见度下稀疏点云航空器快速捕获方法技术

本发明专利技术公开了一种低能见度下稀疏点云航空器快速捕获方法，步骤如下：获取机场泊位位置的激光雷达的连续点云数据，对每帧点云进行预处理；密度聚类，获取机头簇集；根据机头簇集中心点坐标与滑行线和停止线的交点相对位置变换，来判断进入激光雷达覆盖范围的航空器目标位置；针对泊位到指定位置的航空器，进行密度聚类，获取完整的航空器点云集，并进行入侵检测。本发明专利技术基于实际应用场景自适应确定点云分辨率，选择距离最小邻居数位置的点计算两点间距，统计平均值作为簇集点云分辨率。统计平均值作为簇集点云分辨率。统计平均值作为簇集点云分辨率。

【技术实现步骤摘要】
低能见度下稀疏点云航空器快速捕获方法


[0001]本专利技术属于激光雷达扫描
，具体涉及一种低能见度下稀疏点云航空器快速捕获方法。

技术介绍

[0002]激光雷达受环境因素影响小，但扫描到的点云相较于图像稀疏，雨雾等能见度低的天气下，遮挡严重，航空器自身也会相互遮挡，点云中同一目标的点云簇集离散。因此激光雷达输出的是稀疏的三维点云数据，扫描物体残缺不全，这对目标检测识别和目标尺寸测量造成很大障碍。
[0003]机场每个机位上可以停靠多种机型，飞机准确停靠位置是前轮压在对应机型的停机线与滑行线路的交点上。飞机进入泊位时需要经过3个阶段：转弯摆正(机身与滑行线平行)/减速/停止。如果飞机在停机线之前停下/滑过或者未滑到停止线，超出一定范围，廊桥无法停靠，必须用牵引车拖至正确位置，泊位引导初始时辅助人工提高泊位准确率，降低人力消耗，成熟的泊位引导取代人工引导，减少人力消耗，降低机场运营成本。
[0004]实现泊位引导的前期工作是快速获取即将泊位的航空器的完整簇集，判断场面动目标是否为即将泊位到指定位置的航空器，防止误启动泊位引导系统，将航空器与场面其他动目标区分出来，实现入侵检测。
[0005]现有技术利用激光雷达扫描到的点云，进行平面分割，背景剔除，利用点云深度学习网络，训练模型，进行目标检测和识别，雷达获取的点云稀疏，点云目标检测识别的深度学习网络，网络模型复杂，此方法适用于点云各种目标检测，但是没有针对机场这种相对特殊且具体的情况，此方法对机器要求高，训练模型耗时，检测识别耗时，准确率不高，尤其出现雨雾天气等时，准确率降低。
[0006]现有技术的缺点是激光雷达扫描获取物体外形残缺不全且点云稀疏，深度学习模型体量大，训练标定数据集耗时，对机器要求高且耗时，且准确率不高。此方法虽然具有广泛使用性，但是与实际场景结合较少，实用性不强。

技术实现思路

[0007]针对于上述现有技术的不足，本专利技术的目的在于提供一种低能见度下稀疏点云航空器快速捕获方法，以解决现有技术中机场场景下激光雷达扫描获取的物体外形残缺不全且点云稀疏，深度学习模型体量大，训练标定数据集复杂且耗时，检测准确率不高的问题。
[0008]为达到上述目的，本专利技术采用的技术方案如下：
[0009]本专利技术的一种低能见度下稀疏点云航空器快速捕获方法，步骤如下：
[0010]1)获取机场泊位位置的激光雷达的连续点云数据，对每帧点云进行预处理；
[0011]2)密度聚类，获取机头簇集；
[0012]3)根据机头簇集中心点坐标与滑行线和停止线的交点相对位置变换，来判断进入激光雷达覆盖范围的航空器目标位置；
[0013]4)针对泊位到指定位置的航空器，进行密度聚类，获取完整的航空器点云集，并进行入侵检测。
[0014]进一步地，所述步骤1)具体包括：
[0015]11)基于航空器运动的区域，剔除不属于该区域的点；
[0016]12)对获取的原始点云逐帧进行平面分割，以剔除点云中地面点集；
[0017]13)基于连续帧进行动目标检测，获取仅包含运动目标的点云。
[0018]进一步地，所述步骤13)具体为：基于八叉树算法(Octree)进行前后帧空间变换检测，航空器进入雷达扫射范围，空间变换检测获取动目标点云，当连续多帧的动目标点云中点数量为增长趋势，簇最高点坐标满足航空器运动轨迹，基于时间保存空间变换点云，得到基于时间序列的动目标点云集。
[0019]进一步地，所述步骤2)具体包括：
[0020]基于所述得到动目标点云集，初始化邻域半径为引导车簇集最远两点间的距离，最小邻居数为点云中人物簇集点的点数目，对空间变换点云进行密度聚类，获取连续两帧中簇集最高点及最高点附近最小邻居数的点高度平均值，其运动方向和所在位置满足航空器运动轨迹，则判断属于航空器，簇集中包含航空器机头，计算簇集点云的分辨率，作为下一帧的邻域半径，保存簇集点，计算簇集中心点坐标并保存；根据航空器运动轨迹和卡尔曼滤波预测下一次聚类簇集中心点坐标，划出包括所述簇集(上一帧聚类得到的簇集，包含航空器机头等的簇集)所在区域，在该区域进行下一帧聚类，保存航空器簇集，计算簇集中心点坐标并保存，计算簇集点云分辨率，基于航空器运动轨迹更新区域，依照上述的方法迭代聚类(在新区域聚类，根据聚类结果簇集计算簇集分辨率，作为下一帧的聚类参数，更新聚类区域)；当前后两帧得到簇集分辨率差小于阈值，或基于簇集中心点坐标变化判断出航空器运动方向时，则不再迭代计算簇集分辨率，以固定的邻域半径继续聚类，此时确定聚类到的簇集为航空器机头簇集。
[0021]进一步地，所述卡尔曼滤波的系统状态方程为：x
k
＝Ax
k
‑1+Bu
k
‑1+w
k
‑1；
[0022]卡尔曼滤波的观测方程：z
k
＝Hx
k
+v
k
；
[0023]其中，X
k
表示系统状态矩阵，A是状态转移矩阵，B是控制输入矩阵，W
k
‑1是过程噪声；Z
k
表示状态矩阵的观测值，即实测值，H是状态观测矩阵，V
k
是测量噪声；
[0024]卡尔曼滤波器时间更新方程：
[0025][0026]其中，和分别表示k时刻和k
‑
1时刻的后验状态估计值，也即最有估计；表示k时刻的先验状态估计值，根据上一时刻的最有估计预测k时刻的结果；u
k
是k时刻外界对系统的作用；Q是预测噪声噪声协方差矩阵；P
k
‑1和P
k
分别表示k
‑
1时刻和k时刻的后验估计协方差，即和的协方差，为k时刻的先验估计协方差，即的协方差；
[0027]卡尔曼滤波器状态更新方程：
[0028][0029][0030][0031]其中，K
k
是滤波增益矩阵，R是测量噪声协方差矩阵，表示实际观测和预测的残差，其和卡尔曼增益修正先验，得到后验，I为1的矩阵，对于单模型单测量，I＝1。
[0032]进一步地，所述步骤3)具体包括：
[0033]设滑行线方向为y轴方向，停止线方向为x轴方向，滑行线与停止线的交点设为D(x0，y0)，机头簇集中心点坐标为(x，y)，若机头集簇中心点坐标与滑行线和停止线的交点坐标D，在滑行线方向上的距离(abs(y，y0))不断减小，则表示航空器向所监控的泊位驶入；若在滑行线方向上的距离(abs(y，y0))增大或丢失满足航空器高度特征的簇集，则表示航空器停止在所监控的泊位的对面泊位；若机头簇集中心点与交点坐标D在滑行线方向上的距离基本不变，机头簇集中心点与交点坐标D在停止线方向上距离先减小后增大，则表示路过。
[0034]进一步地，所述步骤4)具体包括：
[0035]航空器转弯时开始聚类其完整机身，初始化邻域半径为本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种低能见度下稀疏点云航空器快速捕获方法，其特征在于，步骤如下：1)获取机场泊位位置的激光雷达的连续点云数据，对每帧点云进行预处理；2)密度聚类，获取机头簇集；3)根据机头簇集中心点坐标与滑行线和停止线的交点相对位置变换，来判断进入激光雷达覆盖范围的航空器目标位置；4)针对泊位到指定位置的航空器，进行密度聚类，获取完整的航空器点云集，并进行入侵检测。2.根据权利要求1所述的低能见度下稀疏点云航空器快速捕获方法，其特征在于，所述步骤1)具体包括：11)基于航空器运动的区域，剔除不属于该区域的点；12)对获取的原始点云逐帧进行平面分割，以剔除点云中地面点集；13)基于连续帧进行动目标检测，获取仅包含运动目标的点云。3.根据权利要求2所述的低能见度下稀疏点云航空器快速捕获方法，其特征在于，所述步骤13)具体为：基于八叉树算法进行前后帧空间变换检测，航空器进入雷达扫射范围，空间变换检测获取动目标点云，当连续多帧的动目标点云中点数量为增长趋势，簇最高点坐标满足航空器运动轨迹，基于时间保存空间变换点云，得到基于时间序列的动目标点云集。4.根据权利要求1所述的低能见度下稀疏点云航空器快速捕获方法，其特征在于，所述步骤2)具体包括：基于所述得到动目标点云集，初始化邻域半径为引导车簇集最远两点间的距离，最小邻居数为点云中人物簇集点的点数目，对空间变换点云进行密度聚类，获取连续两帧中簇集最高点及最高点附近最小邻居数的点高度平均值，其运动方向和所在位置满足航空器运动轨迹，则判断属于航空器，簇集中包含航空器机头，计算簇集点云的分辨率，作为下一帧的邻域半径，保存簇集点，计算簇集中心点坐标并保存；根据航空器运动轨迹和卡尔曼滤波预测下一次聚类簇集中心点坐标，划出包括所述簇集所在区域，在该区域进行下一帧聚类，保存航空器簇集，计算簇集中心点坐标并保存，计算簇集点云分辨率，基于航空器运动轨迹更新区域，依照上述的方法迭代聚类；当前后两帧得到簇集分辨率差小于阈值，或基于簇集中心点坐标变化判断出航空器运动方向时，则不再迭代计算簇集分辨率，以固定的邻域半径继续聚类，此时确定聚类到的簇集为航空器机头簇集。5.根据权利要求4所述的低能见度下稀疏点云航空器快速捕获方法，其特征在于，所述卡尔曼滤波的系统状态方程为：x
k
＝Ax
k
‑1+Bu
k
‑1+w
k
‑1；卡尔曼滤波的观测方程：z
k
＝Hx
k
+v
k
；其中，X
k
表示系统状态矩阵，A是状态转移矩阵，B是控制输入矩阵，W
k
‑1是过程噪声；Z
k
表示状态矩阵的观测值，即实测值，H是状态观测矩阵，V
k
是测量噪声；卡尔曼滤波器时间更新方程：卡尔曼滤波器时间更新方程：其中，和分别表示k时刻和k
‑
1时刻的后验状态估计值，也即最有估计；表示k时刻的先验状态估计值，根据上一时刻的最有估计预测k时刻的结果；u
k
是k时刻外界对系统的作用；Q是预测噪声噪声协方差矩阵；P
k
‑1和P
k
分别表示k
‑
1时刻和k时刻的后验估计协
方差，即和的协方差，为k时刻的先验估计协方差，即的协方差；卡尔曼滤波器状态更新方程：卡尔曼滤波器状态更新方程：卡尔曼滤波器状态更新方程：其中，K
k
...

【专利技术属性】
技术研发人员：严璐，刘云，顾昕，刘科岢，李平伟，叶鑫鑫，
申请(专利权)人：南京莱斯信息技术股份有限公司，
类型：发明
国别省市：