点击下方卡片，关注「3D视觉工坊」公众号
选择星标，干货第一时间送达
来源：3D视觉工坊
星球内有20多门3D视觉系统课程、3DGS独家系列视频教程、顶会论文最新解读、海量3D视觉行业源码、项目承接、求职招聘等。想要入门3D视觉、做项目、搞科研，欢迎加入！
为提升滚动直线导轨副中含微小圆弧特征复杂曲面的制造精度，保障滚动功能部件的刚度、摩擦力矩等服役性能，本文提出一种基于线结构光的滚动直线导轨型面精密高效测量方法。针对导轨型面的独特场景，构建了一套场景化改进算法的型面点云拟合体系，包括滚道区域加权的自适应体素滤波算法、多特征融合的PCA-K邻域拐点辨识算法及自适应半径约束的改进RANSAC圆弧拟合算法等。通过与三坐标测量机、影像测量仪的对比实验，验证该测量系统的测量精度优于 ，达到SP级测量标准，能够实现导轨型面全轮廓检测，测量效率较三坐标测量机提升15倍以上，为滚动直线导轨型面的高精度、高效检测提供了理论支撑与技术方案。
1 引言
以滚珠丝杠副、滚动直线导轨副为代表的滚动功能部件，是保障数控机床系统刚性和加工精度保持性的核心部件，也是军民两用高端装备的关键基础部件。对于滚动导轨而言，滑块与导轨间常利用多组微小型的球、柱滚动体支撑与传递载荷，使得导轨表面形成了含有微小圆弧特征的复杂曲面。此类复杂型面一般通过精密铣磨工艺制造形成，其制造精度对滚动功能部件的刚度、摩擦力矩等服役性能有着严重影响。因此，精密、高效的导轨型面测量技术对提升导轨服役性能具有重要意义。
导轨型面精度测量可分为接触与非接触式两类。接触式测量技术以千分表测量法为代表，通过触头直接接触待测表面获取数据，虽具备较高精度，但易受人为操作、环境振动等因素干扰，且检测流程繁琐，难以实现导轨型面的全面指标检测。相比之下，非接触式检测方法如激光干涉法、激光准直法、激光三角法及电磁测量法等，凭借多参量同步高效测量、无磨损和广泛的尺寸适用范围等优势，成为近年来的研究热点。国立虎尾科技大学与国立中兴大学基于电涡流传感器合作研发出一台静/动态测量直线导轨副运行平行度误差的检测装置，以滑块为间接测量对象，利用两个电涡流传感器测量直线导轨的水平和垂直两个方向的平行度误差，但仅针对单一参数，无法覆盖型面整体特征。广东凯特基于非接触式气电传感器探测头与导轨底面距离位移信号的转换，实现直线导轨行走平行度的误差测量，同样存在检测维度有限的问题。南理工冯虎田等结合国内外大型筒类工件半径测量技术原理，提出的点激光型面检测方法，虽能获取三维空间误差，但单点扫描效率较低，难以满足批量检测需求。
近年来，线结构光技术日益进步，其非接触、多点同步测量特点极大提高了现有专用量仪的测量效率，已在齿轮、轴承、航空薄壁件等领域取得了诸多成功应用。北京理工大学刘少丽等设计了一套用于大型构件的主动式线结构光旋转扫描检测系统，系统有效扫描范围大于 ，构件表面栅格宽度、栅格筋厚度等尺寸的检测误差小于0.05mm，但所能检测的特征较为简单，对于复杂特征的评估及相关算法的开发存在不足。湖南大学马子骥等提出一种基于多线结构光视觉的钢轨波磨精确测量方法，利用平行多线结构光来校准失真钢轨断面轮廓，60轨1/3轨头宽度处的单点测量精度低于 ，但评估指标较为单一，仍难以满足微米级的测量要求。西南交通大学李文涛等为解决钢轨断面全轮廓测量和线结构光测量系统现场标定困难的问题，提出了一种基于自由平面靶标的线结构光双目现场标定方法，总体测量误差约为0.542mm，但标定过程较为繁琐。北京交通大学祝祥等建立了一种基于线结构光传感器的轨道板测量系统，利用轨道板中特殊位置对整个点云数据进行位置和姿态的校正，实现轨道板四项关键几何参数的快速非接触测量，钳口高度和钳口距离测量精度 ，但对于复杂微小特征，检测精度有待提高。北京工业大学石照耀等基于激光三角测量原理，研制出齿轮线激光三维测量仪，可实现齿轮齿廓偏差、齿距偏差、拓扑偏差等项目的测量与评定，能够满足5级精度（ ）齿轮的检测要求，但仍存在很多未解构的齿面数据，如齿面螺旋线偏差的测量。线结构光技术在齿轮、钢轨等工业检测中展现出良好的应用价值，但仍存在检测指标不完全、微米级精度难以保证、算法开发不足等问题。
针对上述问题，本文提出一种基于线结构光点云技术的非接触式导轨型面检测方案，实现滚动导轨的高效测量与评价。首先，开展基于线结构光的非接触多点可移动测量方法研究，设计导轨型面自动检测装置方案。随后，运用点云处理技术，构建了一套融合场景化改进算法的型面点云拟合体系。最后，进行非接触式线结构光检测方法与接触式三坐标导轨检测方法的对比实验，验证测量系统的精度与稳定性。
2 基于线结构光的导轨型面测量方法
本文提出一种非接触式线结构光的滚动直线导轨型面精密高效测量方法，通过线结构光拼接扫描技术实现导轨型面的全面测量，并对采集点云利用点云滤波、簇类分割、特征拟合等技术实现评价。
2.1 测量系统构架与检测指标
根据国内外高认可度滚动直线导轨生产厂家制定的选型手册，调研滚动直线导轨各规格使用情况，制定四列式单圆弧直线导轨型面检测指标，如图1 所示，待测型面参数主要包括滚道圆弧半径 、滚道中心距 、底面基准和侧面基准的平面度 和平行度 等。
检测装置采用“三向移动+多传感器协同”架构，如图1所示，其中，Y向移动系统由直线电机驱动，气浮导轨导向，定位精度 m，是导轨型面扫描任务的主要运动轴； 向移动系统由伺服电机驱动，滚珠丝杠传动，直线导轨导向，用于线激光测量仪与不同规格导轨的基准配准；测量平台布置三台线激光测量仪（量程 ，分辨率 ），分别扫描导轨的底面、两侧滚道面，获取导轨型面完整的高精度三维数据。
2.2 多传感器坐标对齐与点云配准实现
三台线激光测量仪独立安装，其自身坐标系存在差异，在导轨型面拟合之前，首先要对初始点云实现坐标对齐。
首先进行坐标相对齐，以高精度标准量块（粗糙度 ）为基准，通过千分表手动定位确保量块长度方向与进给轴方向平行。数控系统控制三轴联动，三台传感器同步采集量块点云。以数控系统输出的三轴定位坐标作为基准，建立传感器与世界坐标系的对应关系，求解旋转矩阵 和平移向量 ，确定初始变换矩阵完成相对齐。
其次进行点云精配准，针对机械安装误差导致的重叠区域微小偏差，采用ICP算法优化。读取粗对齐后的源/参考点云并计算法向量，配置核心配准参数，迭代最小化欧式距离误差直至收敛，输出最终变换矩阵完成精配准。其中，变换收敛阈值通常取值范围为 ，取值越小，配准效果越好；最大对应点距离取值 ，通常为点云间距的 倍。点云坐标粗对齐效果如图2所示，精配准效果如图2所示，对三台传感器采集到的数据整体对齐，效果如图2所示。
3 导轨型面点云场景化改进算法
针对导轨型面“微小圆弧密集、基准面与滚道耦合干扰”的独特场景，构建了一套融合场景化改进算法的型面点云拟合体系，如图3所示，为滚动直线导轨型面的高精度、高效检测提供了理论支撑与技术方案。
3.1 滚道区域加权的自适应体素滤波
线激光位移传感器采集到的点云数据高达数百万个点，巨大的数据量会增加后续运算时间，降低计算效率，甚至对运算设备提出较高要求。体素滤波是针对三维点云的滤波手段，可实现点云去噪与质量优化，其原理是在点云模型中构建三维体素栅格，并将每个栅格内的点替换为重心，以此完成密集点云的精简，同时滤除部分漂移噪声。
传统体素滤波采用均匀体素尺寸，在导轨型面检测中存在两大问题：一是曲率较大的滚道圆弧区域需保留更多细节，但均匀体素易导致点云稀疏，影响拟合精度；二是基准面区域形貌平坦且数据冗余度高，增加后续计算量。以 滚道半径为例，导轨滚道段的点云覆盖率仅 ，采用均匀体素滤波，取体素尺寸 时，拟合半径误差达 。
本文提出一种“区域权重因子+动态体素尺寸”的自适应滤波方法。首先进行区域粗分割，基于线激光点云分布特征，将导轨侧面点云划分为滚道区域与基准面区域两部分。其中，滚道区域点云X轴坐标处于 区间，其余部分为基准面区域。此处分割仅用于权重分配，无需精确分割。
分配权重因子时，滚道区域优先保留细节，基准面区域侧重数据精简。定义滚道区域权重因子 ，基准区域权重因子 。滚道区域是导轨承载与传动的核心功能区，其点云细节直接决定滚道圆弧半径、中心距等关键参数的拟合精度，基准面仅作为测量定位的辅助基准，形貌平坦无复杂特征，核心需求是精简数据以提升计算效率。因此权重分配遵循“滚道区域优先保留细节，基准面区域优先做精简”原则，满足 。
具体取值范围通过预实验确定，对 导轨测试，当 时，滚道点云过度稀疏，导致拟合误差显著增大（ ）；当 时，数据精简效果有限，计算效率提升不明显。因此，确定 取值区间为[1.2，1.8]。同理， 取值区间为[0.3，0.8]，确保基准面在精简后仍保有足够的点进行稳定可靠的平面度与平行度评估。
以滚道拟合误差最小为目标，通过正交实验对 和 水平组合进行寻优。对每个因素在区间内设置3个水平， 分别取1.2，1.5，1.8， 分别取0.3，0.5，0.8，采用 正交表，对9组 组合执行自适应体素滤波后，取最优组合。
设基础体素尺寸 ，其中V为点云包围盒体积，N为点云总数，由点云密度确定。计算实际动态体素尺寸，滚道区域 ，基准面区域 。对点云中每个体素，计算内部所有点的重心坐标，保留离重心最近的点作为体素代表，完成下采样。
以某四列式单圆弧导轨为例，对比传统滤波与改进滤波的效果，如图4所示。取最优组合 ，体素尺寸（ 0.05， ），应用改进算法在滚道区域的点云覆盖率由 提升至 ，数据量减少 ，同时兼顾了细节保留与计算效率。
3.2 导轨型面点云切片提取
导轨型面检测中，滚道圆弧的精确拟合直接影响着滚道自身半径参数的计算和基准面平行度的测量精度。若直接将导轨滚道面看作一个不完整的圆柱进行拟合，常见的拟合算法中，最小二乘法需要大量的数学计算，RANSAC算法对外部参数的输入高度依赖。因此在未知半径、未知法向量的情况下，直接对滚道点云进行圆柱拟合的效果不佳。
切片提取是将三维点云转换为二维轮廓的关键环节，本文提出一种对导轨型面点云切片处理，将滚道圆柱面的拟合转换为平面圆的拟合思路。线结构光采集的导轨点云Y轴坐标对应激光器运动位移，连续触发模式下其组内差值固定，基于此特征采用直通滤波器可实现导轨的轮廓切片提取，提取效果见图5所示。
3.3 基于PCA-K邻域的拐点辨识算法
导轨切片提取完成后，需要对导轨轮廓拐点进行辨识，方便滚道圆弧段的分割与拟合。本文根据“主成分分析+法向量突变+轮廓连续性”三个特征，提出一种拐点辨识算法，提取原理及效果如图6所示。
主成分分析（Principal Component Analysis，PCA）是一种常用的数据降维方法，新坐标系下各维度数据量的大小反映在样本特征值的方差上，方差最大的方向作为第一主成分方向。本文采用KD-Tree的方式建立点之间的拓扑关系，对切片点云上每个点的K邻域点集进行PCA分析。设数据集中目标点 ，其中 为邻域内点的总数。
（1）潜在拐点判断 计算主成分向量 ，根据质心和特征向量 拟合一条直线 。计算邻域内所有点到 的距离集合 和距离标准差 ，反映点云偏离局部拟合直线的程度。若：
则将该点标记为潜在拐点。其中， 为距离标准差阈值，统计初始点云的噪声标准差 ，依据3σ原则，取 。
（2）候选拐点判断 对所有潜在拐点进行法向量突变校验，导轨基准面的法向量方向基本一致，而滚道圆弧与基准面过渡处法向量会发生显著偏转。计算潜在拐点与其邻域的法向量夹角 ，若：
则将该点保留为候选拐点。其中， 为法向量夹角阈值，统计标准导轨型面过渡区的法向量夹角均值 ，取 ，确保既能捕捉法向量突变特征，同时避免误判某些微小波动点。
（3）有效拐点判断 判断候选拐点的轮廓连续性，沿导轨截面轮廓方向，基准面区域的坐标变化率近似为常数，滚道圆弧区域的坐标变化率随曲率变化，当变化率 从线性转为非线性，即从轮廓面转为滚道圆弧。以相邻点曲率差与间距的比值定义轮廓坐标变化率，计算候选拐点前后5个点的坐标变化率 ，若：
则确认该点为有效拐点。
3.4 导轨型面点云特征拟合
3.4.1 自适应半径约束的改进RANSAC滚道拟合算法
随机抽样一致性（RANSAC）的基本思想是在假设数据集中包含正确数据（内点）和异常数据（外点）的前提下，通过随机采样和迭代的方式，找到能够最好地描述内点数据的数学模型。为了剔除轮廓连接处倒角等微小干扰区域的影响，本文引入先验半径约束，提出一种改进RANSAC算法。
首先构建先验半径约束，假设圆弧点云集 包含 个点，每个点表示为 ，其中 。基于导轨设计图纸的滚道理论半径 ，结合同类导轨的加工工艺误差 ，构建动态半径区间 。
对点云集随机采样与模型初始化，设迭代计数为 ，从圆周点集 中随机选取 个点，拟合初始圆模型，记为 ；若圆半径不在约束区间内，直接舍弃该模型，减少无效迭代。
筛选内点，设定偏差阈值 ，计算 中每个数据点到初始圆的距离，将距离 的点纳入内点集 。重复迭代上述步骤，不断更新内点集 ，保留内点数量最多的模型作为最终拟合结果，迭代次数 由外点概率 和模型置信度 以及初始点个数 共同决定，如式（4）计算：
对同一导轨的10个滚道截面进行拟合实验，如 时，导轨的工艺误差 ，动态半径区间为 ；设定偏差阈值 ，取置信度 ，外点率 ，采样点数 ，对比传统RANSAC与改进算法的结果，如表1所示，改进算法的拟合半径误差均值降至 以内，稳定性显著提升。
3.4.2 基于最小二乘-RANSAC的基准平面拟合算法
导轨型面检测中，由于基准平面点云形貌结构较为简单，并未对其进行分割、聚类等前置处理，故提出一种改良的平面提取方法，在传统RANSAC方法拟合完成后，重新对RANSAC评估的最大内点模型进行最小二乘拟合。
假设平面点云集 包含 个点，每个点表示为 ，其中 。待拟合平面方程为 ，点集 ，为得到最佳拟合平面，应满足：
使点集中每个点到该平面的距离尽可能最小。分别对 求偏导数，满足 ，整理后求解平面系数方程 ，其中：
由实验可知，改良算法结合了RANSAC剔除外点与最小二乘精修内点的组合优势，实现了精度、效率与稳定性的平衡：既保留了RANSAC法的抗干扰能力，又通过最小二乘提升了拟合精度，平面度降至 的同时，耗时 与基础RANSAC法持平。对比来看，RANSAC拟合法、MILESAC法、Hough变换及改良算法对平面点云拟合，从拟合精度、运行效率、稳定性三方面对比评价，单次实验结果如表2所示，10次重复实验的稳定性指标如图7所示。实验中，测试点云 ，平面方程为 ，外点比例 ，高斯噪声标准差为0.05。
基础RANSAC和MILESAC法精度中等、效率尚可，但稳定性与精度均不及改良算法；最小二乘法适配无外点、低噪声的理想点云环境；Hough变换的参数与真实值偏差极大，稳定性弱，效率较其他方法低2个数量级。因此，改良算法是适配滚动直线导轨基准平面拟合的最优算法。
4 导轨型面点云检测方法验证实验
4.1 有效性验证
为验证导轨型面点云的拟合算法的有效性，设计非接触式线结构光检测系统、接触式三坐标测量机与影像测量仪的对比实验，实验装置如图8所示。其中，三坐标测量机为海克斯康LeitzReference系列，探针直径 ，单点测量精度至亚微米级别，用于提供高精度基准数据；影像测量仪为中旺MVA-4030系列，辅助验证滚道半径、中心距等关键参数，测量精度 ；线激光位移传感器采用基恩士8060系列，点云间隔 ，采集频率 ，扫描进给速度 ，实现全轮廓快速扫描。实验结果对比如表3所示。
对三种测量方式的测量性能综合对比，以三坐标测量值为测量基准，线结构光测量系统滚道半径误差 ，滚道中心距误差 ，核心型面参数测量误差均控制在 以内，已达到SP级测量标准；影像测量仪需将导轨切断后获取截面进行二维成像测量，其测量值与三坐标基准值偏差显著，源于切断过程中没有保证切口与导轨轴线完全垂直，导致截面轮廓产生倾斜畸变，造成测量数据偏离真实值。
另外，三坐标虽能实现全参数测量，但依赖复杂的逐点路径规划与编程，效率极低，且探针与导轨的轻微接触会产生累积性磨损，长期使用可能影响被测件的服役性能，更换不同规格导轨时还需重新调整设备参数与测量路径；影像测量仪样品切断破坏的测量方式完全不适用于批量检测与在役导轨检测场景，并且仅能获取二维截面信息，无法还原导轨三维型面的整体特征；线结构光扫描技术既避免了接触磨损与样品破坏，又能一次性获取导轨全轮廓三维数据，测量效率较三坐标提升数倍，能够实现对导轨型面的高精度快速检测。
4.2 重复性验证
为验证导轨型面点云的拟合算法的稳定性，在相同的实验条件下，对直线导轨型面的每个参数指标进行10次独立测量，测量结果如表4所示。
重复性实验中，测量标准差反映了测量过程中的随机波动。表4中导轨型面各参数的测量标准差均较小，表明线结构光导轨型面检测方法的重复性较好，点云的拟合算法的稳定性较高。
5 结论
为解决滚动直线导轨含微小圆弧特征复杂型面的高精度高效测量问题，本文围绕线结构光非接触测量技术展开研究，完成了测量方案设计、场景化算法体系构建及实验验证，主要结论如下：
（1）提出了“三向移动+多传感器协同”的线结构光测量方案，通过Y向高精度气浮导轨、X/Z向伺服配准系统及高分辨率线激光测量仪，实现导轨型面全轮廓检测，解决传统接触式测量无法覆盖完整轮廓、易磨损、效率低的问题。
（2）构建了场景化点云改进算法体系，滚道加权自适应体素滤波平衡了细节保留和数据精简的需求，PCA-K邻域拐点辨识实现了滚道与基准面精准分割，改进RANSAC圆弧拟合、最小二乘-RANSAC基准平面拟合误差均值 ，显著提升了处理精度与效率。
（3）验证了线结构光导轨型面测量系统具备高精度、高稳定性与强适用性，系统测量误差 ，达SP级标准，测量效率较三坐标测量机提升15倍以上，为滚动直线导轨型面检测提供高精度、高效的技术方案。
本文仅做学术分享，如有侵权，请联系删文。
。
添加微信：cv3d001，备注：姓名+方向+单位，邀请入群。
