CN102445174B - 一种基于支持向量回归的多测点平面度评定方法 - Google Patents
一种基于支持向量回归的多测点平面度评定方法 Download PDFInfo
- Publication number
- CN102445174B CN102445174B CN 201110311170 CN201110311170A CN102445174B CN 102445174 B CN102445174 B CN 102445174B CN 201110311170 CN201110311170 CN 201110311170 CN 201110311170 A CN201110311170 A CN 201110311170A CN 102445174 B CN102445174 B CN 102445174B
- Authority
- CN
- China
- Prior art keywords
- support vector
- measurement point
- point
- plane
- vector regression
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 41
- 239000013598 vector Substances 0.000 title claims abstract description 33
- 238000005259 measurement Methods 0.000 claims description 36
- 238000012886 linear function Methods 0.000 claims description 3
- 238000000605 extraction Methods 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 abstract description 2
- 238000005070 sampling Methods 0.000 abstract description 2
- 238000007781 pre-processing Methods 0.000 abstract 1
- 238000011156 evaluation Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004441 surface measurement Methods 0.000 description 1
Images
Landscapes
- Length Measuring Devices With Unspecified Measuring Means (AREA)
Abstract
本发明公开了一种基于支持向量回归的多测点平面度评定方法,包括:对被测平板表面的点进行采样,获取各点的三坐标测量值;对原始测量点进行预处理,通过计算测量点集的三维凸壳,剔除凸壳内的测量点,保留凸壳上的测量点;采用支持向量回归法——ε-SVR求出测量点集的最小包容区域平面,其对应的包容区域宽度即为所求平面度。通过本发明可以减少多测点平面度评定中的计算数据量,提高计算效率和评定结果的准确度。
Description
技术领域
本发明涉及一种平面度的评定方法,尤其涉及一种基于支持向量回归的多测点平面度评定方法。
背景技术
随着国家经济的发展和先进制造水平的提高,在现代化工业制造中特别是在精密制造等领域,对工件和结构的几何尺寸精度都有较高的精度要求。平面作为工业产品一种最常见的几何尺寸,在实际应用中多作为承压平面和基准平面。如果工件的平面度一旦超出了设计要求,将会导致承压面的受力不均,容易引起承压面局部变形或损坏,或是引起以该平面为基准的部件的几何尺寸失准,最终会影响整个结构的功能、寿命和安全性等。
因此,为了满足工业产品平面设计的精度要求,不仅需要精密的加工手段,同时更需要高准确度的形位误差评定。准确的评定一方面能够准确判断工件是否合格,另一方面可以反过来指导加工,避免不必要的重复加工而带来的资源浪费。因此,平面度的准确评定有着非常重要的意义。但是,为了能够尽可能真实准确地反映表面的情况,必然要对被测表面进行充分采样,尽可能多地获取表面测量点。在这种情况下,如果用传统方法评定平面度,则很难同时达到高准确度和高计算速度,因此,如何能够快速准确地对多测量点的平面度进行评定成为一个重要研究问题。
目前,平面度评定的方法有最小二乘法,搜索法、计算几何法等。最小二乘法计算速度较快,但使用最小二乘平面代替最小区域平面,评定准确度难以保证;搜索法计算准确度较高,但在迭代搜索中容易陷入局部最优;计算几何法计算准确度高,但是需要通过全部枚举发求得最小区域平面,计算复杂。
发明内容
为解决上述中存在的问题与缺陷,本发明提供了一种基于支持向量回归的多测点平面度评定方法,降低计算复杂度、提高计算效率和评定结果的准确度。
本发明是通过以下技术方案实现的:
本发明所涉及的一种基于支持向量回归的多测点平面度评定方法,该方法包括以下步骤:
对被测平板表面的点进行采样,获取各点的三维坐标测量值;
对测量点进行预处理,计算所有测量点对应的三维凸壳,剔除凸壳内的测量点,保留凸壳上的测量点;
采用支持向量回归法——ε-SVR求出测量点集的最小包容区域平面;
计算各测量点到最小包容区域平面的距离,最大距离与最小距离之差为所求平面度。
本发明提供的技术方案的有益效果是:
利用凸壳算法对原始测量点集进行预处理,剔除与最小区域平面无关的测量点,保留与最小区域平面相关的测量点,从而减少了计算数据量,提高了计算效率;支持向量回归方法与最小包容区域的平面度定义在机理上一致,利用该方法求测量点集的平面度,完全符合最小条件的平面度的定义,因此提高了评定结果的准确度。
附图说明
图1为基于支持向量回归的多测点平面度评定方法流程图;
图2为测量点集预处理方法示意图;
图3为ε-SVR求最小包容区域平面的流程图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明实施方式作进一步地详细描述:
本实施例提供了一种基于支持向量回归的多测点平面度评定方法。
参见图1,该方法主要包括以下步骤:
步骤10对被测对象表面采样,获取所有测量点的三维坐标;
步骤20求测量点集的凸壳,剔除凸壳内的测量点,保留凸壳上的测量点;
步骤30采用支持向量回归法ε-SVR求出测量点集的最小包容区域平面;
步骤40计算各测量点到最小包容区域平面的距离,最大距离与最小距离之差即为所求平面度。
参见图2所示:假设原始测量点集为D0={pi(xi,yi,zi),i=1,…,l}(如图2-a),利用快速凸壳算法求出D0凸壳面CH(如图2-b),提取位于CH表面的测量点pi(xi,yi,zi)∈CH,剔除位于CH内的测量点(如图2-c),组成了新的测量点集D1={pi(xi,yi,zi),pi∈CH,i=1,…,m}(如图2-d)。凸壳CH是原始测量点集D0的一个最小的凸集和,它所构成的包容区域即是包容D0的一个最小区域,因此凸壳CH上的测量点集D1所对应的最小包容区域即是原始测量点集D0的最小包容区域。
上述利用ε-SVR求测量点集的最小包容区域平面的方法,是建立在支持向量回归方法和基于最小条件的平面度定义的机理一致基础之上,二者都是寻求包容区域最小,因此可以直接应用支持向量回归法求解平面度。又因为平面度的最小包容区域面是一个平面,因此需要将支持向量回归的核函数设为线性函数,另外,根据平面度的1-3和2-2判定原则,即最小包容平面是由4个点决定,因此其支持向量测量点的个数为4。此方法的具体步骤为(参见图3所示):将新测量点集D1={(xi,yi,zi),i=1,…,m},首先初始化不敏感函数参数ε0,选用线性函数K(Vi,Vj)=Vi·Vj作为支持向量回归法的核函数,并采用支持向量回归法 ——ε-SVR求测量点集D1的回归超平面设所求得的支持向量元素的个数为如果则该支持向量元素对应的这四个测量点即为D1中构成最小包容区域的四个极值点,即为所有测量点集的最小区域平面ΠMZ。
如果nSV≠4,则采用下山单纯形法进行迭代搜索εi。在每次迭代中都利用支持向量回归法ε-SVR求在不敏感函数参数εi的条件下,测量点集D1的回归超平面Πεi和支持向量个数并将作为目标函数值。最后,通过m次迭代搜索,得到时的不敏感函数参数εm,此时求得的回归平面即为最小区域平面ΠMZ。
求出D1中各测量点到最小区域平面ΠMZ的距离为di,则被测对象表面的平面度值fMZ=dmax-dmin。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。
Claims (1)
1.一种基于支持向量回归的多测点平面度评定方法,其特征在于,该方法包括以下步骤:
对被测平板表面的点进行采样,获取各点的三维坐标测量值;
对测量点进行预处理,计算所有测量点对应的三维凸壳,剔除凸壳内的测量点,保留凸壳上的测量点;
采用支持向量回归法——ε-SVR求出测量点集的最小包容区域平面;
计算各测量点到最小包容区域平面的距离,最大距离与最小距离之差为所求平面度;
所述对测量点预处理的具体步骤为:
设原始测量点集为D0={pi(xi,yi,zi),i=1,…,l},利用快速凸壳算法求出D0凸壳面CH,剔除位于CH内的测量点提取位于CH表面的测量点pi(xi,yi,zi)∈CH,从而组成新测量点集D1={pi(xi,yi,zi),pi∈CH,i=1,…,m};
所述ε-SVR求出测量点集的最小包容区域平面具体步骤为:
初始化不敏感函数参数ε0,选用线性函数K(Vi,Vj)=Vi·Vj作为支持向量回归法的核函数,利用支持向量回归法——ε-SVR求测量点集D1的回归超平面其中Vi、Vj表示两个不同的样本向量;设所求得的支持向量元素的个数为如果则该支持向量元素对应的测量点即为D1的极值点,即为所求的最小区域平面ΠMZ;
如果nSV≠4,则采用下山单纯形法进行迭代搜索εi,在每次迭代中利用支持向量回归法ε-SVR求在不同的不敏感函数参数εi的条件下,测量点集D1的回归超平面Πεi和支持向量个数并将作为作为目标函数值;通过m次迭代搜索,得到时的不敏感函数参数εm,此时的回归平面即为所求最小区域平面ΠMZ;
求出D1中各测量点到最小区域平面ΠMZ的距离为di,则被测对象表面的平面度值fMZ=dmax-dmin;所述nSV表示采用ε-SVR方法求出的支持向量的个数;dmax和dmin表示测量点到最小区域平面ΠMZ距离di的最大距离和最小距离。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110311170 CN102445174B (zh) | 2011-10-14 | 2011-10-14 | 一种基于支持向量回归的多测点平面度评定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110311170 CN102445174B (zh) | 2011-10-14 | 2011-10-14 | 一种基于支持向量回归的多测点平面度评定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102445174A CN102445174A (zh) | 2012-05-09 |
CN102445174B true CN102445174B (zh) | 2013-10-02 |
Family
ID=46007904
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110311170 Expired - Fee Related CN102445174B (zh) | 2011-10-14 | 2011-10-14 | 一种基于支持向量回归的多测点平面度评定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102445174B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102749064B (zh) * | 2012-06-19 | 2015-10-14 | 歌尔声学股份有限公司 | 一种测量液晶显示装置的背板平面度的方法及其装置 |
CN102914270B (zh) * | 2012-10-12 | 2015-12-02 | 西安理工大学 | 基于支持向量机回归的晶体直径测量方法 |
CN103256916B (zh) * | 2013-06-10 | 2015-07-08 | 陈磊磊 | 一种基于最小区域的零件平面度误差评定方法 |
CN105841640B (zh) * | 2016-04-29 | 2018-09-18 | 北京航空航天大学 | 平面度误差评定方法及装置 |
CN107798324B (zh) * | 2016-08-29 | 2020-12-29 | 北京君正集成电路股份有限公司 | 一种车牌图像定位方法及设备 |
CN108645369A (zh) * | 2018-05-16 | 2018-10-12 | 大连大学 | 相对基准要素轴线的平行度基准方向判断方法 |
CN108775883B (zh) * | 2018-06-30 | 2020-06-02 | 北京动力机械研究所 | 一种叶轮类零件快速换装精度在线检测方法 |
CN111177645B (zh) * | 2019-12-26 | 2023-08-29 | 哈尔滨工业大学 | 一种基于大规模点云数据的大型高速回转装备误差混合评定方法 |
CN111721255A (zh) * | 2020-06-30 | 2020-09-29 | 中建材轻工业自动化研究所有限公司 | 一种平面度检测方法及*** |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101561495A (zh) * | 2009-06-01 | 2009-10-21 | 长讯通信服务有限公司 | 一种无线传感器网络节点三维定位方法 |
CN101726251A (zh) * | 2009-11-13 | 2010-06-09 | 江苏大学 | 基于支持向量机的苹果采摘机器人果实自动识别方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002185103A (ja) * | 2000-12-18 | 2002-06-28 | Kyocera Corp | 実装用基板の実装面における電極パッドの平坦性評価方法 |
-
2011
- 2011-10-14 CN CN 201110311170 patent/CN102445174B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101561495A (zh) * | 2009-06-01 | 2009-10-21 | 长讯通信服务有限公司 | 一种无线传感器网络节点三维定位方法 |
CN101726251A (zh) * | 2009-11-13 | 2010-06-09 | 江苏大学 | 基于支持向量机的苹果采摘机器人果实自动识别方法 |
Non-Patent Citations (4)
Title |
---|
JP特开2002-185103A 2002.06.28 |
大尺寸平面度的智能评价;王时礼;《全国优秀硕士学位论文全文数据库 工程科技I辑》;20110630(第6期);正文第1-9页 * |
王时礼.大尺寸平面度的智能评价.《全国优秀硕士学位论文全文数据库 工程科技I辑》.2011,(第6期),全文. |
薛小强.基于近似凸壳的直线度误差评价方法.《南京工程学院学报(自然科学版)》.2009,第7卷(第1期),20-24. * |
Also Published As
Publication number | Publication date |
---|---|
CN102445174A (zh) | 2012-05-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102445174B (zh) | 一种基于支持向量回归的多测点平面度评定方法 | |
CN102785129B (zh) | 复杂零件的曲面加工精度的在线检测方法 | |
CN102506805B (zh) | 一种基于支持向量分类的多测点平面度评定方法 | |
Chen et al. | Spatial statistical analysis and compensation of machining errors for complex surfaces | |
CN101975560B (zh) | 一种面阵ccd靶面与安装定位面平行度的光学检测方法 | |
CN110285781B (zh) | 一种相对于基准面的平面平行度快速评定方法 | |
CN102252638A (zh) | 用于测量超大平面平面度的数据拼接技术 | |
CN101907448A (zh) | 一种基于双目三维视觉的深度测量方法 | |
CN103256916A (zh) | 一种基于最小区域的零件平面度误差评定方法 | |
CN103499318B (zh) | 一种光学元件自重变形量的测定方法 | |
CN103439085B (zh) | 一种接触式测量曲面棱镜参数的方法及装置 | |
CN103292773A (zh) | 一种基于最小区域的对称度误差评定方法 | |
CN104778181A (zh) | 一种进行测量光谱与库光谱匹配的方法及其设备 | |
CN104915473A (zh) | 一种测量晶体表面结构之间距离的方法 | |
CN104655496B (zh) | 自重对离轴反射镜面形影响的测试方法 | |
CN102410822B (zh) | 厚度量测装置及其方法 | |
CN104751005A (zh) | 一种基于正交实验的平面度误差评定方法 | |
Cheung et al. | Measuring ultra-precision freeform surfaces using a robust form characterization method | |
CN104279956A (zh) | 一种岩石结构面表面基准面的确定方法 | |
CN205785048U (zh) | 自动识别半径补偿方向的三坐标测量机专用测头 | |
CN103310106B (zh) | 一种零件孔系作用尺寸的计算方法 | |
Wang et al. | A fiducial-aided data fusion method for the measurement of multiscale complex surfaces | |
CN112833790A (zh) | 减小相位光栅非对称性对位置测量精度影响的方法 | |
Ye et al. | Minimum zone evaluation of flatness error using an adaptive iterative strategy for coordinate measuring machines data | |
Cong et al. | Research on the method of improving the accuracy of CMM (coordinate measuring machine) testing aspheric surface |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20131002 |