CN108230336B - 一种点云杆塔提取方法及装置 - Google Patents

一种点云杆塔提取方法及装置 Download PDF

Info

Publication number
CN108230336B
CN108230336B CN201711481537.4A CN201711481537A CN108230336B CN 108230336 B CN108230336 B CN 108230336B CN 201711481537 A CN201711481537 A CN 201711481537A CN 108230336 B CN108230336 B CN 108230336B
Authority
CN
China
Prior art keywords
tower
point
height
point cloud
minimum
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
Application number
CN201711481537.4A
Other languages
English (en)
Other versions
CN108230336A (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.)
State Grid Power Space Technology Co ltd
State Grid Anhui Electric Power Co Ltd
Original Assignee
Sgcc General Aviation Co ltd
State Grid Anhui Electric Power 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 Sgcc General Aviation Co ltd, State Grid Anhui Electric Power Co Ltd filed Critical Sgcc General Aviation Co ltd
Priority to CN201711481537.4A priority Critical patent/CN108230336B/zh
Publication of CN108230336A publication Critical patent/CN108230336A/zh
Application granted granted Critical
Publication of CN108230336B publication Critical patent/CN108230336B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/11Region-based segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/60Analysis of geometric attributes
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/70Determining position or orientation of objects or cameras
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10028Range image; Depth image; 3D point clouds
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30181Earth observation

Abstract

本发明公开了一种点云杆塔提取方法及装置,方法包括:根据杆塔点云的位置坐标获取杆塔半径、杆塔最小高度值以及杆塔最大高度值;利用杆塔点云的位置坐标、杆塔半径、杆塔最小高度和杆塔最大高度获得圆柱,利用圆柱内的点云确定杆塔塔臂的高度;以杆塔塔臂的高度为起点,往下对杆塔点云依次进行截取,根据截取的杆塔点云的俯视图获得最小外接四边形,根据最小外接四边形确定当前截取的杆塔点云为最后截取的杆塔点云,将当前截取的杆塔点云内点的最小z轴坐标值记录为Z2,将高度大于等于Z2的点设置为杆塔点;并将获得所有最小外接四边形的四个顶点拟合出四条空间直线;利用四条空间直线判断高度小于Z2的点云中的点是否属于杆塔点。

Description

一种点云杆塔提取方法及装置
技术领域
本申请涉及激光点云分类技术领域,特别涉及一种点云杆塔提取方法及装置。
背景技术
机载激光雷达(LIDAR)技术是近年来发展起来的一种可以直接获取地面三维信息的技术,具有传统航空摄影测量无可比拟的优势,已广泛应用于杆塔点云提取和三维重建。但是,LIDAR技术获取的点云数据除包含杆塔信息外,也包含了其他非杆塔信息。比如:地面信息、植被信息及导线等信息。要准确获取杆塔的三维信息,必须对LIDAR技术获得的数据进行处理。
目前,对点云杆塔提取没有权威的方法可以使用,很多公司在提取杆塔点云时,将离杆塔中心一定距离内的点云都归为杆塔点。这种方法有很大的缺陷,由于各个杆塔半径不一,如果杆塔半径选择过小,则杆塔外部点会丢失,若杆塔半径选择过大,将把杆塔附近点云错化为杆塔点。另外,使用该方法时,杆塔下面四个角支架很多时候分辨不清,影响美观。
发明内容
本申请实施方式的目的是提供一种点云杆塔提取方法及装置,能够解决现有技术从点云中提出杆塔的准确率低的技术问题。
为实现上述目的,本申请实施方式提供一种点云杆塔提取方法,包括:
根据杆塔点云的位置坐标获取杆塔半径、杆塔最小高度值以及杆塔最大高度值;
以所述杆塔点云的位置坐标为中心,利用所述杆塔半径、所述杆塔最小高度和所述杆塔最大高度获得第一圆柱体,利用所述第一圆柱体内的点云确定杆塔塔臂的高度;
以所述杆塔塔臂的高度为起点,往下对杆塔点云依次进行截取,根据截取的杆塔点云的俯视图获得最小外接四边形,求取所述最小外接四边形的对角线长度,并比较每个最小外接四边形的对角线长度,如果当前截取的杆塔点云的最小外接四边形的对角线长度与上一个截取的杆塔点云的最小外接四边形的对角线长度之差大于阈值,或当前截取的杆塔点云内的点的最小z轴坐标值小于Zbottom时,将当前截取的杆塔点云内点的最小z轴坐标值记录为Z2,将高度大于等于Z2的点设置为杆塔点;并将获得所有最小外接四边形的四个顶点拟合出四条空间直线;其中,Zbottom=Zmin+经验值;Zmin为杆塔最小高度值;
利用四条空间直线以及Zbottom对高度小于Z2的点云中的点判断,确定高度小于Z2的点是否属于杆塔点,将不属于杆塔点的点过滤掉,高度大于等于Z2的点以及高度小于Z2的杆塔点构成点云杆塔,提取出点云杆塔。
为实现上述目的,本申请实施方式还提供一种点云杆塔提取装置,所述装置包括:存储器和处理器,所述存储器中存储计算机程序,所述计算机程序被所述处理器执行时,实现以下功能:
根据杆塔点云的位置坐标获取杆塔半径、杆塔最小高度值以及杆塔最大高度值;
以所述杆塔点云的位置坐标为中心,利用所述杆塔半径、所述杆塔最小高度和所述杆塔最大高度获得第一圆柱体,利用所述第一圆柱体内的点云确定杆塔塔臂的高度;
以所述杆塔塔臂的高度为起点,往下对杆塔点云依次进行截取,根据截取的杆塔点云的俯视图获得最小外接四边形,求取所述最小外接四边形的对角线长度,并比较每个最小外接四边形的对角线长度,如果当前截取的杆塔点云的最小外接四边形的对角线长度与上一个截取的杆塔点云的最小外接四边形的对角线长度之差大于阈值,或当前截取的杆塔点云内的点的最小z轴坐标值小于Zbottom时,将当前截取的杆塔点云内点的最小z轴坐标值记录为Z2,将高度大于等于Z2的点设置为杆塔点;并将获得所有最小外接四边形的四个顶点拟合出四条空间直线;其中,Zbottom=Zmin+经验值;Zmin为杆塔最小高度值;
利用四条空间直线以及Zbottom对高度小于Z2的点云中的点判断,确定高度小于Z2的点是否属于杆塔点,将不属于杆塔点的点过滤掉,高度大于等于Z2的点以及高度小于Z2的杆塔点构成点云杆塔,提取出点云杆塔。
由上可见,与现有技术相比较,本技术方案从点云中提取杆塔时,无需用户输入参数即可自动化提取杆塔,并且提取杆塔准确率高。
附图说明
为了更清楚地说明本申请实施方式或现有技术中的技术方案,下面将对实施方式或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施方式,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本申请实施例提出的一种点云杆塔提取方法流程图;
图2为本申请实施例提出的一种点云杆塔提取装置的示意图;
图3为本申请实施例的点云原始示意图;
图4为本申请实施例的点云提取结果示意图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案,下面将结合本申请实施方式中的附图,对本申请实施方式中的技术方案进行清楚、完整地描述,显然,所描述的实施方式仅仅是本申请一部分实施方式,而不是全部的实施方式。基于本申请中的实施方式,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施方式,都应当属于本申请保护的范围。
本申请提供一种点云杆塔提取方法,根据图1所示步骤对点云进行处理,起到提精降噪的目的,最终从点云中获得杆塔。所述方法可以应用于具备数据处理功能的终端设备中。所述终端设备例如可以是台式电脑、笔记本电脑、平板电脑、工作站等。所述方法可以包括以下步骤:
步骤101):根据杆塔点云的位置坐标获取杆塔半径、杆塔最小高度值以及杆塔最大高度值。
在本实施例中,以杆塔点云的位置坐标为中心,以T为半径,获得第四圆柱体,根据所述杆塔最大高度选取合理的高度,该高度所在平面上的点云的点投影到杆塔点云的位置坐标所在平面上,获得投影点,将投影点到杆塔点云的位置坐标之间的距离小于阈值的投影点构成投影点集,利用所述投影点集获得投影点集外接四边形,求取所述投影点集外接四边形的对角线长度,所述投影点集外接四边形的对角线长度的一半即为杆塔半径。
在本实施例中,以杆塔点云的位置坐标为中心,以r为半径获得第二圆柱体,获取所述第二圆柱体内的点云中的点的最大z轴坐标值Zmax和最小z轴坐标值Zmin,所述第二圆柱体内的点云中的点的最小z轴坐标值Zmin为杆塔最小高度。再次以杆塔点云的位置坐标为中心,以R为半径获得第三圆柱体,获取所述第三圆柱体内的点云中高度值大于(Zmax-2.0)的点的最大z轴坐标值,该最大z轴坐标值为杆塔最大高度;其中,R远大于r,R、r均通过经验确定。
步骤102):以所述杆塔点云的位置坐标为中心,利用所述杆塔半径、所述杆塔最小高度和所述杆塔最大高度获得第一圆柱体,利用所述第一圆柱体内的点云确定杆塔塔臂的高度。
在本实施例中,根据经验选取高度经验值,以所述高度经验值为起点,向上对杆塔点云依次进行截取,以1米(并不仅仅限于该数值)厚度为例,根据截取的杆塔点云的俯视图获得外接四边形,比较每个外接四边形的边长,如果当前截取的杆塔点云的外接四边形的边长与上一个截取的杆塔点云的外接四边形的边长之差大于阈值,则上一个截取的杆塔点云内的点的最小z轴坐标值为杆塔塔臂的高度。
步骤103):以所述杆塔塔臂的高度为起点,往下对杆塔点云依次进行截取,根据截取的杆塔点云的俯视图获得最小外接四边形,求取所述最小外接四边形的对角线长度,并比较每个最小外接四边形的对角线长度,如果当前截取的杆塔点云的最小外接四边形的对角线长度与上一个截取的杆塔点云的最小外接四边形的对角线长度之差大于阈值,或当前截取的杆塔点云内的点的最小z轴坐标值小于Zbottom时,将当前截取的杆塔点云内点的最小z轴坐标值记录为Z2,把高度大于等于Z2的点设置为杆塔点;并利用所有最小外接四边形的四个顶点拟合出四条空间直线;其中,Zbottom=Zmin+经验值;Zmin为杆塔最小高度值。
步骤104):利用四条空间直线以及Zbottom对高度小于Z2的点云中的点判断,确定高度小于Z2的点是否属于杆塔点,将不属于杆塔点的点过滤掉,高度大于等于Z2的点以及高度小于Z2的杆塔点构成点云杆塔,提取出点云杆塔。
在本实施例中,首先确定高度小于Z2的点云中的每个点的高度值,将高度小于Z2的点云中的每个点所在平面与所述四条空间直线相交,获得四个交点。然后,判断高度小于Z2的点云中的每个点是否在对应的四个交点构成的四边形内;如果不在,则该点为非杆塔点;否则,当前点的高度是否大于Zbottom;如果成立,则该点为杆塔点;否则,计算当前点到对应四个交点的距离,如果距离小于阈值T,则该点为杆塔点,否则,该点为非杆塔点;其中,阈值T=4.0*(p0.z-Zmin)/(5.0)+0.3;p0.z表示当前点的高度值,Zmin为杆塔最小高度值。
在本实施例中,将非杆塔点过滤掉,保留杆塔点,由杆塔点构成点云杆塔。如图3所示,为本申请实施例的点云原始示意图。图中的点有的属于杆塔点,有的属于非杆塔点,有的属于杂质点。先过滤杂质点,然后利用本技术方案判断过滤后的点哪些属于杆塔点,哪些属于非杆塔点。最终获得图4所示的结果图。将图3和图4进行比对可知,本技术方案对杆塔下面的四个塔角分辨的比较清楚,提取杆塔准确率高。
请参阅图2,本申请还提供一种点云杆塔提取装置。所述装置包括:存储器a和处理器b,所述存储器a中存储计算机程序,所述计算机程序被所述处理器b执行时,实现以下功能:
根据杆塔点云的位置坐标获取杆塔半径、杆塔最小高度值以及杆塔最大高度值;
以所述杆塔点云的位置坐标为中心,利用所述杆塔半径、所述杆塔最小高度和所述杆塔最大高度获得第一圆柱体,利用所述第一圆柱体内的点云确定杆塔塔臂的高度;
以所述杆塔塔臂的高度为起点,往下对杆塔点云依次进行截取,根据截取的杆塔点云的俯视图获得最小外接四边形,求取所述最小外接四边形的对角线长度,并比较每个最小外接四边形的对角线长度,如果当前截取的杆塔点云的最小外接四边形的对角线长度与上一个截取的杆塔点云的最小外接四边形的对角线长度之差大于阈值,或当前截取的杆塔点云内的点的最小z轴坐标值小于Zbottom时,将当前截取的杆塔点云内点的最小z轴坐标值记录为Z2,把高度大于等于Z2的点设置为杆塔点;并利用所有最小外接四边形的四个顶点拟合出四条空间直线;其中,Zbottom=Zmin+经验值;Zmin为杆塔最小高度值;
利用四条空间直线以及Zbottom对高度小于Z2的点云中的点判断,确定高度小于Z2的点是否属于杆塔点,将不属于杆塔点的点过滤掉,高度大于等于Z2的点以及高度小于Z2的杆塔点构成点云杆塔,提取出点云杆塔。
在本实施例中,所述计算机程序被所述处理器执行时,还实现以下功能:
过滤杆塔点云上的杂质点。
在本实施例中,确定高度小于Z2的点是否属于杆塔点,所述计算机程序被所述处理器执行时,实现以下功能:
确定高度小于Z2的点云中的每个点的高度值,将高度小于Z2的点云中的每个点的所在平面与所述四条空间直线相交,获得四个交点;
判断高度小于Z2的点云中的每个点是否在对应的四个交点构成的四边形内;如果不在,则该点为非杆塔点;否则,当前点的高度是否大于Zbottom;如果成立,则该点为杆塔点;否则,计算当前点到对应四个交点的距离,如果距离小于阈值T,则该点为杆塔点,否则,该点为非杆塔点;其中,阈值T=4.0*(p0.z-Zmin)/(5.0)+0.3;p0.z表示当前点的高度值,Zmin为杆塔最小高度值。
在本实施例中,利用所述圆柱内的点云确定杆塔塔臂的高度,所述计算机程序被所述处理器执行时,实现以下功能:
根据经验选取高度经验值,以所述高度经验值为起点,向上对杆塔点云依次进行截取,根据截取的杆塔点云的俯视图获得外接四边形,比较每个外接四边形的边长,如果当前截取的杆塔点云的外接四边形的边长与上一个截取的杆塔点云的外接四边形的边长之差大于阈值,则上一个截取的杆塔点云内的点的最小z轴坐标值为杆塔塔臂的高度。
在本实施方式中,所述存储器包括但不限于随机存取存储器(Random AccessMemory,RAM)、只读存储器(Read-Only Memory,ROM)、缓存(Cache)、硬盘(Hard DiskDrive,HDD)或者存储卡(Memory Card)。
在本实施方式中,所述处理器可以按任何适当的方式实现。例如,所述处理器可以采取例如微处理器或处理器以及存储可由该(微)处理器执行的计算机可读程序代码(例如软件或固件)的计算机可读介质、逻辑门、开关、专用集成电路(Application SpecificIntegrated Circuit,ASIC)、可编程逻辑控制器和嵌入微控制器的形式等等。
本说明书实施方式提供的点云杆塔提取装置,其存储器和处理器实现的具体功能,可以与本说明书中的前述实施方式相对照解释,并能够达到前述实施方式的技术效果,这里便不再赘述。
在20世纪90年代,对于一个技术的改进可以很明显地区分是硬件上的改进(例如,对二极管、晶体管、开关等电路结构的改进)还是软件上的改进(对于方法流程的改进)。然而,随着技术的发展,当今的很多方法流程的改进已经可以视为硬件电路结构的直接改进。设计人员几乎都通过将改进的方法流程编程到硬件电路中来得到相应的硬件电路结构。因此,不能说一个方法流程的改进就不能用硬件实体模块来实现。例如,可编程逻辑器件(Programmable Logic Device,PLD)(例如现场可编程门阵列(Field Programmable GateArray,FPGA))就是这样一种集成电路,其逻辑功能由用户对器件编程来确定。由设计人员自行编程来把一个数字系统“集成”在一片PLD上,而不需要请芯片制造厂商来设计和制作专用的集成电路芯片。而且,如今,取代手工地制作集成电路芯片,这种编程也多半改用“逻辑编译器(logic compiler)”软件来实现,它与程序开发撰写时所用的软件编译器相类似,而要编译之前的原始代码也得用特定的编程语言来撰写,此称之为硬件描述语言(Hardware Description Language,HDL),而HDL也并非仅有一种,而是有许多种,如ABEL(Advanced Boolean Expression Language)、AHDL(Altera Hardware DescriptionLanguage)、Confluence、CUPL(Cornell University Programming Language)、HDCal、JHDL(Java Hardware Description Language)、Lava、Lola、MyHDL、PALASM、RHDL(RubyHardware Description Language)等,目前最普遍使用的是VHDL(Very-High-SpeedIntegrated Circuit Hardware Description Language)与Verilog2。本领域技术人员也应该清楚,只需要将方法流程用上述几种硬件描述语言稍作逻辑编程并编程到集成电路中,就可以很容易得到实现该逻辑方法流程的硬件电路。
本领域技术人员也知道,除了以纯计算机可读程序代码方式实现客户端、服务器以外,完全可以通过将方法步骤进行逻辑编程来使得客户端、服务器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种客户端、服务器可以被认为是一种硬件部件,而对其内部包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至,可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
通过以上的实施方式的描述可知,本领域的技术人员可以清楚地了解到本申请可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解,本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本申请各个实施方式或者实施方式的某些部分所述的方法。
本说明书中的各个实施方式均采用递进的方式描述,各个实施方式之间相同相似的部分互相参见即可,每个实施方式重点说明的都是与其他实施方式的不同之处。尤其,针对客户端的实施方式来说,均可以参照前述方法的实施方式的介绍对照解释。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述,例如程序模块。一般地,程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请,在这些分布式计算环境中,由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中,程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
虽然通过实施方式描绘了本申请,本领域普通技术人员知道,本申请有许多变形和变化而不脱离本申请的精神,希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。

Claims (10)

1.一种点云杆塔提取方法,其特征在于,包括:
根据杆塔点云的位置坐标获取杆塔半径、杆塔最小高度值以及杆塔最大高度值;
以所述杆塔点云的位置坐标为中心,利用所述杆塔半径、所述杆塔最小高度和所述杆塔最大高度获得第一圆柱体,利用所述第一圆柱体内的点云确定杆塔塔臂的高度;
以所述杆塔塔臂的高度为起点,往下对杆塔点云依次进行截取,根据截取的杆塔点云的俯视图获得最小外接四边形,求取所述最小外接四边形的对角线长度,并比较每个最小外接四边形的对角线长度,如果当前截取的杆塔点云的最小外接四边形的对角线长度与上一个截取的杆塔点云的最小外接四边形的对角线长度之差大于阈值,或当前截取的杆塔点云内的点的最小z轴坐标值小于Zbottom时,将当前截取的杆塔点云内点的最小z轴坐标值记录为Z2,把高度大于等于Z2的点设置为杆塔点;并利用所有最小外接四边形的四个顶点拟合出四条空间直线;其中,Zbottom=Zmin+经验值;Zmin为杆塔最小高度值;
利用四条空间直线以及Zbottom对高度小于Z2的点云中的点判断,确定高度小于Z2的点是否属于杆塔点,将不属于杆塔点的点过滤掉,高度大于等于Z2的点以及高度小于Z2的杆塔点构成点云杆塔,提取出点云杆塔。
2.如权利要求1所述的方法,其特征在于,还包括:过滤杆塔点云上的杂质点。
3.如权利要求1或2所述的方法,其特征在于,确定高度小于Z2的点是否属于杆塔点的步骤包括:
确定高度小于Z2的点云中的每个点的高度值,将高度小于Z2的点云中的每个点所在平面与所述四条空间直线相交,获得四个交点;
判断高度小于Z2的点云中的每个点是否在对应的四个交点构成的四边形内;如果不在,则当前点为非杆塔点;否则,当前点的高度是否大于Zbottom;如果成立,则当前点为杆塔点;否则,计算当前点到对应四个交点的距离,如果距离小于阈值T,则当前点为杆塔点,否则,当前点为非杆塔点;其中,阈值T=4.0*(p0.z-Zmin)/(5.0)+0.3;p0.z表示当前点的高度值,Zmin为杆塔最小高度值。
4.如权利要求1或2所述的方法,其特征在于,利用所述圆柱内的点云确定杆塔塔臂的高度的步骤包括:
根据经验选取高度经验值,以所述高度经验值为起点,向上对杆塔点云依次进行截取,根据截取的杆塔点云的俯视图获得外接四边形,比较每个外接四边形的边长,如果当前截取的杆塔点云的外接四边形的边长与上一个截取的杆塔点云的外接四边形的边长之差大于阈值,则上一个截取的杆塔点云内的点的最小z轴坐标值为杆塔塔臂的高度。
5.如权利要求1或2所述的方法,其特征在于,获取杆塔最小高度值以及杆塔最大高度值的步骤包括:
以杆塔点云的位置坐标为中心,以r为半径获得第二圆柱体,获取所述第二圆柱体内的点云中的点的最大z轴坐标值Zmax和最小z轴坐标值Zmin,所述第二圆柱体内的点云中的点的最小z轴坐标值Zmin为杆塔最小高度;
再次以杆塔点云的位置坐标为中心,以R为半径获得第三圆柱体,获取所述第三圆柱体内的点云中高度值大于(Zmax-2.0)的点的最大z轴坐标值,该最大z轴坐标值为杆塔最大高度;其中,R远大于r,R、r均通过经验确定。
6.如权利要求5所述的方法,其特征在于,获取杆塔半径的步骤包括:
以杆塔点云的位置坐标为中心,以T为半径,获得第四圆柱体,根据所述杆塔最大高度选取合理的高度,该高度所在平面上的点云的点投影到杆塔点云的位置坐标所在平面上,获得投影点,将投影点到杆塔点云的位置坐标之间的距离小于阈值的投影点构成投影点集,利用所述投影点集获得投影点集外接四边形,求取所述投影点集外接四边形的对角线长度,所述投影点集外接四边形的对角线长度的一半即为杆塔半径。
7.一种点云杆塔提取装置,其特征在于,所述装置包括:存储器和处理器,所述存储器中存储计算机程序,所述计算机程序被所述处理器执行时,实现以下功能:
根据杆塔点云的位置坐标获取杆塔半径、杆塔最小高度值以及杆塔最大高度值;
以所述杆塔点云的位置坐标为中心,利用所述杆塔半径、所述杆塔最小高度和所述杆塔最大高度获得第一圆柱体,利用所述第一圆柱体内的点云确定杆塔塔臂的高度;
以所述杆塔塔臂的高度为起点,往下对杆塔点云依次进行截取,根据截取的杆塔点云的俯视图获得最小外接四边形,求取所述最小外接四边形的对角线长度,并比较每个最小外接四边形的对角线长度,如果当前截取的杆塔点云的最小外接四边形的对角线长度与上一个截取的杆塔点云的最小外接四边形的对角线长度之差大于阈值,或当前截取的杆塔点云内的点的最小z轴坐标值小于Zbottom时,将当前截取的杆塔点云内点的最小z轴坐标值记录为Z2,把高度大于等于Z2的点设置为杆塔点;并利用所有最小外接四边形的四个顶点拟合出四条空间直线;其中,Zbottom=Zmin+经验值;Zmin为杆塔最小高度值;
利用四条空间直线以及Zbottom对高度小于Z2的点云中的点判断,确定高度小于Z2的点是否属于杆塔点,将不属于杆塔点的点过滤掉,高度大于等于Z2的点以及高度小于Z2的杆塔点构成点云杆塔,提取出点云杆塔。
8.如权利要求7所述的装置,其特征在于,所述计算机程序被所述处理器执行时,还实现以下功能:
过滤杆塔点云上的杂质点。
9.如权利要求7或8所述的装置,其特征在于,确定高度小于Z2的点是否属于杆塔点,所述计算机程序被所述处理器执行时,实现以下功能:
确定高度小于Z2的点云中的每个点的高度值,将高度小于Z2的点云中的每个点的所在平面与所述四条空间直线相交,获得四个交点;
判断高度小于Z2的点云中的每个点是否在对应的四个交点构成的四边形内;如果不在,则当前点为非杆塔点;否则,当前点的高度是否大于Zbottom;如果成立,则当前点为杆塔点;否则,计算当前点到对应四个交点的距离,如果距离小于阈值T,则当前点为杆塔点,否则,当前点为非杆塔点;其中,阈值T=4.0*(p0.z-Zmin)/(5.0)+0.3;p0.z表示当前点的高度值,Zmin为杆塔最小高度值。
10.如权利要求7或8所述的装置,其特征在于,利用所述圆柱内的点云确定杆塔塔臂的高度,所述计算机程序被所述处理器执行时,实现以下功能:
根据经验选取高度经验值,以所述高度经验值为起点,向上对杆塔点云依次进行截取,根据截取的杆塔点云的俯视图获得外接四边形,比较每个外接四边形的边长,如果当前截取的杆塔点云的外接四边形的边长与上一个截取的杆塔点云的外接四边形的边长之差大于阈值,则上一个截取的杆塔点云内的点的最小z轴坐标值为杆塔塔臂的高度。
CN201711481537.4A 2017-12-29 2017-12-29 一种点云杆塔提取方法及装置 Active CN108230336B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711481537.4A CN108230336B (zh) 2017-12-29 2017-12-29 一种点云杆塔提取方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711481537.4A CN108230336B (zh) 2017-12-29 2017-12-29 一种点云杆塔提取方法及装置

Publications (2)

Publication Number Publication Date
CN108230336A CN108230336A (zh) 2018-06-29
CN108230336B true CN108230336B (zh) 2021-09-24

Family

ID=62647294

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711481537.4A Active CN108230336B (zh) 2017-12-29 2017-12-29 一种点云杆塔提取方法及装置

Country Status (1)

Country Link
CN (1) CN108230336B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109872384B (zh) * 2018-12-29 2021-03-09 中国科学院遥感与数字地球研究所 一种基于机载lidar点云数据的杆塔自动化建模方法
CN110864646B (zh) * 2019-11-28 2021-11-19 北京百度网讯科技有限公司 可用于自动驾驶车辆的用于检测抬杆的方法和装置
CN113009452B (zh) * 2019-12-20 2024-03-19 广东科诺勘测工程有限公司 一种激光点云电力杆塔提取方法
CN111814691B (zh) * 2020-07-10 2022-01-21 广东电网有限责任公司 一种输电杆塔图像的空间扩展显示方法及装置
CN112945198B (zh) * 2021-02-02 2023-01-31 贵州电网有限责任公司 一种基于激光lidar点云的输电线路铁塔倾斜自动检测方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103412296A (zh) * 2013-06-28 2013-11-27 广东电网公司电力科学研究院 无序激光点云数据中自动提取电力塔方法
CN104732588A (zh) * 2015-03-30 2015-06-24 中国测绘科学研究院 一种基于机载激光雷达点云的电力线三维重建方法
CN106529029A (zh) * 2016-10-25 2017-03-22 北京煜邦电力技术股份有限公司 输电线路杆塔的点云数据提取方法及装置
CN106709946A (zh) * 2016-12-16 2017-05-24 武汉大学 基于LiDAR点云的多分裂导线自动提取与精细建模方法
CN106897686A (zh) * 2017-02-19 2017-06-27 北京林业大学 一种机载lidar电力巡检点云分类方法
CN107273902A (zh) * 2017-05-19 2017-10-20 中国科学院遥感与数字地球研究所 一种从机载LiDAR数据自动提取电塔点云的方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8374431B2 (en) * 2010-07-27 2013-02-12 Aerotec, Llc Method and apparatus for direct detection, location, analysis, identification, and reporting of vegetation clearance violations

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103412296A (zh) * 2013-06-28 2013-11-27 广东电网公司电力科学研究院 无序激光点云数据中自动提取电力塔方法
CN104732588A (zh) * 2015-03-30 2015-06-24 中国测绘科学研究院 一种基于机载激光雷达点云的电力线三维重建方法
CN106529029A (zh) * 2016-10-25 2017-03-22 北京煜邦电力技术股份有限公司 输电线路杆塔的点云数据提取方法及装置
CN106709946A (zh) * 2016-12-16 2017-05-24 武汉大学 基于LiDAR点云的多分裂导线自动提取与精细建模方法
CN106897686A (zh) * 2017-02-19 2017-06-27 北京林业大学 一种机载lidar电力巡检点云分类方法
CN107273902A (zh) * 2017-05-19 2017-10-20 中国科学院遥感与数字地球研究所 一种从机载LiDAR数据自动提取电塔点云的方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
机载激光雷达点云电力线三维重建方法研究;段敏燕 等;《中国博士学位论文全文数据库 信息科技辑》;20170315(第2017年第03期);全文 *
激光点云中输电线拟合与杆塔定位方法研究;游安清 等;《计算机科学》;20130430;第40卷(第4期);第298-300页 *

Also Published As

Publication number Publication date
CN108230336A (zh) 2018-06-29

Similar Documents

Publication Publication Date Title
CN108230336B (zh) 一种点云杆塔提取方法及装置
KR102062680B1 (ko) 레이저 포인트 클라우드 기반의 도시 도로 인식 방법, 장치, 저장 매체 및 기기
TWI709866B (zh) 設備型號識別方法、裝置及處理設備
CN102930277B (zh) 一种基于识别反馈的字符图像验证码识别方法
WO2021051604A1 (zh) Osd的文字区域的识别方法、装置及存储介质
CN110990516B (zh) 地图数据的处理方法、装置和服务器
KR20170137170A (ko) 텍스트 이미지 처리 방법 및 장치
CN112560862B (zh) 文本识别方法、装置及电子设备
CN109598670B (zh) 地图信息采集的内存管理方法、装置、存储介质和系统
CN111311709A (zh) 一种生成高精地图的方法及装置
US20230042968A1 (en) High-definition map creation method and device, and electronic device
CN111402161B (zh) 点云障碍物的去噪方法、装置、设备和存储介质
CN107292002B (zh) 一种数字岩心重构的方法及装置
CN110929664B (zh) 图像识别方法以及装置
CN114283343B (zh) 基于遥感卫星图像的地图更新方法、训练方法和设备
CN110135224B (zh) 一种监控视频的前景目标提取方法及系统、存储介质及终端
CN114283148A (zh) 一种道路元素提取方法、装置、存储介质及电子设备
CN113570727A (zh) 场景文件的生成方法、装置、电子设备及存储介质
CN106407941A (zh) 基于面向对象影像分析的港口围填海监测的方法及系统
CN111104908A (zh) 一种道路边沿确定方法及装置
CN110737652B (zh) 露天矿山三维数字模型的数据清洗方法、系统及存储介质
CN111783651B (zh) 路面元素识别方法、装置、电子设备和存储介质
CN117392631B (zh) 一种道路边界提取方法、装置、电子设备及存储介质
CN117274536B (zh) 一种实景三维模型重构方法及装置
CN113591847B (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
CP01 Change in the name or title of a patent holder

Address after: 102209 B218, No. 10 Zhongxing Road, Changping District Science Park, Beijing

Patentee after: State Grid Power Space Technology Co.,Ltd.

Patentee after: State Grid Anhui Electric Power Company

Address before: 102209 B218, No. 10 Zhongxing Road, Changping District Science Park, Beijing

Patentee before: SGCC GENERAL AVIATION Co.,Ltd.

Patentee before: State Grid Anhui Electric Power Company

CP01 Change in the name or title of a patent holder