1、672023年第66卷第7期航空制造技术数字化测量技术及应用Digital Measurement Technology and Application*基金项目：基础加强计划重大基础研究项目（2019JCJQDA001129）。引文格式：潘新,赵建国,刘标,等.机器人原位铣削离散点云快速精简方法研究及应用J.航空制造技术,2023,66(7):6772.PANXin,ZHAOJianguo,LIUBiao,etal.Methodandapplicationofscatteredpointclouddatafastsimplificationusedtorobot in-situ millin
2、gJ.Aeronautical Manufacturing Technology,2023,66(7):6772.机器人原位铣削离散点云快速精简方法研究及应用*潘 新，赵建国，刘 标，台春雷，刘 哲，黄 宁（航空工业沈阳飞机工业（集团）有限公司，沈阳 110034）摘要 为了满足飞机装配机器人原位铣削对数字化高精度快速测量需求，开展了适应三维点云的双边滤波处理和快速精简算法研究，开发了点云快速处理与分析应用软件，实现了空间点云数据的降噪和精简处理，降低滤波对边界特征的损伤，有效保留了复杂构件圆角、筋条、空间曲线边界等关键特征。以典型样件为测试对象，通过与三坐标测量结果对比，验证了本文研究的有效
3、性，软件接口开放性高，实现了与自动化设备的高度集成应用，避免了大量人机交互操作，缩短了数据处理时间，满足了机器人原位铣削的时效性需求。关键词：飞机装配；机器人原位铣削；离散点云；数据精简；集成应用DOI:10.16080/j.issn1671-833x.2023.07.067潘 新 正高级工程师，研究方向为飞机智能装配与数字化测量。等获得大量应用，测量得到的原始点云数据存在规模大、测点噪音与层叠等问题，严重影响数据分析精度和效率。为此大量学者对点云数据的滤波和精简进行研究，开发了多种处理算法1，国外学者率先开展了散乱点云滤波方法研究，常用方法有高斯滤波1、双边滤波24、拉普拉斯滤波5等，高斯滤
4、波通过加权平均处理抑制噪音影响，但不适用于高曲率和薄壁特征的降噪。拉普拉斯滤波计算简单，适应平滑曲面点云光顺，对尖锐和小曲率半径区域点云有损伤。双边滤波应用较多，采用邻域点的加权平均进行处理，有效保留了几何边缘特征，适用于孔洞、筋条、空间曲线边界的处理。国内学者基于前人成果，开展了大量点云精简算法研究，主要包括包围盒法6、随机采机器人原位铣削技术应用越来越广泛，三维数字化测量在保证飞机装配蒙皮铣削精度方面起到了重要作用。如采用激光测距、相位光栅、三角定位以及立体视觉等原理的扫描设备，与机器人协同完成自动铣削/制孔的部件型面扫描、全机外形测量等工作，实现了点云数据的快速采集，显著提高了飞机装配精
5、度。但对于飞机大型曲面（面积超过 20 m2），一次测量点云数量数以亿计，直接处理点云数据极大增加计算机运算负担，占用大量存储容量和内存空间，而且提取用于质量评价或加工的有效信息需耗费大量时间，显著降低了产品检测和自动化铣削的效率。因此，为提高测量数据处理效率，对点云进行压缩处理，开展面向机器人原位铣削的数据快速精简技术研究势在必行。随着数字化测量技术的发展，非接触光学测量设备，如激光跟踪仪、激光雷达、ATOS、Creaform 扫描仪68航空制造技术2023年第66卷第7期FORUM论坛样法7、角度法8、均匀网格法9和曲率采样法10，包围盒法基于八叉树理论，将所有点云数据包围在内并进行有序分
6、裂，通过采样点的融合进行点云压缩，有着良好的序列性和分层性，便于空间划分，提高了点云轻量化速度。随机采样法易于实现，按设定精简率对点云数据进行随机索引删除，精简速度快，但容易造成细节特征的损失，应用较少。角度法计算采样点与邻近两点连线间的夹角，并与设定的角度阈值比较，控制点云删除与保留，完成点云精简，计算效率较高，但造成曲率变化小区域精度损失。均匀网格法将点云分配给与扫描方向垂直的均匀平面方格，计算同一方格内各点到方格平面的距离并排序，压缩掉除“中值点”以外测点，从而完成压缩，该法对尖锐和小曲率半径区域特征损伤大。曲率采样法对点云数据进行曲面拟合，得到各处曲率估算值，曲率较小区域采用大压缩率，
7、曲率较大区域采用较小压缩率，该方法对曲面特征表述有利。基于上述理论开发的商业软件有 Geomagic、Polyworks 等，基本满足了用户使用需求，但商业软件需要大量的人机交互操作，不便于自动化设备的集成化应用。为此，国内学者针对飞机质量控制和自动加工需求，开发了快速检测分析系统11，将数字化测量技术与自动化加工设备集成12，在蒙皮下陷检测和加工方面取得了较好的应用效果，保证了检测精度，提高了检测分析效率。综上所述，为满足飞机装配机器人原位铣削快速检测分析需求，亟须开展离散点云数据快速精简技术研究，保证关键特征的检测精度，实现大型构件的高效检测与分析，提高机器人铣削自动化程度，促进飞机装配工
8、艺技术水平的不断完善。1 数据精简算法飞机装配原位铣削对关键特征测量的时效性要求高，通过测量点云降噪和压缩处理，保留关键特征信息，降低点云数据规模，为铣削边界信息的快速提取奠定基础。1.1 双边滤波双边滤波降噪效果明显，且能够保持点模型的几何特征，在扫描图像处理中具有良好的鲁棒性，其核心思想是在点云数据上设置多个邻域，并依据判定准则为域内点赋予差异化的滤波权值，进而实现点云的滤波降噪。飞机原位铣削关键特征包括圆角、筋条、加工边界等，以下采用双边滤波法处理点云数据，可有效保留特征边缘特征，降低滤波对边界特征的损伤。二维双边滤波根据图像领域内灰度值进行处理，但不适于点云数据滤波。本研究对二维滤波进
9、行改进，首先根据点云数据拟合确定表面法向，再以法向为参考计算出领域内各点高度值，然后仿照灰度值处理制定三维点云双边滤波算法，具体为pppp=+n（1）其中，pWGpqGnqpnqppq Spp=-1dr(|)(|,|),pWGpqGnqpnqppq Spp=-1dr(|)(|,|),（2）式（1）和（2）中，p 代表点 p 经过双边滤波处理后的计算值；p 为处理点；q 为 p 邻域内点；S 为点云域；np为单位法向量；d和 r分别代表定义域和值域中高斯滤波器的滤波器核，决定了滤波过程中邻域内点云值的权重。定义域滤波是空间邻域内处理点的加权平均值，其中 q 和 p 之间距离越大，相应权重系数就越
10、小，反之则越大。类似的，值域滤波也是一个加权平均值，对处理点数值邻域内的点进行处理，其中坐标差值越大，权重系数就越小，反之则越大。Wp是一个标准量，其值为WGpqGnqppq Sp=-dr(|)(|,|)WGpqGnqppq Sp=-dr(|)(|,|)（3）Gd是空间函数，定义为Gd p qd=-exp(,)122d（4）Gr是范围函数，定义为Gp qrr=-exp(,)122（5）式（4）中，d（p，q）为两点间距离差；式（5）中，（p，q）为两点间沿法向高度差值。1.2 点云精简点云精简的方法较多，工程上常用的有均匀精简法和曲率适应性精简法，均匀精简主要用于相对平缓（曲率变化小）的曲面，
11、曲率适应性精简主要用于曲率变化大具有高曲率特征和薄壁特征的曲面。点云精简算法大都基于八叉树理论制定，采用八叉树划分，可处理任意拓扑结构点云数据，便于数据的组织和压缩处理。八叉树划分的压缩算法主要流程分八叉树构建和数据编码两个阶段。首先根据点云数据极限值建立包围盒，保证所有数据均处于包围盒内。已知点云数据集 P=p1，p2，pm，计算点云集最大坐标值（xmax，ymax，zmax）和最小坐标值（xmin，ymin，zmin），建立边长为 lx、ly、lz的包围盒，有lxxlyylzzxyz=-=-=-maxminmaxminmaxmin（6）然后对包围盒进行均匀划分，一个栅格被均分为 8 个子栅
12、格。图1 为栅格划分和八叉树结构的示例，栅格细化的过程中，判断栅格是否为空，只有非空栅格才会进一步细化，细化程度根据精度要求进行。如图 2 所示，按照文献 1，为了避免部分测点处于包围盒边界，适当增加包围盒的 3 条边 lx、ly、lz边长，以保证所有点云数据点均位于包围盒内。根据点云规模和精简比例要求，将包围盒划分为 s 个边长为 a 的小栅格。slalalaxyz=（7）式中，表示向上取整运算。令点云692023年第66卷第7期航空制造技术数字化测量技术及应用Digital Measurement Technology and Application包含测量点数为 m，定义每个栅格中含有测
13、量点数均值为 m-grid，则有m-grid=m/s，可以得出每个细化栅格边长估计值为aml l lmx y z=grid1 3/（8）定义点云集为 P，通过八叉树划分，将所有点云观测点 piP 投射到相应的栅格中，每个栅格其内的点被压缩，最后观测点 pi所在包围盒栅格空间位置 gix giy gizT表示为gxxagyyagzzaixiiyiizi=-=-=-()()()minminmin （9）由式（9）即得压缩后采样点。通过如此划分，八叉树中每一个层次上的非空栅格集合就成为点云的一个近似表示，各层的近似模型就构成了原始点云的一个多细节层次的表示。栅格的每次划分均对栅格进行位置编码，非空为
14、 1，反之为 0，这样对于每次栅格细分，用一个固定的顺序遍历所有子栅格，得到占位编码信息，而且与之关联的位置和法线信息均被编码。为了降低占位编码的熵，将占位码重新排序，把非空集都推到前面，便于后续遍历索引。兼顾点云压缩规模和模型标定精度，在曲面平坦区域按上述算法保留均匀测点，在曲面高曲率区域采用基于局部法矢量变化和平均曲率的适应性精简算法，保留相对密集测点。假设某一栅格内各观测点为pi（i=1，m），则栅格内对应观测点平均曲率为 i（i=1，m），可得该栅格平均曲率为=11miim。根据精度要求制定曲率精简阈值，若-，则仅保留栅格内曲率最接近-的观测点；若-，则对该栅格再行细分，细分栅格边长
15、a 为aa=-1212cos(/)cos(/)max（10）式中，a 为原栅格长度；max为原栅格内任意两观测点法矢最大夹角；为给定的法矢变化阈值；代表向下取整运算符号。由于各栅格最大法矢夹角 max不尽相同，细分栅格长度与平均曲率-亦不相同，若-则接着细分，直至满足-，保留栅格内曲率最接近-的采样点。2 点云精简验证与集成应用基于上述算法，通过科研课题合作方式，与华中科技大学联合开发应用程序 iPoint3D-AeroSpect，如图 3 所示，采用典型样件对点云数据精简效果进行验证，并与自动铣削设备进行集成应用，证明本文方法的有效性。2.1 测量点云精简验证飞机表面对缝处多采用锯齿结构以提
16、高隐身性，该类零件具有典图 1 八叉树构建Fig.1 Octree building（a）栅格划分（b）八叉树结构图 2 点云栅格精简Fig.2 Voxel-grid simplificationlzlxlya采样点图 3 iPoint3D-AeroSpect 软件界面Fig.3 Software interface of iPoint3D-AeroSpect70航空制造技术2023年第66卷第7期FORUM论坛型边界特征，扫描测量与特征提取直接影响自动化铣削精度。仿照锯齿口盖类零件（长 1102 mm，宽 695 mm），采用高精度机床加工试验样件，模拟实际装机零件加工后喷漆处理，漆层厚度不超过 10 m，同时喷漆可降低金属零件反光性，便于点云数据采集。试验样件准备完成后，采用三坐标测量机（测量精度 3 m+4L/1000 m，L 为长度）进行精确标定，作为验证检测标准。采用三维扫描仪进行数据采集（单幅测量误差不超过 0.03 mm），能够获得精细特征细节信息，原始点云数据如图 4（a）所示（其中圆孔为夹持模块，用于自动铣削夹紧和减振），点云量为 119970312 个，包括大量噪点
