CN113804115B - 一种基于聚类的混凝土浇筑坯层厚度监控方法 - Google Patents

一种基于聚类的混凝土浇筑坯层厚度监控方法 Download PDF

Info

Publication number
CN113804115B
CN113804115B CN202110912492.1A CN202110912492A CN113804115B CN 113804115 B CN113804115 B CN 113804115B CN 202110912492 A CN202110912492 A CN 202110912492A CN 113804115 B CN113804115 B CN 113804115B
Authority
CN
China
Prior art keywords
point
blank layer
real
thickness
layer
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
CN202110912492.1A
Other languages
English (en)
Other versions
CN113804115A (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.)
Tianjin University
PowerChina Huadong Engineering Corp Ltd
Original Assignee
Tianjin University
PowerChina Huadong Engineering Corp 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 Tianjin University, PowerChina Huadong Engineering Corp Ltd filed Critical Tianjin University
Priority to CN202110912492.1A priority Critical patent/CN113804115B/zh
Publication of CN113804115A publication Critical patent/CN113804115A/zh
Application granted granted Critical
Publication of CN113804115B publication Critical patent/CN113804115B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • 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
    • G06F18/23213—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions with fixed number of clusters, e.g. K-means clustering
    • 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
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00—Energy generation through renewable energy sources
    • Y02E10/20—Hydro energy

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Artificial Intelligence (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Biology (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Probability & Statistics with Applications (AREA)
  • On-Site Construction Work That Accompanies The Preparation And Application Of Concrete (AREA)

Abstract

本发明公开了一种基于聚类的混凝土浇筑坯层厚度监控方法，采用聚类算法实时监测混凝土浇筑坯层厚度，结合增强现实技术实时可视化展示并预警。本发明通过监控设备获取坯层表面点云数据；通过K‑means聚类算法划分历史坯层点云数据，得到历史坯层点云区域划分；通过k‑近邻算法实时计算区域混凝土浇筑坯层厚度；通过增强现实技术可视化混凝土实际浇筑坯层厚度，实时评价坯层浇筑厚度质量，并对坯层厚度不合格区域进行预警。本发明能够实时监测评价仓面局部区域浇筑坯层厚度，并可视化标记不合格区域，实现混凝土浇筑坯层厚度精准控制。

Description

一种基于聚类的混凝土浇筑坯层厚度监控方法
技术领域
本发明涉及水利水电工程混凝土浇筑技术，特别涉及一种基于聚类的混凝土浇筑坯层厚度监控方法。
背景技术
混凝土坝多采用分层浇筑和间歇上升的方法进行大坝的浇筑施工。在坝块浇筑过程中考虑到混凝土水化放热和振捣施工作业，一般在浇筑坝块时采用分坯层浇筑和坯层间间歇控制坯层混凝土水化放热量。已有研究发现浇筑坯层厚度对坝体温度场有较大的影响，若水化热发展较快，坯层厚度减小，坯层间歇增大，则坝体受坯层水化放热的影响，其最高温度会有所提高。因此，为减小浇筑坯层划分对坝体温度应力场的影响，在坝块浇筑时需要严格控制浇筑坯层的厚度。通常在仓面浇筑设计时会合理划分坯层浇筑厚度来指导施工现场的混凝土浇筑。但现场影响混凝土浇筑施工因素较多。根据混凝土下料量多少，平仓机驾驶人员凭经验进行平仓作业，坯层厚度受主观因素影响，不确定性较大，难以保证浇筑坯层厚度满足设计标准要求，对坝体安全存在不利影响。
随着科技的进步和水利工程事业的发展，监控技术广泛应用于水利工程施工领域中，监控施工各项指标是否满足设计要求(包括振捣、温度、应力等)，防止发生较大的失误并及时对施工过程的问题进行纠偏，保证工程建设安全的可靠性。现阶段监控技术已能实现仓面机械作业的实时监控，实现了施工机械作业轨迹监控、混凝土振捣质量信息监控以及温度监控等，并能满足监控数据误差在可接受范围内。但未能实现仓面坯层厚度的有效监控并及时反馈施工现场人员，未能有效控制混凝土浇筑坯层厚度并及时纠偏。
因此在仓面混凝土浇筑监控过程中，通过浇筑坯层厚度自动监控的方法，自动监测浇筑坯层厚度，结合仓面设计标准，并通过增强现实技术将分析结果与预警信息真实展示给平仓机驾驶员，实现坯层厚度实时分析、判断和预警，帮助驾驶员针对仓面区域浇筑厚度问题实行有效纠偏操作。
发明内容
本发明的目的是克服现有技术中的不足，提供一种基于聚类的混凝土坯层厚度监控方法，实时监测坝块坯层的浇筑厚度，并将浇筑信息可视化展示，实现浇筑坯层厚度实时监测和预警。
本发明所采用的技术方案是：一种基于聚类的混凝土浇筑坯层厚度监控方法，其特征在于，包括以下步骤：
步骤1，根据仓面坯层设计厚度(h0,h1,h2…hj…hm-1,hm)和坝块起始浇筑高程，计算各混凝土浇筑坯层表面设计高程(H0,H1,H2…Hj…Hm-1,Hm)，其中，j＝0,1,2…(m-1)，m表示混凝土浇筑坯层表面设计高程的总层数，hj表示仓面第j坯层设计厚度，Hj表示第j坯层混凝土浇筑坯层表面设计高程，第0坯层表示第1坯层以下部位；
步骤2，浇筑当前坯层，即j+1坯层时，通过仓面施工监控设备，获取当前坯层施工机械的轨迹点(p1,p2…pi…pn)，其中，i＝1,2,…,n，pi表示施工机械的第i个实时监控点，n表示当前坯层施工机械的实时监控点总数，根据实时监控点pi的真实坐标与浇筑混凝土面空间坐标点的相互关系，得到实时监控点pi在浇筑混凝土表面上的实时投影点Pi，其中，Pi点的空间坐标为(Pix，Piy，Piz)；
步骤3，根据实时投影点Pi的高程信息Piz赋予实时投影点Pi坯层编码属性：若|Piz-Hj+1|＜0.5(Hj+2-Hj+1)，则实时投影点Pi属于j+1坯层，赋予j+1坯层编码属性的点Pi记为P(j+1)i，P(j+1)i即为j+1坯层第i个实时投影点的浇筑混凝土表面点云数据，P(j+1)i点的空间坐标为(P(j+1)ix，P(j+1)iy，P(j+1)iz)；
步骤4，对上一坯层，即j坯层的浇筑混凝土表面点云数据Pja进行分组，其中，Pja为赋予j坯层属性的实时投影点Pa，Pa为浇筑j坯层时所获取的实时监控点pa在浇筑混凝土表面上的实时投影点，a＝1,2,…,A，Pja在浇筑j坯层时根据步骤2和步骤3获得：若j≠0时，采用K-means算法将分散密集的该浇筑混凝土表面点云数据Pja分成K组，并计算每个分组的特征中心点Sjak，k＝1,2,…,K，Sjak点的空间坐标为(Sjakx，Sjaky，Sjakz)；若j＝0，则将0坯层均匀划分K个点云区域，第k个点云区域的中心点则为S0ak的取值；
步骤5，采用k-近邻算法计算当前坯层，即j+1坯层实际浇筑坯层厚度：以j坯层分组的点云数据的各簇类特征中心点Sjak为聚类中心，通过k-近邻算法判断当前坯层实时监控点pi所对应的P(j+1)i的所属簇类，设P(j+1)i的所属簇类即为k组，计算P(j+1)i的高程与P(j+1)i所述簇类的特征中心点Sjak的高程之差Δh＝P(j+1)iz-Sjakz；若下一实时监控点pi+1所对应的P(j+1)(i+1)的所属簇类发生变化，则按照“Δh＝P(j+1)iz-Sjakz”计算k组中所有实时监控点高程与特征中心点Sjak的高程之差为(Δh1,Δh2…Δhb…ΔhB)，b＝1,2,…,B，B为k组中实时监控点的个数，Δhb为k组中第b个实时监控点高程与特征中心点Sjak的高程之差，并计算(Δh1,Δh2…Δhb…ΔhB)的平均值
Figure GDA0003650718620000032
将平均值
Figure GDA0003650718620000033
作为该区域的实时浇筑坯层厚度，并与坯层设计厚度hj+1对比分析；
步骤6、通过相机三维注册技术实现真实场景与虚拟场景的融合，根据实时计算分析获得浇筑坯层厚度
Figure GDA0003650718620000034
建立实际坯层厚度虚拟模型，并展示在真实场景的相应位置；进行浇筑厚度质量评价，若发生坯层实际厚度与设计厚度的偏差大于设定阈值，则发出报警信息；对坯层浇筑厚度不合格区域进行标记，若不合格区域离开相机视野，则采用虚拟指针指示不合格区域位置，提醒平仓机驾驶员及时纠偏。
进一步地，步骤4中，所述的采用K-means算法将分散密集的该浇筑混凝土表面点云数据Pja分成K组，并计算每个分组的特征中心点Sjak，包括：
步骤4-1，将j坯层仓面区域划分为1m×1m网格，获得仓面区域中大小为c×d的点阵K作为初始聚类中心Sjak0，其中，c为点阵列数，d为点阵行数，采用欧式距离法判定Pja所属分类，即，采用公式(2)计算j坯层点Pja与各个初始聚类中心Sjak0的欧式距离dk：
Figure GDA0003650718620000031
式中，(Sjak0x，Sjak0y，Sjak0z)为Sjak0点的空间坐标；(Pjax，Pjay，Pjaz)为Pja点的空间坐标；
若dk最小，则点Pja属于第k簇类，从而获得了各簇类所包含的实时投影点；
步骤4-2，计算各个簇类所包含的实时投影点的Pja均值作为该簇类的新特征中心点Sjak’；
步骤4-3，将步骤4-2获得的Sjak’代入公式(2)，重复步骤4-1和步骤4-2，直到收敛，则j坯层数据分组完成，并获得最终各簇类的特征中心点Sjak。
进一步地，步骤5中，所述的以上一坯层分组的点云数据的各特征中心点Sjak为聚类中心，通过k-近邻算法判断当前坯层实时监控点pi所对应的P(j+1)i的所属簇类，包括：
采用公式(3)计算j+1坯层点P(j+1)i与j坯层各簇类特征中心点Sjak的欧式距离rk：
Figure GDA0003650718620000041
若rk最小，则P(j+1)i属于第k簇类。
本发明的有益效果是：本发明提出的仓面混凝土浇筑厚度监测方法以仓面施工监控系统为基础，分析监控系统获取的混凝土坯层的点云数据，利用聚类算法实时检测混凝土浇筑的厚度，通过增强现实技术直观展现在车载显示屏上，并对厚度偏差较大的区域进行有效预警，帮助施工现场及时精准解决坯层浇筑问题，改善了施工现场的决策条件并极大提高了仓面混凝土浇筑的质量。
附图说明
图1：本发明基于聚类的混凝土浇筑坯层厚度监控方法流程图；
图2：本发明的点云数据坯层编号示意图。
具体实施方式
为能进一步了解本发明的发明内容、特点及功效，兹例举以下实施例，并配合附图详细说明如下：
本发明结合仓面施工监控设备和实现手段，提供了一种基于聚类的混凝土浇筑坯层自动监测和可视化的方法。主要利用仓面施工监控设备获取实时定位数据和处理历史数据，采用聚类算法划分历史点云数据，并通过近邻算法实时分析浇筑坯层的厚度，结合增强现实技术，将实际浇筑坯层厚度在机载显示屏中可视化展示，进而及时处理混凝土浇筑坯层厚度偏差，提高坯层浇筑质量，增强坝体安全稳定性。本发明的优点在于能够实时监测评价仓面局部区域浇筑坯层厚度，并可视化标记不合格区域，实现混凝土浇筑坯层厚度精准控制。
如附图1所示，一种基于聚类的混凝土浇筑坯层厚度监控方法，采用聚类的方法实时分析点云数据得到浇筑坯层厚度信息，并结合增强现实技术，利用可视化的手段实现信息展示和预警，具体包括以下步骤：
步骤1，根据仓面坯层设计厚度(h0,h1,h2…hj…hm-1,hm)和坝块起始浇筑高程，计算各混凝土浇筑坯层表面设计高程(H0,H1,H2…Hj…Hm-1,Hm)，其中，j＝0,1,2…(m-1)，m表示混凝土浇筑坯层表面设计高程的总层数，hj表示仓面第j坯层设计厚度，Hj表示第j坯层混凝土浇筑坯层表面设计高程，第0坯层表示第1坯层以下部位，将第一坯层开始浇筑高程位置设为编号0，第一坯层顶高程位置编号设为1，以此类推，如图2所示。
步骤2，浇筑当前坯层，即j+1坯层时，通过仓面施工监控设备，获取当前坯层施工机械上A点的实时轨迹(p1,p2…pi…pn)，其中，i＝1,2,…,n，pi表示施工机械的第i个实时监控点，n表示当前坯层施工机械的实时监控点总数，根据实时监控点pi的真实坐标与浇筑混凝土面空间坐标点的相互关系，根据公式(1)得到实时监控点pi在浇筑混凝土表面上的实时投影点Pi，其中，Pi点的空间坐标为(Pix，Piy，Piz)。
Pi＝f(pi) (1)
步骤3，坝块浇筑过程中采用分坯层浇筑的方式进行，系统根据仓面设计标准自动识别点云数据的坯层编号。根据实时投影点Pi的高程信息Piz赋予实时投影点Pi坯层编码属性：若|Piz-Hj+1|＜0.5(Hj+2-Hj+1)，则实时投影点Pi属于j+1坯层，赋予j+1坯层编码属性的点Pi记为P(j+1)i，P(j+1)i即为j+1坯层第i个实时投影点的浇筑混凝土表面点云数据，P(j+1)i点的空间坐标为(P(j+1)ix，P(j+1)iy，P(j+1)iz)。
步骤4，通过聚类算法实监测计算坯层厚度，以第j+1坯层厚度实时监测为例进行说明。进行第j+1坯层厚度实时分析之前需要划分j坯层的点云数据。同一坯层历史监控数据点通过K-means算法进行分组，实现仓面区域划分，进而完成混凝土起伏区域划分，具体包括：
将j坯层的历史数据从数据库中取出。对上一坯层，即j坯层的浇筑混凝土表面点云数据Pja进行分组，其中，Pja为赋予j坯层属性的实时投影点Pa，Pa为浇筑j坯层时所获取的实时监控点pa在浇筑混凝土表面上的实时投影点，a＝1,2,…,A，Pja在浇筑j坯层时根据步骤2和步骤3获得：若j≠0时，采用K-means算法将分散密集的该浇筑混凝土表面点云数据Pja分成K组，并计算每个分组的特征中心点Sjak，k＝1,2,…,K，Sjak点的空间坐标为(Sjakx，Sjaky，Sjakz)。
步骤4-1，将j坯层仓面区域划分为1m×1m网格，获得仓面区域中大小为c×d的点阵K作为初始聚类中心Sjak0，其中，c为点阵列数，d为点阵行数，采用欧式距离法判定Pja所属分类，即，采用公式(2)计算j坯层点Pja与各个初始聚类中心Sjak0的欧式距离dk：
Figure GDA0003650718620000061
式中，(Sjak0x，Sjak0y，Sjak0z)为Sjak0点的空间坐标；(Pjax，Pjay，Pjaz)为Pja点的空间坐标；
若dk最小，则点Pja属于第k簇类，从而获得了各簇类所包含的实时投影点；
步骤4-2，计算各个簇类所包含的实时投影点的Pja均值作为该簇类的新特征中心点Sjak’；
步骤4-3，将步骤4-2获得的Sjak’代入公式(2)，重复步骤4-1和步骤4-2，直到收敛，则j坯层数据分组完成，并获得最终各簇类的特征中心点Sjak。
若j＝0，0坯层的分组与其他坯层的分组方法不同。0坯层高程为第一坯层浇筑起始高程，坯层各点高程相同，将仓面区域划分为1m×1m网格，获得仓面区域中大小为c×d的点阵K，第k个点云区域的中心点则为S0ak的取值。
坯层厚度通过k-近邻算法实现实时数据点分组，完成局部区域平均浇筑厚度监测，具体包括：
分析j+1坯层在浇筑施工过程实时数据时，采用k-近邻算法计算当前坯层，即j+1坯层实际浇筑坯层厚度：以j坯层分组的点云数据的各簇类特征中心点Sjak为聚类中心，通过k-近邻算法判断当前坯层实时监控点pi所对应的P(j+1)i的所属簇类，即，采用公式(3)计算j+1坯层点P(j+1)i与j坯层各簇类特征中心点Sjak的欧式距离rk：
Figure GDA0003650718620000062
若rk最小，则P(j+1)i属于第k簇类。聚类中心为Sjak。
计算P(j+1)i的高程与P(j+1)i所述簇类的特征中心点Sjak的高程之差Δh＝P(j+1)iz-Sjakz。P(j+1)iz为j+1坯层Pi点高程，Sjakz为j坯层第k簇内中心点高程。
若下一实时监控点pi+1所对应的P(j+1)(i+1)的所属簇类发生变化，则按照“Δh＝P(j+1)iz-Sjakz”计算k组中所有实时监控点高程与特征中心点Sjak的高程之差为(Δh1,Δh2…Δhb…ΔhB)，b＝1,2,…,B，B为k组中实时监控点的个数，Δhb为k组中第b个实时监控点高程与特征中心点Sjak的高程之差，并计算(Δh1,Δh2…Δhb…ΔhB)的平均值
Figure GDA0003650718620000072
将平均值
Figure GDA0003650718620000073
作为该区域的实时浇筑坯层厚度，并与坯层设计厚度hj+1对比分析。
步骤5，混凝土浇筑坯层厚度自动检测后评价浇筑厚度质量，采用增强现实技术可视化展示混凝土浇筑坯层厚度，建立虚拟混凝土坯层模型，并实现厚度不合格的区域精准预警。
车载主机通过相机三维注册技术实现真实场景与虚拟场景的融合，根据实时计算分析获得浇筑坯层厚度
Figure GDA0003650718620000071
建立实际坯层厚度虚拟模型，并添加到真实场景中；进行浇筑厚度质量评价，若发生坯层实际厚度与设计厚度的偏差大于设定阈值，则在虚拟模型中进行标识；对坯层浇筑厚度不合格区域进行标记，当厚度缺陷位置离开相机视野范围，则在相机视野区域生成指针指向厚度缺陷位置，帮助及时修正混凝土浇筑的缺陷。
在仓面施工过程中持续进行上述过程，则完成全仓面的浇筑厚度监测和预警，可控制仓面混凝土浇筑保持在合理区间内，保证仓面混凝土浇筑的质量，提高混凝土坝安全稳定性。
尽管上面结合附图对本发明的优选实施例进行了描述，但是本发明并不局限于上述的具体实施方式，上述的具体实施方式仅仅是示意性的，并不是限制性的，本领域的普通技术人员在本发明的启示下，在不脱离本发明宗旨和权利要求所保护的范围情况下，还可以做出很多形式，这些均属于本发明的保护范围之内。

Claims (3)

1.一种基于聚类的混凝土浇筑坯层厚度监控方法，其特征在于，包括以下步骤：
步骤1，根据仓面坯层设计厚度(h0,h1,h2…hj…hm-1,hm)和坝块起始浇筑高程，计算各混凝土浇筑坯层表面设计高程(H0,H1,H2…Hj…Hm-1,Hm)，其中，j＝0,1,2…(m-1)，m表示混凝土浇筑坯层表面设计高程的总层数，hj表示仓面第j坯层设计厚度，Hj表示第j坯层混凝土浇筑坯层表面设计高程，第0坯层表示第1坯层以下部位；
步骤2，浇筑当前坯层，即j+1坯层时，通过仓面施工监控设备，获取当前坯层施工机械的轨迹点(p1,p2…pi…pn)，其中，i＝1,2,…,n，pi表示施工机械的第i个实时监控点，n表示当前坯层施工机械的实时监控点总数，根据实时监控点pi的真实坐标与浇筑混凝土面空间坐标点的相互关系，得到实时监控点pi在浇筑混凝土表面上的实时投影点Pi，其中，Pi点的空间坐标为(Pix，Piy，Piz)；
步骤3，根据实时投影点Pi的高程信息Piz赋予实时投影点Pi坯层编码属性：若|Piz-Hj+1|＜0.5(Hj+2-Hj+1)，则实时投影点Pi属于j+1坯层，赋予j+1坯层编码属性的点Pi记为P(j+1)i，P(j+1)i即为j+1坯层第i个实时投影点的浇筑混凝土表面点云数据，P(j+1)i点的空间坐标为(P(j+1)ix，P(j+1)iy，P(j+1)iz)；
步骤4，对上一坯层，即j坯层的浇筑混凝土表面点云数据Pja进行分组，其中，Pja为赋予j坯层属性的实时投影点Pa，Pa为浇筑j坯层时所获取的实时监控点pa在浇筑混凝土表面上的实时投影点，a＝1,2,…,A，Pja在浇筑j坯层时根据步骤2和步骤3获得：若j≠0时，采用K-means算法将分散密集的该浇筑混凝土表面点云数据Pja分成K组，并计算每个分组的特征中心点Sjak，k＝1,2,…,K，Sjak点的空间坐标为(Sjakx，Sjaky，Sjakz)；若j＝0，则将0坯层均匀划分K个点云区域，第k个点云区域的中心点则为S0ak的取值；
步骤5，采用k-近邻算法计算当前坯层，即j+1坯层实际浇筑坯层厚度：以j坯层分组的点云数据的各簇类特征中心点Sjak为聚类中心，通过k-近邻算法判断当前坯层实时监控点pi所对应的P(j+1)i的所属簇类，设P(j+1)i的所属簇类即为k组，计算P(j+1)i的高程与P(j+1)i所述簇类的特征中心点Sjak的高程之差Δh＝P(j+1)iz-Sjakz；若下一实时监控点pi+1所对应的P(j+1)(i+1)的所属簇类发生变化，则按照“Δh＝P(j+1)iz-Sjakz”计算k组中所有实时监控点高程与特征中心点Sjak的高程之差为(Δh1,Δh2…Δhb…ΔhB)，b＝1,2,…,B，B为k组中实时监控点的个数，Δhb为k组中第b个实时监控点高程与特征中心点Sjak的高程之差，并计算(Δh1,Δh2…Δhb…ΔhB)的平均值
Figure FDA0003650718610000021
将平均值
Figure FDA0003650718610000022
作为该区域的实时浇筑坯层厚度，并与坯层设计厚度hj+1对比分析；
步骤6、通过相机三维注册技术实现真实场景与虚拟场景的融合，根据实时计算分析获得浇筑坯层厚度
Figure FDA0003650718610000023
建立实际坯层厚度虚拟模型，并展示在真实场景的相应位置；进行浇筑厚度质量评价，若发生坯层实际厚度与设计厚度的偏差大于设定阈值，则发出报警信息；对坯层浇筑厚度不合格区域进行标记，若不合格区域离开相机视野，则采用虚拟指针指示不合格区域位置，提醒平仓机驾驶员及时纠偏。
2.根据权利要求1所述的一种基于聚类的混凝土浇筑坯层厚度监控方法，其特征在于，步骤4中，所述的采用K-means算法将分散密集的该浇筑混凝土表面点云数据Pja分成K组，并计算每个分组的特征中心点Sjak，包括：
步骤4-1，将j坯层仓面区域划分为1m×1m网格，获得仓面区域中大小为c×d的点阵K作为初始聚类中心Sjak0，其中，c为点阵列数，d为点阵行数，采用欧式距离法判定Pja所属分类，即，采用公式(2)计算j坯层点Pja与各个初始聚类中心Sjak0的欧式距离dk：
Figure FDA0003650718610000024
式中，(Sjak0x，Sjak0y，Sjak0z)为Sjak0点的空间坐标；(Pjax，Pjay，Pjaz)为Pja点的空间坐标；
若dk最小，则点Pja属于第k簇类，从而获得了各簇类所包含的实时投影点；
步骤4-2，计算各个簇类所包含的实时投影点的Pja均值作为该簇类的新特征中心点Sjak’；
步骤4-3，将步骤4-2获得的Sjak’代入公式(2)，重复步骤4-1和步骤4-2，直到收敛，则j坯层数据分组完成，并获得最终各簇类的特征中心点Sjak。
3.根据权利要求1所述的一种基于聚类的混凝土浇筑坯层厚度监控方法，其特征在于，步骤5中，所述的以上一坯层分组的点云数据的各特征中心点Sjak为聚类中心，通过k-近邻算法判断当前坯层实时监控点pi所对应的P(j+1)i的所属簇类，包括：
采用公式(3)计算j+1坯层点P(j+1)i与j坯层各簇类特征中心点Sjak的欧式距离rk：
Figure FDA0003650718610000031
若rk最小，则P(j+1)i属于第k簇类。
CN202110912492.1A 2021-08-10 2021-08-10 一种基于聚类的混凝土浇筑坯层厚度监控方法 Active CN113804115B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110912492.1A CN113804115B (zh) 2021-08-10 2021-08-10 一种基于聚类的混凝土浇筑坯层厚度监控方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110912492.1A CN113804115B (zh) 2021-08-10 2021-08-10 一种基于聚类的混凝土浇筑坯层厚度监控方法

Publications (2)

Publication Number Publication Date
CN113804115A CN113804115A (zh) 2021-12-17
CN113804115B true CN113804115B (zh) 2022-08-30

Family

ID=78942944

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110912492.1A Active CN113804115B (zh) 2021-08-10 2021-08-10 一种基于聚类的混凝土浇筑坯层厚度监控方法

Country Status (1)

Country Link
CN (1) CN113804115B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117888717A (zh) * 2024-03-14 2024-04-16 武汉博雅宏科技有限公司 一种混凝土浇注流量的补差控制方法及系统

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103195254A (zh) * 2013-02-21 2013-07-10 中国水电顾问集团成都勘测设计研究院 混凝土浇筑振捣坯层自动监控方法
CN104266625A (zh) * 2014-10-09 2015-01-07 中国电建集团成都勘测设计研究院有限公司 混凝土坯层厚度监测方法
WO2016016957A1 (ja) * 2014-07-30 2016-02-04 株式会社 日立ハイテクノロジーズ パターン高さ測定装置、及び荷電粒子線装置
CN107218894A (zh) * 2017-04-28 2017-09-29 合肥雅视智能科技有限公司 一种快速稳定的亚像素级精度器件厚度检测方法
CN108871149A (zh) * 2018-05-08 2018-11-23 青岛理工大学 钢筋混凝土保护层厚度测量装置及测量方法
CN111492200A (zh) * 2020-03-17 2020-08-04 长江存储科技有限责任公司 用于半导体结构厚度测量的方法和系统
CN113155043A (zh) * 2021-04-26 2021-07-23 深圳市鹰眼在线电子科技有限公司 透明材质屏幕厚度测量系统、方法及装置

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2527954B (en) * 2013-01-30 2019-01-16 Univ Cincinnati Measuring wall thickness loss for a structure
JP7244336B2 (ja) * 2019-04-10 2023-03-22 清水建設株式会社 山岳トンネルコンクリート厚測定方法および測定装置

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103195254A (zh) * 2013-02-21 2013-07-10 中国水电顾问集团成都勘测设计研究院 混凝土浇筑振捣坯层自动监控方法
WO2016016957A1 (ja) * 2014-07-30 2016-02-04 株式会社 日立ハイテクノロジーズ パターン高さ測定装置、及び荷電粒子線装置
CN104266625A (zh) * 2014-10-09 2015-01-07 中国电建集团成都勘测设计研究院有限公司 混凝土坯层厚度监测方法
CN107218894A (zh) * 2017-04-28 2017-09-29 合肥雅视智能科技有限公司 一种快速稳定的亚像素级精度器件厚度检测方法
CN108871149A (zh) * 2018-05-08 2018-11-23 青岛理工大学 钢筋混凝土保护层厚度测量装置及测量方法
CN111492200A (zh) * 2020-03-17 2020-08-04 长江存储科技有限责任公司 用于半导体结构厚度测量的方法和系统
CN113155043A (zh) * 2021-04-26 2021-07-23 深圳市鹰眼在线电子科技有限公司 透明材质屏幕厚度测量系统、方法及装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Process Evaluation for Smart Concrete Road Construction:Road Surface and Thickness Evaluation Using High-Speed LiDAR Technology;Patric Skalecki等;《Automation》;20210303;第31-47页 *
基于因果模糊聚类法预测覆岩裂隙带高度研究;李波等;《中国安全生产科学技术》;20131130;第9卷(第11期);第25-29页 *

Also Published As

Publication number Publication date
CN113804115A (zh) 2021-12-17

Similar Documents

Publication Publication Date Title
CN106182765B (zh) 基于支持向量机的3d打印模型尺寸误差预测方法
CN103136393B (zh) 一种基于网格划分的区域覆盖率计算方法
CN106339707B (zh) 一种基于对称性特征的仪表指针图像识别方法
CN108257114A (zh) 一种基于深度学习的输电设备缺陷自动识别方法
WO2023045299A1 (zh) 基于三维轮廓的路面技术状况检测方法及设备
CN103400373B (zh) 相机检校控制场人工标志像点坐标自动识别定位方法
CN107818430A (zh) 一种基于bim的建筑工程施工进度信息管理方法及系统
CN113804115B (zh) 一种基于聚类的混凝土浇筑坯层厚度监控方法
CN105034366A (zh) 一种3d打印机及其打印还原度的检测方法
CN102147982B (zh) 一种扇区动态容量预测的方法
CN106204560B (zh) 菌落挑选仪自动标定方法
CN109857097A (zh) 一种路面碾压过程的监控方法及系统
CN104713491A (zh) 可获取边坡变形三维数据的边坡监测系统及其获取边坡变形三维数据的方法
CN108197283A (zh) 利用河道实测水位进行洪水二、三维动态展示的方法
CN109814536A (zh) 一种基于机器视觉套准控制系统的套准偏差检测方法
CN101968877A (zh) 一种双层圆弧的编码标记点设计方法
CN110929710A (zh) 一种基于视觉的仪表指针读数自动识别方法与系统
CN112381944B (zh) 一种基于点云数据的输电线路杆塔形变4d观测方法
CN104899852A (zh) 一种基于遗传算法及支持向量机的微小颗粒检测系统标定方法
CN103605680A (zh) 一种数据库地图的功能实现方法
CN114463932B (zh) 非接触式施工安全距离主动动态识别预警系统及方法
CN105180886A (zh) 一种测量冷轧薄钢板应变分布的方法
CN109584208B (zh) 一种针对工业结构缺陷智能识别模型的检验方法
CN107341589A (zh) 精工生产制造过程精细化管理方法
CN114932341B (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