CN113513978A - 高低温环境下的端面位姿相对变化高精度测量方法和系统 - Google Patents
高低温环境下的端面位姿相对变化高精度测量方法和系统 Download PDFInfo
- Publication number
- CN113513978A CN113513978A CN202110612495.3A CN202110612495A CN113513978A CN 113513978 A CN113513978 A CN 113513978A CN 202110612495 A CN202110612495 A CN 202110612495A CN 113513978 A CN113513978 A CN 113513978A
- Authority
- CN
- China
- Prior art keywords
- face
- measurement
- pose
- angle change
- temperature working
- 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.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
一种高低温环境下的端面位姿相对变化高精度测量方法,包括步骤如下:步骤1,计算被测端面每次测量的位姿初始函数;步骤2,进行被测端面采集点云噪点剔除,筛选确认被测端面每次测量的最终位姿状态函数;步骤3,计算被测端面各温度工况下位姿变化结果的初始结果;步骤4,进行端面位姿变化测量数据筛选,确认各温度工况下的端面相对位姿变化最终测量结果。本发明的方法和系统能够有效的剔除测量中的端面噪点及奇异数据,进一步减小测量误差提高测量精度。
Description
技术领域
本发明涉及一种端面位姿相对变化高精度测量方法和系统。
背景技术
目前航空、航天等诸多高精尖工业领域内,受外界高低温变化环境影响许多设备产品等材料件结构件的各种端面均会因热变形及自身因素发生变形情况使端面的位姿发生变化,如何准确合理的测量被测端面在外界环境加载条件下的端面位姿变化是全面了解评估被测件性能的关键。目前常见的端面热变形测量算法是利用端面微变形测量系统例如:数字图像相关测量系统、摄影测量系统、激光多普勒测量系统、三坐标机等获取被测端面上多点的位置坐标信息,进一步利用多点坐标信息通过数据处理及相关端面形状拟合方法直接分析出端面的函数模型再根据每次测量的函数模型按温度工况区分后计算变形均值作为最终测量结果。由于现有的测量算法是对直接测量数据进行全体处理获得的端面函数模型,仅能对所有输入点信息进行统一分析处理,无法避免因测量系统误差及端面自身缺陷带来的测量噪点对端面函数拟合结果的影响,同时各温度工况下所有测量数据直接取均值的方式无法消除测量奇异数据对温度工况最终测量结果的影响,这些影响将导致测量方法精度表现不佳、测量误差较大。
发明内容
本发明要解决的技术问题是:克服现有技术的不足,本发明提供了一种高低温环境下的端面位姿相对变化高精度测量方法和系统,能够有效的剔除测量中的端面噪点及奇异数据,进一步减小测量误差提高测量精度。
本发明所采用的技术方案是:一种高低温环境下的端面位姿相对变化高精度测量方法,包括步骤如下:
步骤一,根据被测端面的结构,确定被测端面函数形式f(r0,r1,r2,...,rN),和函数相关未知参数r0,r1,r2,...,rN,N为正整数;
步骤二,导出每次测量获得的被测端面的点云图片,点云图片上包含被测端面的点云坐标点集{Ij1(Xj1,Yj1,Zj1),Ij2(Xj2,Yj2,Zj2),...,Iji(Xji,Yji,Zji),...,Inp(Xnp,Ynp,Znp)};Iji(Xji,Yji,Zji)表示第j次测量被测端面上的第i个点的坐标;i=1,2,3,...,p,p为被测端面上的点的总数;j=1,2,3,...,n,n为对被测端面进行点云测量的总数,p,n均为正整数;
步骤三,利用端面拟合算法分别对每次测量采集的点云进行端面相关函数拟合,确定每次采集的端面函数的未知参数;
步骤四,将步骤三中用于计算的每次测量获得的所有点云坐标Iji(Xji,Yji,Zji)与相应的拟合出的端面函数进行距离计算,公式如下:
ΔLi=f(Xji,Yji,r0,r1,r2...rN)-Zji
ΔLi表示第i个点距拟合出的平面的距离;
步骤五,进行循环判断:如果每次测量得到的所有点坐标到拟合函数的距离都小于α则拟合的函数为测量得到的端面函数,进入步骤七;如果存在点坐标到拟合函数的距离大于α的点,则剔除相应的点形成新的点集,进入步骤六;
α为设定的理论测量精度;
步骤六,将每次测量筛选后的剩余点重新作为新点云,返回步骤三,对每次测量采集的剩余点云进行重新函数拟合,获取新的拟合函数,再次计算每次采集的剩余点云相对新拟合函数的距离,重新进行点的筛选;
步骤七,得到每次测量对应的被测端面拟合函数后,以第一次测量对应的拟合出的函数的特征向量f0为基准,对比分析出其他次测量对应的拟合出的函数的特征向量fj相对特征向量f0的角度变化值A{A1,A2,A3,...,Aj,...,An};
步骤八,按照各温度工况条件给所有测量得到的点云图片进行分组,进行各温度工况下平均位姿角度变化值计算:
将各温度工况下的k个端面角度变化值分别进行取平均处理,得到所有m个温度工况下的平均角度变化值a{a1,a2,a3,...,am}:
第m个温度工况下的平均角度变化值am=(A1+A2+....+Ak)/k;
a={a1,a2,a3,...,am};
k、m均为正整数,n=k·m;
步骤九,将每个温度工况下的所有角度变化值与该温度工况对应的平均角度变化值进行对比,获取每个温度工况下k张点云图片的角度变化值相对平均角度变化值的差值绝对值ΔA{ΔA1,ΔA2,ΔA3,...,ΔAk};
第m个温度工况对应的平均角度变化差值点集为ΔA{ΔA1,ΔA2,ΔA3,...,ΔAk},其中,ΔA1=A1-am、ΔA2=A2-am、......、ΔAk=Ak-am;
步骤十,计算各温度工况下所有角度变化值的标准差A#,将各温度工况下所有角度变化值相对平均角度变化值的差值绝对值ΔA与相应的标准差A#进行对比,筛除各温度工况下角度变化值相对平均角度变化值的差值绝对值ΔA大于标准差A#的所有奇异数据;
步骤十一,进行循环判断:如果各温度工况剩余数据数量均小于初始数量k的60%,则进入步骤十二;如果存在剩余数据数量大于初始数量k的60%的温度工况,则进入步骤十三;
步骤十二,对各温度工况剩余数据进行重新计算平均角度变化值,将重新计算的均值作为各温度工况下的端面相对位姿变化最终测量值;
步骤十三,将剩余数据数量大于初始数量k的60%的所有温度工况的剩余数据按温度工况区分,重新进行平均处理获得角度变化值的新平均值a’,返回步骤九,对比剩余点云图片的角度变化值与相应各角度变化值的新平均值a’的差值绝对值,继续剔除新角度变化值相对平均角度变化值的差值绝对值大于标准差A#的所有数据,再进入步骤十一进行判断,直至所有温度工况下剩余数据均满足要求。
一种高低温环境下的端面位姿相对变化高精度测量系统,包括:
第一模块,用于计算被测端面每次测量的位姿初始函数;
第二模块,用于进行被测端面采集点云噪点剔除,筛选确认被测端面每次测量的最终位姿状态函数;
第三模块,用于计算被测端面各温度工况下位姿变化结果的初始结果;
第四模块,用于进行端面位姿变化测量数据筛选,确认各温度工况下的端面相对位姿变化最终测量结果。
计算被测端面每次测量的位姿初始函数的具体步骤如下:
步骤1.1,根据被测端面的结构,确定被测端面函数形式f(r0,r1,r2,...,rN),和函数相关未知参数r0,r1,r2,...,rN,N为正整数;
步骤1.2,导出每次测量获得的被测端面的点云图片,点云图片上包含被测端面的点云坐标点集{Ij1(Xj1,Yj1,Zj1),Ij2(Xj2,Yj2,Zj2),...,Iji(Xji,Yji,Zji),...,Inp(Xnp,Ynp,Znp)};Iji(Xji,Yji,Zji)表示第j次测量被测端面上的第i个点的坐标;i=1,2,3,...,p,p为被测端面上的点的总数;j=1,2,3,...,n,n为对被测端面进行点云测量的总数,p,n均为正整数;
步骤1.3,利用端面拟合算法分别对每次测量采集的点云进行端面相关函数拟合,确定每次采集的端面函数的未知参数,得到被测端面每次测量的位姿初始函数。
进行被测端面采集点云噪点剔除,筛选确认被测端面每次测量的最终位姿状态函数的具体步骤如下:
步骤2.1,将步骤1.3中用于计算的每次测量获得的所有点云坐标Iji(Xji,Yji,Zji)与相应的拟合出的端面函数进行距离计算,公式如下:
ΔLi=f(Xji,Yji,r0,r1,r2...rN)-Zji
ΔLi表示第i个点距拟合出的平面的距离;
步骤2.2,进行循环判断:如果每次测量得到的所有点坐标到拟合函数的距离都小于α则拟合的函数为测量得到的端面的最终位姿状态函数;如果存在点坐标到拟合函数的距离大于α的点,则剔除相应的点形成新的点集,进入步骤2.3;α为设定的理论测量精度;
步骤2.3,将每次测量筛选后的剩余点重新作为新点云,返回步骤1.3,对每次测量采集的剩余点云进行重新函数拟合,获取新的拟合函数,再次计算每次采集的剩余点云相对新拟合函数的距离,重新进行点的筛选。
计算被测端面各温度工况下位姿变化结果的初始结果的具体步骤如下:
步骤3.1,根据步骤2得到的每次测量对应的被测端面的最终位姿状态函数,以第一次测量对应的拟合出的函数的特征向量f0为基准,对比分析出其他次测量对应的拟合出的函数的特征向量fj相对特征向量f0的角度变化值A{A1,A2,A3,...,Aj,...,An};
步骤3.2,按照各温度工况条件给所有测量得到的点云图片进行分组,进行各温度工况下平均位姿角度变化值计算:
将各温度工况下的k个端面角度变化值分别进行取平均处理,得到所有m个温度工况下的平均角度变化值a{a1,a2,a3,...,am}:
第m个温度工况下的平均角度变化值am=(A1+A2+....+Ak)/k;
a={a1,a2,a3,...,am};
k、m均为正整数,n=k·m。
进行端面位姿变化测量数据筛选,确认各温度工况下的端面相对位姿变化最终测量结果的具体步骤如下:
步骤4.1,将每个温度工况下的所有角度变化值与该温度工况对应的平均角度变化值进行对比,获取每个温度工况下k张点云图片的角度变化值相对平均角度变化值的差值绝对值ΔA{ΔA1,ΔA2,ΔA3,...,ΔAk};
第m个温度工况对应的平均角度变化差值点集为ΔA{ΔA1,ΔA2,ΔA3,...,ΔAk},其中,ΔA1=A1-am、ΔA2=A2-am、......、ΔAk=Ak-am;
步骤4.2,计算各温度工况下所有角度变化值的标准差A#,将各温度工况下所有角度变化值相对平均角度变化值的差值绝对值ΔA与相应的标准差A#进行对比,筛除各温度工况下角度变化值相对平均角度变化值的差值绝对值ΔA大于标准差A#的所有奇异数据;
步骤4.3,进行循环判断:如果各温度工况剩余数据数量均小于初始数量k的60%,则进入步骤4.4;如果存在剩余数据数量大于初始数量k的60%的温度工况,则进入步骤4.5;
步骤4.4,对各温度工况剩余数据进行重新计算平均角度变化值,将重新计算的均值作为各温度工况下的端面相对位姿变化最终测量值;
步骤4.5,将剩余数据数量大于初始数量k的60%的所有温度工况的剩余数据按温度工况区分,重新进行平均处理获得角度变化值的新平均值a’,返回步骤4.1,对比剩余点云图片的角度变化值与相应各角度变化值的新平均值a’的差值绝对值,继续剔除新角度变化值相对平均角度变化值的差值绝对值大于标准差A#的所有数据,再进入步骤4.3进行判断,直至所有温度工况下剩余数据均满足要求。
本发明与现有技术相比的优点在于:
(1)本发明适用于各种端面变形测量系统对各种被测端面在高低温环境甚至其他复杂工况条件下进行变形测量时测量数据处理及位姿变形计算。本发明旨在针对现有高低温等其他工况条件下端面位姿测量过程中测量数据处理方法无法消除测量噪点及较差数据影响造成测量误差大精度不高等不足,发明一种高精度端面位姿测量数据处理方法,主要解决现有端面位姿测量数据处理中噪点消除及数据筛选等问题。
(2)本发明的方法可解决各种微变形测量系统对被测端面在高低温环境下的端面相对变形进行测量时因测量系统自身因素及外界影响产生的测量噪点及测量奇异数据对测量结果的影响,有效的减小了端面变形测量中的测量误差,提高了测量精度。该测量算法还可适用于高低温环境外的其他工况加载端面变形情况测量,具有较广的适用性。
(3)本发明提出了一种端面测量采集点云噪点筛除方法,实现了对端面进行测量过程中对采集到的点云进行噪点筛选,并进一步利用剩余较好的测量点云拟合出最佳被测端面的函数模型,有效地实现了端面位姿高精度拟合测量。
(4)本发明提出了一种测量数据高精度筛选方法,针对各工况下的所有初始测量数据进行标准差筛选提出测量奇异数据,并且为保证计算采样次数足够进行相关筛选中止设定确保最终处理数据量足够。该方法可解决端面位姿测量过程中每次测量间数据波动大、测量精度不稳定的问题,有效的避免了测量过程因随机误差带较大而来的测量奇异数据对测试结果的影响,有效的保证了测量过程中的全程测量高精度。
附图说明
图1为高低温环境下的端面位姿相对变化高精度测量方法流程图。
图2为端面测量噪点筛选示意图。
图3为测量数据筛选示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明作进一步详细描述。
一种高低温环境下的端面位姿相对变化高精度测量方法,主要按以下四阶段完成:
1.被测端面每次采集位姿初始函数计算。首先针对采集到的各温度工况下的端面点云坐标信息进行端面形态拟合处理获取端面的初始函数,主要过程是:根据端面的形态设定函数形式建立函数未知参数,根据每次采集的所有端面点云坐标解算出每次采集的函数未知参数,建立端面位姿的初始函数。
2.端面采集点云噪点剔除筛选确认被测端面每次测量的最终位姿状态函数。主要过程是:对比每次采集得到的所有测量点相对拟合函数端面的位置情况剔除位置偏移大于测量系统精度范围的噪点,对剩余点重新拟合,反复迭代筛选最后得到合理准确的端面位姿拟合函数。
3.端面各工况位姿变化结果初始结果计算。主要过程是:针对各工况下多次采集拟合出的端面函数,将第一次采集端面作为基准计算出每次采集时相对第一次采集的角度变化情况,按工况拆分后,进行均值处理,将均值作为端面各工况位姿变化初始测量结果。
4.端面位姿变化测量数据筛选确认该工况的端面角度最终变化测量结果。主要过程是:对比剔除各工况下每次采集的端面角度变化测量值与相应工况位姿变化均值计算结果偏差较大的奇异数据,对剩余数据重新取均值,反复迭代筛选,直到所有数据均在筛选范围内,则剩余数据再取均值为该工况的端面角度最终变化测量结果。
如图1所示,一种高低温环境下的端面位姿相对变化高精度测量方法,包括步骤如下:
步骤一,根据被测端面的结构,确定被测端面函数形式f(r0,r1,r2,...,rN),和函数相关未知参数r0,r1,r2,...,rN,如端面为平面则设定为平面函数,曲面为相应曲面函数,N为正整数;
步骤二,导出每次测量获得的被测端面点云图片,点云图片上包含被测端面的点云坐标点集{Ij1(Xj1,Yj1,Zj1),Ij2(Xj2,Yj2,Zj2),...,Iji(Xji,Yji,Zji)},...,
Inp(Xnp,Ynp,Znp);Iji(Xji,Yji,Zji)表示第j次测量被测端面上的第i个点的坐标;i=1,2,3,...,p,p为被测端面上的点的总数;j=1,2,3,...,n,n为对被测端面进行点云测量的总数,p,n均为正整数;
步骤三,利用端面拟合算法如最小二乘法等分别对每次测量采集的点云进行端面相关函数拟合,确认每次采集的端面函数的未知参数;
例如利用最小二乘法,即需拟合的点集中的所有点Iji(Xji,Yji,Zji)距拟合的端面函数距离之和应为最小:
∑Li2=∑{f(Xji,Yji,r0,r1,r2...rN)-Zji}2=min,
其中Li为点距拟合的端面函数距离,根据对应的偏导关系:
计算出以当前点云拟合出的端面函数模型参数:r0,r1,r2...rN。
步骤四,将上一步骤中用于计算的每次测量获得的所有点云坐标Iji(Xji,Yji,Zji)与相应的拟合出的端面函数进行距离计算。
ΔLi=f(Xji,Yji,r0,r1,r2...rN)-Zji
ΔLi表示第i个点距拟合出的平面的距离;
步骤五,进行循环判断,设设定的理论测量精度为α,如果所有点坐标到拟合函数的距离都小于α则拟合的函数即为测量得到的端面函数进入步骤七,如果存在距离大于α的点如图2所示,则剔除这些点形成新的点集,进入步骤六。
步骤六,将每次测量筛选后的剩余点重新作为新点云,重复步骤三对每次测量采集的剩余点云进行重新函数拟合,获取新的拟合函数,重复步骤四再次计算每次采集的剩余点云相对新拟合函数的距离,重新筛选。
步骤七,所有测量次数均筛选并完成函数拟合后,以第一次采集拟合出的函数的特征向量(例如:平面的法向量,曲面的轴向量等)f0为基准,对比分析出其他次采集拟合函数的特征向量fj相对f0的角度变化值A{A1,A2,A3,...,Aj,...,An}(假设共采集n次)。
步骤八,按照各温度工况条件给所有测量得到的点云图片进行分组,进行各温度工况平均位姿角度变化值进行计算:
设各温度工况下分别有k个测量数据,一共m个温度工况,将各温度工况下的k个端面角度变化值A分别进行取平均处理得到所有m个温度工况下各温度工况的平均角度变化值a{a1,a2,a3,...,am}。
a1=(A1+A2+....+Ak)/k,
a2=(A1+A2+....+Ak)/k,
......
am=(A1+A2+....+Ak)/k;
a={a1,a2,a3,...,am};
k、m均为正整数,n=k·m;
步骤九,将每个温度工况下的所有角度变化值与该温度工况的平均角度变化值进行对比,获取每个温度工况下k张点云图片的角度变化值相对平均角度变化值的差值绝对值ΔA{ΔA1,ΔA2,ΔA3,...,ΔAk};
以第一个温度工况为例,a1为第一温度工况的平均角度变化值,则第一温度工况的平均角度变化差值点集ΔA1{ΔA1,ΔA2,ΔA3,...,ΔAk}:ΔA1=A1-a1、ΔA2=A2-a1、......、ΔAk=Ak-a1;
a2为第二温度工况的平均角度变化值,则第二温度工况的平均角度变化差值点集ΔA2{ΔA1,ΔA2,ΔA3,...,ΔAk}:ΔA1=A1-a2、ΔA2=A2-a2、......、ΔAk=Ak-a2;
......
am为第m个温度工况的平均角度变化值,第m个温度工况对应的平均角度变化差值点集为ΔAm{ΔA1,ΔA2,ΔA3,...,ΔAk},其中,ΔA1=A1-am、ΔA2=A2-am、......、ΔAk=Ak-am;
步骤十,计算各温度工况下所有角度变化值的标准差A#,将各温度工况下所有角度变化值相对平均角度变化值的差值绝对值ΔA与相应的标准差A#进行对比,筛除各温度工况下角度变化值相对平均角度变化值的差值绝对值ΔA大于标准差A#的所有奇异数据,如图3。
其中标准差为:
步骤十一,进行循环判断,各温度工况测量数据筛选完成后,如果各温度工况剩余数据数量均小于初始数量k的60%,则进入步骤十二,如果存在剩余数据数量大于初始数量k的60%的温度工况,则进入步骤十三。
步骤十二,对各温度工况剩余数据进行重新计算平均角度变化值,将重新计算的均值作为各温度工况下的端面相对位姿变化最终测量值。
步骤十三,将剩余数据数量大于初始数量k的60%的所有温度工况的剩余数据按温度工况区分,重新进行平均处理获得角度变化值的新平均值a’,返回步骤九,对比剩余点云图片的角度变化值与相应各角度变化值的新平均值a’的差值绝对值,继续剔除新角度变化值相对平均角度变化值的差值绝对值大于标准差A#的所有数据,再进入步骤十一进行判断,直至所有温度工况下剩余数据均满足要求。
一种高低温环境下的端面位姿相对变化高精度测量系统,包括:
第一模块,用于计算被测端面每次测量的位姿初始函数;
第二模块,用于进行被测端面采集点云噪点剔除,筛选确认被测端面每次测量的最终位姿状态函数;
第三模块,用于计算被测端面各温度工况下位姿变化结果的初始结果;
第四模块,用于进行端面位姿变化测量数据筛选,确认各温度工况下的端面相对位姿变化最终测量结果。
计算被测端面每次测量的位姿初始函数的具体步骤如下:
步骤1.1,根据被测端面的结构,确定被测端面函数形式f(r0,r1,r2,...,rN),和函数相关未知参数r0,r1,r2,...,rN,N为正整数;
步骤1.2,导出每次测量获得的被测端面的点云图片,点云图片上包含被测端面的点云坐标点集{Ij1(Xj1,Yj1,Zj1),Ij2(Xj2,Yj2,Zj2),...,Iji(Xji,Yji,Zji),...,Inp(Xnp,Ynp,Znp)};Iji(Xji,Yji,Zji)表示第j次测量被测端面上的第i个点的坐标;i=1,2,3,...,p,p为被测端面上的点的总数;j=1,2,3,...,n,n为对被测端面进行点云测量的总数,p,n均为正整数;
步骤1.3,利用端面拟合算法分别对每次测量采集的点云进行端面相关函数拟合,确定每次采集的端面函数的未知参数,得到被测端面每次测量的位姿初始函数。
进行被测端面采集点云噪点剔除,筛选确认被测端面每次测量的最终位姿状态函数的具体步骤如下:
步骤2.1,将步骤1.3中用于计算的每次测量获得的所有点云坐标Iji(Xji,Yji,Zji)与相应的拟合出的端面函数进行距离计算,公式如下:
ΔLi=f(Xji,Yji,r0,r1,r2...rN)-Zji
ΔLi表示第i个点距拟合出的平面的距离;
步骤2.2,进行循环判断:如果每次测量得到的所有点坐标到拟合函数的距离都小于α则拟合的函数为测量得到的端面的最终位姿状态函数;如果存在点坐标到拟合函数的距离大于α的点,则剔除相应的点形成新的点集,进入步骤2.3;α为设定的理论测量精度;
步骤2.3,将每次测量筛选后的剩余点重新作为新点云,返回步骤1.3,对每次测量采集的剩余点云进行重新函数拟合,获取新的拟合函数,再次计算每次采集的剩余点云相对新拟合函数的距离,重新进行点的筛选。
计算被测端面各温度工况下位姿变化结果的初始结果的具体步骤如下:
步骤3.1,根据步骤2得到的每次测量对应的被测端面的最终位姿状态函数,以第一次测量对应的拟合出的函数的特征向量f0为基准,对比分析出其他次测量对应的拟合出的函数的特征向量fj相对特征向量f0的角度变化值A{A1,A2,A3,...,Aj,...,An};
步骤3.2,按照各温度工况条件给所有测量得到的点云图片进行分组,进行各温度工况下平均位姿角度变化值计算:
将各温度工况下的k个端面角度变化值分别进行取平均处理,得到所有m个温度工况下的平均角度变化值a{a1,a2,a3,...,am}:
第1个温度工况下的平均角度变化值a1=(A1+A2+....+Ak)/k,
第2个温度工况下的平均角度变化值a2=(A1+A2+....+Ak)/k,
......
第m个温度工况下的平均角度变化值am=(A1+A2+....+Ak)/k;
a={a1,a2,a3,...,am};
k、m均为正整数,n=k·m。
进行端面位姿变化测量数据筛选,确认各温度工况下的端面相对位姿变化最终测量结果的具体步骤如下:
步骤4.1,将每个温度工况下的所有角度变化值与该温度工况对应的平均角度变化值进行对比,获取每个温度工况下k张点云图片的角度变化值相对平均角度变化值的差值绝对值ΔA{ΔA1,ΔA2,ΔA3,...,ΔAk};
a1为第一温度工况的平均角度变化值,则第一温度工况的平均角度变化差值点集ΔA{ΔA1,ΔA2,ΔA3,...,ΔAk}:ΔA1=A1-a1、ΔA2=A2-a1、......、ΔAk=Ak-a1;
a2为第二温度工况的平均角度变化值,则第二温度工况的平均角度变化差值点集ΔA{ΔA1,ΔA2,ΔA3,...,ΔAk}:ΔA1=A1-a2、ΔA2=A2-a2、......、ΔAk=Ak-a2;
......
am为第m个温度工况的平均角度变化值,第m个温度工况对应的平均角度变化差值点集为ΔA{ΔA1,ΔA2,ΔA3,...,ΔAk},其中,ΔA1=A1-am、ΔA2=A2-am、......、ΔAk=Ak-am;
步骤4.2,计算各温度工况下所有角度变化值的标准差A#,将各温度工况下所有角度变化值相对平均角度变化值的差值绝对值ΔA与相应的标准差A#进行对比,筛除各温度工况下角度变化值相对平均角度变化值的差值绝对值ΔA大于标准差A#的所有奇异数据;
步骤4.3,进行循环判断:如果各温度工况剩余数据数量均小于初始数量k的60%,则进入步骤4.4;如果存在剩余数据数量大于初始数量k的60%的温度工况,则进入步骤4.5;
步骤4.4,对各温度工况剩余数据进行重新计算平均角度变化值,将重新计算的均值作为各温度工况下的端面相对位姿变化最终测量值;
步骤4.5,将剩余数据数量大于初始数量k的60%的所有温度工况的剩余数据按温度工况区分,重新进行平均处理获得角度变化值的新平均值a’,返回步骤4.1,对比剩余点云图片的角度变化值与相应各角度变化值的新平均值a’的差值绝对值,继续剔除新角度变化值相对平均角度变化值的差值绝对值大于标准差A#的所有数据,再进入步骤4.3进行判断,直至所有温度工况下剩余数据均满足要求。
本发明未详细说明部分属于本领域技术人员公知技术。
Claims (10)
1.一种高低温环境下的端面位姿相对变化高精度测量方法,其特征在于,包括步骤如下:
步骤1,计算被测端面每次测量的位姿初始函数;
步骤2,进行被测端面采集点云噪点剔除,筛选确认被测端面每次测量的最终位姿状态函数;
步骤3,计算被测端面各温度工况下位姿变化结果的初始结果;
步骤4,进行端面位姿变化测量数据筛选,确认各温度工况下的端面相对位姿变化最终测量结果。
2.根据权利要求1所述的一种高低温环境下的端面位姿相对变化高精度测量方法,其特征在于,步骤1的具体步骤如下:
步骤1.1,根据被测端面的结构,确定被测端面函数形式f(r0,r1,r2,...,rN),和函数相关未知参数r0,r1,r2,...,rN,N为正整数;
步骤1.2,导出每次测量获得的被测端面的点云图片,点云图片上包含被测端面的点云坐标点集{Ij1(Xj1,Yj1,Zj1),Ij2(Xj2,Yj2,Zj2),…,Iji(Xji,Yji,Zji),…,Inp(Xnp,Ynp,Znp)};Iji(Xji,Yji,Zji)表示第j次测量被测端面上的第i个点的坐标;i=1,2,3,...,p,p为被测端面上的点的总数;j=1,2,3,...,n,n为对被测端面进行点云测量的总数,p,n均为正整数;
步骤1.3,利用端面拟合算法分别对每次测量采集的点云进行端面相关函数拟合,确定每次采集的端面函数的未知参数,得到被测端面每次测量的位姿初始函数。
3.根据权利要求2所述的一种高低温环境下的端面位姿相对变化高精度测量方法,其特征在于,步骤2的具体步骤如下:
步骤2.1,将步骤1.3中用于计算的每次测量获得的所有点云坐标Iji(Xji,Yji,Zji)与相应的拟合出的端面函数进行距离计算,公式如下:
ΔLi=f(Xji,Yji,r0,r1,r2...rN)-Zji
ΔLi表示第i个点距拟合出的平面的距离;
步骤2.2,进行循环判断:如果每次测量得到的所有点坐标到拟合函数的距离都小于α则拟合的函数为测量得到的端面的最终位姿状态函数;如果存在点坐标到拟合函数的距离大于α的点,则剔除相应的点形成新的点集,进入步骤2.3;α为设定的理论测量精度;
步骤2.3,将每次测量筛选后的剩余点重新作为新点云,返回步骤1.3,对每次测量采集的剩余点云进行重新函数拟合,获取新的拟合函数,再次计算每次采集的剩余点云相对新拟合函数的距离,重新进行点的筛选。
4.根据权利要求3所述的一种高低温环境下的端面位姿相对变化高精度测量方法,其特征在于,步骤3的具体步骤如下:
步骤3.1,根据步骤2得到的每次测量对应的被测端面的最终位姿状态函数,以第一次测量对应的拟合出的函数的特征向量f0为基准,对比分析出其他次测量对应的拟合出的函数的特征向量fj相对特征向量f0的角度变化值A{A1,A2,A3,...,Aj,...,An};
步骤3.2,按照各温度工况条件给所有测量得到的点云图片进行分组,进行各温度工况下平均位姿角度变化值计算:
将各温度工况下的k个端面角度变化值分别进行取平均处理,得到所有m个温度工况下的平均角度变化值a{a1,a2,a3,...,am}:
第m个温度工况下的平均角度变化值am=(A1+A2+....+Ak)/k;
a={a1,a2,a3,...,am};
k、m均为正整数,n=k·m。
5.根据权利要求4所述的一种高低温环境下的端面位姿相对变化高精度测量方法,其特征在于,步骤4的具体步骤如下:
步骤4.1,将每个温度工况下的所有角度变化值与该温度工况对应的平均角度变化值进行对比,获取每个温度工况下k张点云图片的角度变化值相对平均角度变化值的差值绝对值ΔA{ΔA1,ΔA2,ΔA3,...,ΔAk};
第m个温度工况对应的平均角度变化差值点集为ΔA{ΔA1,ΔA2,ΔA3,...,ΔAk},其中,ΔA1=A1-am、ΔA2=A2-am、……、ΔAk=Ak-am;
步骤4.2,计算各温度工况下所有角度变化值的标准差A#,将各温度工况下所有角度变化值相对平均角度变化值的差值绝对值ΔA与相应的标准差A#进行对比,筛除各温度工况下角度变化值相对平均角度变化值的差值绝对值ΔA大于标准差A#的所有奇异数据;
步骤4.3,进行循环判断:如果各温度工况剩余数据数量均小于初始数量k的60%,则进入步骤4.4;如果存在剩余数据数量大于初始数量k的60%的温度工况,则进入步骤4.5;
步骤4.4,对各温度工况剩余数据进行重新计算平均角度变化值,将重新计算的均值作为各温度工况下的端面相对位姿变化最终测量值;
步骤4.5,将剩余数据数量大于初始数量k的60%的所有温度工况的剩余数据按温度工况区分,重新进行平均处理获得角度变化值的新平均值a’,返回步骤4.1,对比剩余点云图片的角度变化值与相应各角度变化值的新平均值a’的差值绝对值,继续剔除新角度变化值相对平均角度变化值的差值绝对值大于标准差A#的所有数据,再进入步骤4.3进行判断,直至所有温度工况下剩余数据均满足要求。
6.一种高低温环境下的端面位姿相对变化高精度测量系统,其特征在于,包括:
第一模块,用于计算被测端面每次测量的位姿初始函数;
第二模块,用于进行被测端面采集点云噪点剔除,筛选确认被测端面每次测量的最终位姿状态函数;
第三模块,用于计算被测端面各温度工况下位姿变化结果的初始结果;
第四模块,用于进行端面位姿变化测量数据筛选,确认各温度工况下的端面相对位姿变化最终测量结果。
7.根据权利要求6所述的一种高低温环境下的端面位姿相对变化高精度测量系统,其特征在于,计算被测端面每次测量的位姿初始函数的具体步骤如下:
步骤1.1,根据被测端面的结构,确定被测端面函数形式f(r0,r1,r2,...,rN),和函数相关未知参数r0,r1,r2,...,rN,N为正整数;
步骤1.2,导出每次测量获得的被测端面的点云图片,点云图片上包含被测端面的点云坐标点集{Ij1(Xj1,Yj1,Zj1),Ij2(Xj2,Yj2,Zj2),…,Iji(Xji,Yji,Zji),…,Inp(Xnp,Ynp,Znp)};Iji(Xji,Yji,Zji)表示第j次测量被测端面上的第i个点的坐标;i=1,2,3,...,p,p为被测端面上的点的总数;j=1,2,3,...,n,n为对被测端面进行点云测量的总数,p,n均为正整数;
步骤1.3,利用端面拟合算法分别对每次测量采集的点云进行端面相关函数拟合,确定每次采集的端面函数的未知参数,得到被测端面每次测量的位姿初始函数。
8.根据权利要求7所述的一种高低温环境下的端面位姿相对变化高精度测量系统,其特征在于,进行被测端面采集点云噪点剔除,筛选确认被测端面每次测量的最终位姿状态函数的具体步骤如下:
步骤2.1,将步骤1.3中用于计算的每次测量获得的所有点云坐标Iji(Xji,Yji,Zji)与相应的拟合出的端面函数进行距离计算,公式如下:
ΔLi=f(Xji,Yji,r0,r1,r2...rN)-Zji
ΔLi表示第i个点距拟合出的平面的距离;
步骤2.2,进行循环判断:如果每次测量得到的所有点坐标到拟合函数的距离都小于α则拟合的函数为测量得到的端面的最终位姿状态函数;如果存在点坐标到拟合函数的距离大于α的点,则剔除相应的点形成新的点集,进入步骤2.3;α为设定的理论测量精度;
步骤2.3,将每次测量筛选后的剩余点重新作为新点云,返回步骤1.3,对每次测量采集的剩余点云进行重新函数拟合,获取新的拟合函数,再次计算每次采集的剩余点云相对新拟合函数的距离,重新进行点的筛选。
9.根据权利要求8所述的一种高低温环境下的端面位姿相对变化高精度测量系统,其特征在于,计算被测端面各温度工况下位姿变化结果的初始结果的具体步骤如下:
步骤3.1,根据步骤2得到的每次测量对应的被测端面的最终位姿状态函数,以第一次测量对应的拟合出的函数的特征向量f0为基准,对比分析出其他次测量对应的拟合出的函数的特征向量fj相对特征向量f0的角度变化值A{A1,A2,A3,...,Aj,...,An};
步骤3.2,按照各温度工况条件给所有测量得到的点云图片进行分组,进行各温度工况下平均位姿角度变化值计算:
将各温度工况下的k个端面角度变化值分别进行取平均处理,得到所有m个温度工况下的平均角度变化值a{a1,a2,a3,...,am}:
第m个温度工况下的平均角度变化值am=(A1+A2+....+Ak)/k;
a={a1,a2,a3,...,am};
k、m均为正整数,n=k·m。
10.根据权利要求9所述的一种高低温环境下的端面位姿相对变化高精度测量系统,其特征在于,进行端面位姿变化测量数据筛选,确认各温度工况下的端面相对位姿变化最终测量结果的具体步骤如下:
步骤4.1,将每个温度工况下的所有角度变化值与该温度工况对应的平均角度变化值进行对比,获取每个温度工况下k张点云图片的角度变化值相对平均角度变化值的差值绝对值ΔA{ΔA1,ΔA2,ΔA3,...,ΔAk};
第m个温度工况对应的平均角度变化差值点集为ΔA{ΔA1,ΔA2,ΔA3,...,ΔAk},其中,ΔA1=A1-am、ΔA2=A2-am、……、ΔAk=Ak-am;
步骤4.2,计算各温度工况下所有角度变化值的标准差A#,将各温度工况下所有角度变化值相对平均角度变化值的差值绝对值ΔA与相应的标准差A#进行对比,筛除各温度工况下角度变化值相对平均角度变化值的差值绝对值ΔA大于标准差A#的所有奇异数据;
步骤4.3,进行循环判断:如果各温度工况剩余数据数量均小于初始数量k的60%,则进入步骤4.4;如果存在剩余数据数量大于初始数量k的60%的温度工况,则进入步骤4.5;
步骤4.4,对各温度工况剩余数据进行重新计算平均角度变化值,将重新计算的均值作为各温度工况下的端面相对位姿变化最终测量值;
步骤4.5,将剩余数据数量大于初始数量k的60%的所有温度工况的剩余数据按温度工况区分,重新进行平均处理获得角度变化值的新平均值a’,返回步骤4.1,对比剩余点云图片的角度变化值与相应各角度变化值的新平均值a’的差值绝对值,继续剔除新角度变化值相对平均角度变化值的差值绝对值大于标准差A#的所有数据,再进入步骤4.3进行判断,直至所有温度工况下剩余数据均满足要求。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110612495.3A CN113513978B (zh) | 2021-06-02 | 2021-06-02 | 高低温环境下的端面位姿相对变化高精度测量方法和系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110612495.3A CN113513978B (zh) | 2021-06-02 | 2021-06-02 | 高低温环境下的端面位姿相对变化高精度测量方法和系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113513978A true CN113513978A (zh) | 2021-10-19 |
CN113513978B CN113513978B (zh) | 2023-04-14 |
Family
ID=78065197
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110612495.3A Active CN113513978B (zh) | 2021-06-02 | 2021-06-02 | 高低温环境下的端面位姿相对变化高精度测量方法和系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113513978B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114001654A (zh) * | 2021-11-01 | 2022-02-01 | 北京卫星制造厂有限公司 | 工件端面位姿评价方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20160171735A1 (en) * | 2014-12-16 | 2016-06-16 | Fu Tai Hua Industry (Shenzhen) Co., Ltd. | Computing device and method for joining point clouds |
CN106600639A (zh) * | 2016-12-09 | 2017-04-26 | 江南大学 | 遗传算法结合自适应阈值约束的icp位姿定位技术 |
CN110119438A (zh) * | 2019-04-23 | 2019-08-13 | 东华理工大学 | 基于主动学习的机载LiDAR点云滤波方法 |
CN110967206A (zh) * | 2019-10-31 | 2020-04-07 | 北京卫星制造厂有限公司 | 高低温环境下细长件端面相对位姿关系测量装置和方法 |
CN111059998A (zh) * | 2019-12-31 | 2020-04-24 | 中国地质大学(北京) | 一种基于高分辨率的时序InSAR形变监测方法及系统 |
CN111055293A (zh) * | 2019-11-21 | 2020-04-24 | 东南大学 | 一种基于曲面自适应的工业机器人高精度恒力打磨方法 |
CN111127658A (zh) * | 2019-12-23 | 2020-05-08 | 北京工商大学 | 一种基于点云重建三角网格曲面的保特征曲面重建方法 |
CN111243008A (zh) * | 2020-01-19 | 2020-06-05 | 广西师范大学 | 一种用于高精度工件的圆弧数据拟合方法 |
CN111964575A (zh) * | 2020-07-06 | 2020-11-20 | 北京卫星制造厂有限公司 | 一种面向移动机器人铣削加工的数字孪生建模方法 |
CN112257722A (zh) * | 2020-11-11 | 2021-01-22 | 南京工业大学 | 基于抗差非线性高斯-赫尔默特模型的点云拟合方法 |
-
2021
- 2021-06-02 CN CN202110612495.3A patent/CN113513978B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20160171735A1 (en) * | 2014-12-16 | 2016-06-16 | Fu Tai Hua Industry (Shenzhen) Co., Ltd. | Computing device and method for joining point clouds |
CN106600639A (zh) * | 2016-12-09 | 2017-04-26 | 江南大学 | 遗传算法结合自适应阈值约束的icp位姿定位技术 |
CN110119438A (zh) * | 2019-04-23 | 2019-08-13 | 东华理工大学 | 基于主动学习的机载LiDAR点云滤波方法 |
CN110967206A (zh) * | 2019-10-31 | 2020-04-07 | 北京卫星制造厂有限公司 | 高低温环境下细长件端面相对位姿关系测量装置和方法 |
CN111055293A (zh) * | 2019-11-21 | 2020-04-24 | 东南大学 | 一种基于曲面自适应的工业机器人高精度恒力打磨方法 |
CN111127658A (zh) * | 2019-12-23 | 2020-05-08 | 北京工商大学 | 一种基于点云重建三角网格曲面的保特征曲面重建方法 |
CN111059998A (zh) * | 2019-12-31 | 2020-04-24 | 中国地质大学(北京) | 一种基于高分辨率的时序InSAR形变监测方法及系统 |
CN111243008A (zh) * | 2020-01-19 | 2020-06-05 | 广西师范大学 | 一种用于高精度工件的圆弧数据拟合方法 |
CN111964575A (zh) * | 2020-07-06 | 2020-11-20 | 北京卫星制造厂有限公司 | 一种面向移动机器人铣削加工的数字孪生建模方法 |
CN112257722A (zh) * | 2020-11-11 | 2021-01-22 | 南京工业大学 | 基于抗差非线性高斯-赫尔默特模型的点云拟合方法 |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114001654A (zh) * | 2021-11-01 | 2022-02-01 | 北京卫星制造厂有限公司 | 工件端面位姿评价方法 |
WO2023072086A1 (zh) * | 2021-11-01 | 2023-05-04 | 北京卫星制造厂有限公司 | 工件端面位姿评价方法 |
CN114001654B (zh) * | 2021-11-01 | 2024-03-26 | 北京卫星制造厂有限公司 | 工件端面位姿评价方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113513978B (zh) | 2023-04-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110543618A (zh) | 基于概率密度函数估计的圆度不确定度评定方法 | |
CN112082491A (zh) | 基于点云的高度检测方法 | |
CN101718528A (zh) | 基于数字图像的快速求解圆参数方法 | |
CN113513978B (zh) | 高低温环境下的端面位姿相对变化高精度测量方法和系统 | |
CN113269673B (zh) | 一种基于标准球架的三维点云拼接方法 | |
CN114001654B (zh) | 工件端面位姿评价方法 | |
CN111540001A (zh) | 航空发动机涡轮叶片气膜孔轴线方向检测方法 | |
CN111765902A (zh) | 一种基于多棱锥标靶的激光点云精度评估方法 | |
CN110017797B (zh) | 一种基于图像等值面分割法的尺寸测量结果不确定度评价方法 | |
CN113702994B (zh) | 一种基于刚性约束的激光跟踪仪测量精度提升方法 | |
CN110864671B (zh) | 基于线结构光拟合平面的机器人重复定位精度测量方法 | |
CN112344863B (zh) | 一种自由曲面回转体工件壁厚检测方法 | |
CN112097692B (zh) | 一种轮对垂直度动态测量的方法及系统 | |
CN111060050B (zh) | 一种综合验证三坐标测量机检测软件的方法 | |
CN110222382B (zh) | 一种飞机坐标系优化拟合方法 | |
CN108375337B (zh) | 一种机器人及其工艺设备相对位姿的测定方法及测定装置 | |
CN112747773A (zh) | 基于Allan方差和随机多项式提高陀螺仪精度的方法 | |
CN110021027A (zh) | 一种基于双目视觉的切边点计算方法 | |
Shammari et al. | A Study to Evaluate the Calibration of Optical Three-Dimensional Scanner | |
Zhao et al. | Position measurement of the hole group based on point cloud registration algorithm | |
CN114719800B (zh) | 一种小球冠中心点和曲率半径参数的精确评定方法 | |
CN114571506B (zh) | 工业机器人性能测量的姿态对齐方法 | |
Ferreira et al. | Performance of articulated arm CMM using virtual spheres gauge and geometry deviation analysis | |
Babadzhanov et al. | A Topographic Method of Assessing Coating Thickness Measures | |
CN113409367A (zh) | 条纹投影测量点云逐点加权配准方法、设备和介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |