一种机载LIDAR航带重叠区冗余辨识及消除方法技术

一种机载LIDAR航带重叠区冗余辨识及消除方法，属于遥感数据处理技术领域；该方法包括：判别机载LIDAR点云中各个激光点所属航带，提取重叠区域；重叠区域中提取同名点；根据最小二乘法则计算转换参数，变换重叠区域激光点坐标，实现同名特征重合；计算距离阈值，根据激光点距离进行消冗；利用信息熵计算消冗前、后重叠区域内LIDAR点云数据的信息量，对消冗结果进行定量评价；本发明专利技术可有效去除重叠区冗余信息，降低数据处理压力、提高效率，可直接应用于地面、车载LIDAR相邻测站数据间的冗余消除；自适应得到距离阈值，判断激光点是否冗余，制定数据点的选取规则，实现冗余激光点的剔除；提出了评价消冗有效性的信息熵测度。

Redundancy identification and elimination method of airborne LIDAR overlap zone

An airborne LIDAR strip overlap redundancy identification and elimination method, which belongs to the technical field of remote sensing data processing; the method comprises the following steps: identification of airborne LIDAR point clouds of each laser point is strip extraction, overlapping area; extraction of corresponding points in the overlapping region; calculating the conversion parameters based on the least squares principle, coordinate transform overlap area laser, implementation the same feature of coincidence; calculate the distance threshold, according to the distance of the laser point to eliminate redundant computation; eliminate redundant overlap information before and after the LIDAR point cloud data in the region by using the information entropy, to eliminate redundant results quantitatively evaluation; the invention can effectively remove the overlap of redundant information, reduce pressure and improve the efficiency of data processing, can be direct application in ground vehicle LIDAR adjacent station redundancy elimination; adaptive distance threshold, to judge whether the laser point redundancy, setting the number of The selection rules of the base points are used to eliminate the redundant laser points, and the entropy measure is proposed to evaluate the effectiveness of the redundancy elimination.

【技术实现步骤摘要】
一种机载LIDAR航带重叠区冗余辨识及消除方法
本专利技术属于遥感数据处理
，具体涉及一种机载LIDAR航带重叠区冗余辨识及消除方法。
技术介绍
受航高及扫描视场角的限制，机载激光雷达(LightDetectionandRanging，LIDAR)系统每条航带扫描只能覆盖一定宽度的区域。为了实现大面积数据采集，必须飞行多条航带且相邻航带间需有10％-20％的旁向重叠以完全覆盖待测区域。对航带重叠区的二次扫描或重复采样原有激光点或采样新的激光点。其中，前者产生冗余数据，导致数据量增加，进而降低数据处理效率。为了获得有关扫描区域分布相对均匀的点云数据并提高数据处理的效率，必须消除原始点云中存在的大量冗余数据。目前，针对机载LIDAR航带重叠区域内冗余数据消除方法主要可分为：基于点云密度的方法和基于测量精度的数据消冗方法。前者本质是基于局部点云数据密度的简化抽样，可获得相对均匀的消冗数据，但消冗的随机性较大、点云去除过于盲目，缺乏测量误差精度理论的支撑。后者又可细分为两类：航迹信息辅助下的消冗和无需航迹信息的消冗。航迹信息辅助下的消冗方法理论严密、精度较高，但必须依赖航迹信息。而在实际数据处理时，航迹信息缺失的情况时有发生，因此，研究无需航迹信息的消冗方法具有很大的现实意义。但现有的无需航迹信息的消冗方法均未考虑消除原始点云中存在的系统误差，由此导致消冗结果的准确性降低。并且，现有方法均未对消冗算法的有效性进行定量评价。
技术实现思路
针对上述现有技术存在的不足，本专利技术提供一种机载LIDAR航带重叠区冗余辨识及消除方法。本专利技术的技术方案：一种机载LIDAR航带重叠区冗余辨识及消除方法，包括如下步骤：步骤1：读取原始机载LIDAR点云数据；步骤2：判别原始机载LIDAR点云中各个激光点所属航带；步骤3：消除相邻两条航带重叠区域内同名特征间的空间偏移，实现航带重叠区域各同名特征的重合：步骤3-1：计算相邻两条航带重叠区域；步骤3-2：分别提取两条航带中属于重叠区域的激光点；步骤3-3：转换参数赋初始值；步骤3-4：将重叠区域中平面坐标相同的激光点作为同名点；步骤3-5：根据最小二乘法则计算使同名点高程差的平方和最小的转换参数，所述转换参数包括：平移参数、旋转参数和缩放系数；步骤3-6：根据转换参数变换重叠区域激光点坐标，实现航带重叠区域各同名特征的重合。步骤4：对同名特征重合处理后的重叠区域内激光点进行冗余辨识及消除，得到消除冗余后的机载LIDAR点云数据：步骤4-1：比较航带重叠区域内的两次扫描激光点间的距离，判断是否存在冗余：步骤4-1-1：统计单条航带内最邻近两激光点的欧式距离频率分布得到距离阈值；步骤4-1-2：比较航带重叠区域内的两次扫描激光点间的距离，若距离小于等于距离阈值，则两次扫描激光点是对同名特征点的重复采样，得到冗余激光点，否则，不存在冗余。步骤4-2：分别计算存在冗余的激光点到自身所在的航带中心距离，若距离相等则任意删除其中一个激光点，否则，删除距离较大的激光点，保留距离较小的激光点，得到消除冗余后的机载LIDAR点云数据。步骤5：利用信息熵计算消除冗余前和消除冗余后重叠区域内LIDAR点云数据的信息量，对消冗结果进行定量评价。有益效果：一种机载LIDAR航带重叠区冗余辨识及消除方法与现有技术相比，具有如下优势：(1)可有效地去除机载LIDAR重叠区冗余信息，面对海量、高密度的LIDAR点云数据，可显著降低其数据处理压力、提高数据处理效率，同时，该成果可直接应用于地面、车载LIDAR相邻测站数据间的冗余消除；(2)通过对航带重叠区内激光点的统计分析自适应地得到距离阈值，并以此阈值为判据，判定航带重叠区的两次扫描激光点是否为同名特征的重复采样，定义重复采样激光点之一为冗余数据点的选取规则，实现冗余激光点的剔除；(3)提出了一种用于评价消冗算法有效性的信息熵测度，以评价提出的消冗算法的有效性。附图说明图1为本专利技术具体实施方式的机载LIDAR航带重叠区冗余辨识及消除方法流程图；图2为本专利技术具体实施方式中步骤4的具体流程图；图3为本专利技术具体实施方式中步骤4中的基于实验数据统计所得最邻近两点间的欧式距离频率直方图；图4为本专利技术具体实施方式中的实验数据图像，其中，(a)为航带1点云顶视图，(b)为航带2点云顶视图，灰度代表高程值。具体实施方式下面结合附图对本专利技术的一种实施方式作详细说明。本专利技术具体实施方式：如图1所示，一种机载LIDAR航带重叠区冗余辨识及消除方法，包括如下步骤：步骤1：读取原始机载LIDAR点云数据；步骤2：对原始机载LIDAR点云中各个激光点的航带归属进行判别，并确定航带数目；各激光点具有“GPSTime”属性：1)每个激光点的“GPSTime”值唯一；2)同一航带相邻两点间的“GPSTime”值是均匀且连续的；3)外业实测时，遥感平台飞完一条航带后需要飞出测区，然后再进行下一条航带，因而不同航带点云之间的GPSTime值相差较大。本实施例中，基于激光点“GPSTime”属性，对相邻的两激光点的“GPSTime”值做差，将“GPSTime”差值大于阈值的两个激光点分别作为两个航带的分界点。依据相邻航带的“GPSTime”分界点确定每条航带内点云的GPSTime的起止范围，由此可得每一个激光点的航带归属并进而确定航带数目。本实施例中，计算得到相邻两条航带：第1条航带和第2条航带的激光点集分别为P＝{pi(xi，yi，zi)}和Q＝{qj(xj，yj，zj)}，其中，i、j是激光点的索引，i＝1，2，…，n，j＝1，2，…，m，n、m是激光点数，(xi，yi，zi)为第1条航带第i个激光点坐标，(xj，yj，zj)为第2条航带第j个激光点坐标。步骤3：消除相邻两条航带重叠区域内同名特征间的空间偏移，实现航带重叠区域各同名特征的重合；本实施例中，采用最小高程差算法消除相邻两条航带同名特征间的空间偏移，实现航带重叠区域各同名特征的重合，最小高程差算法的任务在于找到相邻两条航带之间的变换7参数(3个平移参数tx、ty、tz，3个旋转参数rx、ry、rz，1个缩放系数S)，对第2条航带进行转换，使第2条航带和第1条航带同名点之间的Z坐标差的平方和最小。具体为：步骤3-1：第1条航带和第2条航带激光点在XOY平面的映射坐标集分别为Pxy＝{pxyi(xi，yi)}和Qxy＝{qxyj(xj，yj)}，分别做Pxy和Qxy的最小外接矩形Rp和Rq，提取Rp和Rq的交集得到重叠区域Ω；步骤3-2：分别提取第1条航带和第2条航带中属于重叠区域的激光点集Ps＝{pi'(xi'，yi'，zi')}和Qs＝{qj'(xj'，yj'，zj')}，其中，xi'、yi'∈Ω，xj'、yj'∈Ω；步骤3-3：转换参数赋初始值：3个旋转参数均赋值0，缩放系数S赋值1，3个平移参数可目视某一同名特征在两个航带的激光点坐标相减获得；步骤3-4：对于在Ps中搜索与qj'具有相同平面坐标的点作为qj'的同名点，如果qj'没有同名点，则在Ps中内插出与qj'平面相同坐标点即同名点的高程；步骤3-5：根据转换参数建立同名点坐标间转换关系：P＝S·R·P′+T(1)其中，P(X，Y，Z)和P'(x，y，z)为第1条航带本文档来自技高网...

【技术保护点】
一种机载LIDAR航带重叠区冗余辨识及消除方法，其特征在于，包括以下步骤：步骤1：读取原始机载LIDAR点云数据；步骤2：判别原始机载LIDAR点云中各个激光点所属航带；步骤3：消除相邻两条航带重叠区域内同名特征间的空间偏移，实现航带重叠区域各同名特征的重合；步骤4：对同名特征重合处理后的重叠区域内激光点进行冗余辨识及消除，得到消除冗后的机载LIDAR点云数据：步骤4‑1：比较航带重叠区域内的两次扫描激光点间的距离，判断是否存在冗余；步骤4‑2：对存在冗余的激光点进行消冗，得到消除冗余后的机载LIDAR点云数据；步骤5：利用信息熵计算消除冗余前和消除冗余后重叠区域中LIDAR点云数据的信息量，对消冗结果进行定量评价。

【技术特征摘要】
1.一种机载LIDAR航带重叠区冗余辨识及消除方法，其特征在于，包括以下步骤：步骤1：读取原始机载LIDAR点云数据；步骤2：判别原始机载LIDAR点云中各个激光点所属航带；步骤3：消除相邻两条航带重叠区域内同名特征间的空间偏移，实现航带重叠区域各同名特征的重合；步骤4：对同名特征重合处理后的重叠区域内激光点进行冗余辨识及消除，得到消除冗后的机载LIDAR点云数据：步骤4-1：比较航带重叠区域内的两次扫描激光点间的距离，判断是否存在冗余；步骤4-2：对存在冗余的激光点进行消冗，得到消除冗余后的机载LIDAR点云数据；步骤5：利用信息熵计算消除冗余前和消除冗余后重叠区域中LIDAR点云数据的信息量，对消冗结果进行定量评价。2.根据权利要求1所述的一种机载LIDAR航带重叠区冗余辨识及消除方法，其特征在于，所述步骤3具体包含如下步骤：步骤3-1：计算相邻两条航带重叠区域；步骤3-2：分别提取两条航带中属于重叠区域的激光点；步骤3-3：转换参数赋初始值；步骤3-4：将重叠...

【专利技术属性】
技术研发人员：王丽英，徐辛超，张正鹏，
申请(专利权)人：辽宁工程技术大学，
类型：发明
国别省市：辽宁,21