一种基于CFT的牛顿环参数估计方法技术

本发明专利技术公开了一种基于CFT的牛顿环参数估计方法，属于光学测量领域，该方法包括：将牛顿环的灰度值变换到[0,1]，计算相邻两行的差，从上述结果中随机选取

【技术实现步骤摘要】
一种基于CFT的牛顿环参数估计方法
本专利技术涉及一种基于CFT的牛顿环参数估计方法，属于光学测量领域。
技术介绍
干涉测量是通过干涉条纹获得物理参数的光学测量技术，在工程计量、光学计量、地震、光谱学等领域具有重要的应用价值。牛顿环作为经典的二次相位干涉条纹之一，在干涉测量中有着广泛的应用，如判断透镜表面的凹凸、测量透镜表面曲率半径和液体折射率等。随着光电技术的出现，现代信号处理理论可以应用到干涉测量中，形成了信号处理与传统干涉测量交叉的新的研究领域。这一新的交叉研究领域应用了一些现代信号处理工具来估计干涉条纹的参数，如傅里叶变换、加窗傅里叶变换、调频傅里叶变换、S变换、Kalmana滤波器、小波变换、极坐标变换和分数傅里叶变换。这些方法改善了传统干涉测量中的噪声和障碍物鲁棒性，但这些方法仍然存在一些缺点。例如，在精细确定条纹参数时，基于傅里叶变换的方法需要三次迭代来分离一个边带，迭代过程耗时较长；基于小波变换的方法可以提取任意局部范围内干涉条纹的频率特性，但仍存在阈值选择等问题，窗口大小的选择和耗时；基于小波变换的方法具有较好的时频局部分析能力，但小波函数的选取仍是一个难点；基于FRFT的方法采用二维干涉条纹和遍历搜索来获得匹配角，速度和精度不能同时得到保证。作为傅立叶变换的推广，Chirp-Fourier变换(CFT)将时间信号变换到频率和调频率平面上的二元函数，可以用来匹配具有多个分量的chirp型信号，这些良好的性质使得CFT成为信号处理和光学测量领域中的一个重要工具。专利技术内容本专利技术的目的是提供一种基于CFT的牛顿环参数估计方法；该方法不需要相位解缠绕，可以直接准确地估计牛顿环的相位导数，进而可以应用到球面曲率半径的测量；同时本专利技术比传统的牛顿环测量方法具有更好的噪声和遮挡鲁棒性，且具有计算复杂度低、精度高、速度快等优点。本专利技术的目的是通过下述技术方案实现的：一种基于CFT的牛顿环参数估计方法，包含如下步骤：步骤一、将牛顿环的灰度值变换到,记为,其中；步骤二、计算中相邻两行的差，记为，其中；步骤三、从中随机选取K行，分别利用DTCFT计算其幅度谱,即其中表示DTCFT算子；步骤四、令为数值优化算法ALO的代价函数,优化变量为和,设置ALO的初始参数并运行得到全局最优值,即步骤五、利用离群点检测算法去掉其中的离群点，并基于剩余结果计算其平均值得到；步骤六、牛顿环的相位导数可由下述方程给出，即其中表示牛顿环的相位；步骤七、利用相位导数中所蕴含的信息，计算球面的曲率半径，即其中表示入射光的波长，表示牛顿环的物理尺寸，表示球面曲率半径的估计值。进一步，牛顿环的每一行或列可以看作一个包含3个分量的chirp信号，其中常数分量对参数估计是无用的，但是它的存在会影响估计的稳定性。因此，本专利技术提出了一种简单实用的常数分量去除办法，通过计算牛顿环相邻两行或两列的差，有效降低了常数分量对估计结果的影响。另外，由于相邻两行或两列的差仍是chirp信号，且保持初始频率和调频率不变，同样可用于参数估计。进一步，本专利技术使用了异常点检测算法，仅需使用牛顿环的某些行来估计参数，大大降低了计算量。本专利技术的有益效果在于：1.本专利技术使用离散时间CFT来估计牛顿环的参数，突破了离散CFT对时频分辨率的限制，提高了参数估计的精度。2.由于异常点检测算法的应用，该方法仅使用牛顿环的某些行来估计参数，大大降低了计算量。3.建立了牛顿环参数估计的优化模型，并进一步应用ALO算法搜索最优参数，有效地改善了现有一些方法的不足。与传统的牛顿环干涉法相比，基于CFT的测量方法具有更好的噪声和遮挡鲁棒性。与基于FRFT的方法相比，该测量方法具有计算复杂度低、精度高、速度快等优点。附图说明图1—512x512的仿真牛顿环。图2—仿真牛顿环的相位导数及其误差。图3—720x720的实际牛顿环。图4—实际牛顿环的相位导数及其误差。图5—图1的和相应的CFT幅度谱。图6—高斯白噪声鲁棒性分析。图7—遮挡鲁棒性分析。图8—实际牛顿环。图9—基于720x720实际牛顿环的性能分析。图10—2幅中心不可见的实际牛顿环及其CFT幅度谱。具体实施方式下面结合具体的实例和附图对本专利技术做详细说明：为了验证所提方法在参数估计中的性能，基于仿真和实际牛顿环进行了实验分析。数值优化算法ALO的初始化参数为其中表示种群数目，表示最大迭代次数，LB和UB表示优化变量的下界和上界，表示优化变量的个数。对于如图1中所示的仿真牛顿环，其相位参数为，估计得到的相位参数为。为了更直观地显示估计结果，估计的相位导数及其误差如图2(a)和(b)所示。对于实际的牛顿环（见图3），实际牛顿环的相位参数是，对应的估计相位参数为。估计的相位导数及其误差如图4(a)和(b)所示。这些数据表明，仿真和实际牛顿环的相位导数可以在不使用解缠绕方法的情况下直接准确地估计出来。基于牛顿环参数所包含的信息，基于CFT的牛顿环参数估计方法可以进一步应用于测量球面的曲率半径。对于图1中的纯净仿真牛顿环，其和相应的CFT幅度谱如图5所示，从图中可以看出明显的冲激特性。该方法得到的曲率半径估计值分别为0.8606m和0.8595m，相对误差分别为0.0668%和0.0591%，处理时间约为3.8s，基于FRFT方法得到的曲率半径相对误差为0.22%，处理时间约为990s。这些数据说明了，基于CFT的测量方法的效率远高于基于FRFT的方法。由于在实际操作中，光学仪器和人工操作不可避免地会造成污染，为了更好地模拟实际情况，有必要在牛顿环受到污染的情况下对该方法的性能进行分析。主要从以下两点分析：1）污染源为高斯白噪声，信噪比在-10db到10db之间变化；2）污染源为较大的噪声点，由半径为100像素、中心位置不同的圆形遮挡物模拟。对于以上两点，所得结果列于表1。通过分析这些数据，可以得到以下结论：表1基于污染牛顿环的性能分析1）该方法对噪声具有较强的鲁棒性，当时，噪声对测量结果的影响较小。为了将上述结果可视化，将图1中的牛顿环添加了的高斯白噪声(见图6(a))，未添加噪声和添加噪声之后的如图6(b)所示。应用CFT到添加噪声之后的，得到的幅度谱及其投影如图6(c)和(d)所示。从图中可得，虽然牛顿环几乎被噪声湮没，但在图6(d)中可以看到明显的冲激。因此，基于CFT的测量方法对高斯白噪声具有较好的鲁棒性。2）该方法对遮挡具有很强的鲁棒性。与纯净牛顿环的结果相比，遮挡物对估计结果的影响很小。为了将上述结果可视化，牛顿环被圆形障碍物遮挡，其中心位于（256,256）像素(见图7(a))。图7(b)中展示了被遮挡牛顿环的和。由此图可以看出，由于障碍物的存在，不再是c本文档来自技高网...

【技术保护点】
1.一种基于CFT的牛顿环参数估计方法，其特征在于：/n步骤一、将

【技术特征摘要】
1.一种基于CFT的牛顿环参数估计方法，其特征在于：
步骤一、将牛顿环的灰度值变换到,记为,其中；
步骤二、计算中相邻两行的差，记为，其中；
步骤三、从中随机选取K行，分别利用DTCFT计算其幅度谱,即



其中表示DTCFT算子；
步骤四、令为数值优化算法ALO的代价函数,优化变量为和,设置ALO的初始参数并运行得到全局最优值,即



步骤五、利用离群点检测算法去掉其中的离群点，并基于剩余结果计算其平均值得到；
步...

【专利技术属性】
技术研发人员：郭勇，杨立东，唐俊，曹富军，
申请(专利权)人：内蒙古科技大学，
类型：发明
国别省市：内蒙;15