《基于机器视觉的大长径比轴类零件直线度误差检测方法.pdf》由会员分享，可在线阅读，更多相关《基于机器视觉的大长径比轴类零件直线度误差检测方法.pdf（14页完整版）》请在专利查询网上搜索。
1、(10)申请公布号 CN 102519400 A (43)申请公布日 2012.06.27 CN 102519400 A *CN102519400A* (21)申请号 201110419626.2 (22)申请日 2011.12.15 G01B 11/27(2006.01) (71)申请人 东南大学 地址 210096 江苏省南京市玄武区新街口街 道四牌楼 2 号 申请人 南京工程学院 (72)发明人 郝飞 史金飞 朱松青 张志胜 陈茹雯 韩延祥 (74)专利代理机构 南京知识律师事务所 32207 代理人 蒋海军 (54) 发明名称 基于机器视觉的大长径比轴类零件直线度误 差检测方法 (57。
2、) 摘要 本发明公开了一种基于机器视觉的大长径比 轴类零件直线度误差检测方法， 属于轴类零件直 线度误差检测领域。本发明的步骤为 ：（1） 标准件 安装 ；（2）相机安装 ；（3） 工件安装 ；（4）图像采 集 ；（5） 图像处理 ；（6） 建立坐标系统 ；（7） 计算得 出圆心位置坐标 ；（8） 统一坐标 ；（9） 根据截面的 圆心坐标， 利用现有的直线度评定方法， 计算直线 度。 本发明以标准件相互间位置关系为纽带， 建立 各图像间的坐标位置关系， 从而避免使用耗时较 长的图像拼接技术， 同时， 不需要严格保证相机到 零件的距离和相机到标准件之间距离相等， 降低 了测量难度， 提高了测量精。
3、度。 (51)Int.Cl. 权利要求书 2 页 说明书 9 页 附图 2 页 (19)中华人民共和国国家知识产权局 (12)发明专利申请 权利要求书 2 页 说明书 9 页 附图 2 页 1/2 页 2 1. 一种基于机器视觉的大长径比轴类零件直线度误差检测方法， 其步骤为 ： (1) 标准件安装 ： 标准件等高安装多组， 所有标准件的对应边相互平行， 两标准件的几何尺寸以及标准 件之间的位置关系采用常规的高精度检测工具测量获得 ； (2) 相机安装 ： 确保成像平面与标准件的上表面平行， 并且相机能沿被测轴的轴线方向和垂直于工作 台面方向移动 ； (3) 工件安装 ： 将被测轴置于视场内，。
4、 调整位置， 使被测轴的轴线与标准件边平行或垂直 ； (4) 图像采集 ： 当相机处于第 i 个采样点处， 调整相机位置， 使得相机到被测轴的距离为 D1， 采集 1 幅图像， 记为 image1i， 调整相机位置， 使得相机到被测轴距离为 D2， 再采集 1 幅图像， 记为 image2i； 沿被测轴轴线方向移动相机到下一个采样点 i+1 处， 调整相机位置， 使得相机到被 测轴的距离为 D1， 采集 1 幅图像， 记为 image1i+1， 调整相机位置， 使得相机到被测轴距离为 D2， 再采集 1 幅图像， 记为 image2i+1， 如此沿轴线方向移动相机， 直至到最后一个采样点， 完。
5、 成整个轴线方向上图像序列的采集， 其中， D1和 D2是两个未知的变量， 但是 D1和 D2之间的 差值是已知量 ； (5) 图像处理 ： image1序列图像中第 i 幅图像处理， 采用 Harris 角点提取算法提取两标准件的 8 个 角点， 用 Canny 边缘提取算法提取轴的轮廓边缘 ； image2序列图像中第 i 幅图像处理， 采用 Harris 角点提取算法提取两标准件的 8 个角点， 用 Canny 边缘提取算法提取轴的轮廓边 缘 ； (6) 建立坐标系统 ： image1图像序列中标准件上最左下角的角点为坐标原点， x 轴平行于轴线， y 轴垂直于 轴线， z 轴垂直于工作。
6、台面 ； (7) 计算圆心位置坐标 ： 在序列图像image1i的被测轴轮廓上取2个点， 保证2个点在垂直于轴线方向上图像坐 标相同， 在序列图像 image2i的被测轴轮廓上取 2 个与 image1i中的点相对应的点 ； 根据小 孔成像原理、 上述所取的 4 个点共圆、 光心与轮廓点所确定的直线和轮廓点与圆心所确定 的直线垂直、 相机移动的距离 |D1-D2|4 个条件， 计算得出圆心位置坐标 ； (8) 统一坐标 ： 根据标准件之间的安装尺寸关系， 将上述计算得到的坐标统一到同一个坐标系下 ； (9) 根据截面的圆心坐标， 利用直线度评定方法， 计算直线度。 2. 根据权利要求 1 所述。
7、的基于机器视觉的大长径比轴类零件直线度误差检测方法， 其 特征在于 ： 步骤(7)中计算圆心位置坐标时建立系统数学模型， 当相机处于1#位置处时， 被 测轴轮廓上的对称点 p9和 p11成像 ； 当相机处于 2# 位置处时， 被测轴轮廓上的对称点 p8和 p10成像， 标准件上坐标原点及 pi(i 1， 2， 3， ， 7) 的 8 个角点间的距离是已知的参数， 再 根据小孔成像原理、 pi(i 8， 9， 10， 11)4 点共圆、 Ocp8 p8o、 Ocp9 p9o、 相机从 1# 位置到 2# 位置移动的距离， 已知这些条件建立系统数学模型如下 ： 权 利 要 求 书 CN 10251。
8、9400 A 2 2/2 页 3 Ow-XwYwZw为世界坐标系， Oc-XcYcZc为相机坐标系， 1#、 2# 为相机所处两个位置， o 为轴 某处截面圆的圆心， pi(i 1， 2， 3， 11) 坐标分别为 ： p1(x1， 0， 0)， p2(x1， y2， 0)， p3(0， y2， 0) ； p4(x4， y4， 0)， p5(x5， y4， 0)， p6(x4， y6， 0)， p7(x5， y6， 0) ； p8(x8， y， z8)， p9(x9， y， z9)， p10(x10， y， z10)， p11(x11， y， z11) ； 当相机处于 1# 位置时， 相机光心。
9、 Oc到圆心 o 距离为 D1， 当相机处于 2# 位置时， 相机光 心 Oc到圆心 o 距离为 D2， 并假设圆心在世界坐标系下坐标为 o(x， y， z)， 则有 ： x9+x11 2x (1.1) z9 z11 (1.2) x8+x10 2x (1.7) z8 z10 (1.8) (x8-x)2+(z8-z)2 (x9-x)2+(z9-z)2 (1.11) D1-D2 C (1.12) 其中， 为直线Ocp9方向，为直线p9o方向， 为直线Ocp8方向，为直线p8o方向 ； C1为 p9点的像到 p11点像之间的距离 ； C2为 p8点的像到 p10点像之间的距离 ； 为当相机处于 1#。
10、 位置处时， 沿 Xc轴方向测量比， 由两标准件平行于 Xc轴的 4 条边的长度与其对应像的长度 的比值所确定， 为当相机处于 1# 位置处时， 沿 Yc轴方向测量比， 由两标准件平行于 Yc轴 的 4 条边的长度与其对应像的长度的比值所确定， 为当相机处于 2# 位置处时， 沿 Xc轴方 向测量比， 由两标准件平行于 Xc轴的 4 条边的长度与其对应像的长度的比值所确定， 为 当相机处于 2# 位置处时， 沿 Yc轴方向测量比， 由两标准件平行于 Yc轴的 4 条边的长度与其 对应像的长度的比值所确定 ； C 为相机从 1# 位置到 2# 位置所移动过的距离 ； f 为相机的焦 距 ； dx。
11、 是像素在 X 方向上的尺寸， dy 是像素在 Y 方向上的尺寸 ； 联立式 (1.1) 式 (1.14)， 即为本测量系统的数学模型， 根据该测量系统的数学模型， 计算得出圆心位置坐标。 权 利 要 求 书 CN 102519400 A 3 1/9 页 4 基于机器视觉的大长径比轴类零件直线度误差检测方法 技术领域 0001 本发明涉及一种轴类零件直线度误差检测方法， 更具体地说， 涉及一种基于机器 视觉的大长径比轴类零件直线度误差检测方法。 背景技术 0002 视觉是人类感知和理解周边的世界， 从大自然中获取知识的最主要的手段。在人 类获取的信息中 80来自视觉， 其余的 20来自听觉、 。
12、味觉、 触觉等。机器视觉， 是指利用 计算机和一些辅助设备来实现人的视觉功能， 从而实现通过二维图像感知外界事物和客观 三维世界， 是一门新兴的学科， 是令人感兴趣的前沿研究领域之一。 0003 随着机器视觉技术的发展， 已经出现了一种把图像当作检测和传递信息的手段或 载体加以利用的新的检测方法， 即基于机器视觉的图像检测技术。 它是以现代光学为基础， 融光电子学、 计算机图像学、 信息处理、 机器视觉等科学技术为一体的现代检测技术。基于 机器视觉的图像检测技术的仪器设备能够实现智能化、 数字化、 小型化、 网络化和多功能 化， 具备在线检测、 实时分析、 实时控制的能力， 在军事、 工业、 。
13、医学等领域得到广泛关注和 应用。 0004 基于机器视觉的图像检测方法在机械零件几何尺寸检测中已经有应用， 主要集中 在基于单幅图像的尺寸测量和基于图像拼接的尺寸测量。 单幅图像检测方法是用一幅图像 来展示整个零件的全貌， 通过图像运算得到零件的几何尺寸。但是由于相机的视野范围受 到约束， 因此方法只适用于较小零件的测量。 对于大型零件， 图像检测技术可以对零件的不 同部位进行拍摄， 得到多幅局部重叠的高精度图像， 然后利用图像之间的信息冗余进行图 像拼接得到零件的完整图像， 对拼接后的图像进行分析可以得到零件的完整结构尺寸。大 量的文献资料已经指出， 这一检测方法需要进行大量的数据存储和运算。
14、， 势必影响检测方 法的实时性。 0005 现有的机械零件直线度测量方法， 主要集中在基于单幅图像的直线度测量。基于 单幅图像的直线度测量方法是用一幅图像展示被测轴的全貌， 通过对这一幅图像的运算处 理得到零件的直线度误差。但是由于相机的视野范围与精度是反比例关系， 为了获得较高 的检测精度， 被测零件的几何尺寸必须约束在一定范围内。天津大学博士后研究工作报告 视觉准则在线测量技术研究 ( 卢荣胜 ) 对直径为 140mm， 长度为 1.5m 的无缝钢管的直线 度检测方法进行了研究， 是国内开展大长径比轴类零件直线度检测方法研究较早的学者。 其报告中采用的方法是一种激光视觉方法， 检测设备精度。
15、较高， 但设备比较昂贵。此外， 需 要对图像传感器进行全局标定， 算法实时性受到一定的限制。 0006 1. 与本发明相关的现有技术一 0007 现有技术一的技术方案 ： 某轨道车辆门用变导程丝杠长约 2m， 直径约为 2cm， 长径 比为100， 属于大长径比轴类零件。 其直线度是质量控制的重要指标， 也是其功能性指标， 从 客户反馈的信息来看， 即使在生产环节进行了检测， 使用一段时间后， 其直线度仍会出现较 大的误差。研究发现， 主要是因为检测手段的精度不够高。 说 明 书 CN 102519400 A 4 2/9 页 5 0008 目前， 生产线在用的检测方法是将丝杠置于工作台面上， 。
16、用眼睛去观察丝杠与工 作台面的接触情况， 透光地方用标准的塞规去检测， 某处能放进塞规， 说明该处的直线度误 差已经 0009 现有技术一的缺点 ： 上述检测方法存在以下几个缺点 ： 0010 (1) 不能准确的给出直线度误差的数据 ； 0011 (2) 不能根据数据制定合理的校正方案 ； 0012 (3) 该岗位对工人的技术熟练程度依赖性比较大 ； 0013 (4) 工人的劳动强度比较大 ； 0014 (5) 检测速度比较慢， 无法实现生产线自动化。 0015 2. 与本发明相关的现有技术二 0016 现有技术二的技术方案 ： 图 1 摘自天津大学卢荣胜博士后工作报告 视觉准则在 线测量技术。
17、研究 ， 该测试系统由半导体激光器 LD 发出一束光平面， 它与无缝钢管相截后， 在钢管的外圆表面形成一条椭圆弧， 由 CCD 摄象机接收各椭圆弧的图像， 经计算机实时图 像处理求得各椭圆弧中心在世界坐标系中的坐标， 运用相应的算法， 求出钢管的直线度误 差。 0017 现有技术二的缺点 ： 图 1 中采用了半导体激光投射器， 该元件属于高精度产品， 价 格比较昂贵。方法中采用多个 CCD 摄像机采集图像， 需要采集的图像统一到同一个坐标系 下方可进行直线度的计算， 因此需要较为复杂的标定算法。 标定是机器视觉中的一个难点， 现有的各种标定技术都有一定的局限性， 在一定精度范围内可用。 如线性。
18、标定计算精度低， 非线性标定需要多次非线性迭代， 计算速度低， 因此标定技术的成为一个瓶颈。因此， 上述 方法又借助了其他高精度设备进行了全局标定， 这就给方法的推广带来了困难。 0018 3. 与本发明相关的现有技术三 0019 现有技术三的技术方案 ： 山东大学学位论文 基于机器视觉的零件平面尺寸自动 测量 ( 马文娟 ) 中写道用机器视觉方法检测零件的直线度， 详细过程如下 ： 图像采集设备 获得零件的全貌， 运用图像处理方法得到灰度图， 对灰度图进行二值化处理， 再用边缘检测 算子， 得到边缘像素级点集， 如果需要进一步提高精度， 再使用亚像素边缘定位技术， 获得 亚像素级边缘点集， 。
19、对边缘点集进行最小二乘拟合得到理想直线， 边缘点集中距离理想直 线距离最大值即被认为是直线度误差。 0020 现有技术三的缺点 ： 基于单幅图像的直线度测量方法是用一幅图像展示被测零件 的全貌， 通过对这一幅图像的运算处理得到零件的直线度误差。但是由于相机的视野范围 与精度是反比例关系， 为了获得零件的全貌， 必须要牺牲一定的精度 ； 又或者为了获得较高 的检测精度， 被测零件的几何尺寸必须约束在一定范围内。 发明内容 0021 1. 发明要解决的技术问题 0022 本发明在于克服现有技术中的不足， 提供一种基于机器视觉的大长径比轴类零件 直线度误差检测方法， 通过快速获得被测轴各个局部区域的。
20、高精度图像， 采用简单的图像 处理算法得到多个位置处的中心坐标， 再将相互独立的中心点统一到一个坐标系下进行直 线度误差的评价， 以得到整个轴直线度误差值， 并对超差部分进行记录， 为后续校正提供数 说 明 书 CN 102519400 A 5 3/9 页 6 据。 0023 2. 技术方案 0024 为达到上述目的， 本发明提供的技术方案为 ： 0025 本发明的一种基于机器视觉的大长径比轴类零件直线度误差检测方法， 其步骤 为 ： 0026 (1) 标准件安装 ： 0027 标准件等高安装多组， 所有标准件的对应边相互平行， 两标准件的几何尺寸以及 标准件之间的位置关系采用常规的高精度检测。
21、工具测量获得 ； 0028 (2) 相机安装 ： 0029 确保成像平面与标准件的上表面平行， 并且相机能沿被测轴的轴线方向和垂直于 工作台面方向移动 ； 0030 (3) 工件安装 ： 0031 将被测轴置于视场内， 调整位置， 使被测轴的轴线与标准件边平行或垂直 ； 0032 (4) 图像采集 ： 0033 当相机处于第 i 个采样点处， 调整相机位置， 使得相机到被测轴的距离为 D1， 采集 1 幅图像， 记为 image1i， 调整相机位置， 使得相机到被测轴距离为 D2， 再采集 1 幅图像， 记为 image2i； 沿被测轴轴线方向移动相机到下一个采样点 i+1 处， 调整相机位置。
22、， 使得相机到被 测轴的距离为 D1， 采集 1 幅图像， 记为 image1i+1， 调整相机位置， 使得相机到被测轴距离为 D2， 再采集 1 幅图像， 记为 image2i+1， 如此沿轴线方向移动相机， 直至到最后一个采样点， 完 成整个轴线方向上图像序列的采集， 其中， D1和 D2是两个未知的变量， 但是 D1和 D2之间的 差值是已知量 ； 0034 (5) 图像处理 ： 0035 image1序列图像中第 i 幅图像处理， 采用 Harris 角点提取算法提取两标准件的 8 个角点， 用 Canny 边缘提取算法提取轴的轮廓边缘 ； image2序列图像中第 i 幅图像处理， 。
23、采 用 Harris 角点提取算法提取两标准件的 8 个角点， 用 Canny 边缘提取算法提取轴的轮廓边 缘 ； 0036 (6) 建立坐标系统 ： 0037 image1图像序列中标准件上最左下角的角点为坐标原点， x 轴平行于轴线， y 轴垂 直于轴线， z 轴垂直于工作台面 ； 0038 (7) 计算圆心位置坐标 ： 0039 在序列图像image1i的被测轴轮廓上取2个点， 保证2个点在垂直于轴线方向上图 像坐标相同， 在序列图像 image2i的被测轴轮廓上取 2 个与 image1i中的点相对应的点 ； 根 据小孔成像原理、 上述所取的 4 个点共圆、 光心与轮廓点所确定的直线和。
24、轮廓点与圆心所 确定的直线垂直、 相机移动的距离 (|D1-D2|)4 个条件， 计算得出圆心位置坐标 ； 0040 其中， 计算圆心位置坐标时建立系统数学模型， 当相机处于 1# 位置处时， 被测轴 轮廓上的对称点p9和p11成像 ； 当相机处于2#位置处时， 被测轴轮廓上的对称点p8和p10成 像， 标准件上的 8 个角点 ( 坐标原点及 pi(i 1， 2， 3， ， 7) 间的距离是已知的参数， 再根 据小孔成像原理、 pi(i 8， 9， 10， 11)4 点共圆、 Ocp8 p8o、 Ocp9 p9o、 相机移动的距离 (1# 位置到 2# 位置距离 ) 已知这些条件建立系统数学模。
25、型如下 ： 说 明 书 CN 102519400 A 6 4/9 页 7 0041 Ow-XwYwZw为世界坐标系， Oc-XcYcZc为相机坐标系， 1#、 2# 为相机所处两个位置， o 为 轴某处截面圆的圆心， pi(i 1， 2， 3， ， 11) 坐标分别为 ： p1(x1， 0， 0)， p2(x1， y2， 0)， p3(0， y2， 0) ； p4(x4， y4， 0)， p5(x5， y4， 0)， p6(x4， y6， 0)， p7(x5， y6， 0) ； p8(x8， y， z8)， p9(x9， y， z9)， p10(x10， y， z10)， p11(x11， y。
26、， z11) ； 0042 当相机处于 1# 位置时， 相机光心 Oc到圆心 o 距离为 D1， 当相机处于 2# 位置时， 相 机光心 Oc到圆心 o 距离为 D2， 并假设圆心在世界坐标系下坐标为 o(x， y， z)， 则有 ： 0043 x9+x11 2x (1.1) 0044 z9 z11 (1.2) 0045 0046 0047 0048 0049 x8+x10 2x (1.7) 0050 z8 z10 (1.8) 0051 0052 0053 (x8-x)2+(z8-z)2 (x9-x)2+(z9-z)2 (1.11) 0054 D1-D2 C (1.12) 0055 0056 。
27、0057 其中， 为直线Ocp9方向，为直线p9o方向， 为直线Ocp8方向， 为直线p8o方向 ； C1为 p9点的像到 p11点像之间的距离 ； C2为 p8点的像到 p10点像之间的距离 ； 为当相机处 于 1# 位置处时， 沿 Xc轴方向测量比， 由两标准件平行于 Xc轴的 4 条边的长度与其对应像的 长度的比值所确定， 为当相机处于1#位置处时， 沿Yc轴方向测量比， 由两标准件平行于Yc 轴的 4 条边的长度与其对应像的长度的比值所确定， 为当相机处于 2# 位置处时， 沿 Xc轴 方向测量比， 由两标准件平行于 Xc轴的 4 条边的长度与其对应像的长度的比值所确定， 为当相机处于。
28、 2# 位置处时， 沿 Yc轴方向测量比， 由两标准件平行于 Yc轴的 4 条边的长度与 其对应像的长度的比值所确定 ； C 为相机从 1# 位置到 2# 位置所移动过的距离 ； f 为相机的 焦距 ； dx 是像素在 X 方向上的尺寸， dy 是像素在 Y 方向上的尺寸 ； 0058 联立式 (1.1) 式 (1.14)， 即为本测量系统的数学模型， 根据该测量系统的数学 说 明 书 CN 102519400 A 7 5/9 页 8 模型， 计算得出圆心位置坐标 ； 0059 (8) 统一坐标 ： 0060 根据标准件之间的安装尺寸关系， 将上述计算得到的坐标统一到同一个坐标系 下 ； 00。
29、61 (9) 根据截面的圆心坐标， 利用直线度评定方法， 计算直线度。 0062 3. 有益效果 0063 采用本发明提供的技术方案， 与已有的公知技术相比， 具有如下显著效果 ： 0064 (1) 本发明的一种基于机器视觉的大长径比轴类零件直线度误差检测方法， 对照 现有技术一， 本发明用图像作为检测和传递信息的载体， 存储于计算机内， 可以获得准确的 直线度误差数据， 根据获得的直线度误差， 制定合理的校正方案， 机器视觉长时间稳定工 作， 人类难以长时间对同一对象进行观察， 而机器视觉则可以长时间地作测量、 分析和识别 任务 ； 0065 (2) 本发明的一种基于机器视觉的大长径比轴类零。
30、件直线度误差检测方法， 对照 现有技术二， 本发明属于普通视觉， 不需要采用昂贵的高精度激光投射器， 机器视觉尺寸测 量系统的标定方法大都采用标准件法， 即把标准件的精确尺寸传递给图像， 这种方法的主 要优点在于不需要标定摄像机的内外参数， 只需要标定出摄像机的物面分辨率， 可以简化 标定过程， 该方法只适用于测量平面物体尺寸， 只要保证每次采集图像时， 视场内同时有两 个标准件， 即可根据两标准件间相互位置关系将全部坐标统一起来， 这种标定算法不需要 进行多次非线性迭代计算， 计算速度较高， 精度主要取决标准件精度和相机的成像精度， 可 根据精度等级要求配置标准件和摄像机 ； 0066 (3。
31、) 本发明以图像为检测和传递信息载体， 可以获得用数据准确描述的直线度误 差数据， 数据存储于计算上， 可以和设计、 生产以及产品质量控制等多个部门适时沟通， 以 提升产品的质量 ； 0067 (4) 本发明的一种基于机器视觉的大长径比轴类零件直线度误差检测方法， 不需 要耗时较长的图像拼接技术， 得到零件全貌， 也不要进行十分复杂且需要高精度设备支持 的全局坐标标定， 得到轴线各点绝对坐标， 进行直线度的计算 ； 0068 (5) 本发明的检测精度只与相机的成像精度和标准件精度有关， 相机沿被测轴轴 线方向移动与定位不需要做特别的处理 ； 0069 (6) 本发明可以实现任意长径比轴类零件直。
32、线度检测任务， 只需要适当延长工作 台的尺寸， 增加一定数量的标准件即可。 附图说明 0070 图 1 为天津大学卢荣胜博士后工作报告 视觉准则在线测量技术研究 中的测试 系统 ； 0071 图 2 是本发明的测量系统安装示意图 ； 0072 图 3 为本发明的测量系统模型图。 0073 示意图中的标号说明 ： 0074 1、 3- 标准件 ； 2- 被测轴。 说 明 书 CN 102519400 A 8 6/9 页 9 具体实施方式 0075 下面结合附图对本发明作进一步的描述。 0076 结合图 2 和图 3 对本发明技术方案进行阐述。图 2 中的 “1” 和 “3” 是 2 个标准 件，。
33、 其尺寸不相同， 图 2 中的 “2” 为某被测轴。测量前要进行相机、 被测轴 2 以及标准件 1、 3 的安装。标准件 1 和标准件 3 等高安装， 且保证边的平行 ( 或垂直 )， 两标准件 1、 3 的几 何尺寸以及标准件 1、 3 之间的位置关系用常规的高精度检测工具测量获得。相机安装需要 确保成像平面与标准件 1、 3 的上表面平行， 并且相机能沿被测轴 2 的轴线方向和垂直于工 作台面方向移动。被测轴 2 安装于视场内， 调整其位置， 保证被测轴 2 的轴线与标准件 1、 3 的边平行 ( 或垂直 )。当相机处于第 i 个采样点处， 调整相机位置， 使得相机到被测轴 2 的 距离为。
34、 D1， 采集 1 幅图像， 记为 image1i； 调整相机位置， 使得相机到被测轴 2 距离为 D2， 再采 集1幅图像， 记为image2i。 沿被测轴2轴线方向移动相机到下一个采样点i+1处， 调整相机 位置， 使得相机到被测轴 2 的距离为 D1， 采集 1 幅图像， 记为 image1i+1； 调整相机位置， 使得 相机到被测轴 2 距离为 D2， 再采集 1 幅图像， 记为 image2i+1。如此沿轴线方向移动相机， 直 至到最后一个采样点， 完成整个轴线方向上图像序列的采集。D1和 D2是两个未知的变量， 但是 D1和 D2之间的差值是已知量。对 image1序列图像中第 i。
35、 幅图像处理， 用 Harris 角 点 (C.Harris， etc.1998) 提取算法提取两标准件 1、 3 的 8 个角点， 用 Canny(John Harris， 1986) 边缘提取算法提取轴的轮廓边缘。对 image2序列图像中第 i 幅图像处理， 用 Harris 角点提取算法提取两标准件1、 3的8个角点， 用Canny边缘提取算法提取轴的轮廓边缘。 坐 标系统为 ： image1图像序列中标准件 1、 3 上最左下角的角点为坐标原点， x 轴垂直于轴线， y 轴平行于轴线， z 轴垂直于工作台面。在序列图像 image1i的被测轴 2 轮廓上取 2 个点， 保证 2 个点。
36、在垂直于轴线方向上图像坐标相同。在序列图像 image2i的被测轴 2 轮廓上取 2 个与 image1i中点相对应的点。 0077 图 3 为测量系统模型图， 当相机处于 1# 位置处时， 被测轴 2 轮廓上的点 p9和 p11 成像。当相机处于 2# 位置处时， 被测轴 2 轮廓上的点 p8和 p10成像。标准件 1、 3 上的 8 个 角点 ( 坐标原点及 pi(i 1， 2， 3， 7) 间的距离是已知的参数。再根据小孔成像原理、 pi(i 8， 9， 10， 11)4 点共圆、 直线 Ocp8与 p8o 垂直和直线 Ocp9与 p9o 也垂直、 相机移动的距 离 (1# 位置到 2#。
37、 位置距离 ) 已知， 计算得出圆心位置坐标。利用标准件 1、 3 之间的安装尺 寸关系， 将计算得到的圆心坐标统一到同一个坐标系下。 根据每个截面的圆心坐标， 运用现 有的直线度评定方法， 计算得出该被测轴 2 的直线度。 0078 下面结合实施例对本发明作进一步的描述。 0079 实施例 ： 0080 本实施例的一种基于机器视觉的大长径比轴类零件直线度误差检测方法， 其步骤 为 ： 0081 (1) 标准件安装 ： 0082 标准件等高安装多组， 图 2 所示，“1” 和 “3” 是 2 个标准件， 其尺寸不相同， 图 2 中 的 “2” 为某被测轴， 所有标准件的对应边相互平行， 两标准。
38、件的几何尺寸以及标准件之间的 位置关系采用常规的高精度检测工具测量获得 ； 0083 (2) 相机安装 ： 0084 确保成像平面与标准件的上表面平行， 并且相机能沿被测轴的轴线方向和垂直于 说 明 书 CN 102519400 A 9 7/9 页 10 工作台面方向移动 ； 0085 (3) 工件安装 ： 0086 将被测轴置于视场内， 调整位置， 使被测轴的轴线与标准件边平行或垂直 ； 0087 (4) 图像采集 ： 0088 当相机处于第 i 个采样点处， 调整相机位置， 使得相机到被测轴的距离为 D1， 采集 1 幅图像， 记为 image1i， 调整相机位置， 使得相机到被测轴距离为。
39、 D2， 再采集 1 幅图像， 记为 image2i； 沿被测轴轴线方向移动相机到下一个采样点 i+1 处， 调整相机位置， 使得相机到被 测轴的距离为 D1， 采集 1 幅图像， 记为 image1i+1， 调整相机位置， 使得相机到被测轴距离为 D2， 再采集 1 幅图像， 记为 image2i+1， 如此沿轴线方向移动相机， 直至到最后一个采样点， 完 成整个轴线方向上图像序列的采集， 其中， D1和 D2是两个未知的变量， 但是 D1和 D2之间的 差值是已知量 ； 0089 (5) 图像处理 ： 0090 image1序列图像中第 i 幅图像处理， 采用 Harris 角点提取算法提。
40、取两标准件的 8 个角点， 用 Canny 边缘提取算法提取轴的轮廓边缘 ； image2序列图像中第 i 幅图像处理， 采 用 Harris 角点提取算法提取两标准件的 8 个角点， 用 Canny 边缘提取算法提取轴的轮廓边 缘 ； 0091 (6) 建立坐标系统 ： 0092 image1图像序列中标准件上最左下角的角点为坐标原点， x 轴平行于轴线， y 轴垂 直于轴线， z 轴垂直于工作台面 ； 0093 (7) 计算圆心位置坐标 ： 0094 在序列图像image1i的被测轴轮廓上取2个点， 保证2个点在垂直于轴线方向上图 像坐标相同， 在序列图像 image2i的被测轴轮廓上取 。
41、2 个与 image1i中的点相对应的点 ； 根 据小孔成像原理、 上述所取的 4 个点共圆、 光心与轮廓点所确定的直线和轮廓点与圆心所 确定的直线垂直、 相机移动的距离 (|D1-D2|)4 个条件， 计算得出圆心位置坐标 ； 0095 其中， 计算圆心位置坐标时建立如图 3 所示的系统数学模型， 当相机处于 1# 位置 处时， 被测轴轮廓上的对称点 p9和 p11成像 ； 当相机处于 2# 位置处时， 被测轴轮廓上的对称 点 p8和 p10成像， 标准件上的 8 个角点 ( 坐标原点及 pi(i 1， 2， 3， 7) 间的距离是已 知的参数， 再根据小孔成像原理、 pi(i 8， 9， 。
42、10， 11)4 点共圆、 Ocp8 p8o、 Ocp9 p9o、 相机 移动的距离 (1# 位置到 2# 位置距离 ) 已知这些条件建立系统数学模型如下 ： 0096 Ow-XwYwZw为世界坐标系， Oc-XcYcZc为相机坐标系， 1#、 2# 为相机所处两个位置， o 为 轴某处截面圆的圆心， pi(i 1， 2， 3， ， 11) 坐标分别为 ： p1(x1， 0， 0)， p2(x1， y2， 0)， p3(0， y2， 0) ； p4(x4， y4， 0)， p5(x5， y4， 0)， p6(x4， y6， 0)， p7(x5， y6， 0) ； p8(x8， y， z8)， 。
43、p9(x9， y， z9)， p10(x10， y， z10)， p11(x11， y， z11) ； 0097 当相机处于 1# 位置时， 相机光心 Oc到圆心 o 距离为 D1， 当相机处于 2# 位置时， 相 机光心 Oc到圆心 o 距离为 D2， 并假设圆心在世界坐标系下坐标为 o(x， y， z)， 则有 ： 0098 x9+x11 2x (1.1) 0099 z9 z11 (1.2) 0100 说 明 书 CN 102519400 A 10 8/9 页 11 0101 0102 0103 0104 x8+x10 2x (1.7) 0105 z8 z10 (1.8) 0106 010。
44、7 0108 (x8-x)2+(z8-z)2 (x9-x)2+(z9-z)2 (1.11) 0109 D1-D2 C (1.12) 0110 0111 0112 其中， 为直线Ocp9方向，为直线p9o方向， 为直线Ocp8方向， 为直线p8o方向 ； C1为 p9点的像到 p11点像之间的距离 ； C2为 p8点的像到 p10点像之间的距离 ； 为当相机处 于 1# 位置处时， 沿 Xc轴方向测量比， 由两标准件平行于 Xc轴的 4 条边的长度与其对应像的 长度的比值所确定， 为当相机处于1#位置处时， 沿Yc轴方向测量比， 由两标准件平行于Yc 轴的 4 条边的长度与其对应像的长度的比值所。
45、确定， 为当相机处于 2# 位置处时， 沿 Xc轴 方向测量比， 由两标准件平行于 Xc轴的 4 条边的长度与其对应像的长度的比值所确定， 为当相机处于 2# 位置处时， 沿 Yc轴方向测量比， 由两标准件平行于 Yc轴的 4 条边的长度与 其对应像的长度的比值所确定 ； C 为相机从 1# 位置到 2# 位置所移动过的距离 ； f 为相机的 焦距 ； dx 是像素在 X 方向上的尺寸， dy 是像素在 Y 方向上的尺寸 ； 0113 联立式 (1.1) 式 (1.14)， 即为本测量系统的数学模型， 根据该测量系统的数学 模型， 计算得出圆心位置坐标 ； 0114 (8) 统一坐标 ： 01。
46、15 根据标准件之间的安装尺寸关系， 将上述计算得到的坐标统一到同一个坐标系 下 ； 0116 (9) 根据截面的圆心坐标， 运用现有的直线度评定方法， 本实施例中采用运用 3DLSA( 胡仲勋等 ) 直线度评定算法， 计算得出该被测轴的直线度。 0117 本发明的基于机器视觉的大长径比轴类零件直线度误差检测方法， 是一种全新的 基于机器视觉实现大长径比轴类零件直线度误差在线测量方法， 其独特之处在于 ： 以图像 为检测和传递信息载体， 可以获得用数据准确描述的直线度误差数据， 数据存储于计算上， 可以和设计、 生产以及产品质量控制等多个部门适时沟通， 以提升产品的质量。 本发明的方 说 明 书 CN 102519400 A 11 9/9 页 12 法中不使用昂贵的高精度的激光投射器。以标准件相互间位置关系为纽带， 建立各图像间 的坐标位置关系， 从而避免使用耗时较长的图像拼接技术。在测量系统数学模型的建立过 程中， 不直接使用测量比法 ( 或称之为当量像素法 )， 因此不需要严格保证相机到零件的距 离和相机到标准件之间距离相等， 降低测量难度， 提高测量精度。 说 明 书 CN 102519400 A 12 1/2 页 13 图 1 图 2 说 明 书 附 图 CN 102519400 A 13 2/2 页 14 图 3 说 明 书 附 图 CN 102519400 A 14 。