一种基于勒让德傅立叶多项式的圆柱干涉拼接方法和系统技术方案

本发明专利技术提供一种基于勒让德傅立叶多项式的圆柱干涉拼接方法和系统，其中，所述方法包括：利用柱面干涉测量系统依次获得待测圆柱的子孔径测量数据，相邻测量区域间没有重叠但紧密相连；根据所述的子孔径测量数据，利用Legendre多项式分离失调像差，获得子孔径面形数据；建立局部坐标系和全局坐标系，将所述的子孔径面形数据转换成全局三维坐标数据；利用勒让德傅立叶多项式拟合所述的全局三维坐标数据，即可获得待测圆柱的360度面形分布。本发明专利技术提供的技术方案无需重叠区，因此能够减少拼接测量所需的子孔径数目，缩短测量时间；无需计算重叠区的对应点，能够降低拼接计算的复杂度；采用最小二乘拟合，能够准确地获得待测圆柱的360度面形误差分布。

Cylindrical interference splicing method and system based on Legendre Fu Liye polynomial

The present invention provides a method and system for stitching interference, Legendre Fu Liye polynomial based on the cylinder, the method comprises the following steps: using cylindrical interference measurement system respectively and obtain sub aperture measurement data to be measured cylindrical, no overlap between the adjacent measurement area but closely linked; according to the sub aperture measurement data, separation of offset aberration by Legendre polynomial, obtain sub aperture surface data; local coordinate system and global coordinate system, the sub aperture surface data into the global three-dimensional coordinate data; using the global 3D coordinates of the Legendre Fu Liye polynomial fitting, can obtain 360 degree cylindrical surface distribution to be measured. The technical scheme provided by the invention without overlap, so it can reduce the required number of sub aperture stitching, shorten the measurement time; no need to calculate the corresponding point overlap, can reduce the complexity of computing the mosaic; the least square fitting to 360 degrees, the surface error of accurately measured cylinder.

【技术实现步骤摘要】

本专利技术涉及光学三维测量领域，尤其涉及一种基于勒让德傅立叶多项式的圆柱干涉拼接方法。
技术介绍
在现代工业中，精密轴类零件是机器设备中的重要组成部分，如精密机床的回转轴及空气静压轴承，精密陀螺中的气浮轴承、动压马达等。精密轴类零件的面形误差直接影响到精密设备的使用性能。因此，需要对精密轴类零件检测，获取它们的面形误差，为零件的加工误差分析及工艺提供指导，保证轴类零件的质量。目前，针对精密轴类零件的面形误差检测，主要采用接触式测量方法，如三坐标测量机和圆柱度仪。然而，接触式测量方法的采样率低，难以获得高分辨率的测量结果，用于表征待测圆柱的360度面形误差分布。此外，接触式探针容易磨损，会产生测量误差，影响最终的测量结果。干涉测量术，作为光学三维测量技术的分支，由于具有非接触性、全场性、高精度、高分辨率等特点，使其在精密零件、光学元件的面形误差检测等领域得到了日益广泛的应用。由于大口径、低F/数的参考镜制造存在困难，使得干涉术还存在解决大空间尺寸物体、大数值孔径(或低F/数)、高分辨率及特殊形状物体(如球面、柱面等)的面形测量等方面的难题。干涉拼接术为解决这类问题提供了新思路。它的基本原理是将被测物体划分为若干个小尺寸的子孔径，相邻的子孔径间有局部重叠区，每次用小口径的干涉系统测量待测物体的局面面形，通过移动、转动待测物体或干涉系统，测得全部子孔径面形，然后采用拼接术得到全孔径的测量结果。文献(J.Peng,H.Xu,Y.Yu,andM.Chen,“Stitchinginterferometryforcylindricalopticswithlargeangularaperture,”Meas.Sci.Technol.,vol.26,no.2,p.25204,2015.)借鉴了干涉拼接术的基本思想，获得了大数值孔径柱面光学元件的全口径面形误差。然而，现有干涉拼接算法要求相邻的子孔径间有局部重叠区，因此子孔径数目较多，测量时间长。此外，为了实现拼接，需要查找重叠区内的对应点对，而调整机构的运动误差会引起点匹配误差，拼接过程需要采用迭代法消除点的匹配误差，导致拼接算法复杂。因此，如何提供一种高效、便捷的圆柱干涉拼接方法成为精密检测领域中的一项技术难题，亟待解决。
技术实现思路
有鉴于此，本专利技术的目的在于提供一种基于勒让德傅立叶(Legendre-Fourier)多项式的圆柱干涉拼接方法，旨在解决现有圆柱干涉拼接测量中需要重叠区、测量时间长；需要查找重叠区内的对应点，拼接算法复杂的问题。本专利技术提出一种基于勒让德傅立叶多项式的圆柱干涉拼接方法，主要包括：步骤一、利用柱面干涉测量系统获取待测圆柱的子孔径测量数据；步骤二、根据所述子孔径测量数据，分离测量数据中的失调像差，获得子孔径面形数据；步骤三、定义局部坐标系和全局坐标系，将所述的子孔径面形数据变换成全局三维坐标。步骤四、利用勒让德傅立叶多项式拟合所述的全局三维坐标，所得拟合结果即为待测圆柱的360度面形误差分布。另一方面，本专利技术还提供一种基于勒让德傅立叶多项式的圆柱干涉拼接系统，所述系统包括：子孔径数据获取模块，用于利用柱面干涉测量系统获取待测圆柱的子孔径测量数据；失调像差分离模块，用于通过二维Legendre多项式(V.N.Mahajan,“Orthonormalaberrationpolynomialsforanamorphicopticalimagingsystemswithrectangularpupils,”Appl.Opt.,vol.49,no.36,pp.6924–6929,2010.)分离所述子孔径数据中的失调像差，获得子孔径面形数据；坐标转换模块，用于定义局部坐标系和全局坐标系，将所述的子孔径面形数据转换成全局三维坐标；拼接模块，用于通过勒让德傅立叶多项式同时拟合所述的全局三维坐标数据，获得待测圆柱的360度面形误差分布。本专利技术提供的技术方案与现有干涉拼接技术相比，其显著优点包括：(1)本专利技术无需相邻的子孔径有局部重叠区，能够减少测量所需子孔径的数目，缩短测量时间；(2)本专利技术无需计算重叠区内的对应点，能够降低拼接算法的复杂性。附图说明图1为本专利技术一实施方式中基于勒让德傅立叶多项式的圆柱干涉拼接方法流程图；图2为本专利技术一实施方式中柱面干涉测量系统的示意图；图3为本专利技术一实施方式中基于勒让德傅立叶多项式的圆柱干涉拼接系统30的内部结构示意图；图4为本专利技术一实施方式中相邻子孔径间没有重叠区，利用柱面干涉系统对金属轴检测时获取的子孔径的测量结果；图5为本专利技术一实施方式中利用基于勒让德傅立叶多项式的拼接算法获得待测金属轴的360度面形误差分布图；图6为本专利技术一实施方式中相邻子孔径间具有局部重叠区，利用柱面干涉系统对金属轴检测时获取的子孔径的测量结果；图7为本专利技术一实施方式中利用重叠干涉拼接算法获得待测金属轴的360度面形误差分布图。具体实施方式为了使本专利技术的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本
技术实现思路
进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本专利技术，并不用于限定本专利技术。本专利技术采用如下技术方案：首先利用柱面干涉测量系统获取待测圆柱的子孔径测量数据；其次利用二维Legendre多项式分离测量结果中的失调像差，获得子孔径面形数据；接着定义局部坐标系和全局坐标系，将所述的子孔径面形数据转换成全局三维坐标数据；最后利用勒让德傅立叶多项式拟合全局三维坐标数据，即可获得待测圆柱的360度面形误差分布。这种方法无需要求相邻的子孔径间局部重叠，只需紧密连接，能够减少子孔径的数目，缩短测量时间；无需查找重叠点对，能够降低拼接算法的复杂性。以下将对本专利技术所提供的一种基于勒让德傅立叶多项式的圆柱干涉拼接方法进行详细说明。请参阅图1，为本专利技术一实施方式中基于勒让德傅立叶多项式的圆柱干涉拼接测量方法流程图。在步骤S1中，步骤一、利用柱面干涉测量系统获取待测圆柱的子孔径测量数据。在本实施方式中，所述的柱面干涉测量系统，请参阅图2，由平面干涉仪21，柱面波转换器22，多自由度的调整工作台23和待测圆柱24组成。它的工作原理是：柱面波转换器将干涉仪发出的平面波转换成理想柱面波前，经待测圆柱面反射后携带有面形误差信息再次经过柱面波转换器变成近似平面波，最后返回到干涉仪内部与参考波前干涉。通过解调干涉仪记录的条纹图即可获得待测柱面的面形误差。在本实施方式中，所述步骤一具体包括：根据柱面波转换器的F/数和通光孔径将待测圆柱划分成若干子孔径，并规划测量路径。根据测量路径调整待测圆柱，使每个子孔径在测量时的条纹数最少，即满足零位干涉测量条件，并利用柱面干涉测量系统获取每个子孔径测量数据。在步骤S2中，步骤二、分离子孔径测量数据中的失调像差，获得子孔径面形数据。在本实施方式中，所述步骤二具体包括：利用二维Legendre多项式拟合子孔径测量数据，其中W表示拟合的相位数据，ai表示多项式的系数，Li表示二维Legendre多项式的基函数，i＝1,…,N表示多项式的项数，ii为多项式的阶数，ii∈[7,15]。对于一般的光学系统，面形误差的高阶分量不会太大，使用7阶Legendre多项式拟合的误差小于0.01。二维Legendre多项式的阶数在理本文档来自技高网...

【技术保护点】
一种基于勒让德傅立叶多项式的圆柱拼接方法，其特征在于，所述方法包括：步骤一、利用柱面干涉测量系统获取待测圆柱的子孔径测量数据；步骤二、根据所述子孔径测量数据，利用Legendre多项式分离失调像差，获得子孔径面形数据；步骤三、定义局部坐标系和全局坐标系，将所述的子孔径面形数据变换到全局坐标系下，获得子孔径的全局三维坐标；步骤四、利用勒让德傅立叶多项式拟合所述的全局三维坐标，所得拟合结果即为待测圆柱的360度面形误差分布。

【技术特征摘要】
1.一种基于勒让德傅立叶多项式的圆柱拼接方法，其特征在于，所述方法包括：步骤一、利用柱面干涉测量系统获取待测圆柱的子孔径测量数据；步骤二、根据所述子孔径测量数据，利用Legendre多项式分离失调像差，获得子孔径面形数据；步骤三、定义局部坐标系和全局坐标系，将所述的子孔径面形数据变换到全局坐标系下，获得子孔径的全局三维坐标；步骤四、利用勒让德傅立叶多项式拟合所述的全局三维坐标，所得拟合结果即为待测圆柱的360度面形误差分布。2.如权利要求1所述的基于勒让德傅立叶多项式的圆柱干涉拼接方法，其特征在于，所述步骤一具体包括：⑴根据柱面波转换器的F/数和通光孔径将待测圆柱划分成若干子孔径，并规划测量路径；⑵根据测量路径调整待测圆柱，使每个子孔径在测量时的条纹数最少，即满足零位干涉测量条件，并利用柱面干涉测量系统获取每个子孔径测量数据。3.如权利要求2所述的基于勒让德傅立叶多项式的圆柱干涉拼接方法，其特征在于，所述步骤二具体包括：利用二维Legendre多项式拟合子孔径测量数据，其中W表示拟合的相位数据，ai表示多项式的系数，Li表示二维Legendre多项式的基函数，i＝1,…,N表示多项式的项数，ii为二维Legendre多项式的阶数；然后利用公式其中Wm表示子孔径测量数据，分离测量数据中的失调像差，获得子孔径面形数据Wr。4.如权利要求3所述的基于勒让德傅立叶多项式的圆柱干涉拼接方法，其特征在于，所述步骤三具体包括：⑴定义局部坐标系，坐标原点位于柱面波转换器的中心，利用公式将子孔径面形数据(u,v,Wr)转换成局部三维坐标(z,θ,Δr)，其中，s为比例缩放因子，d表示柱面波转换器上下两个标志点之间的距离，dp表示这两个标志点的像素坐标差；rbf表示柱面波转换器的后焦距，u和v表示像素坐标，z表示沿着圆柱轴线方向的坐标，θ表示角度坐标，Δr表示待测圆柱的径向误差；⑵定义全局坐标系，Z轴与待测圆柱的中心线重合，根据待测圆柱的名义运动参数，包括转动参数θ0和移动参数z0，用公式将局部三维坐标(z,θ,Δr)转换成全局三维坐标(Z,Θ,Δr)。5.如权利要求4所述的基于勒让德傅立叶多项式的圆柱干涉拼接方法，其特征在于，所述步骤四具体包括：⑴建立两个行数相等的矩阵Mk和Mc，分别用于表示各子孔径相对基准子孔径的位置偏差和待测圆柱的360度面形误差，其中Npx表示子孔径测量数据的像素数，Ns表示完成圆柱的360度面形测量所需要的子孔径数目，(Npx,Ns)表示第Ns个子孔径的第Npx个像素点，LFn,m表示勒让德傅立叶多项式的基函数：LFn,m=2(2n+1)Ln(Z)cos(mΘ),m=-jj,-jj+1,...,-1,02(2n+1)Ln(Z)sin(mΘ),m=1,...,jj]]>其中，n＝0,…,jj，jj表示勒让德傅立叶多项式的阶数；⑵建立如下方程ΔR＝MA其中ΔR＝[Δr(1,1),…,Δr(Npx,Ns)]T，M＝[Mc,Mk]，A为系数矩阵，将所有像素点对应的全局三维坐标(Z,Θ,Δr)代入到上述公式中构建线性方程组，然后求解上述方程组获得系数矩阵A的值，最后将矩阵Mc与其对应的系数相乘，获得待测圆柱的360度面形误差分布。6.一种基于勒让德傅立叶多项式的圆柱干涉拼接系统，其特征在于...

【专利技术属性】
技术研发人员：彭军政，钟金钢，
申请(专利权)人：暨南大学，
类型：发明
国别省市：广东;44