CN112560573B - 一种建筑物震害信息检测提取方法 - Google Patents
一种建筑物震害信息检测提取方法 Download PDFInfo
- Publication number
- CN112560573B CN112560573B CN202011179191.4A CN202011179191A CN112560573B CN 112560573 B CN112560573 B CN 112560573B CN 202011179191 A CN202011179191 A CN 202011179191A CN 112560573 B CN112560573 B CN 112560573B
- Authority
- CN
- China
- Prior art keywords
- point
- building
- point cloud
- wall
- cloud data
- 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 45
- 230000006378 damage Effects 0.000 title claims abstract description 37
- 238000001514 detection method Methods 0.000 title claims abstract description 9
- 238000000513 principal component analysis Methods 0.000 claims abstract description 12
- 238000007619 statistical method Methods 0.000 claims abstract description 12
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 10
- 238000006243 chemical reaction Methods 0.000 claims abstract description 10
- 239000013598 vector Substances 0.000 claims abstract description 6
- 238000000034 method Methods 0.000 claims description 27
- 238000004364 calculation method Methods 0.000 claims description 12
- 230000000694 effects Effects 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 6
- 235000007164 Oryza sativa Nutrition 0.000 claims description 4
- 238000006073 displacement reaction Methods 0.000 claims description 4
- 238000001914 filtration Methods 0.000 claims description 4
- 235000009566 rice Nutrition 0.000 claims description 4
- 238000007781 pre-processing Methods 0.000 claims description 3
- 230000009467 reduction Effects 0.000 claims description 3
- 230000001131 transforming effect Effects 0.000 claims description 2
- 240000007594 Oryza sativa Species 0.000 claims 1
- 238000011156 evaluation Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 14
- 239000004567 concrete Substances 0.000 description 12
- 239000002023 wood Substances 0.000 description 11
- 239000011449 brick Substances 0.000 description 9
- 238000010276 construction Methods 0.000 description 4
- 238000011835 investigation Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 241000209094 Oryza Species 0.000 description 3
- 208000027418 Wounds and injury Diseases 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000036541 health Effects 0.000 description 2
- 208000014674 injury Diseases 0.000 description 2
- 239000011150 reinforced concrete Substances 0.000 description 2
- 239000004568 cement Substances 0.000 description 1
- 239000000919 ceramic Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 230000002787 reinforcement Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/10—Terrestrial scenes
- G06V20/176—Urban or other man-made structures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
- G06F18/232—Non-hierarchical techniques
- G06F18/2321—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0635—Risk analysis of enterprise or organisation activities
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
- G06Q50/265—Personal security, identity or safety
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/40—Extraction of image or video features
- G06V10/44—Local feature extraction by analysis of parts of the pattern, e.g. by detecting edges, contours, loops, corners, strokes or intersections; Connectivity analysis, e.g. of connected components
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- Tourism & Hospitality (AREA)
- Data Mining & Analysis (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Development Economics (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Multimedia (AREA)
- Entrepreneurship & Innovation (AREA)
- Educational Administration (AREA)
- Game Theory and Decision Science (AREA)
- Probability & Statistics with Applications (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- General Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Evolutionary Biology (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Computation (AREA)
- Computer Security & Cryptography (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种建筑物震害信息检测提取方法，采用如下步骤：S1，点云数据预处理；S2，将坐标转换后的建筑物每个墙面的原始数据进行主成分分析，计算墙面点云数据的曲率值，提取曲率正态分布期望值大于阈值的点；S3，计算墙面点云数据的临近点夹角值，统计分析震后建筑物墙体裂缝临近点云夹角特征值的分布特点，提取临近点夹角值大于阈值的点；S4，对点云数据的法向量进行统计分析，得到建筑物破坏区域所对应的判别阈值；S5，采用密度聚类算法对S4计算的初步提取结果进行去噪处理，得到建筑物墙体裂缝点云的分布区域；本发明通过三维激光扫描仪扫描震后建筑物获取三维点云数据，快速提取建筑物震害信息，能够为震后的应急救援和损失评估提供技术支撑和科学保障。
Description
技术领域
本发明涉及一种建筑物震害信息检测提取方法，属于地震灾害评估技术领域。
背景技术
地震灾害会对人们的生命财产造成严重影响，其中建筑物的损坏是造成人员伤亡的主要原因之一。地震发生后快速评估损失对于应急响应处置、救援行动和灾后重建至关重要。破坏性地震发生后，建筑物结构类型的抗震设防等级在人员伤亡中起着至关重要的作用。震后不同结构类型建筑物震害信息的提取分析研究，能够为震后的应急救援、建筑物的损坏程度评估、其他建构筑物(桥梁、隧道、大坝等)的变形特征提取及建筑物抗震加固等方面提供科学依据和技术支撑。
震后建筑物损失评估判定的方法是由地震应急抗震救灾指挥部派出的相关专家组成震后现场科考小组，通过测量和目视解读(解译)判定的方法对震后建筑物破坏情况进行评估分析。震后建筑物的损失评估工作，时间紧任务重，现场工作人员在调查的过程中用尺子量测建筑物裂缝信息，存在人为误差。建筑物墙体发生轻微倾斜变形的，在工作人员现场调查中存在不易被识别的缺点，疏漏建筑物变形信息。破坏程度较高的建筑物存在二次倒塌的风险，为现场科考工作人员的人身安全带来较大的危险。
发明内容
本发明所要解决的技术问题是提供一种利用地面三维激光扫描数据,基于对点云中每个点的局部表面曲率值的分析,进行点云震害信息识别,针对初步提取结果，结合密度聚类算法来定义提取信息程度的建筑物震害信息检测提取方法。
本发明采用如下技术方案：
本发明建筑物震害信息检测提取方法，其采用如下步骤：
S1，点云数据预处理
采用三维激光扫描仪采集建筑物表面的点云数据，所述点云数据包括：X/Y/Z空间位置信息、反射强度、RGB信息以及X/Y/Z的法线信息数据；
利用三维激光扫描仪软件realworks的功能进行点云数据去噪、滤波处理和坐标转换；取坐标转换后的建筑物每个墙面的原始数据，并将其保存为独立的一个点云数据文件，定义墙面点云数据的左下角为原点坐标X/Y/Z的起点(0.0.0)，对处理后的数据按照单位为米的规范进行转存，XYZ为点云数据的三维坐标，其中X为横坐标，Y为纵坐标，Z为垂直于XY平面的坐标；
S2，将坐标转换后的建筑物每个墙面的原始数据进行主成分分析，计算墙面点云数据的曲率值；统计分析震后建筑物墙体裂缝特征值的分布特点，并设定曲率阈值，提取曲率正态分布期望值大于阈值的点；
通过主成分分析降低具有大量相互关联变量数据的维度；缩减包括将原始数据转换为一组新的不相关的变量，这些变量的数量等于所分析数据的维数，在离散点云数据集P中选中任意一点a，其在离散点云数据集P中的位置为i，那么点a曲率Ki的表达式为：
其中x、y为点a的横、纵坐标值，x、y的单位为：米；
S3，计算墙面点云数据的临近点夹角值，统计分析震后建筑物墙体裂缝临近点云夹角特征值的分布特点，并设定夹角阈值，提取临近点夹角值大于阈值的点；
对离散点云数据集P中任意一点，搜索其半径r邻域内的所有临近点，r的单位为：米，设临近点数目为g，点Pi的g临近点对应法线分别为{ni1,ni2,ni3,……nig}，点Pi与r邻域内临近点法线夹角值α的计算公式为：
S4，对点云数据的法向量进行统计分析，得到建筑物破坏区域所对应的判别阈值，通过设定与阈值的关系，剔除临近平面的点，将计算结果与步骤S2和步骤S3的结果取交集，得到墙体裂缝和墙皮脱落区域的初步提取结果，所述交集为具有相同坐标位置的点；
S5，采用密度聚类算法对S4计算的初步提取结果进行去噪处理，降低墙面抹灰不均匀和边界效应对震害特征的影响，得到建筑物墙体裂缝点云的分布区域。
本发明步骤S5中，在密度聚类算法中将数据点分为三类：
(1)核心点，若样本Pi的半径r邻域内至少包含了指定数目MinPts个样本，即Nr(Pi)≥MinPts，则称样本点Pi为核心点；
(2)边界点，若样本Pi的r邻域内包含的样本数目小于MinPts，但是它在其他核心点的邻域内，则称样本点Pi为边界点；
(3)噪点，既不是核心点也不是边界点的点。
本发明积极效果如下：
(1)本发明通过主成分分析计算点云数据的曲率值和临近点夹角值，发挥了三维激光扫描仪(地基LiDAR)数据高精度的优势；
(2)自动化程度高，运算速度快，适合数万点云数据的加载计算。
(3)具有一定的灵活性，通过调整曲率和夹角的阈值，可以降低人为干扰因素。
(4)解决了建筑物震害信息三维激光扫描仪(地基LiDAR)数据的自动检测提取。
采用本发明方法可减少不必要的劳动，降低现场调查人员二次伤害的风险，消除人为主观干扰因素，快速、准确地评价建筑物的健康状况。
本发明通过三维激光扫描仪扫描震后建筑物获取三维点云数据，快速提取建筑物震害信息，能够为震后的应急救援和损失评估提供技术支撑和科学保障。
附图说明
图1为本发明数据处理流程图；
图2为砖木结构建筑现场图；
图3为砖木结构建筑原始点云数据；
图4为砖木结构建筑初步提取结果图；
图5为砖结构建筑现场图；
图6为砖结构建筑原始点云数据图；
图7为砖结构建筑初步提取结果图；
图8为砖混结构建筑现场图；
图9为砖混结构建筑原始点云数据图；
图10为砖混结构建筑初步提取结果图；
图11为框架结构建筑现场照片图；
图12为框架结构原始点云数据图；
图13为墙面脱落区域的信息初步提取结果图；
图14为砖木结构建筑墙体右下角倾斜变形提取结果图；
图15为砖结构建筑墙面震害信息提取结果图；
图16为砖混结构建筑墙面震害信息提取结果图；
图17为框架结构建筑震害信息提取结果图；
图18调查点及实例样点分布图。
具体实施方式
如附图1所示，本发明建筑物震害信息检测提取方法采用如下步骤实现：
S1，点云数据预处理
三维激光扫描仪采集建筑物表面的原始数据主要包括：X/Y/Z空间位置信息(1-3列)、反射强度(4列)、RGB信息(5-7列)以及X/Y/Z的法线信息(8-10列)，如下表1所示：
表1：建筑物墙面原始三维激光扫描仪(地基LiDAR)数据
利用三维激光扫描仪软件realworks的功能进行点云数据去噪、滤波处理和坐标转换；取坐标转换后的建筑物每个墙面的原始数据，并将其保存为独立的一个点云数据文件，定义墙面点云数据的左下角为原点坐标X/Y/Z的起点(0.0.0)，对处理后的数据按照单位为米的规范进行转存，XYZ为点云数据的三维坐标，其中X为横坐标，Y为纵坐标，Z为垂直于XY平面的坐标，如下表2所示：
表2坐标转换后的数据
S2，将坐标转换后的建筑物每个墙面的原始数据进行主成分分析，计算墙面点云数据的曲率值；统计分析震后建筑物墙体裂缝特征值的分布特点，并设定曲率阈值，提取曲率正态分布期望值大于阈值的点；
主成分分析(PCA)是在保留数据变异性的前提下，降低具有大量相互关联变量数据的维度；缩减包括将原始数据转换为一组新的不相关的变量，这些变量的数量等于所分析数据的维数；主成分分析可以从海量数据集中提取特征，目前已成功应用于点云处理(法向量、局部曲率和夹角计算)，用于分割和检测TLS数据中的特征或对象；在离散点云数据集P中选中任意一点a，其在离散点云数据集P中的位置为i，那么点a曲率Ki的表达式为：
其中x、y为点a的横、纵坐标值，x、y的单位为：米；
S3，计算墙面点云数据的临近点夹角值，统计分析震后建筑物墙体裂缝临近点云夹角特征值的分布特点，并设定夹角阈值，提取临近点夹角值大于阈值的点；
对离散点云数据集P中任意一点，搜索其半径r邻域内的所有临近点，r的单位为：米，设临近点数目为g，点Pi的g临近点对应法线分别为{ni1,ni2,ni3,……nig}，点Pi与r邻域内临近点法线夹角值α的计算公式为：
S4，对点云数据的法向量进行统计分析，得到建筑物破坏区域所对应的判别阈值，通过设定与阈值的关系，剔除临近平面的点，将计算结果与步骤S2和步骤S3的结果取交集，交集为具有相同坐标位置的点；得到墙体裂缝和墙皮脱落区域的初步提取结果；
地震现场的建筑物结构类型主要分为砖结构、砖木结构、砖混结构和框架结构；在曲率估计和临近点夹角计算过程中，需充分考虑建筑物单体墙面的曲率分布形态，排除建筑物墙面平整度不均匀性等人为干扰因素。
对震区建筑物单体墙面曲率统计分析，本发明中将曲率阈值设定为0.0015，当墙面曲率正态分布期望值小于0.0015时，建筑物墙面平整度较高，对临近点夹角的阈值设定为大于1.5度，该方法适用于墙体表面平整度较高，墙体裂缝位移错动量较小的建筑；当墙面曲率正态分布期望值大于0.0015时，建筑物墙面平整度较低，对临近点夹角的阈值设定为大于6度，该方法适用于墙体表面平整度较差，墙体裂缝位移错动量较大的建筑如下表3所示：
表3阈值设定对照表
通过上述方法对四种不同建筑结构类型墙体裂缝的提取如图2-13所示：
图2为砖木结构建筑现场照片，该建筑位于此次地震的Ⅶ度区(图18)，属于当地居民自建房屋，其建筑特点：墙面存在抹灰不平整的人为因素，砖墙砌体，没有抗震构造措施，屋顶为木梁和瓦片，墙面曲率正态分布期望值为0.0462；图3为砖木结构建筑原始点云数据；图4为砖木结构建筑提取结果，墙面抹灰不平整掩盖了墙面右下角细微裂缝的变形信息，在曲率计算的初步识别中不能够提取到裂缝的相关信息。
图5为砖结构建筑现场照片，该建筑位于此次地震的Ⅸ度区(图18)，属于当地居民自建二层房屋，其建筑特点：建筑成本低廉，不具备抗震构造措施，承重墙为砖砌体，没有钢筋水泥混凝土柱和梁，屋顶为预制板，墙面曲率正态分布期望值为0.0012；图6为砖结构建筑原始点云数据；图7为砖结构建筑提取结果，墙体裂缝的提取信息在分布形态、长度和宽度等方面与实际值吻合，由于受自建房屋墙面抹灰不均匀及边界效应的影响，在初步提取结果中存在一些噪点，但不影响墙体震害信息的提取。
图8为砖混结构建筑现场照片，该建筑位于此次地震的IX度区(图18)，属于当地居民的自建三层房屋，具备抗震设防措施，其建筑特点为：砖混砌体结构，在墙角有构造柱与墙面共同承重，一楼、二楼屋顶为钢筋水泥浇筑，同时配有承重梁，开间较大；其震害特征较为明显，墙面发生倾斜变形和X裂缝，属于严重破坏；墙面曲率正态分布期望值为0.0213；图9为砖混结构建筑原始点云数据；图10为砖混结构建筑初步提取结果，墙面左侧和顶部的裂缝能够完全提取，在该墙体右侧受瓷砖装饰面不平的影响，产生了一些噪点，掩盖了墙体右下角连接墙皮脱落区域的细微裂缝信息。
图11为框架结构建筑现场照片，该建筑位于此次地震的Ⅷ度区(图18)，属于钢筋混凝土框架结构，具备抗震构造措施；受地震的影响，该建筑存在墙皮脱落的现象，墙面曲率正态分布期望值为0.0017；图12为框架结构原始点云数据；图13为墙面脱落区域的信息提取，其宽度为8m，高度为12m，由于受混凝土抹灰不均匀、墙面悬挂条幅及边界效应的影响，存在一定的噪声点，但不影响墙体脱落范围信息的提取。
S5，采用密度聚类算法对S4计算的初步提取结果进行去噪处理，降低墙面抹灰不均匀和边界效应对震害特征的影响，得到建筑物墙体裂缝点云的分布区域。在密度聚类算法中将数据点分为三类：
(1)核心点，若样本Pi的半径r邻域内至少包含了指定数目(该指定数目为人工输入参数，即设定聚类点的个数，若想提取半径r邻域内≥n的点云，则该数目为n。)MinPts个样本，即Nr(Pi)≥MinPts，则称样本点Pi为核心点；
(2)边界点，若样本Pi的r邻域内包含的样本数目小于MinPts，但是它在其他核心点的邻域内，则称样本点Pi为边界点；
(3)噪点，既不是核心点也不是边界点的点；
对计算过程中出现的边界效应、人为因素造成的墙面点云噪点进行二次过滤去噪；如图14-17所示，在图中：*为噪点，o为边界点，·为核心点；
图14为砖木结构建筑提取结果，由于细微裂缝的宽度小于LiDAR点云数据的点间距，因此在震害信息提取过程中，只检测到了墙体右下角的倾斜变形，未提取到相关的裂缝信息；该墙面的震害信息提取度为80％。图15为砖结构建筑提取结果，墙体裂缝的分布形态、纹理形状等信息与实际情况吻合度较高，该墙面的震害信息提取度约为98％。图16为砖混结构建筑提取结果，墙体裂缝、墙皮脱落及不平整墙面的相关信息均被提取，该墙面的震害信息提取度约为95％。图17为框架结构建筑提取结果，墙皮脱落区域的边界提取度较高，该墙面的震害信息提取度约为95％。
实施例：
如图18所示2013年8月震区调查，使用Trimble GX200三维激光扫描仪(技术指标见表4)在Ⅶ、Ⅷ、Ⅸ度区进行震后建筑物点云数据采集(如图1)，共涉及四类建筑结构类型房屋，为砖木结构、砖结构、砖混结构、框架结构。在数据采集过程中，为提高工作效率，每个建筑物采用夹角式扫描，即同时采集建筑物的正面和侧面；为降低坐标误差，采用相对坐标的形式，即三维激光扫描仪所在点为控制点坐标。
表4：Trimble GX200技术指标
本发明利用地面三维激光扫描数据，检测和测量建筑物的损伤，基于对点云中每个点的局部表面曲率值的分析，该值由主成分分析确定，然后利用临近点夹角算法，结合实际测量值的数据特征进行点云震害信息识别，基于每个点的法向量分析，对原始点云数据去噪处理，针对初步提取结果，结合密度聚类算法来定义提取信息的程度。
(1)本发明通过主成分分析计算点云数据的曲率值和临近点夹角值，发挥了三维激光扫描仪数据高精度的优势；
(2)自动化程度高，运算速度快，适合数万点云数据的加载计算。
(3)具有一定的灵活性，通过调整曲率和夹角的阈值，可以降低人为干扰因素。
(4)解决了建筑物震害信息三维激光扫描仪(LiDAR)数据的自动检测提取。
采用本发明方法可减少不必要的劳动，降低现场调查人员二次伤害的风险，消除人为主观干扰因素，快速、准确地评价建筑物的健康状况。
最后说明的是，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换或改进等，均应包含在本发明的保护范围之内。
Claims (2)
1.一种建筑物震害信息检测提取方法，其特征在于，其采用如下步骤：
S1，点云数据预处理
采用三维激光扫描仪采集建筑物表面的点云数据，所述点云数据包括：X/Y/Z空间位置信息、反射强度、RGB信息以及X/Y/Z的法线信息数据；
利用三维激光扫描仪软件realworks的功能进行点云数据去噪、滤波处理和坐标转换；取坐标转换后的建筑物每个墙面的原始数据，并将其保存为独立的一个点云数据文件，定义墙面点云数据的左下角为原点坐标X/Y/Z的起点(0.0.0)，对处理后的数据按照单位为米的规范进行转存，XYZ为点云数据的三维坐标，其中X为横坐标，Y为纵坐标，Z为垂直于XY平面的坐标；
S2，将坐标转换后的建筑物每个墙面的原始数据进行主成分分析，计算墙面点云数据的曲率值；统计分析震后建筑物墙体裂缝特征值的分布特点，并设定曲率阈值，提取曲率正态分布期望值大于阈值的点；
通过主成分分析降低具有大量相互关联变量数据的维度；缩减包括将原始数据转换为一组新的不相关的变量，这些变量的数量等于所分析数据的维数，在离散点云数据集P中选中任意一点a，其在离散点云数据集P中的位置为i，那么点a曲率Ki的表达式为：
其中x、y为点a的横、纵坐标值，x、y的单位为：米；
S3，计算墙面点云数据的临近点夹角值，统计分析震后建筑物墙体裂缝临近点云夹角特征值的分布特点，并设定夹角阈值，提取临近点夹角值大于阈值的点；
对离散点云数据集P中任意一点，搜索其半径r邻域内的所有临近点，r的单位为：米，设临近点数目为g，点Pi的g临近点对应法线分别为{ni1,ni2,ni3,……nig}，点Pi与r邻域内临近点法线夹角值α的计算公式为：
S4，对点云数据的法向量进行统计分析，得到建筑物破坏区域所对应的判别阈值，通过设定与阈值的关系，剔除临近平面的点，将计算结果与步骤S2和步骤S3的结果取交集，得到墙体裂缝和墙皮脱落区域的初步提取结果，所述交集为具有相同坐标位置的点；
对震区建筑物单体墙面曲率统计分析，将曲率阈值设定为0.0015，当墙面曲率正态分布期望值小于0.0015时，建筑物墙面平整度较高，对临近点夹角的阈值设定为大于1.5度，其适用于墙体表面平整度较高，墙体裂缝位移错动量较小的建筑；当墙面曲率正态分布期望值大于0.0015时，建筑物墙面平整度较低，对临近点夹角的阈值设定为大于6度，其适用于墙体表面平整度较差，墙体裂缝位移错动量较大的建筑；
S5，采用密度聚类算法对S4计算的初步提取结果进行去噪处理，降低墙面抹灰不均匀和边界效应对震害特征的影响，得到建筑物墙体裂缝点云的分布区域。
2.根据权利要求1所述的一种建筑物震害信息检测提取方法，其特征在于，步骤S5中，在密度聚类算法中将数据点分为三类：
(1)核心点，若样本Pi的半径r邻域内至少包含了指定数目MinPts个样本，即Nr(Pi)≥MinPts，则称样本点Pi为核心点；
(2)边界点，若样本Pi的r邻域内包含的样本数目小于MinPts，但是它在其他核心点的邻域内，则称样本点Pi为边界点；
(3)噪点，既不是核心点也不是边界点的点。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011179191.4A CN112560573B (zh) | 2020-10-29 | 2020-10-29 | 一种建筑物震害信息检测提取方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011179191.4A CN112560573B (zh) | 2020-10-29 | 2020-10-29 | 一种建筑物震害信息检测提取方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112560573A CN112560573A (zh) | 2021-03-26 |
| CN112560573B true CN112560573B (zh) | 2023-03-28 |
Family
ID=75042684
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011179191.4A Active CN112560573B (zh) | 2020-10-29 | 2020-10-29 | 一种建筑物震害信息检测提取方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112560573B (zh) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114596313B (zh) * | 2022-05-09 | 2022-10-04 | 深圳大学 | 一种基于室内点云的建筑物构件损伤检测方法及相关设备 |
| CN116310791B (zh) * | 2023-01-19 | 2023-09-05 | 中国地震台网中心 | 基于建筑物震害检测的极灾区快速判定方法及电子设备 |
| CN116030103B (zh) * | 2023-03-07 | 2024-02-27 | 北京百度网讯科技有限公司 | 确定砌筑质量的方法、装置、设备和介质 |
| CN117557971B (zh) * | 2024-01-12 | 2024-03-26 | 湖北神龙工程测试技术有限公司 | 基于物联网的危房监测方法 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104634785A (zh) * | 2013-11-08 | 2015-05-20 | 中冶建筑研究总院有限公司 | 一种全新的混凝土平面裂缝预测评估和图形化输出方法 |
| CN104899883A (zh) * | 2015-05-29 | 2015-09-09 | 北京航空航天大学 | 一种深度图像场景的室内对象立方体检测方法 |
| CN106780509A (zh) * | 2016-12-01 | 2017-05-31 | 山东交通学院 | 融合多维特征的建筑物点云层次聚类分割方法 |
| CN107644452A (zh) * | 2017-09-15 | 2018-01-30 | 武汉大学 | 机载LiDAR点云屋顶面片分割方法及系统 |
| CN109509256A (zh) * | 2018-06-21 | 2019-03-22 | 华南理工大学 | 基于激光雷达的建筑结构自动测量及3d模型生成方法 |
| CN110390683A (zh) * | 2018-04-17 | 2019-10-29 | 河海大学 | 一种基于密集点云的古城墙三维裂缝检测方法 |
-
2020
- 2020-10-29 CN CN202011179191.4A patent/CN112560573B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104634785A (zh) * | 2013-11-08 | 2015-05-20 | 中冶建筑研究总院有限公司 | 一种全新的混凝土平面裂缝预测评估和图形化输出方法 |
| CN104899883A (zh) * | 2015-05-29 | 2015-09-09 | 北京航空航天大学 | 一种深度图像场景的室内对象立方体检测方法 |
| CN106780509A (zh) * | 2016-12-01 | 2017-05-31 | 山东交通学院 | 融合多维特征的建筑物点云层次聚类分割方法 |
| CN107644452A (zh) * | 2017-09-15 | 2018-01-30 | 武汉大学 | 机载LiDAR点云屋顶面片分割方法及系统 |
| CN110390683A (zh) * | 2018-04-17 | 2019-10-29 | 河海大学 | 一种基于密集点云的古城墙三维裂缝检测方法 |
| CN109509256A (zh) * | 2018-06-21 | 2019-03-22 | 华南理工大学 | 基于激光雷达的建筑结构自动测量及3d模型生成方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112560573A (zh) | 2021-03-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112560573B (zh) | 一种建筑物震害信息检测提取方法 | |
| Korumaz et al. | An integrated Terrestrial Laser Scanner (TLS), Deviation Analysis (DA) and Finite Element (FE) approach for health assessment of historical structures. A minaret case study | |
| CN111561884B (zh) | 一种用于预制混凝土叠合板表面粗糙度检测的方法 | |
| CN108824816B (zh) | 一种高空大跨网架滑移定位安装及监测方法 | |
| CN108445188B (zh) | 基于中智区间函数的岩体结构面粗糙度系数尺寸效应下边坡稳定性表达方法 | |
| CN115147408B (zh) | 危房智能检测鉴定方法 | |
| CN108154498A (zh) | 一种裂缝缺陷检测系统及其实现方法 | |
| CN111022119B (zh) | 一种基于模式识别的盾构隧道结构抗震安全评价方法 | |
| Karataş et al. | Analytical Documentation of Stone Material Deteriorations on Facades with Terrestrial Laser Scanning and Photogrammetric Methods: Case Study of Şanlıurfa Kışla Mosque | |
| CN111507972A (zh) | 结合卷积神经网络和支持向量机的隧道表面缺陷检测方法 | |
| Yang et al. | A Model Study of Building Seismic Damage Information Extraction and Analysis on Ground‐Based LiDAR Data | |
| Pöchtrager et al. | Automated reconstruction of historic roof structures from point clouds–Development and examples | |
| Jiao et al. | Nondestructive inspection of a brick–timber structure in a modern architectural heritage building: Lecture hall of the Anyuan Miners' Club, China | |
| Olivito et al. | Inventory and monitoring of historical cultural heritage buildings on a territorial scale: a preliminary study of structural health monitoring based on the CARTIS approach | |
| Hermida et al. | Two-dimensional models of variable inertia from LiDAR data for structural analysis of timber trusses | |
| Yuan et al. | Automatic evaluation of rebar spacing and quality using LiDAR data: Field application for bridge structural assessment | |
| Truong-Hong et al. | Documentation of bridges by terrestrial laser scanner | |
| CN111255220B (zh) | 超高层建筑施工钢平台模架装备的安全监控方法 | |
| Mirzaei et al. | Automatic compliance inspection and monitoring of building structural members using multi-temporal point clouds | |
| Olivito et al. | SHM systems applied to the built heritage inventory at the territorial scale. A preliminary study based on CARTIS approach | |
| Shin | A Study on Predicting the Movement of Columns in Hanok Architecture Using the UMAP and DBSCAN Algorithms | |
| Yang et al. | Three-dimensional point cloud analysis for building seismic damage information | |
| Baraccani et al. | An assessment of the structural behaviour of the Garisenda tower in bologna through finite element modelling and structural health monitoring | |
| Abraldes et al. | Use of polar coordinates for improving the measurement of resistant cross-sections of existing timber elements combining laser scanner and drilling resistance tests | |
| CN111380475A (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 |