CN115508806A - 基于激光雷达的输电线覆冰在线监测方法、介质及设备 - Google Patents

基于激光雷达的输电线覆冰在线监测方法、介质及设备 Download PDF

Info

Publication number
CN115508806A
CN115508806A CN202211373417.3A CN202211373417A CN115508806A CN 115508806 A CN115508806 A CN 115508806A CN 202211373417 A CN202211373417 A CN 202211373417A CN 115508806 A CN115508806 A CN 115508806A
Authority
CN
China
Prior art keywords
line type
laser
icing
point cloud
laser point
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
CN202211373417.3A
Other languages
English (en)
Other versions
CN115508806B (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.)
Guangdong Anheng Power Technology Co ltd
Original Assignee
Guangdong Anheng Power Technology Co ltd
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 Guangdong Anheng Power Technology Co ltd filed Critical Guangdong Anheng Power Technology Co ltd
Priority to CN202211373417.3A priority Critical patent/CN115508806B/zh
Publication of CN115508806A publication Critical patent/CN115508806A/zh
Application granted granted Critical
Publication of CN115508806B publication Critical patent/CN115508806B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/4802—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00—Measuring arrangements characterised by the use of optical techniques
    • G01B11/02—Measuring arrangements characterised by the use of optical techniques for measuring length, width or thickness
    • G01B11/06—Measuring arrangements characterised by the use of optical techniques for measuring length, width or thickness for measuring thickness ; e.g. of sheet material
    • G01B11/0616—Measuring arrangements characterised by the use of optical techniques for measuring length, width or thickness for measuring thickness ; e.g. of sheet material of coating
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00—Measuring arrangements characterised by the use of optical techniques
    • G01B11/24—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01G—WEIGHING
    • G01G19/00—Weighing apparatus or methods adapted for special purposes not provided for in the preceding groups
    • 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
    • G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
    • G01S17/88—Lidar systems specially adapted for specific applications
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Remote Sensing (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Algebra (AREA)
  • Pure & Applied Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Electromagnetism (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明公开了一种基于激光雷达的输电线覆冰在线监测方法，包括：获取激光雷达对输电线进行扫描得到的激光点云数据，所述激光雷达架设在输电铁塔上；对所述激光点云数据进行分割，得到至少一种线路类型对应的激光点云；对于每一种线路类型，根据所述线路类型对应的激光点云计算所述线路类型在每一数据段中的每一片层对应的外轮廓多边形；根据所述线路类型在每一数据段中的每一片层对应的外轮廓多边形，计算所述线路类型对应的覆冰重量和覆冰厚度，从而有效地提高了输电线在覆冰情况下的覆冰重量和覆冰厚度的计算精度，提升了输电线覆冰在线监测的效率。

Description

基于激光雷达的输电线覆冰在线监测方法、介质及设备
技术领域
本发明涉及机器人技术领域，尤其涉及一种基于激光雷达的输电线覆冰在线监测方法、介质及设备。
背景技术
随着全球气候变化加剧，极端天气现象频发，输配电线路结冰的发生概率和极端性进一步提升。在特定地区，输配电线路结冰是造成电网停电的主要原因之一。虽然这种结冰事件依然属于低频事件，但是代价却十分高昂。因此，需要对线路覆冰进行监测并准确评估其状态，以采取相应手段进行治理，防止线路跳闸甚至倒塌的事故发生。
现有技术假定导体周围的覆冰均沿线路分布，且成规则的圆柱体形状，其几何中心与导体一致，在此假定条件下依靠拉力传感器、风速传感器、温度传感器采集的数据，结合理论推导的解析模型计算得出线路在覆冰情况下的覆冰重量和覆冰厚度。然而，实际覆冰情况受多种情况的影响，除了线路的应力值和图片成像识别的综合方式外，在形成覆冰的气象条件下覆冰监测摄像头会出现起雾、结冰现象，从而导致监控图像模糊、甚至完全是一片空白；而拉力传感器本身故障率较高，现场使用不便，通过拉力传感器测量再经过推导的覆冰厚度及覆冰重量精度不高。
发明内容
本发明实施例提供了一种基于激光雷达的输电线覆冰在线监测方法、装置、介质及设备，以解决现有输电线在覆冰情况下的覆冰重量和覆冰厚度的计算精度不高的问题。
一种基于激光雷达的输电线覆冰在线监测方法，所述方法包括：
获取激光雷达对输电线进行扫描得到的激光点云数据，所述激光雷达架设在输电铁塔上；
对所述激光点云数据进行分割，得到至少一种线路类型对应的激光点云；
对于每一种线路类型，根据所述线路类型对应的激光点云计算所述线路类型在每一数据段中的每一片层对应的外轮廓多边形；
根据所述线路类型在每一数据段中的每一片层对应的外轮廓多边形，计算所述线路类型对应的覆冰重量和覆冰厚度。
可选地，所述对所述激光点云数据进行分割，得到至少一种线路类型对应的激光点云包括：
从所述激光点云数据中获取初始激光点，所述初始激光点为与坐标原点的相对距离最小的激光点；
从所述初始激光点开始，计算所述激光点云数据中的两两激光点之间的相对距离，构建距离比邻矩阵；
构建原始激光点云，所述原始激光点云中包括所述激光点云数据中除所述初始激光点外其他的激光点；
以所述初始激光点为搜索点，按预设方式搜索所述距离比邻矩阵，获取与所述搜索点的相对距离满足预设条件的激光点构成链接矩阵，其中，一个链接矩阵表示一种线路类型对应的激光点云。
可选地，所述按预设方式搜索所述距离比邻矩阵，获取与所述搜索点的相对距离满足预设条件的激光点构成链接矩阵包括：
构建链接矩阵，将所述搜索点添加到所述链接矩阵；
在所述距离比邻矩阵中进行搜索，获取所述搜索点对应的最小的n个相对距离，n>0；
若所述n个相对距离均小于预设的线路直径阈值，则获取所述n个相对距离对应的激光点，将所获得的激光点从原始激光点云中删除，并添加到所述链接矩阵；
以所述n个相对距离中的最小值对应的激光点作为新的搜索点，继续在所述距离比邻矩阵中进行搜索；
若所述n个相对距离中存在至少一个相对距离大于或等于预设的线路直径阈值，则获取小于所述线路直径阈值的相对距离对应的激光点，将所获得的激光点从原始激光点云中删除，并添加到所述链接矩阵，以所述链接矩阵为一种线路类型对应的激光点云；
构建下一个链接矩阵，在大于或等于预设的线路直径阈值的相对距离中获取最小值，以所述最小值对应的激光点作为新的搜索点，添加到所述下一个链接矩阵，在所述距离比邻矩阵中进行下一轮搜索，以生成下一种线路类型对应的激光点云；
直至所述原始激光点云中的所述激光点纳入不同的链接矩阵中。
可选地，所述获取所述搜索点对应的最小的n个相对距离为：
获取所述搜索点对应的最小的6个相对距离。
可选地，所述根据所述线路类型对应的激光点云计算所述线路类型在每一数据段中的每一片层对应的外轮廓多边形包括：
获取一种线路类型对应的激光点云，按预设分割距离将所述激光点云沿X轴方向 分割为若干个数据段，根据每一个数据段中激光点云在Z轴方向的坐标最大值和最小值进 行等间距切割，得到
Figure 697109DEST_PATH_IMAGE001
个片层及对应的点云片层数据，
Figure 83091DEST_PATH_IMAGE001
为大于或等于1的正整数，每一点 云片层数据中包括若干个激光点，每一激光点的Z轴坐标均落在所述点云片层对应的Z轴坐 标间隔内；
遍历每一个点云片层数据，采用步进法对所述点云片层数据中的激光点进行排序，生成所述线路类型在每一个数据段中的每一片层对应的外轮廓多边形。
可选地，所述根据所述线路类型在每一数据段中的每一片层对应的外轮廓多边形，计算所述线路类型对应的覆冰重量包括：
根据点云片层数据计算所述线路类型在每一数据段中的每一个片层对应的外轮廓多边形的面积；
对每一个点云片层数据的高度和外轮廓多边形的面积之积进行求和处理，得到所述线路类型中每一数据段对应的体积；
获取气象特征信息，所述气象特征信息包括温度信息、湿度信息、风速信息、时间信息；
根据所述气象特征信息对冰密度进行修正；
根据修正后的冰密度、所述线路类型在每一数据段对应的体积计算所述线路类型在每一数据段的覆冰重量。
可选地，所述根据所述线路类型在每一数据段中的每一片层对应的外轮廓多边形，计算所述线路类型对应的覆冰厚度包括：
获取所述线路类型中每一个数据段对应的导体中心点；
对于所述线路类型每一数据段中的每一片层，计算所述片层对应的外轮廓多边形上的激光点与所述导体中心点的相对距离，对所述相对距离进行加权求和，得到所述线路类型在每一数据段中的片层对应的覆冰厚度；
计算所述线路类型在所有数据段的所有片层对应的覆冰厚度的平均值，得到所述线路类型对应的覆冰平均厚度。
可选地，所述方法还包括：
根据所述线路类型对应的导线密度和截面积计算所述线路类型对应的导线自重比载；
根据所述线路类型的覆冰重量获取单位距离覆冰重量，根据所述单位距离覆冰重量计算冰重比载；
根据所述导线自重比载和冰重比载，计算所述线路类型对应的冰重总比载；
获取风载体型系数，根据所述风载体型系数和所述线路类型对应的覆冰厚度计算覆冰风压比载；
根据所述线路类型对应的冰重总比载和覆冰风压比载计算所述线路类型在覆冰时的综合比载。
一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现如上所述的基于激光雷达的输电线覆冰在线监测方法。
一种计算机设备，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现如上所述的基于激光雷达的输电线覆冰在线监测方法。
本发明实施例通过获取激光雷达对输电线进行扫描得到的激光点云数据；然后对所述激光点云数据进行分割，得到至少一种线路类型对应的激光点云；对于每一种线路类型，对不规则覆冰切分片层，根据所述线路类型对应的激光点云计算所述线路类型在每一数据段中的每一片层对应的外轮廓多边形，实现轮廓边界确定；最后根据所述线路类型在每一数据段中的每一片层对应的外轮廓多边形，计算所述线路类型对应的覆冰重量和覆冰厚度，从而有效地提高了输电线在覆冰情况下的覆冰重量和覆冰厚度的计算精度，提升了输电线覆冰在线监测的效率。
附图说明
为了更清楚地说明本发明实施例的技术方案，下面将对本发明实施例的描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1是本发明一实施例提供的基于激光雷达的输电线覆冰在线监测方法的实现流程图；
图2是本发明一实施例提供的基于激光雷达的输电线覆冰在线监测方法中步骤S102的实现流程图；
图3是本发明一实施例提供的基于激光雷达的输电线覆冰在线监测方法中步骤S204的实现流程图；
图4是本发明一实施例提供的基于激光雷达的输电线覆冰在线监测装置的示意图；
图5是本发明一实施例中计算机设备的一示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明实施例提供了一种基于激光雷达的输电线覆冰在线监测方法，采用点云切片法，对不规则覆冰切分片层，根据所述线路类型对应的激光点云计算所述线路类型在每一数据段中的每一片层对应的外轮廓多边形，实现对每一片层的轮廓边界确定，最后根据所述线路类型在每一数据段中的每一片层对应的外轮廓多边形，计算所述线路类型对应的覆冰重量和覆冰厚度，从而有效地提高了输电线路在覆冰情况下的覆冰重量和覆冰厚度的计算精度，提升了输电线覆冰在线监测的效率。
以下对本发明实施例提供的基于激光雷达的输电线覆冰在线监测方法进行详细的描述。图1为本发明实施例提供的基于激光雷达的输电线覆冰在线监测方法。如图1所示，所述基于激光雷达的输电线覆冰在线监测方法包括：
在步骤S101中，获取激光雷达对输电线进行扫描得到的激光点云数据。
其中，所述激光雷达架设在输电铁塔上。本发明实施例通过所述激光雷达对输电线进行扫描，从而可以得到所述输电线对应的激光点云数据。所述激光点云数据中包括若干个激光点，每一个激光点均对应一个三维坐标，该三维坐标以激光雷达的坐标系原点作为坐标原点。
在步骤S102中，对所述激光点云数据进行分割，得到至少一种线路类型对应的激光点云。
在实际应用中，同时架设在输电铁塔上的线路有多种类型，激光雷达扫描得到的激光点云数据中包括若干种类型的输电线路导线和地线，本发明实施例首先对所述线路类型进行辨识，分割激光点云数据，得到每一种线路类型对应的激光点云。可选地，作为本发明的一个优选示例，如图2所示，步骤S102所述的对所述激光点云数据进行分割，得到至少一种线路类型对应的激光点云包括：
在步骤S201中，从所述激光点云数据中获取初始激光点。
其中，所述初始激光点为与坐标原点的相对距离最小的激光点。所述坐标原点为激光雷达的坐标系原点。本发明实施例遍历所述激光点云数据，获取初始激光点，为了便于描述，这里记为q0=（x0,y0,z0）。
在步骤S202中，从所述初始激光点开始，计算所述激光点云数据中的两两激光点之间的相对距离，构建距离比邻矩阵。
其中，所述距离比邻矩阵为以激光点的相对距离为元素构成的矩阵，其元素为所述激光点云数据中的两两激光点之间的相对距离。在得到初始激光点之后，以所述初始激光点为起始点，遍历所述激光点云数据中的每一激光点与其他激光点，计算两两激光点之间的相对距离，并根据所述相对距离构建距离比邻矩阵。
可选地，为了便于理解，假设所述激光点云数据中包括J个激光点，标记为激光点q0、q1、q2、q3……、qJ-1，所述距离比邻矩阵A所包括的元素如下表1所示：
Figure 580938DEST_PATH_IMAGE003
表1
在步骤S203中，构建原始激光点云，所述原始激光点云中包括所述激光点云数据中除所述初始激光点外其他的激光点。
可选地，为了便于理解，假设所述激光点云数据中包括J个激光点，标记为激光点q0、q1、q2、q3……、qJ-1，所述原始激光点云B为{ q1、q2、q3……、qJ-1 },将每一个激光点按照序号放置于所述原始激光点云B中。应当理解，每一个激光点均为一个三维坐标（xj,yj,zj）,j=1、2、3、…、J-1。
在步骤S204中，以所述初始激光点为搜索点，按预设方式搜索所述距离比邻矩阵，获取与所述搜索点的相对距离满足预设条件的激光点构成链接矩阵。
其中，一个链接矩阵表示一种线路类型对应的激光点云，实现对激光点云数据的分割。可选地，本发明实施例通过在所述距离比邻矩阵A中搜索相对距离满足预设条件的目标激光点，来实现对所述激光点云数据的分割。在搜索的时候，首先以所述初始激光点为搜索点进行搜索。作为本发明的一个优选示例，如图3所示，步骤S204所述的按预设方式搜索所述距离比邻矩阵，获取与所述搜索点的相对距离满足预设条件的激光点构成链接矩阵包括：
在步骤S301中，构建链接矩阵，将所述搜索点添加到所述链接矩阵。
在确定搜索点之后，建立链接矩阵，为了便于描述，这里用符号P表示链接矩阵。所述链接矩阵P为待生成的链接矩阵，初始创建时其元素个数为0。在开始搜索前，将搜索点添加到链接矩阵P，此时所述链接矩阵P包括一个元素，以初始化所述链接矩阵P，然后以所述搜索点开始搜索所述距离比邻矩阵A。示例性地，以所述初始激光点q0为搜索点为例，将所述初始激光点添加到所述链接矩阵P，即P为{ q0}，以所述搜索点q0开始搜索所述距离比邻矩阵A。
在步骤S302中，在所述距离比邻矩阵中进行搜索，获取所述搜索点对应的最小的n个相对距离，n>0。
在所述距离比邻矩阵A中搜索出所述搜索点与其他激光点的相对距离，将所搜索出来的相对距离进行排序，根据排序结果获取距离值最小的n个相对距离。可选地，作为本发明的一个优选示例，在实际应用中所述n取6，即获取所述搜索点对应的最小的6个相对距离。
在步骤S303中，若所述n个相对距离均小于预设的线路直径阈值，则获取所述n个相对距离对应的激光点，将所获得的激光点从原始激光点云中删除，并添加到所述链接矩阵。
在步骤S304中，以所述n个相对距离中的最小值对应的激光点作为新的搜索点，继续在所述距离比邻矩阵中进行搜索。
其中，所述线路直径阈值为不同输电线类型或地线的分割标准，作为一轮搜索的结束标志。可选地，所述线路直径阈值可以为10倍的导体直径。
在搜索得到n个相对距离之后，本发明实施例将所搜索到的n个相对距离分别与预设的线路直径阈值进行比较。
若所述n个相对距离均小于预设的线路直径阈值，表明所获得的n个相对距离对应的激光点均属于同一线路类型的扫描激光点，则获取所述n个相对距离对应的激光点，将所获得的激光点从原始激光点云B中删除，与此同时将所获得的激光点添加到所述链接矩阵P。示例性地，假设所述n个相对距离对应的激光点的序号分别为2、4、7、9、10、21，将所获得的激光点添加到所述链接矩阵P，所述链接矩阵P变为{q0、q2、q4、q7、q9、q10、q21}，所述原始激光点云B变为{q1、q3、q5、q6、q8、q11、……q20、q22、……、qJ-1}。
从所述n个相对距离中获取最小值及其对应的激光点，以所述相对距离最小值的激光点作为新的搜索点，返回步骤S302继续在所述距离比邻矩阵A中进行搜索。承接前文示例，假设在所述n个相对距离对应的激光点q2、q4、q7、q9、q10、q21中激光点q2与初始激光点q0的相对距离最小，则以所述激光点q2作为新的搜索点，返回所述距离比邻矩阵A中搜索与所述激光点q2之间最小的n个相对距离，重复执行步骤S302至步骤S304。
在步骤S305中，若所述n个相对距离中存在至少一个相对距离大于或等于预设的线路直径阈值，则获取小于所述线路直径阈值的相对距离对应的激光点，将所获得的激光点从原始激光点云中删除，并添加到所述链接矩阵，以所述链接矩阵为一种线路类型对应的激光点云。
若所述n个相对距离中存在至少一个相对距离，其大于或等于预设的线路直径阈值，表明检测到属于不同线路类型的扫描激光点。其中，对于小于所述线路直径阈值的相对距离，获取对应的激光点，将所获得的激光点从原始激光点云B中删除，与此同时将所获得的激光点添加到所述链接矩阵P。示例性地，所述n个相对距离中，假设有1个相对距离大于或等于所述预设的线路直径阈值，其余的n-1个相对距离小于所述预设的线路直径阈值，则将所述n-1个相对距离的激光点添加到所述链接矩阵P中；假设有2个相对距离大于或等于所述预设的线路直径阈值，其余的n-2个相对距离小于所述预设的线路直径阈值，则将所述n-2个相对距离小于所述预设的线路直径阈值的激光点添加到所述链接矩阵P中，以此类推。
此时所述链接矩阵P表示一个线路类型的全部数据，一个线路类型分割完成，将所 述链接矩阵P进行存储，并不再更改。在实际应用中可以根据链接矩阵P的生成顺序进行标 号命名，比如第一链接矩阵
Figure 710568DEST_PATH_IMAGE004
、第二连接矩阵
Figure 259361DEST_PATH_IMAGE005
、……。应当理解，一种线路类型包括但 不限于一种独立输电线路导线或地线。
在步骤S306中，构建下一个链接矩阵，在大于或等于预设的线路直径阈值的相对距离中获取最小值，以所述最小值对应的激光点作为新的搜索点，添加到所述下一个链接矩阵，在所述距离比邻矩阵中进行下一轮搜索，以生成下一种线路类型对应的激光点云。
在通过步骤S305完成对一种线路类型分割之后，构建下一个链接矩阵，从所述大 于或等于所述预设的线路直径阈值的相对距离中获取最小值及其对应的激光点，以所获得 的激光点作为新的搜索点，返回步骤S301继续在所述距离比邻矩阵A中进行搜索。示例性 地，假设所述大于或等于所述预设的线路直径阈值的相对距离中，最小值对应的激光点为 k，则将所述激光点k添加到新的链接矩阵，得到链接矩阵
Figure 816244DEST_PATH_IMAGE006
，以所述激光点
Figure 552119DEST_PATH_IMAGE007
作为新 的搜索点，继续在所述距离比邻矩阵A中进行搜索，重复执行步骤S301至步骤S306。
在步骤S307中，直至所述原始激光点云中的所述激光点纳入不同的链接矩阵中。
通过上述步骤S301至步骤S306，可以生成至少一个链接矩阵，若所述原始激光点云B不为空，则继续构建下一个链接矩阵，并以大于或等于所述预设的线路直径阈值的相对距离中的最小值对应的激光点作为搜索点进行搜索。直至所述原始激光点云B中的激光点均被纳入不同的链接矩阵中，从而完成对激光点云数据的分割，得到至少一种线路类型及其对应的激光点云，实现对不同的线路类型的辨识，提高了线路类型分割的精度。
在步骤S103中，对于每一种线路类型，根据所述线路类型对应的激光点云计算所述线路类型在每一数据段中的每一片层对应的外轮廓多边形。
其中，所述片层是指所述线路类型在Z轴方向上所划分区，所述片层对应的外轮廓 多边形是指根据所述片层上的激光点云构成的激光点外延。在得到至少一种线路类型之 后，分别就每一种线路类型
Figure 485440DEST_PATH_IMAGE008
分割数据段以及切割片层，并基于凸包理念 计算每一片层对应的外轮廓多边形。可选地，作为本发明的一个优选示例，步骤S103所述的 根据所述线路类型对应的激光点云计算所述线路类型在每一数据段中的每一片层对应的 外轮廓多边形包括：
在步骤S401中，获取一种线路类型对应的激光点云，按预设分割距离将所述激光 点云沿X轴方向分割为若干个数据段，根据每一个数据段中激光点云在Z轴方向的坐标最大 值和最小值进行等间距切割，得到
Figure 888739DEST_PATH_IMAGE009
个片层及对应的点云片层数据。
其中，所述线路沿X轴方向延伸。所述Z轴坐标是指纵轴方向上的坐标，所述Z轴是 指三维坐标系中垂直于横轴X和竖轴Y所成平面的坐标轴。本发明实施例以首先以X轴方向 为基准，对每一种线路类型对应的所有激光点云，按照预设分割距离进行切分，得到若干个 数据段。可选地所述预设分割距离为20厘米，最后不能整除的算作独立的一个数据段。然后 以Z轴方向为基准，对于一种线路类型对应的所有激光点云，获取其中的Z轴坐标的最大值 和最小值，对所述最大值和最小值之间的激光点云进行等间距切割，得到
Figure 367256DEST_PATH_IMAGE010
个片层及点云 片层数据。每一点云片层数据中包括若干个激光点，每一激光点的Z轴坐标均落在所述点云 片层对应的Z轴坐标间隔内。其中，
Figure 590427DEST_PATH_IMAGE009
为大于或等于1的正整数。
在步骤S402中，遍历每一个点云片层数据，采用步进法对所述点云片层数据中的激光点进行排序，生成所述线路类型在每一个数据段中的每一片层对应的外轮廓多边形。
其中，所述外轮廓多边形表示点云片层数据中所包括的激光点的轮廓边界。在所 得到的每一个数据段中的每一个片层中，激光点的分布是混乱的。本发明实施例采用步进 法对乱序的点云片层数据进行排序，从而得到每一个片层对应的外轮廓多边形。具体地，本 发明实施例运用凸包理念，通过调用一个convexhull函数包来实现，最后留在栈里面的激 光点是构成外轮廓多边形的所有凸包外围点，所述外轮廓多边形为由所有凸包外围点按照 逆时针回路首尾相接构成。为了便于描述，承接前文示例，假设第
Figure 796281DEST_PATH_IMAGE011
种线路类型划分出了
Figure 788507DEST_PATH_IMAGE012
个数据段、每个数据段中划分出了
Figure 952772DEST_PATH_IMAGE009
个片层及点云片层数据，所述外轮廓多边形记为
Figure 646928DEST_PATH_IMAGE013
, 其中，
Figure 656472DEST_PATH_IMAGE011
表示第
Figure 34364DEST_PATH_IMAGE011
种线路类型，
Figure 838372DEST_PATH_IMAGE014
表示所述第
Figure 301714DEST_PATH_IMAGE011
种线路类型的第
Figure 865682DEST_PATH_IMAGE014
个数据段，
Figure 98080DEST_PATH_IMAGE015
，
Figure 338569DEST_PATH_IMAGE016
表 示所述第
Figure 758049DEST_PATH_IMAGE011
种线路类型的第
Figure 374975DEST_PATH_IMAGE014
个数据段中的第
Figure 445568DEST_PATH_IMAGE017
个片层，
Figure 856958DEST_PATH_IMAGE018
，外轮廓多边形
Figure 29313DEST_PATH_IMAGE013
包 括
Figure 918772DEST_PATH_IMAGE019
个顶点。
在步骤S104中，根据所述线路类型在每一数据段中的每一片层对应的外轮廓多边形，计算所述线路类型对应的覆冰重量和覆冰厚度。
在得到线路类型在每一个数据段中的每一片层对应的外轮廓多边形
Figure 125762DEST_PATH_IMAGE013
之后，所 述外轮廓多边形
Figure 973632DEST_PATH_IMAGE013
将用于计算所述线路类型对应的覆冰重量和覆冰厚度。
可选地，作为本发明的一个优选示例，步骤S104所述的根据所述线路类型在每一数据段中的每一片层对应的外轮廓多边形，计算所述线路类型对应的覆冰重量包括：
在步骤S501中，根据点云片层数据计算所述线路类型在每一数据段中的每一个片层对应的外轮廓多边形的面积。
承接前文示例，假设第
Figure 384016DEST_PATH_IMAGE020
种线路类型的第
Figure 342745DEST_PATH_IMAGE021
个数据段中的第
Figure 873083DEST_PATH_IMAGE022
个片层的外轮廓多 边
Figure 891855DEST_PATH_IMAGE013
包括
Figure 38802DEST_PATH_IMAGE023
个顶点，所述外轮廓多边形的面积
Figure 535643DEST_PATH_IMAGE024
为所述
Figure 435335DEST_PATH_IMAGE023
个顶点按照逆时针回路首 尾相接构成的多边形面积。所述外轮廓多边形
Figure 93849DEST_PATH_IMAGE013
对应的面积
Figure 462514DEST_PATH_IMAGE024
为：
Figure 497466DEST_PATH_IMAGE025
其中，
Figure 284287DEST_PATH_IMAGE026
表示顶点个数，
Figure 379282DEST_PATH_IMAGE027
表示第
Figure 500822DEST_PATH_IMAGE028
个顶点的横轴坐标，
Figure 605044DEST_PATH_IMAGE029
表示第
Figure 230061DEST_PATH_IMAGE028
个顶点 的竖轴坐标，
Figure 495957DEST_PATH_IMAGE030
。
在步骤S502中，对每一个点云片层数据的高度和外轮廓多边形的面积之积进行求和处理，得到所述线路类型中每一数据段对应的体积。
遍历所述第
Figure 822902DEST_PATH_IMAGE031
种线路类型的所有数据段中的所有点云片层数据，计算得到每一点 云片层数据对应的外轮廓多边形的面积
Figure 730815DEST_PATH_IMAGE032
。然后针对所述第
Figure 210338DEST_PATH_IMAGE033
个 数据段，求取每一点云片层数据的高度
Figure 912715DEST_PATH_IMAGE034
,再根据高度
Figure 477688DEST_PATH_IMAGE034
和每一点云片层数据的外轮廓多 边形的面积
Figure 923713DEST_PATH_IMAGE035
之积得到每一片层对应的体积，对每一片层对应的体积进行求和处理，从 而得到所述第
Figure 274054DEST_PATH_IMAGE031
种线路类型对应的第
Figure 881753DEST_PATH_IMAGE033
个数据段的体积。计算公式为：
Figure 199602DEST_PATH_IMAGE036
其中，
Figure 183738DEST_PATH_IMAGE037
表示所述第
Figure 903433DEST_PATH_IMAGE031
种线路类型中的第
Figure 947612DEST_PATH_IMAGE033
个数据段对应的点云数据在Z轴上的 高度，通过求取点云数据中Z轴坐标最大值和最小值之差得到；
Figure 736445DEST_PATH_IMAGE038
表示所述第
Figure 524273DEST_PATH_IMAGE031
种线路类 型的第
Figure 98473DEST_PATH_IMAGE033
个数据段所切分出来的点云片层个数，
Figure 47975DEST_PATH_IMAGE039
表示一个点云片层数据对应的高度；
Figure 809258DEST_PATH_IMAGE040
表示第
Figure 885929DEST_PATH_IMAGE031
种线路类型的第
Figure 314636DEST_PATH_IMAGE033
个数据段对应的体积。
遍历所述第
Figure 435039DEST_PATH_IMAGE031
种线路类型的每一个数据段，得到所述线路类型中每一个数据段对 应的体积。
在步骤S503中，获取气象特征信息，所述气象特征信息包括温度信息、湿度信息、风速信息、时间信息。
在步骤S504中，根据所述气象特征信息对冰密度进行修正。
其中，冰密度是指输电线上所覆冰的单位质量。在对冰密度进行修正时可由当前气象台检测数据构建感知机去预测冰的密度。修正方法如下：
Figure 683618DEST_PATH_IMAGE041
Figure 796936DEST_PATH_IMAGE042
其中，
Figure 80150DEST_PATH_IMAGE043
表示气象特征信息，
Figure 902613DEST_PATH_IMAGE044
表示感知机的隐藏层输出，
Figure 904067DEST_PATH_IMAGE045
表示修正后的冰密 度，
Figure 571808DEST_PATH_IMAGE046
表示权重系数矩阵，
Figure 709529DEST_PATH_IMAGE047
表示激活函数，
Figure 176327DEST_PATH_IMAGE048
，
Figure 665077DEST_PATH_IMAGE049
表示隐藏层输出 的变量维度，
Figure 870930DEST_PATH_IMAGE050
表示第
Figure 128736DEST_PATH_IMAGE051
维隐藏层输出，
Figure 27422DEST_PATH_IMAGE052
，
Figure 721578DEST_PATH_IMAGE053
为1.7，为个人通过数据调参后的结 果。
在步骤S505中，根据修正后的冰密度、所述线路类型在每一数据段对应的体积计算所述线路类型在每一数据段的覆冰重量。
本发明实施例预先根据现有线路类型构建线路数据库，存储现有线路类型对应的导线密度及截面积，导线密度为单位长度的线路具有的质量。在得到修正后的冰密度和线路类型对应的体积之后，从所述线路数据库中获得所述线路类型对应的导线密度。示例性地，类型为LGJ-400/35的架空线路，其导线密度为1.349kg/m。
然后按照以下公式计算：
Figure 996701DEST_PATH_IMAGE054
其中，
Figure 374593DEST_PATH_IMAGE055
表示第
Figure 178601DEST_PATH_IMAGE056
种线路类型在第
Figure 641943DEST_PATH_IMAGE057
个数据段的单位覆冰重量，
Figure 455178DEST_PATH_IMAGE058
表示第
Figure 687577DEST_PATH_IMAGE056
种 线路类型在第
Figure 678798DEST_PATH_IMAGE057
个数据段对应的体积，
Figure 629436DEST_PATH_IMAGE059
表示所述第
Figure 246362DEST_PATH_IMAGE056
种线路类型在第
Figure 333267DEST_PATH_IMAGE057
个数据段对应 的点云数据在Z轴上的高度，通过求取点云数据中Z轴坐标最大值和最小值之差得到，
Figure 744657DEST_PATH_IMAGE060
表示修正后的冰密度，
Figure 651433DEST_PATH_IMAGE061
表示线路类型对应的导线密度。
可选地，作为本发明的一个优选示例，步骤S104中所述的根据所述线路类型在每一数据段中的每一片层对应的外轮廓多边形，计算所述线路类型对应的覆冰厚度包括：
在步骤S601中，获取所述线路类型中每一个数据段对应的导体中心点。
本发明实施例获取所述线路类型对应的导体中心点，记为
Figure 55738DEST_PATH_IMAGE062
。
在步骤S602中，对于所述线路类型每一数据段中的每一片层，计算所述片层对应的外轮廓多边形上的激光点与所述导体中心点的相对距离，对所述相对距离进行加权求和，得到所述线路类型在每一数据段中的片层对应的覆冰厚度。
在这里，本发明实施例针对每个点云片层数据形成的外轮廓多边形，求取外轮廓多边形上的激光点与导体中心点之间的相对距离并进行加权平均，公式如下：
Figure 997149DEST_PATH_IMAGE063
其中，
Figure 579440DEST_PATH_IMAGE064
表示第
Figure 239092DEST_PATH_IMAGE065
种线路类型在第
Figure 197821DEST_PATH_IMAGE066
个数据段的第
Figure 259318DEST_PATH_IMAGE067
片层上对应的覆冰厚度，
Figure 763242DEST_PATH_IMAGE068
表示第
Figure 910190DEST_PATH_IMAGE065
种线路类型在第
Figure 407030DEST_PATH_IMAGE066
个数据段的第
Figure 323034DEST_PATH_IMAGE067
片层上对应的外轮廓多边形上的第
Figure 512707DEST_PATH_IMAGE069
个激光点，
Figure 146950DEST_PATH_IMAGE070
表示第
Figure 447482DEST_PATH_IMAGE065
种线路类型在第
Figure 201680DEST_PATH_IMAGE066
个数据段的第
Figure 562254DEST_PATH_IMAGE067
片层上对应的外轮廓多边形上 的激光点个数，
Figure 683794DEST_PATH_IMAGE071
表示第
Figure 522437DEST_PATH_IMAGE066
个数据段的导体中心点。
在步骤S603中，计算所述线路类型在所有数据段的所有片层对应的覆冰厚度的平均值，得到所述线路类型对应的覆冰厚度。
遍历所述线路类型的每一个数据段的每一个片层，计算得到每一个数据段的每一片层对应的覆冰厚度，对所有数据段的所有片层对应的覆冰厚度计算平均值，从而得到所述线路类型对应的覆冰厚度。
其中，计算公式为
Figure 413032DEST_PATH_IMAGE072
Figure 678929DEST_PATH_IMAGE073
表示第
Figure 22185DEST_PATH_IMAGE074
种线路类型对应的覆冰厚度，
Figure 680831DEST_PATH_IMAGE075
表示所述第
Figure 425933DEST_PATH_IMAGE076
种线路类型所切分出 来的数据段个数，
Figure 862731DEST_PATH_IMAGE077
表示所述第
Figure 427704DEST_PATH_IMAGE078
种线路类型的第
Figure 608150DEST_PATH_IMAGE079
个数据段所切分出来的点云片层个 数，
Figure 473338DEST_PATH_IMAGE080
表示第
Figure 330304DEST_PATH_IMAGE081
种线路类型在第
Figure 648153DEST_PATH_IMAGE082
个数据段的第
Figure 632289DEST_PATH_IMAGE083
片层上对应的覆冰厚度。
至此，本发明实施例通过采用点云切片法，对不规则覆冰进行分层，并通过对每一数据段的每一层确定外轮廓多边形，基于外轮廓多边形实现对线路的覆冰重量和覆冰厚度计算，有效地提高了线路在覆冰情况下的覆冰重量和覆冰厚度的计算精度，提升了线路覆冰在线监测的效率。且通过分割数据段，可以有效地解决覆冰不连续不均匀，进一步增加计算覆冰重量和覆冰厚度的精度。
可选地，作为本发明的一个优选示例，本发明实施例还可以根据覆冰重量和覆冰厚度计算覆冰情况下线路的综合比载，所述综合比载是指单位长度架空线上所受的载荷折算到单位截面积上的数值，所述载荷受到覆冰厚度和覆冰重量的影响，因此叫做综合比载。所述方法还包括：
在步骤S701中，根据所述线路类型对应的导线密度和截面积计算所述线路类型对应的导线自重比载。
其中，导线比载是指单位长度架空线上所受的载荷折算到单位截面积上的数值，由于导体本身重量导致的比载叫自重比载，由于覆冰导致的比载叫做冰重比载。本发明实施例首先从所述线路数据库获取线路类型对应的密度和截面积，然后按照公式计算线路对应的导线自重比载：
Figure 351984DEST_PATH_IMAGE084
其中，
Figure DEST_PATH_IMAGE085
表示所述线路类型对应的导线自重比载，单位为N/m.mm2,
Figure 865005DEST_PATH_IMAGE086
表示所述 线路类型对应的密度,为每公里导线的质量，单位为kg/km，
Figure DEST_PATH_IMAGE087
表示导线的截面积，单位为 mm2。
在步骤S702中，根据所述线路类型的覆冰重量获取单位距离覆冰重量，根据所述单位距离覆冰重量计算冰重比载。
冰重比载的计算公式为：
Figure 155303DEST_PATH_IMAGE088
在上式中，
Figure 943130DEST_PATH_IMAGE089
表示冰重比载，
Figure 251752DEST_PATH_IMAGE090
表示第种线路类型在第
Figure 466833DEST_PATH_IMAGE091
个数据段的单位覆冰 重量。
在步骤S703中，根据所述导线自重比载和冰重比载，计算所述线路类型对应的冰重总比载。
其中，所述冰重总比载的计算公式为：
Figure 493694DEST_PATH_IMAGE092
，
Figure 85213DEST_PATH_IMAGE093
表示冰重总比载。
在步骤S704中，获取风载体型系数，根据所述风载体型系数和所述线路类型对应的覆冰厚度计算覆冰风压比载。
其中，所述覆冰风压比载是指线路在覆冰之后每平方毫米的风压载荷。所述覆冰 风压比载的计算公式为
Figure 763188DEST_PATH_IMAGE094
。
在上式中，
Figure 149170DEST_PATH_IMAGE095
表示风载体型系数，当线路直径小于17毫米时，
Figure DEST_PATH_IMAGE096
，否则
Figure 397748DEST_PATH_IMAGE097
，
Figure 527378DEST_PATH_IMAGE098
表示线路的截面积，单位为mm2，
Figure 810592DEST_PATH_IMAGE099
表示导体直径，单位为mm，
Figure 633055DEST_PATH_IMAGE100
表示设计风 速，单位为m/s，
Figure 385241DEST_PATH_IMAGE101
表示风速不均匀系数，可通过表2查询得到：
设计风速（m/s） 20以下 20-30 30-35 35以上
Figure 52983DEST_PATH_IMAGE102
 1.0 0.85 0.75 0.70
表2
在步骤S705中，根据所述线路类型对应的冰重总比载和覆冰风压比载计算所述线路类型在覆冰时的综合比载。
其中，综合比载的计算公式为
Figure 456282DEST_PATH_IMAGE103
，
综上所述，本发明实施例提供的基于激光雷达的输电线覆冰在线监测方法，既提高了线路在覆冰情况下的覆冰重量和覆冰厚度的计算精度，而且提高了对线路的综合比载的计算精度，有效地提升了线路覆冰在线监测的效率。
应理解，上述实施例中各步骤的序号的大小并不意味着执行顺序的先后，各过程的执行顺序应以其功能和内在逻辑确定，而不应对本发明实施例的实施过程构成任何限定。
在一实施例中，本发明还提供一种基于激光雷达的输电线覆冰在线监测装置，该基于激光雷达的输电线覆冰在线监测装置与上述实施例中基于激光雷达的输电线覆冰在线监测方法一一对应。如图4所示，该基于激光雷达的输电线覆冰在线监测包括获取模块41、分割模块42、多边形获取模块43、覆冰监测模块44。各功能模块详细说明如下：
获取模块41，用于获取激光雷达对输电线进行扫描得到的激光点云数据，所述激光雷达架设在输电铁塔上；
分割模块42，用于对所述激光点云数据进行分割，得到至少一种线路类型对应的激光点云；
多边形获取模块43，用于对于每一种线路类型，根据所述线路类型对应的激光点云计算所述线路类型在每一数据段中的每一片层对应的外轮廓多边形；
覆冰监测模块44，用于根据所述线路类型在每一数据段中的每一片层对应的外轮廓多边形，计算所述线路类型对应的覆冰重量和覆冰厚度。
可选地，所述分割模块42包括：
初始激光点获取单元，用于从所述激光点云数据中获取初始激光点，所述初始激光点为与坐标原点的相对距离最小的激光点；
第一构建单元，用于从所述初始激光点开始，计算所述激光点云数据中的两两激光点之间的相对距离，构建距离比邻矩阵；
第二构建单元，用于构建原始激光点云，所述原始激光点云中包括所述激光点云数据中除所述初始激光点外其他的激光点；
搜索分割单元，用于以所述初始激光点为搜索点，按预设方式搜索所述距离比邻矩阵，获取与所述搜索点的相对距离满足预设条件的激光点构成链接矩阵，其中，一个链接矩阵表示一种线路类型对应的激光点云。
可选地，所述搜索分割单元包括：
连接矩阵构建子单元，用于构建链接矩阵，将所述搜索点添加到所述链接矩阵；
搜索子单元，用于在所述距离比邻矩阵中进行搜索，获取所述搜索点对应的最小的n个相对距离，n>0；
删除添加子单元，用于若所述n个相对距离均小于预设的线路直径阈值，则获取所述n个相对距离对应的激光点，将所获得的激光点从原始激光点云中删除，并添加到所述链接矩阵；
所述搜索子单元还用于，以所述n个相对距离中的最小值对应的激光点作为新的搜索点，继续在所述距离比邻矩阵中进行搜索；
所述删除添加子单元还用于，若所述n个相对距离中存在至少一个相对距离大于或等于预设的线路直径阈值，则获取小于所述线路直径阈值的相对距离对应的激光点，将所获得的激光点从原始激光点云中删除，并添加到所述链接矩阵，以所述链接矩阵为一种线路类型对应的激光点云；
所述连接矩阵构建子单元还用于，构建下一个链接矩阵，在大于或等于预设的线路直径阈值的相对距离中获取最小值，以所述最小值对应的激光点作为新的搜索点，添加到所述下一个链接矩阵，在所述距离比邻矩阵中进行下一轮搜索，以生成下一种线路类型对应的激光点云；直至所述原始激光点云中的所述激光点纳入不同的链接矩阵中。
可选地，所述获取所述搜索点对应的最小的n个相对距离为：
获取所述搜索点对应的最小的6个相对距离。
可选地，所述多边形获取模块43包括：
片层切割单元，用于获取一种线路类型对应的激光点云，按预设分割距离将所述 激光点云沿X轴方向分割为若干个数据段，根据每一个数据段中激光点云在Z轴方向的坐标 最大值和最小值进行等间距切割，得到
Figure 184067DEST_PATH_IMAGE104
个片层及对应的点云片层数据，
Figure 407238DEST_PATH_IMAGE104
为大于或等 于1的正整数，每一点云片层数据中包括若干个激光点，每一激光点的Z轴坐标均落在所述 点云片层对应的Z轴坐标间隔内；
多边形生成单元，用于遍历每一个点云片层数据，采用步进法对所述点云片层数据中的激光点进行排序，生成所述线路类型在每一个数据段中的每一片层对应的外轮廓多边形。
可选地，所述覆冰监测模块44包括：
面积计算单元，用于根据点云片层数据计算所述线路类型在每一数据段中的每一个片层对应的外轮廓多边形的面积；
体积计算单元，用于对每一个点云片层数据的高度和外轮廓多边形的面积之积进行求和处理，得到所述线路类型中每一数据段对应的体积；
气象信息获取单元，用于获取气象特征信息，所述气象特征信息包括温度信息、湿度信息、风速信息、时间信息。
修正单元，用于根据所述气象特征信息对冰密度进行修正；
覆冰重量计算单元，用于根据修正后的冰密度、所述线路类型在每一数据段对应的体积计算所述线路类型在每一数据段的覆冰重量。
可选地，所述覆冰监测模块44包括：
导体中心点获取单元，用于获取所述线路类型中每一个数据段对应的导体中心点；
第一覆冰厚度计算单元，用于对于所述线路类型每一数据段中的每一片层，计算所述片层对应的外轮廓多边形上的激光点与所述导体中心点的相对距离，对所述相对距离进行加权求和，得到所述线路类型在每一数据段中的片层对应的覆冰厚度；
第二覆冰厚度计算单元，用于计算所述线路类型在所有数据段的所有片层对应的覆冰厚度的平均值，得到所述线路类型对应的覆冰平均厚度。
可选地，所述装置还包括：
导线自重比载计算模块，用于根据所述线路类型对应的导线密度和截面积计算所述线路类型对应的导线自重比载；
冰重比载计算模块，用于根据所述线路类型的覆冰重量获取单位距离覆冰重量，根据所述单位距离覆冰重量计算冰重比载；
冰重总比载计算模块，用于根据所述导线自重比载和冰重比载，计算所述线路类型对应的冰重总比载；
覆冰风压比载计算模块，用于获取风载体型系数，根据所述风载体型系数和所述线路类型对应的覆冰厚度计算覆冰风压比载；
综合比载计算模块，用于根据所述线路类型对应的冰重总比载和覆冰风压比载计算所述线路类型在覆冰时的综合比载。
关于基于激光雷达的输电线覆冰在线监测装置的具体限定可以参见上文中对于基于激光雷达的输电线覆冰在线监测方法的限定，在此不再赘述。上述基于激光雷达的输电线覆冰在线监测装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中，也可以以软件形式存储于存储器中，以便于处理器调用执行以上各个模块对应的操作。
在一个实施例中，提供了一种计算机设备，该计算机设备可以是服务器，其内部结构图可以如图5所示。该计算机设备包括通过系统总线连接的处理器、存储器、网络接口和数据库。其中，该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作系统、计算机程序和数据库。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该计算机设备的网络接口用于与外部的终端通过网络连接通信。该计算机程序被处理器执行时以实现一种基于激光雷达的输电线覆冰在线监测方法。
在一个实施例中，提供了一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，处理器执行计算机程序时实现以下步骤：
获取激光雷达对输电线进行扫描得到的激光点云数据，所述激光雷达架设在输电铁塔上；
对所述激光点云数据进行分割，得到至少一种线路类型对应的激光点云；
对于每一种线路类型，根据所述线路类型对应的激光点云计算所述线路类型在每一数据段中的每一片层对应的外轮廓多边形；
根据所述线路类型在每一数据段中的每一片层对应的外轮廓多边形，计算所述线路类型对应的覆冰重量和覆冰厚度。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一非易失性计算机可读取存储介质中，该计算机程序在执行时，可包括如上述各方法的实施例的流程。其中，本发明所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用，均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器（ROM）、可编程ROM（PROM）、电可编程ROM（EPROM）、电可擦除可编程ROM（EEPROM）或闪存。易失性存储器可包括随机存取存储器（RAM）或者外部高速缓冲存储器。作为说明而非局限，RAM以多种形式可得，诸如静态RAM（SRAM）、动态RAM（DRAM）、同步DRAM（SDRAM）、双数据率SDRAM（DDRSDRAM）、增强型SDRAM（ESDRAM）、同步链路（Synchlink） DRAM（SLDRAM）、存储器总线（Rambus）直接RAM（RDRAM）、直接存储器总线动态RAM（DRDRAM）、以及存储器总线动态RAM（RDRAM）等。
所属领域的技术人员可以清楚地了解到，为了描述的方便和简洁，仅以上述各功能单元、模块的划分进行举例说明，实际应用中，可以根据需要而将上述功能分配由不同的功能单元、模块完成，即将所述装置的内部结构划分成不同的功能单元或模块，以完成以上描述的全部或者部分功能。
以上所述实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围，均应包含在本发明的保护范围之内。

Claims (10)

1.一种基于激光雷达的输电线覆冰在线监测方法，其特征在于，所述方法包括：
获取激光雷达对输电线进行扫描得到的激光点云数据，所述激光雷达架设在输电铁塔上；
对所述激光点云数据进行分割，得到至少一种线路类型对应的激光点云；
对于每一种线路类型，根据所述线路类型对应的激光点云计算所述线路类型在每一数据段中的每一片层对应的外轮廓多边形；
根据所述线路类型在每一数据段中的每一片层对应的外轮廓多边形，计算所述线路类型对应的覆冰重量和覆冰厚度。
2.如权利要求1所述的基于激光雷达的输电线覆冰在线监测方法，其特征在于，所述对所述激光点云数据进行分割，得到至少一种线路类型对应的激光点云包括：
从所述激光点云数据中获取初始激光点，所述初始激光点为与坐标原点的相对距离最小的激光点；
从所述初始激光点开始，计算所述激光点云数据中的两两激光点之间的相对距离，构建距离比邻矩阵；
构建原始激光点云，所述原始激光点云中包括所述激光点云数据中除所述初始激光点外其他的激光点；
以所述初始激光点为搜索点，按预设方式搜索所述距离比邻矩阵，获取与所述搜索点的相对距离满足预设条件的激光点构成链接矩阵，其中，一个链接矩阵表示一种线路类型对应的激光点云。
3.如权利要求2所述的基于激光雷达的输电线覆冰在线监测方法，其特征在于，所述按预设方式搜索所述距离比邻矩阵，获取与所述搜索点的相对距离满足预设条件的激光点构成链接矩阵包括：
构建链接矩阵，将所述搜索点添加到所述链接矩阵；
在所述距离比邻矩阵中进行搜索，获取所述搜索点对应的最小的n个相对距离，n>0；
若所述n个相对距离均小于预设的线路直径阈值，则获取所述n个相对距离对应的激光点，将所获得的激光点从原始激光点云中删除，并添加到所述链接矩阵；
以所述n个相对距离中的最小值对应的激光点作为新的搜索点，继续在所述距离比邻矩阵中进行搜索；
若所述n个相对距离中存在至少一个相对距离大于或等于预设的线路直径阈值，则获取小于所述线路直径阈值的相对距离对应的激光点，将所获得的激光点从原始激光点云中删除，并添加到所述链接矩阵，以所述链接矩阵为一种线路类型对应的激光点云；
构建下一个链接矩阵，在大于或等于预设的线路直径阈值的相对距离中获取最小值，以所述最小值对应的激光点作为新的搜索点，添加到所述下一个链接矩阵，在所述距离比邻矩阵中进行下一轮搜索，以生成下一种线路类型对应的激光点云；
直至所述原始激光点云中的所述激光点纳入不同的链接矩阵中。
4.如权利要求3所述的基于激光雷达的输电线覆冰在线监测方法，其特征在于，所述获取所述搜索点对应的最小的n个相对距离为：
获取所述搜索点对应的最小的6个相对距离。
5.如权利要求1至4任一项所述的基于激光雷达的输电线覆冰在线监测方法，其特征在于，所述根据所述线路类型对应的激光点云计算所述线路类型在每一数据段中的每一片层对应的外轮廓多边形包括：
获取一种线路类型对应的激光点云，按预设分割距离将所述激光点云沿X轴方向分割为若干个数据段，根据每一个数据段中激光点云在Z轴方向的坐标最大值和最小值进行等间距切割，得到个片层及对应的点云片层数据，为大于或等于1的正整数，每一点云片层数据中包括若干个激光点，每一激光点的Z轴坐标均落在所述点云片层对应的Z轴坐标间隔内；
遍历每一个点云片层数据，采用步进法对所述点云片层数据中的激光点进行排序，生成所述线路类型在每一个数据段中的每一片层对应的外轮廓多边形。
6.如权利要求5所述的基于激光雷达的输电线覆冰在线监测方法，其特征在于，所述根据所述线路类型在每一数据段中的每一片层对应的外轮廓多边形，计算所述线路类型对应的覆冰重量包括：
根据点云片层数据计算所述线路类型在每一数据段中的每一个片层对应的外轮廓多边形的面积；
对每一个点云片层数据的高度和外轮廓多边形的面积之积进行求和处理，得到所述线路类型中每一数据段对应的体积；
获取气象特征信息，所述气象特征信息包括温度信息、湿度信息、风速信息、时间信息；
根据所述气象特征信息对冰密度进行修正；
根据修正后的冰密度、所述线路类型在每一数据段对应的体积计算所述线路类型在每一数据段的覆冰重量。
7.如权利要求5所述的基于激光雷达的输电线覆冰在线监测方法，其特征在于，所述根据所述线路类型在每一数据段中的每一片层对应的外轮廓多边形，计算所述线路类型对应的覆冰厚度包括：
获取所述线路类型中每一个数据段对应的导体中心点；
对于所述线路类型每一数据段中的每一片层，计算所述片层对应的外轮廓多边形上的激光点与所述导体中心点的相对距离，对所述相对距离进行加权求和，得到所述线路类型在每一数据段中的片层对应的覆冰厚度；
计算所述线路类型在所有数据段的所有片层对应的覆冰厚度的平均值，得到所述线路类型对应的覆冰平均厚度。
8.如权利要求1所述的基于激光雷达的输电线覆冰在线监测方法，其特征在于，所述方法还包括：
根据所述线路类型对应的导线密度和截面积计算所述线路类型对应的导线自重比载；
根据所述线路类型的覆冰重量获取单位距离覆冰重量，根据所述单位距离覆冰重量计算冰重比载；
根据所述导线自重比载和冰重比载，计算所述线路类型对应的冰重总比载；
获取风载体型系数，根据所述风载体型系数和所述线路类型对应的覆冰厚度计算覆冰风压比载；
根据所述线路类型对应的冰重总比载和覆冰风压比载计算所述线路类型在覆冰时的综合比载。
9.一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现如权利要求1至8任一项所述的基于激光雷达的输电线覆冰在线监测方法。
10.一种计算机设备，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现如权利要求1至8任一项所述的基于激光雷达的输电线覆冰在线监测方法。
CN202211373417.3A 2022-11-04 2022-11-04 基于激光雷达的输电线覆冰在线监测方法、介质及设备 Active CN115508806B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211373417.3A CN115508806B (zh) 2022-11-04 2022-11-04 基于激光雷达的输电线覆冰在线监测方法、介质及设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211373417.3A CN115508806B (zh) 2022-11-04 2022-11-04 基于激光雷达的输电线覆冰在线监测方法、介质及设备

Publications (2)

Publication Number Publication Date
CN115508806A true CN115508806A (zh) 2022-12-23
CN115508806B CN115508806B (zh) 2023-03-14

Family

ID=84512148

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211373417.3A Active CN115508806B (zh) 2022-11-04 2022-11-04 基于激光雷达的输电线覆冰在线监测方法、介质及设备

Country Status (1)

Country Link
CN (1) CN115508806B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104122560A (zh) * 2014-07-11 2014-10-29 国家电网公司 一种输电线路广域冰情监测方法
US20190235011A1 (en) * 2018-01-26 2019-08-01 LineVision, Inc. System and method for power transmission line monitoring
CN112329221A (zh) * 2020-11-02 2021-02-05 云南师范大学 输电线覆冰厚度确定方法
CN113983987A (zh) * 2021-09-26 2022-01-28 国电南瑞科技股份有限公司 一种输电线路覆冰动态监测方法、系统、装置及存储介质
CN114241292A (zh) * 2021-11-05 2022-03-25 南方电网数字电网研究院有限公司 基于智能网关的输电线路覆冰监测方法和装置
CN115014216A (zh) * 2022-06-27 2022-09-06 山东信通电子股份有限公司 一种输电线路覆冰检测方法及装置
CN115143916A (zh) * 2022-06-14 2022-10-04 山东信通电子股份有限公司 一种输电线路覆冰状态监测方法、设备及介质

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104122560A (zh) * 2014-07-11 2014-10-29 国家电网公司 一种输电线路广域冰情监测方法
US20190235011A1 (en) * 2018-01-26 2019-08-01 LineVision, Inc. System and method for power transmission line monitoring
CN112329221A (zh) * 2020-11-02 2021-02-05 云南师范大学 输电线覆冰厚度确定方法
CN113983987A (zh) * 2021-09-26 2022-01-28 国电南瑞科技股份有限公司 一种输电线路覆冰动态监测方法、系统、装置及存储介质
CN114241292A (zh) * 2021-11-05 2022-03-25 南方电网数字电网研究院有限公司 基于智能网关的输电线路覆冰监测方法和装置
CN115143916A (zh) * 2022-06-14 2022-10-04 山东信通电子股份有限公司 一种输电线路覆冰状态监测方法、设备及介质
CN115014216A (zh) * 2022-06-27 2022-09-06 山东信通电子股份有限公司 一种输电线路覆冰检测方法及装置

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
彭业昌等: "基于激光测量的输电线路覆冰监测系统", 《国外电子测量技术》 *
林先堪等: "基于点云数据的电力线三维重建技术研究", 《机电信息》 *
麻卫峰等: "基于激光点云的高压输电线覆冰厚度反演", 《电力系统保护与控制》 *

Also Published As

Publication number Publication date
CN115508806B (zh) 2023-03-14

Similar Documents

Publication Publication Date Title
CN110908014B (zh) 舞动精细化订正预报方法及系统
CN108169751B (zh) 一种天气雷达基数据三维栅格化方法、计算机可读存储介质及电子设备
CN108228722B (zh) 破碎化区域采样点的地理空间分布均匀度检测方法
CN110784883B (zh) 基站建设评估方法、装置、设备和存储介质
CN110909260B (zh) 车辆显示方法、装置、计算机可读存储介质和计算机设备
CN111460237B (zh) 一种数据查询的方法、装置、可读存储介质和电子设备
CN116182859A (zh) 路径规划方法及装置、设备、存储介质
CN112926465A (zh) 基于点云类型的海岸线性质识别方法及装置
JP2011033987A (ja) 地図生成装置及び地図生成方法
CN116540775A (zh) 基于通信盲点的无人飞行基站自动巡航方法、设备
CN112101229A (zh) 点云数据特征点提取方法、装置、计算机设备及存储介质
CN115508806B (zh) 基于激光雷达的输电线覆冰在线监测方法、介质及设备
CN106919783A (zh) 一种浮标数据的多目标关联度划分处理方法
CN113747342A (zh) 网络监控方法、装置及计算设备
JP7193384B2 (ja) 残差特性推定モデル作成方法および残差特性推定モデル作成システム
US8428878B2 (en) Lightning strike simulation apparatus, method thereof, and program
CN109360215B (zh) 三维模型的外轮廓查找方法、装置、设备和存储介质
CN116127661A (zh) 终端接入变电站的路径确定方法、装置、设备和存储介质
CN111898687B (zh) 一种基于狄隆涅三角剖分的雷达反射率数据融合方法
CN115061222A (zh) 一种气候预测方法、计算机设备及存储介质
CN113689526A (zh) 地图中无效区域的划分方法及装置、电子设备
CN110413656B (zh) 一种基于高空间分辨率数值模式的急流轴自动识别方法
CN113704309A (zh) 图数据处理方法、装置、计算机设备和存储介质
CN114511571A (zh) 一种点云数据语义分割方法、系统及相关组件
CN113840299A (zh) 广播天线权值配置方法、装置、设备及存储介质

Legal Events

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