CN110068279B - 一种基于点云数据的预制构件平面圆孔提取方法 - Google Patents

一种基于点云数据的预制构件平面圆孔提取方法 Download PDF

Info

Publication number
CN110068279B
CN110068279B CN201910341004.9A CN201910341004A CN110068279B CN 110068279 B CN110068279 B CN 110068279B CN 201910341004 A CN201910341004 A CN 201910341004A CN 110068279 B CN110068279 B CN 110068279B
Authority
CN
China
Prior art keywords
point
edge
circular hole
cloud data
list
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
CN201910341004.9A
Other languages
English (en)
Other versions
CN110068279A (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.)
Chongqing University Industrial Technology Research Institute
Original Assignee
Chongqing University Industrial Technology Research Institute
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 Chongqing University Industrial Technology Research Institute filed Critical Chongqing University Industrial Technology Research Institute
Priority to CN201910341004.9A priority Critical patent/CN110068279B/zh
Publication of CN110068279A publication Critical patent/CN110068279A/zh
Application granted granted Critical
Publication of CN110068279B publication Critical patent/CN110068279B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/08Measuring arrangements characterised by the use of optical techniques for measuring diameters
    • G01B11/12Measuring arrangements characterised by the use of optical techniques for measuring diameters internal diameters
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0004Industrial image inspection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/13Edge detection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/168Segmentation; Edge detection involving transform domain methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/60Analysis of geometric attributes
    • G06T7/62Analysis of geometric attributes of area, perimeter, diameter or volume
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/60Analysis of geometric attributes
    • G06T7/64Analysis of geometric attributes of convexity or concavity
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20048Transform domain processing
    • G06T2207/20061Hough transform
    • 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/30108Industrial image inspection
    • G06T2207/30164Workpiece; Machine component

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • Quality & Reliability (AREA)
  • Image Analysis (AREA)

Abstract

发明提供一种基于点云数据的预制构件平面圆孔提取方法。该方法包括输入点云数据、获取边缘点列表、拟合圆孔和修正圆孔半径等步骤。该方法使得预制构件圆孔逆向建模有效地减轻了在传统逆向建模工作中的大量人为干预,为预制构件的非接触式检测与圆孔逆向建模提供了技术支持。

Description

一种基于点云数据的预制构件平面圆孔提取方法
技术领域
本发明涉及建筑管理与信息技术领域,特别涉及一种基于点云数据的预制构件平面圆孔提取方法。
背景技术
三维激光扫描技术是上世纪九十年代中期开始出现的一项高新技术,是继GPS空间定位***之后又一项测绘技术上的突破。它通过高速激光扫描测量的方法,快速大面积高分辨率地获取被测对象表面的三维坐标以及被测物体表面反射率。可以快速、大量的采集空间点位信息,为建立物体的三维影像模型提供了一种全新的技术手段。由于其具有快速性,不接触性,实时、动态、主动性,高密度、高精度,数字化、自动化等特性,其应用很有可能会像GPS一样又一次引起测量技术的革命。利用三维激光扫描技术,可以对预制构件的进行轮廓扫描,获得预制构件的表面点云数据。高密度的点云可以对预制构件的进行非接触式检测,并通过处理点云数据完成对被测量构件高效逆向建模。
目前,技术人员在对预制构件的孔洞检测时,都需要利用卷尺或钢尺等人工测量手段,效率低下,将消耗大量的劳动力与时间。另外,根据现有的技术手段,预制构件点云数据进行BIM逆向建模时,其圆孔通常根据建筑图纸进行尺寸复原的翻模,或根据实际手工测量构件圆孔位置进行翻模,或按照点云数据的大致轮廓进行粗略的人工建模。然而,以上方法均需要通过大量的人工操作才可以确定圆孔位置与大小。
因此,为了避免在数据规模较大的情况下,人工孔洞检测与逆向建模易出现的主观性误差,以及劳动力消耗过大的问题,亟需开发一种平面圆孔提取方法。
发明内容
本发明的目的是提供一种基于点云数据的预制构件平面圆孔提取方法,以解决现有技术中存在的问题。
为实现本发明目的而采用的技术方案是这样的,一种基于点云数据的预制构件平面圆孔提取方法,包括以下步骤:
1)输入预制构件具有圆孔的平面的点云数据。
2)利用边缘提取算法对步骤1)输入的点云数据进行处理,获得包含内边缘和外边缘的边缘点列表。
3)增大最近邻点个数,再次利用边缘提取算法对步骤1)输入的点云数据进行处理,获得只含外边缘的边缘点列表。
4)将步骤2)获得的边缘点列表和步骤3)获得的边缘点列表做差,获得只包含内边缘的边缘点列表。
5)对步骤4)获得的只含内边缘的边缘点列表进行拟合,输出拟合的圆孔信息,并将圆孔信息添加到圆孔列表中。
6)对步骤5)获得的圆孔半径进行修正,并替换圆孔列表中对应的圆孔半径。
进一步,步骤2)具体包括以下步骤:
2.1)计算每个点的M个最近邻点。
2.2)以每个点为中心,将其邻域按角度八等分。判断这八个区域中是否都落入该点的最近邻点。若否,将该点标记为边缘点。
进一步,步骤3)具体包括以下步骤:
3.1)计算每个点的N个最近邻点。其中,N=5M。
3.2)以每个点为中心,将其邻域按角度八等分,判断这八个区域中是否都落入该点的最近邻点。若否,将该点标记为边缘点。
进一步,步骤5)中采用霍夫圆形检测算法求解圆的圆心坐标和半径大小。
进一步,步骤6)中,以每个圆背离圆心的径向为搜索方向,以步骤5)获得的圆孔列表中每个圆的圆心作为起点,按照搜索方向进行搜索。当满足终止条件条件时,搜索停止。停止点距圆心距离即为修正后的圆孔半径。使用修正后的圆孔半径替换步骤5)中圆孔列表中对应的圆孔半径。
本发明的技术效果是毋庸置疑的:
A.通过对预制构件的点云数据进行边缘提取,可得到内边缘点的集合;
B.有效解决了预制构件的圆孔提取问题,使得预制构件圆孔提取变得自动化,极大的减轻了大规模预制构件的圆孔逆向建模过程;
C.有效拟合与修正后得到圆孔信息的集合。
附图说明
图1为方法流程图;
图2为边缘提取算法示意图;
图3为原始预制构件点云数据;
图4为包含内边缘和外边缘的点云数据;
图5为只包含外边缘的点云数据;
图6为只包含内边缘的点云数据;
图7为只包含内边缘的点云数据拟合与修正后得到的圆孔。
具体实施方式
下面结合实施例对本发明作进一步说明,但不应该理解为本发明上述主题范围仅限于下述实施例。在不脱离本发明上述技术思想的情况下,根据本领域普通技术知识和惯用手段,做出各种替换和变更,均应包括在本发明的保护范围内。
实施例1:
参见图1,本实施例公开一种基于点云数据的预制构件平面圆孔提取方法,包括以下步骤:
1)输入预制构件存在圆孔的平面的点云数据。
2)对步骤1)输入的点云数据利用边缘提取算法进行处理,获得包含内边缘和外边缘的边缘点列表。
2.1)计算每个点一定数量的最近邻点。
2.2)如图2所示,以每个点为中心,将其邻域按角度八等分,判断这八个区域中是否都落入该点的最近邻点,若否,将该点标记为边缘点。特别地,通过控制最近邻点的数量,可以得到不同类别的边缘点。参见图4,根据具体实施实例的点云个数,在采用边缘提取时,其最近邻点的个数M设置为点云个数/100,用边缘提取算法进行处理,得到包含内边缘和外边缘的边缘点列表。
3)最近邻点个数增大至步骤2)中的5倍,对步骤1)输入的点云数据再次利用边缘提取算法进行处理,获得只含外边缘的边缘点列表。
3.1)计算每个点的N个最近邻点;其中,N=5M;
3.2)以每个点为中心,将其邻域按角度八等分,判断这八个区域中是否都落入该点的最近邻点;若否,将该点标记为边缘点。
参见图5,根据具体实施实例的点云个数,在采用边缘提取时,其最近邻点的个数N设置为点云个数/20,用边缘提取算法进行处理,得到包含外边缘的边缘点列表。
4)将步骤3)中得到的只含外边缘的边缘点列表从步骤2)中得到的包含内边缘和外边缘的边缘点列表中减去,留下如图6所示的只含内边缘的边缘点列表。
5)对步骤4)中获得的只含内边缘的边缘点列表进行拟合,输出拟合的圆孔信息,将其添加到圆孔列表中。本实施例中,采用霍夫圆形检测(Hough圆形检测)算法求解圆孔的圆心坐标和半径大小。
所述霍夫圆形检测算法原理为:给定圆形的函数关系:
(x-a)2+(y-b)2=r2
在参数空间中,x与y为已知量,a、b与r为未知量,图像空间上的一点(x,y),在参数空间中对应着一个圆锥,而图像空间的一个圆就对应着这一簇圆锥相交的一个点,这个特定点在参数空间的三维参数一定,就表示一定半径一定圆心坐标的图像空间的那个圆,将参数空间栅格化,利用图像空间上的每个点执行投票,由参数空间里累加票数的局部最大值来决定输出参数,从而检测到圆形并得到圆心与半径。
进行霍夫圆形检测时,在参数空间中以每个内边缘上的点(x,y)为圆心,遍历一定范围内的r,仅选取圆的交点参与投票以减小计算量,选出投票结果大于10的圆孔参数,将重复的参数只保留一组,得到圆孔的坐标与半径列表。
6)对步骤5)中获得的圆孔半径进行修正,替换步骤5)的圆孔列表中对应的圆孔半径。本实施例中,以每个圆背向圆心的径向为搜索方向,以步骤5)获得的圆孔列表中每个圆的圆心作为起点,按照搜索方向进行搜索,设置相应终止条件,当满足条件时,搜索停止,停止点距圆心距离即为修正后的圆孔半径。将所得圆孔半径替换步骤5)的圆孔列表中对应的圆孔半径。如图7所示,背向圆心沿径向搜索时,当以点距圆心距离为半径的1mm宽圆带内的点超过5个时,搜索停止,此时的点距圆心距离即为修正后的圆孔半径,将其替换步骤5)的圆孔列表中对应的圆孔半径。
值得说明的是,本实施例可以应用于任何具有孔洞的预制构件平面,对预制构件点云数据中边缘进行提取,进一步处理后可以获得构件上的圆孔。本实施例解决了预制构件孔洞质量检测过程中的操作繁琐,劳动成本消耗大的问题,以及对预制构件点云数据进行BIM常规逆向建模时,圆孔建模需要进行反复测量的问题,为大规模的预制构件点云数据的圆孔逆向建模提供技术支持,是切实有效的。

Claims (5)

1.一种基于点云数据的预制构件平面圆孔提取方法,其特征在于,包括以下步骤:
1)输入预制构件具有圆孔的平面的点云数据;
2)利用边缘提取算法对步骤1)输入的点云数据进行处理,获得包含内边缘和外边缘的边缘点列表;
3)增大最近邻点个数,再次利用边缘提取算法对步骤1)输入的点云数据进行处理,获得只含外边缘的边缘点列表;
4)将步骤2)获得的边缘点列表和步骤3)获得的边缘点列表做差,获得只包含内边缘的边缘点列表;
5)对步骤4)获得的只含内边缘的边缘点列表进行拟合,输出拟合的圆孔信息,并将圆孔信息添加到圆孔列表中;
6)对步骤5)获得的圆孔半径进行修正,并替换圆孔列表中对应的圆孔半径。
2.根据权利要求1所述的基于点云数据的预制构件平面圆孔提取方法,其特征在于,步骤2)具体包括以下步骤:
2.1)计算每个点的M个最近邻点;
2.2)以每个点为中心,将其邻域按角度八等分;判断这八个区域中是否都落入该点的最近邻点;若否,将该点标记为边缘点。
3.根据权利要求2所述的基于点云数据的预制构件平面圆孔提取方法,其特征在于,步骤3)具体包括以下步骤:
3.1)计算每个点的N个最近邻点;其中,N=5M;
3.2)以每个点为中心,将其邻域按角度八等分,判断这八个区域中是否都落入该点的最近邻点;若否,将该点标记为边缘点。
4.根据权利要求1所述的基于点云数据的预制构件平面圆孔提取方法,其特征在于,步骤5)中采用霍夫圆形检测算法求解圆的圆心坐标和半径大小。
5.根据权利要求1所述的基于点云数据的预制构件平面圆孔提取方法,其特征在于,步骤6)中,以每个圆背离圆心的径向为搜索方向,以步骤5)获得的圆孔列表中每个圆的圆心作为起点,按照搜索方向进行搜索;当满足终止条件条件时,搜索停止;停止点距圆心距离即为修正后的圆孔半径;使用修正后的圆孔半径替换步骤5)中圆孔列表中对应的圆孔半径。
CN201910341004.9A 2019-04-25 2019-04-25 一种基于点云数据的预制构件平面圆孔提取方法 Active CN110068279B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910341004.9A CN110068279B (zh) 2019-04-25 2019-04-25 一种基于点云数据的预制构件平面圆孔提取方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910341004.9A CN110068279B (zh) 2019-04-25 2019-04-25 一种基于点云数据的预制构件平面圆孔提取方法

Publications (2)

Publication Number Publication Date
CN110068279A CN110068279A (zh) 2019-07-30
CN110068279B true CN110068279B (zh) 2021-02-02

Family

ID=67368882

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910341004.9A Active CN110068279B (zh) 2019-04-25 2019-04-25 一种基于点云数据的预制构件平面圆孔提取方法

Country Status (1)

Country Link
CN (1) CN110068279B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114087989B (zh) * 2021-11-19 2023-09-22 江苏理工学院 汽车缸体工件定位孔圆心三维坐标测量方法及***

Family Cites Families (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101477691B (zh) * 2008-12-26 2011-04-20 武汉大学 基于边长比约束的离散点区域拓扑边界追踪方法
CN101750015B (zh) * 2009-12-11 2011-04-20 东南大学 基于数字图像技术的沙石坑土方量测量方法
US10019812B2 (en) * 2011-03-04 2018-07-10 General Electric Company Graphic overlay for measuring dimensions of features using a video inspection device
CN103226833B (zh) * 2013-05-08 2015-08-05 清华大学 一种基于三维激光雷达的点云数据分割方法
US9547901B2 (en) * 2013-11-05 2017-01-17 Samsung Electronics Co., Ltd. Method and apparatus for detecting point of interest (POI) in three-dimensional (3D) point clouds
CN103886593B (zh) * 2014-03-07 2016-08-17 华侨大学 一种基于三维点云曲面圆孔检测方法
CN104463851B (zh) * 2014-11-19 2018-05-22 哈尔滨工业大学深圳研究生院 一种基于机器人的鞋底边缘线自动跟踪方法
US20160379366A1 (en) * 2015-06-25 2016-12-29 Microsoft Technology Licensing, Llc Aligning 3d point clouds using loop closures
CN105261047B (zh) * 2015-09-08 2019-04-09 北京控制工程研究所 一种基于近距离短弧段图像的对接环圆心提取方法
CN106845535B (zh) * 2017-01-06 2019-08-02 沈阳工业大学 基于点云的典型零部件识别方法
CN106887003B (zh) * 2017-01-06 2019-06-07 沈阳工业大学 基于八邻域深度差的点云边缘提取方法
KR101858902B1 (ko) * 2017-06-26 2018-05-16 한국도로공사 컴포넌트를 활용한 점군 데이터의 객체 위치정보 추출 시스템
US10297074B2 (en) * 2017-07-18 2019-05-21 Fuscoe Engineering, Inc. Three-dimensional modeling from optical capture
CN107516098B (zh) * 2017-07-30 2021-08-10 华南理工大学 一种基于边缘曲率角的目标轮廓三维信息提取方法
CN108376266A (zh) * 2018-03-13 2018-08-07 中国电子科技集团公司第二十八研究所 基于样本边缘点内部点的单类支持向量机核参数优化方法
CN109345523B (zh) * 2018-09-21 2022-08-16 中国科学院苏州生物医学工程技术研究所 表面缺陷检测和三维建模方法
CN109410342A (zh) * 2018-09-28 2019-03-01 昆明理工大学 一种保留边界点的点云精简方法
CN109655019B (zh) * 2018-10-29 2021-02-02 北方工业大学 一种基于深度学习和三维重建的货物体积测量方法

Also Published As

Publication number Publication date
CN110068279A (zh) 2019-07-30

Similar Documents

Publication Publication Date Title
CN106969703B (zh) 用于自动成型冷却孔测量的方法和***
CN102589435B (zh) 噪声环境下激光束中心高效精确检测方法
CN110060338B (zh) 一种基于bim模型的预制构件点云识别方法
CN102441581A (zh) 基于机器视觉的型钢断面尺寸在线检测的装置及方法
CN113865508B (zh) 一种蜂窝夹芯复合材料声衬的通孔率自动化检测装置与方法
CN116740060B (zh) 基于点云几何特征提取的装配式预制构件尺寸检测方法
CN113487722A (zh) 一种基于三维激光扫描方法的混凝土构件自动化检测方法
CN115014198B (zh) 一种基于三维激光扫描的钢筋安装检测方法
CN109118476B (zh) 一种零部件边缘轮廓完整性检测方法及装置
CN108917711B (zh) 一种隧道工程三维激光扫描分段测量方法及***
CN111524154B (zh) 一种基于影像的隧道管片自动化分割方法
CN116305436A (zh) 一种基于三维激光扫描结合bim的既有桥梁监测方法
CN110068279B (zh) 一种基于点云数据的预制构件平面圆孔提取方法
CN116309118A (zh) 基于射线法的隧道点云去噪方法及***
CN115235375A (zh) 盖板类工件的多圆特征参数测量方法、检测方法和装置
CN112085708A (zh) 产品外轮廓中的直线边缘的缺陷检测方法及设备
JP5964803B2 (ja) データ処理方法及びデータ処理装置
JP2012037488A (ja) 形状検査装置及び形状検査方法
CN116385356A (zh) 一种基于激光视觉的正六边形孔特征的提取方法和***
CN115545426A (zh) 一种基于bim和三维激光扫描的预制构件生产进度分析***
CN110021027B (zh) 一种基于双目视觉的切边点计算方法
CN115511718A (zh) Pcb图像校正方法、装置、终端设备和存储介质
Wieczorowski et al. Influence of Selected Measurement Conditions on the Reliability of the Representation of Ring and Rim Features
CN112629432B (zh) 交互式孔位多角度扫描控制方法及装置
Kucharski et al. Influence of Selected Measurement Conditions on the Reliability of the Representation of Ring and Rim Features

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