基于点云数据的网架结构提升变形智能监测方法_金钦明.pdf

.，.，工业建筑 年第 卷第 期 基于点云数据的网架结构提升变形智能监测方法金钦明 程国忠 李东声 王 聪 陈莎莎 王瑞荣 毕静刚（中建铁路投资建设集团有限公司，北京；重庆大学土木工程学院，重庆；深圳大学土木与交通工程学院，广东深圳）摘 要：变形监测是保障网架结构提升施工安全的重要手段之一，目前传统方法仅能实现对局部点的变形监测，难以实现对整体结构的变形监测。三维激光扫描技术可全覆盖地得到已完成结构的精准点云数据，这为解决上述难题提供了新思路。为此，以泸州高铁站为工程背景，开展基于点云数据的网架结构提升变形智能监测研究，包括点云数据预处理、点云数据非刚性配准以及提升变形可视化三个方面。针对网架结构提升前后的点云数据非刚性配准问题，基于聚类算法、随机采样一致性算法、图结构方法以及正交普氏分析等提出了集球心智能定位、球心粗匹配、球心非刚性配准于一体的非刚性配准算法。研究结果表明，基于点云数据的网架结构提升变形智能监测方法高效、全面且实用。关键词：网架结构；变形监测；三维激光扫描；点云数据；非刚性配准 ：.（，；，；，）：，（），：；国家自然科学基金项目（）。第一作者：金钦明，男，年出生，工程师。通信作者：程国忠，。收稿日期：引 言网架结构是多根杆件按照一定规律的网格形式并通过螺栓球或焊接球连接形成的高次超静定空间结构，它可充分发挥空间三维捷径传力的优势，满足日益复杂化、美观化的建筑外观设计需求。网架结构具有刚度大、质量轻、受力合理、稳定性高、施工周期短以及工业化程度高等优势，被广泛用作体育场馆、机场航站楼、博览会展厅和高铁站房等大型民用基础设施的屋盖。网架结构的常用施工方法包括：高空散装法、分块吊装法、整体吊装法以及整体提升法等。其中，高空散装法、分块吊装法和整体吊装法存在高空作业多、施工效率低、施工安全性难以保障等限制，常用于小型网架结构的施工。对于大型网架结构而言，应用最广泛的施工方法是整体提升法。整体提升法是指在设计位置就地完成 工业建筑 年第 卷第 期网架结构拼接后，采用液压同步提升技术对网架结构进行整体提升。对于整体提升法而言，影响施工质量和安装精度的关键因素之一是网架结构提升的变形量。目前，网架结构提升变形监测手段主要包括全站仪、光电测距仪等，这些设备只能对局部点进行监测，无法对网架提升变形进行全覆盖的监测。因此，目前亟需一种全面且高效的变形监测方法。目前，三维激光扫描技术因测量精度高、扫描效率高、受外界影响小、可操作性强等优点而受到建筑业研究人员的青睐。三维激光扫描仪（图）通过主动发射激光束的方式来完成对目标点的测量，可以快速获得扫描环境内的全景三维点云数据，具有数据精度高、受外界影响小、可操作性强等优点。三维激光扫描仪应用于施工变形检测的最大优势是可以全覆盖地获得结构变形，提供更加全面、充分的变形信息。为了改善建筑业的传统监测方法和提高监测智能化程度，学者们开展了基于三维激光扫描仪的施工进度智能监测，但研究对象主要集中于刚性的建筑构部件。提升前与提升后的网架结构处于不同受力状态，存在不同程度的变形。基于三维激光扫描技术的结构变形智能监测需要将结构变形前后的点云数据进行配准。目前点云数据的非刚性配准通常需要人为选择配准点，存在人为主观误差。因此，目前相关研究方法不适用于对网架结构提升变形进行智能监测。图 三维激光扫描技术 为此，以泸州高铁站为工程背景，开展基于点云数据的网架结构提升变形智能监测研究，包括点云数据预处理、点云数据非刚性配准、提升变形可视化 个方面。针对提升前与提升后网架结构点云数据的非刚性配准，提出了集球心智能定位、球心粗匹配、球心非刚性配准于一体的综合算法。研究成果以期为网架结构提升变形监测提供理论和算法基础支撑。工程概况泸州高铁站（图）位于四川省泸州市马潭区境内，总建筑面积.万，建筑高度.。泸州高铁站主要包括侧式站房和高架站房两部分，高铁站屋盖均采用大跨网架结构，侧式站房屋盖最大跨度为 ，高架站房屋盖最大跨度为 。大跨网架结构中，圆杆均通过焊接球进行连接（图），圆杆 最 大 直 径 为 ，焊 接 球 最 大 直 径 为。图 泸州高铁站结构示意 图 典型节点 高架站房需进行 次提升：第 次是将网架结构从标高.提升到标高.，提升结构质量为 ；第 次是将网架结构从标高.提升到标高.，提升结构质量为 ；第 次是将网架结构从标高.提升到标高.，提升结构质量 。第 次提升时，网架结构共设置 个吊点，吊点的布置见图。以高架站房屋盖的第 次提升变形监测为例，对基于点云数据的网架结构提升变形智能监测方法进行介绍（图）。点云数据预处理为实现网架结构的全覆盖扫描，选用最大量程为 的 三维激光扫描仪对提升前与提升后的网架结构分别进行多站点扫描（图）。基于重叠区域的焊接球球心，利用点云处理软件对扫描点云数据进行粗拼接，继而采用迭代最近邻算法对扫描点云数据进行精细化拼接。拼接后的点云数基于点云数据的网架结构提升变形智能监测方法 金钦明，等 图 吊点布置 图 网架结构提升变形智能监测技术流程 据中包含防护网、混凝土构部件等大量背景噪声数据，这会显著地增加点云数据处理的难度。为此，拼接后的点云数据须要进行人工去噪。采用基于图结构滤波算法对去噪后的网架结构点云数据进行轻量化（图），轻量化后的点云数据见图。图 数据采集 点云数据非刚性配准非刚性配准问题通常指具有一定差异的两个目标对象通过某个仿射变换实现最优匹配。考虑到完整点云数据的非刚性配准难度大、效率低，提出采用焊接球球心代替完整点云数据进行网架结构变形监测的策略。网架结构点云数据的非刚性配准具体包括球心智能定位、球心粗匹配以及球心非刚性原始点云数据；采样的点云数据。图 点云数据轻量化 提升前；提升后。图 网架结构的部分点云数据 配准。.球心智能定位焊接球球心智能定位包括球点云数据智能识别、球点云数据智能分割以及球拟合。.球点云数据智能识别去噪后的点云数据由杆件点云数据与焊接球点云数据组成。考虑到杆件和焊接球点云数据的曲率特征、线性特征、发散性特征等局部邻域特征存在明显差异，因此可以通过逐点计算局部邻域特征对球点云数据进行智能识别。由于网架结构的点云数据量大且点云数据无拓扑关系，为提高球特征检测效率，需对点云数据建立多维二叉树（树）结构。对于任意计算数据点，采用基于范围搜索的 最近邻算法。可得到 的局部邻域点集合 ，其中 表示预设的邻域距离阈值。可由矩阵，进行表示，其中 为 的局部邻域点数量，的协方差矩阵 通过下式进行计算：()()（）式中：为 的局部邻域点均值。对 采用特征值分解可得矩阵特征值（）。线性特征 与发散性特征 可按下式进行计算：（）工业建筑 年第 卷第 期（）邻域距离阈值和特征筛选阈值是影响球点云数据智能识别的关键参数。经多次测试，计算 和的邻域距离阈值分别取.，.，同时选取.和.作为球点云数据识别准则，可实现球点云数据智能识别。图 为球点云数据智能识别的结果，其中红色数据点为球点云数据，黄色数据点为非球点云数据。图 球点云数据智能识别 .球点云数据智能分割采用聚类算法（）对识别出的球点云数据进行智能分割，从而得到每个球对应的点云数据。球点云数据智能分割时，选取任意点为种子点，以固定搜索半径进行点云数据扩充，直到固定搜索半径无点云数据为止，从而得到一个球的点云数据。根据 聚类算法原理可知，固定搜索半径是智能分割的关键参数，高架站房屋盖的杆件长度均超过 ，故固定搜索半径可取.。球节点云数据智能分割效果见图。图 球节点云数据智能分割 .球拟合采用算法（）对分割后的球点云数据进行球拟合。球拟合时，首先从分割后的球点云数据中随机选取 个点；然后根据选取的 个点计算球心坐标和球半径；接着对拟合的球进行评价，评价函数 按下式计算：（）（）其中（）式中：为第 个点到球面的欧式距离；为球面点的允许误差。重复上述过程 次，最后输出评价函数最大对应的球心坐标和球半径。对于 算法而言，最关键的 个控制参数为、以及，可根据蒙特卡洛法进行计算：（）（）（）式中：为球面点的百分比，设置为；为球拟合正确的概率，设置为.。设置为，设置为.。球拟合结果见图。图 球节点数据拟合 .球心粗匹配提升后的大跨度网架结构通常会发生较大的整体变形且网架结构存在大量的局部四点全等集，因此采用基于欧式距离的四点全等集（）算法难以实现球心匹配。网架结构提升前后相邻焊接球之间的欧式距离存在较小的差异，故可采用图结构特征作为节点球心匹配的依据。.图结构特征计算首先，对球心集建立 树结构；然后，采用 最近邻算法确定每个球心的邻近图：（）（）（）（）（）式中：和 分别为第 个和第 个节点球心；和均为球心 与节点球心 的邻近关系系数，取 或，表示节点球心 与节点球心 不存在邻近关系，表示球心 与球心 存在邻近关系；（）表示 最近邻筛选，取。最后，通过邻近关系系数可得到任意球心 的集聚系数：（），（）图 分别展示了提升前后的球心集聚系数。基于点云数据的网架结构提升变形智能监测方法 金钦明，等 提升前；提升后。图 球心集聚系数 .图结构特征匹配对.节得到的球心集聚系数进行统计，可得到球心集聚系数的分布图（图、）。可以看出：提升前后球心集聚系数的分布基本一致，这为图结构特征匹配提供了良好的理论依据。首先，从提升前的球心集中选取低频球心集聚系数（）对应的球心集；然后，从提升后的节点球心集中选取集聚系数为的球心集；接着，对球心集 和 进行全排列，生成候选的球心匹配对（）；最后，对候选的球心匹配对进行评价，选出最优的球心匹配对。候选的节点球心匹配进行评价时，对每一个候选的球心匹配对按下式求得旋转矩阵 与平动矩阵：（）（）（）式中：和 分别是 中提升前与提升后的球心坐标；和 分别是 中提升前与提升后球心坐标的均值；为 包含球心匹配对的数量；为协方差矩阵；对角矩阵、左奇异向量 以及右奇异向量 均由矩阵 奇异值分解得到。（）（）（）最优变换矩阵 与 通过下式确定：，（）式中：为球的总数量；和 分别为提升前和提升后的球心坐标，对应关系通过最近邻算法确定。最优变换矩阵 与 对应的球心匹配对即为最优的球心匹配对，基于图结构特征的球心最优匹配见图。基于图结构特征的球心匹配确定了提升前后球心的对应关系，为球心非刚性配准提供了良好的初始解。图 提升前的球心集聚系数分布 图 提升后的球心集聚系数分布 图 基于图结构特征的球心最优匹配 .球心非刚性配准由于网架结构的提升变形通常为低阶模态，提升前后的网架结构存在大量相同局部特征。然而，现有的非刚性配准算法是建立在目标对象变形前后存在少量相同局部特征且重叠区域不规则的基础上。为此，亟须提出一种面向网架结构提升变形智能监测的非刚性配准算法。低阶变形具有同号性质（图）。因此，球心非刚性配准问题可以表述为：，（）工业建筑 年第 卷第 期 ，（）式中：表示提升变形在 轴上的分量。图 低阶变形的同号性 为了高效、快速地求解式（），基于贪心策略和局部正交普氏算法提出球心非刚性配准算法（表）。算法具体步骤为：）基于粗匹配后的球心坐标，计算负变形球心的占比；）选出所有负变形球心对，组成集合；）采用正交普氏算法对集合进行匹配，从而获得变换矩阵 与；）基于变换矩阵 与，对提升后的球心进行临时变换；）基于临时变换后的球心，重新计算负变形球心占比；）若 值变小，基于变换矩阵 与 对提升后球心进行正式变换，并重复步骤）；）若 值未变小，基于贪心策略从集合提出负变形最大的球心对，并重复），直至集合的元素小于。表 球心非刚性配准算法 基于贪心策略和局部正交普氏算法的球心非刚性配准算法输入：粗匹配后的球心坐标输出：非刚性配准后的球心坐标）计算负变形球心的占比；）为中间变量，取值为；）选取负变形球心对，组成集合；）同时：）采用正交普氏算法对进行匹配，得到变换矩阵 与；）基于变换矩阵 与，对提升后的球心进行临时坐标变换；）计算负变形球心的占比 ）如果 ）基于变换矩阵 与，对提升后的球心进行坐标变换；）选取负变形球心对，组成集合）结束）此外）基于贫心策略，从中去掉负变形最大的球心对；）如果）断开）结束）结束）结束 从图 可以看出：负