CN110060256B - 一种基于机载LiDAR点云的杆塔提取方法 - Google Patents
一种基于机载LiDAR点云的杆塔提取方法 Download PDFInfo
- Publication number
- CN110060256B CN110060256B CN201910176449.6A CN201910176449A CN110060256B CN 110060256 B CN110060256 B CN 110060256B CN 201910176449 A CN201910176449 A CN 201910176449A CN 110060256 B CN110060256 B CN 110060256B
- Authority
- CN
- China
- Prior art keywords
- point cloud
- point
- pole
- method based
- extraction method
- 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.)
- Active
Links
- 238000000605 extraction Methods 0.000 title claims abstract description 15
- 230000011218 segmentation Effects 0.000 claims abstract description 13
- 230000008859 change Effects 0.000 claims description 9
- 238000000034 method Methods 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000012546 transfer Methods 0.000 claims description 6
- 239000003550 marker Substances 0.000 claims description 4
- 238000012549 training Methods 0.000 claims description 4
- 238000007689 inspection Methods 0.000 description 11
- 230000005540 biological transmission Effects 0.000 description 9
- 230000007547 defect Effects 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T15/00—3D [Three Dimensional] image rendering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10028—Range image; Depth image; 3D point clouds
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20081—Training; Learning
-
- 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
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Graphics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及点云语义分割技术领域,具体涉及一种基于机载LiDAR点云的杆塔提取技术,本发明首先,输入点云预处理,对预处理后的点云数据按给定的格网尺寸划分格网,计算出杆塔点云和其他点云的坡度特征和高差特征,作为权重,利用最小割算法将点云分割成小点云,再使用泰森多边形划分出起小点云之间的边缘,之后,将小点云集作为数据集,使用条件随机场,将小点云集进行分类,本发明相对于PointNet,PointCNN,本算法可以处理大规模点云,相对于传统的点云提取杆塔算法,具有更高的鲁棒性。
Description
技术领域
本发明涉及点云语义分割技术领域,具体涉及一种基于机载LiDAR点云的杆塔提取方法。
背景技术
高压输电线路是电力***的重要组成部分,具有线路距离长、覆盖范围大、安全可靠性要求高等特点。我国电网输电线路巡检目前主要通过运行人员手持仪器或肉眼来观测设备缺陷,人工巡检方式周期长且作业强度大,难以适应现代化电网快速发展和安全运行需要。
无人机技术的发展为输电线路巡检提供了全新的移动平台。利用无人机搭载巡检设备进行输电线路巡检,不受地理条件的限制,且巡检效率高,每小时可达几km至几十km。将无人机应用于输电线路巡检,可以弥补人工巡检的不足,实现对输电线路的快速、安全巡检,具有重要的研究价值和良好的应用前景。
基于无人机载激光雷达的***能以离散采样点云形式,快速、准确采集输电线路走廊电力设备与走廊地物的真三维坐标,为输电线路巡检工作提供了一种新的解决方案,对于丘陵或山区的架空线路来说意义尤为重大,因此正被逐渐应用到输电线路运行维护中,而杆塔的提取对高压输电线路的位置定位具有重要意义。
发明内容
针对现有技术的不足,本发明公开了一种基于机载LiDAR点云的杆塔提取方法,用于解决现阶段,一些点云分割的算法陆续出来,PointNet,PointCNN等神经网络分割点云都取得较好的效果,但不能对大规模的点云进行处理。传统的根据杆塔所具有的特征提取杆塔算法的鲁棒性不是很强的问题。
本发明通过以下技术方案予以实现:
一种基于机载LiDAR点云的杆塔提取方法,其特征在于,包括以下步骤:
S1输入点云预处理;
S2对预处理后的点云数据按给定的格网尺寸划分格网;
S3计算出杆塔点云和其他点云的坡度特征和高差特征并作为权重;
S4利用最小割算法将点云分割成小点云,再使用泰森多边形划分出起小点云之间的边缘;
S5将小点云集作为数据集,使用条件随机场,将小点云集进行分类。
优选的,所述坡度特征中,首先对预处理后的点云数据按给定的格网尺寸划分格网,遍历点云,计算点云所在格网的位置。计算出格网中的最低程HL和最高程HT;高差特征为H=HT-HL。
优选的,所述高差特征中,坡度特征需要反映高程沿任意方向的变化程度,坡度越大的区域高程变化越剧烈;
坡度的计算公式为:
式中gx和gy分别为X方向和Y方向上高程变化率。
优选的,所述gx和gy计算公式为:
式中的ki,j为3×3窗口的中心点的最高高程值,其余为其八邻域点的最高高程;d为格网的宽度。
优选的,所述最小割算法中,最小割分割点云,将每个点云体素作为一个节点,有向图G=(V,E)中,有源点S和汇点T,每条弧有一个非负容量wi(i=1,2,3...n),V为节点集合,E为边缘集合,我们设计n=10,每个体素为一个节点,wi(i=1,2,3...n)为两个节点之间特征,定义:
wi=0.5*(si+Hi)(i=1,2,3...n)
使用分割的依据的函数:
[.]“艾佛森括号”,是一种用方括号记号,如果方括号内的条件满足则为1,不满足则为0,W是w的集合,C是特征值个数的集合。
优选的,将小点云集视为空间中随机点,随机取一个随机点作为一个顶点,绕着该顶点,构成三角网,三角网要求为三角形三个顶点都在一个圆上,且圆内没有其他点,取三角网中,每个圆的圆心,一次连接形成泰森多边形,循环小点云集,形成泰森多边形划分点云的边缘。
优选的,包括条件随机场CRF分类点云集,其条件随机场CRF:根据给定的输出标识序列Y和观察序列X,经过训练后,得出符合模型的条件概率P(Y|X),标记序列Y为杆塔、电线、地面和其他,分别由t,l,g,e表示,观察序列X为分割出来的点云集合。
优选的,在观察序列X其中i位置的标记序列Y的概率可以定义为:
其中tj(yi-1,yi,x,i)是转移函数,sp(yi,x,i)是状态函数,和βp分别是转移函数和状态函数的权重。
本发明的有益效果为:
1、相对于PointNet,PointCNN,本算法可以处理大规模点云。
2、相对于传统的点云提取杆塔算法,具有更高的鲁棒性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明的方法步骤原理图;
图2是本发明的最小割分割点云原理图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例1
本实施例公开如图1所示的一种基于机载LiDAR点云的杆塔提取方法,包括以下步骤:
S1输入点云预处理;
S2对预处理后的点云数据按给定的格网尺寸划分格网;
S3计算出杆塔点云和其他点云的坡度特征和高差特征并作为权重;
S4利用最小割算法将点云分割成小点云,再使用泰森多边形划分出起小点云之间的边缘;
S5将小点云集作为数据集,使用条件随机场,将小点云集进行分类。
所述坡度特征中,首先对预处理后的点云数据按给定的格网尺寸划分格网,遍历点云,计算点云所在格网的位置。计算出格网中的最低程HL和最高程HT;高差特征为H=HT-HL。
所述高差特征中,坡度特征需要反映高程沿任意方向的变化程度,坡度越大的区域高程变化越剧烈;
坡度的计算公式为:
式中gx和gy分别为X方向和Y方向上高程变化率。
所述gx和gy计算公式为:
式中的ki,j为3×3窗口的中心点的最高高程值,其余为其八邻域点的最高高程;d为格网的宽度。
所述最小割算法中,如图2所示的最小割分割点云,将每个点云体素作为一个节点,有向图G=(V,E)中,有源点S和汇点T,每条弧有一个非负容量wi(i=1,2,3...n),V为节点集合,E为边缘集合,我们设计n=10,每个体素为一个节点,wi(i=1,2,3...n)为两个节点之间特征,定义:
wi=0.5*(si+Hi)(i=1,2,3...n)
使用分割的依据的函数:
[.]“艾佛森括号”,是一种用方括号记号,如果方括号内的条件满足则为1,不满足则为0,W是w的集合,C是特征值个数的集合。
将小点云集视为空间中随机点,随机取一个随机点作为一个顶点,绕着该顶点,构成三角网,三角网要求为三角形三个顶点都在一个圆上,且圆内没有其他点,取三角网中,每个圆的圆心,一次连接形成泰森多边形,循环小点云集,形成泰森多边形划分点云的边缘。
包括条件随机场CRF分类点云集,其条件随机场CRF:根据给定的输出标识序列Y和观察序列X,经过训练后,得出符合模型的条件概率P(Y|X),标记序列Y为杆塔、电线、地面和其他,分别由t,l,g,e表示,观察序列X为分割出来的点云集合。
在观察序列X其中i位置的标记序列Y的概率可以定义为:
其中tj(yi-1,yi,x,i)是转移函数,sp(yi,x,i)是状态函数,和βp分别是转移函数和状态函数的权重。
例如
表示第i个观察值为“杆塔点云”时,相对的标记分别是g,t,其状态和特征函数如下:
条件概率定义如下:
其中
Z(x)是规范化因子,求和时在所有可能的输出序列上进行的。j,k=1,2,3...n,n=标记序列Y的元素个数;由此,进行条件随机场的训练。
以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (4)
1.一种基于机载LiDAR点云的杆塔提取方法,其特征在于,包括以下步骤:
S1输入点云预处理;
S2对预处理后的点云数据按给定的格网尺寸划分格网;
S3计算出杆塔点云和其他点云的坡度特征和高差特征并作为权重;
S4利用最小割算法将点云分割成小点云,再使用泰森多边形划分出小点云之间的边缘;
S5将小点云集作为数据集,使用条件随机场,将小点云集进行分类;
所述最小割算法中,最小割分割点云,将每个点云体素作为一个节点,有向图G=(V,E)中,有源点S和汇点T,每条弧有一个非负容量wi,i=1,2,3...n,V为节点集合,E为边缘集合,设计n=10,每个体素为一个节点,wi,i=1,2,3...n,为两个节点之间特征,定义:
wi=0.5*(si+Hi),i=1,2,3...n
使用分割的依据的函数:
[.]为艾佛森括号,是一种用方括号记号,如果方括号内的条件满足则为1,不满足则为0,W是w的集合,C是特征值个数的集合;
将小点云集视为空间中随机点,随机取一个随机点作为一个顶点,绕着该顶点,构成三角网,三角网要求为三角形三个顶点都在一个圆上,且圆内没有其他点,取三角网中,每个圆的圆心,一次连接形成泰森多边形,循环小点云集,形成泰森多边形划分点云的边缘;
所述的基于机载LiDAR点云的杆塔提取方法还包括条件随机场CRF分类点云集,其条件随机场CRF:根据给定的输出标识序列Y和观察序列X,经过训练后,得出符合模型的条件概率P(Y|X),标记序列Y为杆塔、电线、地面和其他,分别由t,l,g,e表示,观察序列X为分割出来的点云集合;
在观察序列X中的i位置发生标记序列Y的概率为:
其中tb(ya-1,ya,xa)是转移函数,sp(ya,xa)是状态函数,αa和βp分别是转移函数和状态函数的权重。
2.根据权利要求1所述的基于机载LiDAR点云的杆塔提取方法,其特征在于,所述高差特征中,首先对预处理后的点云数据按给定的格网尺寸划分格网,遍历点云,计算点云所在格网的位置;计算出格网中的最低程HL和最高程HT;高差特征为H=HT-HL。
3.根据权利要求1所述的基于机载LiDAR点云的杆塔提取方法,其特征在于,所述坡度特征中,坡度特征需要反映高程沿任意方向的变化程度,坡度越大的区域高程变化越剧烈;
坡度的计算公式为:
式中gx和gy分别为X方向和Y方向上高程变化率。
4.根据权利要求3所述的基于机载LiDAR点云的杆塔提取方法,其特征在于,所述gx和gy计算公式为:
式中的ki,j为3×3窗口的中心点的最高高程值,其余为其八邻域点的最高高程;d为格网的宽度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910176449.6A CN110060256B (zh) | 2019-03-08 | 2019-03-08 | 一种基于机载LiDAR点云的杆塔提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910176449.6A CN110060256B (zh) | 2019-03-08 | 2019-03-08 | 一种基于机载LiDAR点云的杆塔提取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110060256A CN110060256A (zh) | 2019-07-26 |
CN110060256B true CN110060256B (zh) | 2023-07-25 |
Family
ID=67316774
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910176449.6A Active CN110060256B (zh) | 2019-03-08 | 2019-03-08 | 一种基于机载LiDAR点云的杆塔提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110060256B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110807774B (zh) * | 2019-09-30 | 2022-07-12 | 九天创新(广东)智能科技有限公司 | 一种点云分类与语义分割方法 |
CN111242799B (zh) * | 2019-12-10 | 2024-01-16 | 国网电力空间技术有限公司 | 基于机载LiDAR点云的高压线塔中心坐标提取编号方法及介质 |
CN111260660B (zh) * | 2020-01-14 | 2023-04-11 | 电子科技大学 | 一种基于元学习的3d点云语义分割迁移方法 |
CN112946689A (zh) * | 2021-03-08 | 2021-06-11 | 苏州岭纬智能科技有限公司 | 一体化激光雷达***及其检测方法 |
CN113554595B (zh) * | 2021-06-21 | 2022-11-25 | 国网山西省电力公司电力科学研究院 | 无人机激光雷达点云杆塔塔头形变检测装置及方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108765446A (zh) * | 2018-05-30 | 2018-11-06 | 湖北大学 | 基于随机场和随机森林的电力线点云分割方法及*** |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130096886A1 (en) * | 2010-03-31 | 2013-04-18 | Borys Vorobyov | System and Method for Extracting Features from Data Having Spatial Coordinates |
CN102254093B (zh) * | 2011-06-24 | 2014-06-18 | 天津大学 | 基于泰森多边形的连通域统计相关算法 |
GB2532948B (en) * | 2014-12-02 | 2021-04-14 | Vivo Mobile Communication Co Ltd | Object Recognition in a 3D scene |
CN106022259B (zh) * | 2016-05-20 | 2019-04-12 | 江苏得得空间信息科技有限公司 | 一种基于激光点云三维特征描述模型的山区道路提取方法 |
CN106526610B (zh) * | 2016-11-04 | 2019-04-09 | 广东电网有限责任公司电力科学研究院 | 一种基于无人机激光点云的电塔自动定位方法及装置 |
CN107025685B (zh) * | 2017-04-11 | 2020-03-17 | 南京林业大学 | 拓扑感知下的机载建筑屋顶点云建模方法 |
CN107832849B (zh) * | 2017-11-01 | 2021-03-12 | 广东电网有限责任公司电力科学研究院 | 一种基于知识库的电力线廊道三维信息提取方法及装置 |
CN108399424B (zh) * | 2018-02-06 | 2022-01-18 | 深圳市建设综合勘察设计院有限公司 | 一种点云分类方法、智能终端及存储介质 |
-
2019
- 2019-03-08 CN CN201910176449.6A patent/CN110060256B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108765446A (zh) * | 2018-05-30 | 2018-11-06 | 湖北大学 | 基于随机场和随机森林的电力线点云分割方法及*** |
Also Published As
Publication number | Publication date |
---|---|
CN110060256A (zh) | 2019-07-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110060256B (zh) | 一种基于机载LiDAR点云的杆塔提取方法 | |
CN106970375B (zh) | 一种机载激光雷达点云中自动提取建筑物信息的方法 | |
CN105023014B (zh) | 一种无人机巡检输电线路图像内的杆塔目标提取方法 | |
CN112034878B (zh) | 基于无人机的电力巡检方法、装置和电子设备 | |
CN107767389A (zh) | 一种基于机载激光点云数据的输电线自动提取方法 | |
CN102879788B (zh) | 一种基于机载LiDAR数据的输电线路的电力线提取方法 | |
CN112435208B (zh) | 显著性边缘引导的绝缘子区域检测方法 | |
CN105045274B (zh) | 一种用于无人机巡检航迹规划的智能杆塔连通图构建方法 | |
CN112883878B (zh) | 一种基于三维格网的变电站场景下的点云自动分类方法 | |
CN111861933B (zh) | 基于空间划分的点云去噪方法及装置 | |
CN109086833A (zh) | 一种基于激光点云雷达数据的输电线路危险点计算方法 | |
CN112767352B (zh) | 面向机载LiDAR电力巡检的高压线塔自动精细化提取方法 | |
CN115586557B (zh) | 一种基于路网数据的车辆行驶轨迹纠偏方法及装置 | |
CN117132915B (zh) | 一种基于点云自动分类的输电线路树障隐患分析方法 | |
CN111062949B (zh) | 一种基于机载激光雷达的电力线点云提取方法 | |
CN115326075A (zh) | 一种基于无人机实现风场全域自动化巡检的路径规划方法 | |
CN115825984A (zh) | 一种基于激光点云的输电线路巡检***及巡检方法 | |
Zhang et al. | Defect detection of anti-vibration hammer based on improved faster R-CNN | |
CN114241158A (zh) | 基于风机点云数据的部件识别方法、装置、设备及介质 | |
CN104049284A (zh) | 一种基于粒子群遗传混合算法的雷电定位方法 | |
CN115423810B (zh) | 一种风力发电机组叶片覆冰形态分析方法 | |
CN110009035B (zh) | 一种基于图像匹配的测风站群空间聚类方法 | |
CN116520881A (zh) | 基于激光点云的无人机连续塔基三维巡检路径规划方法 | |
CN116452604A (zh) | 一种复杂变电站场景分割方法、设备及存储介质 | |
CN114357843B (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 |