基于LiDAR点云的电力线快速分层提取方法与流程

本发明涉及对地观测领域，具体涉及一种基于lidar点云的电力线快速分层提取方法。

背景技术：

电网是高效快捷的能源输送通道和优化配置平台，在我国的能源供应体系中发挥着举足轻重的作用，电网安全对我国的能源安全至关重要。随着电网规模的不断扩大、复杂程度的不断提高，对输电线路安全和可靠性进行保障的工作量和难度也不断增加，当前电网运行管理部门迫切需求一种高效、科学、快速、低成本的输电线路巡检方式。机载激光雷达测量系统是一种集合了激光雷达测距系统(lightdetectionandranging,lidar)能够快速获取具有较高密度且可以精确描述地面物体三维空间信息的激光点云数据。近年来，机载激光雷达测量技术由于其精度高、成本低、适用性广，可以快速获取输电线路的三维空间信息，实现输电线路安全状况的精确定量分析，在电力巡检中得到了广泛应用。现有的电力线提取算法多数使用了二维投影提取直线的方法，包括hough变换、ransac直线拟合等，未考虑电力线垂直排列方式的影响，对于多层排列的电力线提取效果不佳。

技术实现要素：

本发明的目的在于提供一种基于lidar点云的电力线快速分层提取方法，提高电力线提取效率。

为实现上述目的，本发明采用以下技术方案：

基于lidar点云的电力线快速分层提取方法，包括以下步骤：

s1、以二维格网为单位对输电线走廊点云进行局部高程分布分析，对电力线点云进行粗提取；

s2、利用ransac抛物线拟合算法对电力线点云进行分层提取；

s3、利用ransac直线拟合算法对电力线点云进行精提取，将分层提取得到的每层电力线点云数据投影到xoy平面上并使用ransac直线拟合算法提取水平排列的单根电力线点云；

s4、利用电力线直线和抛物线模型进行模型生长，获取完整的电力线点云；

s5、根据电力线的层数和每层电力线的根数获取电力线的排列方式。

优选地，步骤s1包括：

s11、裁去单档输电走廊两端各占档距5％的塔杆的点云数据；

s12、对输电走廊点云数据进行二维格网化划分；确定格网大小；计算输电走廊点云数据x和y坐标的最值；计算格网的数量及确定每个点云所在的格网编号；记录每个点云所在的格网，完成二维格网的构建；

s13、对输电走廊点云数据进行高程分布统计并确定分隔层，计算格网中最低点到最高点高程的差值，设置距离阈值，剔除差值小于距离阈值的格网；对差值大于阈值的格网，将其中所有点云数据从最低点按照间隔高度统计出每个间隔层的点云数量，将相邻没有点云数据的间隔层合并；最后，设置分隔层高度阈值，将高程最低且高度大于分隔层高度阈值的间隔层作为分隔层；

s14、根据每个格网中识别出的分隔层的高程，剔除地表点，得到粗提取的电力线点云。

优选地，步骤s2包括：

s21、将粗提取后的电力线点云投影到平行于电力线走向的竖直平面内；

s22、使用pca算法计算粗提取获得的电力线点云的主方向，根据主方向对电力线点云进行坐标变换，使得电力线走向与x轴方向平行，再将电力线点云投影到xoz平面上；

s23、对粗提取的电力线点云进行二维格网化和格网归并，计算点云数据x和z坐标的最值；计算出格网的数量以及每个点云所在的格网编号；将存在点云的格网赋值为1，不存在点云的格网赋值为0；设置格网归并阈值，逐列对z方向间距小于格网归并阈值的格网进行归并，将电力线上下格网归并阈值范围内z方向中心格网赋值为1，其余格网赋值为0；

s24、基于二维格网生成坐标数据并对其进行ransac抛物线提取；将所有赋值为1的格网的坐标编号记录，获得一组二维坐标数据，对所述坐标数据进行ransac抛物线提取；

s25、根据抛物线模型上各点的格网坐标，在原始格网中，提取出的网格上下方格网归并阈值范围内搜索赋值为1的网格，提取网格中的抛物线原始点云，得到分层之后的电力线点云。

优选地，步骤s24具体包括：

s241、在候选点中随机选择三个点作为局内点，建立抛物线模型y＝aiterationx²+biterationx+c，(iteration＝0，1，2...k)，iteration为已进行的迭代次数，k为计算最优模型要迭代的次数阈值；

s242、计算其它所有网格到该抛物线模型的距离，根据电力线线路规格，设置点到抛物线模型距离阈值，若计算出的距离满足所述点到抛物线模型距离阈值要求，则把此点也当作局内点；

s243、根据输电线路走廊点云密度，设置电力线最小点数阈值，如果局内点个数满足最小点数阈值，则采用最小二乘拟合算法评估该模型，计算中误差；

s244、重复迭代s241-243步骤k次，得到满足精度条件且局内点数最多的模型作为最优模型，将最优模型上的点从候选点中提取出来并记录；

s243、根据输电线路走廊点云密度，设置电力线最小点数阈值，如果局内点个数满足最小点数阈值，则采用最小二乘拟合算法评估该模型，计算中误差；

s244、重复迭代s241-243步骤k次，得到满足精度条件且局内点数最多的模型作为最优模型，将最优模型上的点从候选点中提取出来并记录；

s245、需重复迭代步骤示s241-s244，直至迭代次数达到次数阈值上限或检测不出符合最优模型的点云数据，迭代终止。

优选地，步骤s3包括:

s31、在所有候选点中随机选择两个点作为局内点，建立直线模型y＝aiterationx+biteration，(iteration＝0，1，2...k)，iteration为已进行的迭代次数，k为计算最优模型要迭代的次数阈值；

s32、计算中其它所有网格到该直线模型的距离，根据电力线线路规格，设置点到直线模型距离阈值，若计算出的距离满足点到直线模型距离阈值要求，则此点也作为局内点；

s33、根据输电线路走廊点云密度，设置单根电力线最小点数阈值，如果局内点个数满足最小点数阈值，则采用最小二乘拟合算法评估该模型，计算中误差；

s34、重复迭代s31-s33步骤k次，得到满足精度条件且局内点数最多的直线模型作为最优模型，将最优模型直线模型上的点从候选点中提取出来并记录；

s35、重复迭代步骤s31-s34，直至迭代次数达到次数阈值上限或检测不出符合最优模型的点云数据，迭代终止。

优选地，所述步骤s4，基于精提取电力线点云建立直线和抛物线三维模型，分别计算原始输电走廊点云数据与每条电力线直线模型和抛物线模型的距离，同时满足点到抛物线模型距离阈值和点到直线模型距离阈值要求的输电走廊点云为此电力线点云上的点；通过电力线模型生长，可以获取完整的每根电力线的点云数据，补全点云裁剪和粗提取造成的电力线点云残缺。

优选地，所述步骤s5，根据提取出的电力线层数和每层电力线的根数，获得电力线的大致排列方式，并识别避雷线及导线。

本发明基于机载激光雷达测量系统采集的单档电力走廊点云数据，提出了一种使用ransac抛物线拟合方法和ransac直线拟合方法的全自动电力线分层提取方法，实现电力线点云的自动提取。该方法参数设置简单，适用于不同的电力线排列方式，并且在提取电力线点云同时可以获取电力线的排列信息，具有很高的实用价值。

附图说明

图1为本发明技术路线流程图；

图2为本发明格网归并示意图；

图3为本发明电力线分层结果示意图；

图4为本发明电力线提取结果示意图。

图5为本发明常见架空输电线路导线排列方式。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。

实施例

参考图1所示，本发明公开了一种基于lidar点云的电力线快速分层提取方法，包括以下步骤：

s1、以二维格网为单位对输电线走廊点云进行局部高程分布分析，对电力线点云进行粗提取。

s11、为了适应不同的线路规格，裁去单档输电走廊两端各占档距5％的塔杆点云数据；

在进行电力线点云粗提取之前，为了避免杆塔点的影响，首先裁去单档电力走廊两端各占档距5％的杆塔点云数据。剩余的点云可以分为两大类：1、地表点，即地面、植被、建筑物、道路等混合点；2、噪声点，电力线点混合点云。根据电网相关安全规定，电力线与地面之间需要保留一定的安全距离，通过局部点云数据高程分析的方法可以识别此段安全距离并以此来分割地表点和噪声点。

s12、对输电走廊点云数据进行二维格网化划分。设置格网大小时，需要充分考虑输电线路的规格和点云密度的影响，设置的格网太大将会增加局部数据在高程空间上的连续性，导致过多的电力线点地表被误认为是杂点而剔除，一般格网的大小设置为2-4m。计算输电走廊点云数据x和y坐标的最值xmax，xmin，ymax，ymin，从而可以计算出格网的数量以及每个点所在的格网编号记录每个点所在的格网，完成二维格网的构建。

s13、对输电走廊点云数据进行高程分布统计并确定分隔层。以电力线弧垂最低点到地面的距离的近似值作为阈值，差值小于阈值的格网只包含地面点和植被点，不存在电力线点，予以剔除；差值大于阈值的格网，对其中所有点云数据从最低点按照0.8～1.5m的间隔高度统计每个间隔层的点云数量，若间隔过大易导致统计不完全，若间隔过小，则影响计算速度，随后将相邻没有点云数据的间隔层合并，最后，设置电网运行规范中规定的净空高度值为分隔层高度阈值，将高程最低且高度大于分隔层高度阈值的间隔层作为分隔层。为了更好的统计每个格网中分隔层，步长设置尽可能要小，一般为0.3m左右。

s14、根据每个格网中识别出的分隔层的高程，剔除地表点，得到粗提取的电力线点云。根据输电走廊点云数据的空间分布特征，利用每个格网中识别出的分隔层的高程，可以将绝大部分地表点剔除，得到粗提取的电力线点云。此时，电力线点云中还会包含一些高层植被点等噪声，主要位于电力线正下方；并且由于地表水面、深坑等造成的地面点云空洞以及到电力线点距离小于阈值的危险点、噪声点的影响，电力线会存在残缺。

s2、利用ransac抛物线拟合算法对电力线点云进行分层提取。

通过将电力线向平行于电力线走向的竖直平面内投影，处于同一层的电力线点云会在一个小范围内重合，利用此特征可以实现电力线的分层提取。

s21、将粗提取后的电力线点云投影到平行于电力线走向的竖直平面内。

s22、使用pca算法计算粗提取获得的电力线点云的主方向，根据主方向对电力线点云进行坐标变换，使得电力线走向与x轴方向平行。再将电力线点云投影到xoz平面上，如图4所示，电力线呈现出明显的分层特征。

s23、对粗提取的电力线点云进行二维格网化和格网归并。对于投影后的电力线点云进行格网化处理，计算点云数据x和z坐标的最值xmax，xmin，zmax，zmin，从而可以计算出格网的数量以及每个点所在的格网编号将存在点的格网赋值为1，不存在点的格网赋值为0。设置电力线投影高度为格网归并阈值，取值约0.5～1.5m。由于气温、大风、线路走向等原因，投影后的同一层电力线并不会完全重合，而是存在一个远远小于层间距的偏差。利用这一特点，可以逐列对z方向间距小于网格归并阈值的格网进行归并，将电力线上下网格归并阈值范围内z方向中心格网赋值为1，其余格网赋值为0。完成所有层的格网归并后，电力线呈现出明显的抛物线形态，如图2所示。

s24、基于二维格网生成坐标数据并对其进行ransac抛物线提取。将所有值为1的网格的坐标编号记录下来，获得一组二维坐标数据，对此数据进行ransac抛物线提取。

步骤如下：

s241、在所有候选点中随机选择三个点作为局内点，建立抛物线模型y＝aiterationx²+biterationx+c，(iteration＝0，1，2...k)，iteration为已进行的迭代次数，k为计算最优模型要迭代的次数阈值，为避免出现死循环，通常设置为500～1000次；

s242、计算其它所有网格到该抛物线模型的距离，根据电力线线路规格，设置点到抛物线模型距离阈值为0.5-1.5m，若计算出的距离满足点到抛物线模型距离阈值要求，则把此点也当作局内点；

s243、根据输电线路走廊点云密度，设置电力线最小点数阈值，通常为100～200个点，如果局内点个数满足阈值，则该模型就足够合理，采用最小二乘拟合算法评估该模型，计算中误差；

s244、重复迭代s241-243步骤k次，得到满足精度条件且局内点数最多的模型作为最优模型，符合该模型的所有局内点都为此模型上的点，将此模型上的点从候选点中提取出来并记录；

s245、由于电力线有多层，需重复迭代步骤s241-s244，直至迭代次数达到次数阈值上限或检测不出符合最优模型的点云数据，迭代终止。

s25、提取出的抛物线上的点均为网格坐标，此时需要使用最初生成的格网，通过在提取出的网格上下方格网归并阈值搜索值为1的网格，提取网格中的原始点云，得到分层之后的电力线点，分层结果如图4所示。

s3、利用ransac直线拟合算法对电力线点云进行精提取。

分层获得的电力线点云中，可能包含多条水平排列的电力线点云，此时，通过将每层电力线分别投影到xoy平面上，使用ransac直线拟合方法，可以快速实现单根电力线的提取。具体步骤如下：

s31、在所有候选点中随机选择两个点作为局内点，建立直线模型y＝aiterationx+biteration，(iteration＝0，1，2...k)，iteration为已进行的迭代次数，k为计算最优模型要迭代的次数阈值，通常设置为500-1000次，避免出现死循环；

s32、计算中其它所有网格到该直线模型的距离，根据电力线线路规格，设置点到直线模型距离阈值为0.5-1.5m，若计算出的距离满足点到直线模型距离阈值要求，则此点也作为局内点；

s33、根据输电线路走廊点云密度，设置单根电力线最小点数阈值，通常为50～200个点，如果局内点个数满足最小点数阈值，则采用最小二乘拟合算法评估该模型，计算中误差；

重复迭代s31-s33步骤k次，得到满足精度条件且局内点数最多的直线模型作为最优模型，符合该直线模型的所有局内点都为此模型上的点，将此直线模型上的点从候选点中提取出来并记录；由于每层可能有多条电力线，需重复迭代步骤s31-s34，直至迭代次数达到次数阈值上限或检测不出符合最优模型的点云数据，迭代终止。

s4、利用电力线直线和抛物线模型进行模型生长，获取完整的电力线点云。参考图5所示，基于上述步骤获取的单根电力线点云数据分别建立直线和抛物线数学模型，如公式：

针对原始输电走廊点云数据，分别计算其与每条电力线直线和抛物线的距离，，同时满足点到抛物线模型距离阈值和点到直线模型距离阈值要求的输电走廊点云为此电力线点云上的点。通过电力线模型生长，可以获取完整的每根电力线的点云数据，补全点云裁剪和粗提取造成的电力线点云残缺。

s5、根据电力线的层数和每层电力线的根数获取电力线的排列方式。如图5所示，根据记录的每一档提取出的层数和每层电力线的根数，可以电力线的大致排列方式。由于避雷线一般要远远细于导线，通过对电力线点云密度的判断，可以识别出点云中的避雷线和导线。

以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求的保护范围为准。