一种动态干涉仪干涉图参数失配的相位提取的方法技术

本发明专利技术公开了一种动态干涉仪同步移相干涉图参数失配的相位提取方法，该方法为：首先采集参数失配的空间同步移相干涉图，分割为四幅子干涉图；然后使用条纹细化法计算各干涉图的对比度和载频误差，将对比度、载频误差和预设的移相间隔的初值代入线性方程组计算干涉图的相位、背景以及调制度；接着设定干涉图间的相位、背景、调制度正确，根据方程组计算的对比度、载频误差以及移相间隔的偏差值，更新初值。重复两个过程直至收敛判定函数优于阈值，停止迭代，得到最终的对比度、载频误差以及移相间隔，最后代入最小二乘方程组计算相位。本发明专利技术降低了动态干涉仪相位提取的约束性，干涉图参数失配相位提取的准确性和可靠性高，简单高效、适用范围广。适用范围广。适用范围广。

【技术实现步骤摘要】
一种动态干涉仪干涉图参数失配的相位提取的方法


[0001]本专利技术涉及光干涉计量
，特别是一种动态干涉仪同步移相干涉图参数失配的相位提取方法。

技术介绍

[0002]空间移相干涉可以动态实时捕获多幅移相干涉图，提取相位空间分布的参数，用于测量动态变化的对象，在光学制造、天文、国防等领域具有广泛的应用。解算相位时，若干涉图之间的对比度、载频以及移相间隔不一致，会引入相位误差。为了保证相位的复原精度，需要对失配干涉图进行参数计算，然后再进行相位提取，从而保障干涉仪的测量精度以及准确度。
[0003]其中如AIA[Z.Wang and B.Han,“Advanced iterative algorithm for randomly phase
‑
shifted interferograms with intra
‑
and inter
‑
frame intensity variations,”Opt.Lasers Eng.45(2),274
–
280(2007).]和QCA[J.Vargas and C.O.S.Sorzano,"Quadrature Component Analysis for interferometry,"Opt.Lasers Eng 51,637
–
641(2013)]法，可以快速准确地从包含随机平移相移且载频对比度一致的干涉图中提取相位。Deck[L.L.Deck,“Model
‑/>based phase shifting interferometry,”Appl.Opt.53(21),4628
–
4636(2014).]公开了一种基于多参数物理模型的MPSI方法，采用泰勒级数展开，通过线性回归得到相位，有效抑制了干涉图载频误差。该方法补偿了载频以及移相误差，将背景以及调制度作为时间无关量，但不能补偿对比度。Liu等人[Q.Liu,Y.Wang,J.He,and F.Ji,“Phase shift extraction and wavefront retrieval from interferograms with background and contrast fluctuations,”J.Opt.17(2),025704(2015).]在干涉图参数失配提取相位中做出了大量工作，在MPSI基础上通过条件近似，将非线性最小二乘化简为线性形式计算，得到对比度修正后的相位提取结果，但是不能补偿条纹载频误差。后续他又提出通过三次近似线性的方程组，交替迭代补偿载频误差，收敛速度快，但只适用于高载频干涉图。除此之外，空间载频法[E.McDonnell and L.Deck,“Solutions for environmentally robust interferometric optical testing”in Optical Manufacturing and Testing XIII,R.Rascher and D.W.Kim,eds.(SPIE,2020),p.16.]可以实现单帧干涉图计算，克服了干涉图间干涉参数不一致的影响，但是参考光束与测试光束间的大量倾斜会引入非共光路的回程误差，降低相位提取精度。上述的方法都仅仅补偿三者中部分参数失配，对于干涉图完全失配的相位提取不适用。

技术实现思路

[0004]本专利技术的目的在于提供一种准确性和可靠性高、简单高效、适用范围广的动态干涉仪同步移相干涉图参数失配的相位提取方法。
[0005]实现本专利技术目的的技术解决方案为：一种动态干涉仪同步移相干涉图参数失配的相位提取方法，包括以下步骤：
[0006]步骤1、使用动态干涉仪采集空间同步移相干涉图，并分割为四个子干涉图；
[0007]步骤2、使用等效对比度计算公式计算出每个子干涉图中的等效对比度信息，对各个子干涉图进行条纹细化，并进行插值拟合得到载频量，根据载频量计算各个子干涉图间的载频误差；
[0008]步骤3、将计算的对比度、载频误差以及预设的移相间隔代入干涉图帧之间相关的方程组，计算干涉图表达式中的背景、调制度以及相位分布；
[0009]步骤4、设定计算背景、调制度以及相位分布为实际值，代入干涉图帧内相关的方程组，计算对比度、载频误差以及移相间隔小偏差值，对对比度、载频误差以及移相间隔进行更新；
[0010]步骤5、根据计算的背景、调制度、相位分布、对比度、载频误差以及移相间隔，计算收敛判定函数值，进行判定，如果不满足收敛条件则重复步骤3～步骤5，直至收敛判定函数满足收敛条件；
[0011]步骤6、根据收敛条件下得到的对比度、载频误差以及移相间隔，代入干涉图模型进行最小二乘法计算，确定计算相位。
[0012]本专利技术与现有技术相比，其显著优点在于：(1)使用条纹细化法拟合得到载频误差，更方便简单；(2)通过对对比度、载频误差、移相间隔进行解耦得到干涉图所有的失配参数，准确性和可靠性高；(3)方法简单高效，不需要任何附加的辅助硬件，适用于绝大多数的动态干涉仪。
附图说明
[0013]图1是本专利技术动态干涉仪同步移相干涉图参数失配的相位提取方法的流程图。
[0014]图2是实施例中动态干涉仪采集的空间同步移相干涉图。
[0015]图3是实施例中四幅子干涉图参数失配干涉图的相位解算结果图。
[0016]图4是实施例中点光源异位动态斐索干涉仪采集得到空间同步移相干涉图。
[0017]图5是实施例中参数失配干涉图的相位解算结果图。
具体实施方式
[0018]本专利技术动态干涉仪同步移相干涉图参数失配的相位提取方法，包括以下步骤：
[0019]步骤1、使用动态干涉仪采集空间同步移相干涉图，并分割为四个子干涉图；
[0020]步骤2、使用等效对比度计算公式计算出每个子干涉图中的等效对比度信息，对各个子干涉图进行条纹细化，并进行插值拟合得到载频量，根据载频量计算各个干涉图间的载频误差；
[0021]步骤3、将计算的对比度、载频误差以及预设的移相间隔代入干涉图帧之间相关的方程组，计算干涉图表达式中的背景、调制度以及相位分布；
[0022]步骤4、设定计算背景、调制度以及相位分布为实际值，代入干涉图帧内相关的方程组，计算对比度、载频误差以及移相间隔小偏差值，对对比度、载频误差以及移相间隔进行更新；
[0023]步骤5、根据计算的背景、调制度、相位分布、对比度、载频误差以及移相间隔，计算收敛判定函数值，进行判定，如果不满足收敛条件则重复步骤3、4、5，直至收敛判定函数满
足收敛条件；
[0024]步骤6、根据收敛条件下得到的对比度、载频误差以及移相间隔，代入干涉图模型进行最小二乘法计算，确定计算相位。
[0025]进一步地，步骤1所述的采集空间同步移相干涉图，并分割为四个子干涉图，具体如下：
[0026]在动态干涉仪中采本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种动态干涉仪同步移相干涉图参数失配的相位提取方法，其特征在于，包括以下步骤：步骤1、使用动态干涉仪采集空间同步移相干涉图，并分割为四个子干涉图；步骤2、使用等效对比度计算公式计算出每个子干涉图中的等效对比度信息，对各个子干涉图进行条纹细化，并进行插值拟合得到载频量，根据载频量计算各个子干涉图间的载频误差；步骤3、将计算的对比度、载频误差以及预设的移相间隔代入干涉图帧之间相关的方程组，计算干涉图表达式中的背景、调制度以及相位分布；步骤4、设定计算背景、调制度以及相位分布为实际值，代入干涉图帧内相关的方程组，计算对比度、载频误差以及移相间隔小偏差值，对对比度、载频误差以及移相间隔进行更新；步骤5、根据计算的背景、调制度、相位分布、对比度、载频误差以及移相间隔，计算收敛判定函数值，进行判定，如果不满足收敛条件则重复步骤3～步骤5，直至收敛判定函数满足收敛条件；步骤6、根据收敛条件下得到的对比度、载频误差以及移相间隔，代入干涉图模型进行最小二乘法计算，确定计算相位。2.根据权利要求1中所述的动态干涉仪同步移相干涉图参数失配的相位提取方法，其特征在于，步骤1所述的使用动态干涉仪采集空间同步移相干涉图，并分割为四个子干涉图，具体如下：在动态干涉仪中采集得到空间同步移相干涉图，对其进行分割，得到四个子干涉图I0、I1、I2和I3。3.根据权利要求2中所述的动态干涉仪同步移相干涉图参数失配的相位提取方法，其特征在于，步骤2所述的使用等效对比度计算公式计算出每个子干涉图中的等效对比度信息，对各个子干涉图进行条纹细化，并进行插值拟合得到载频量，根据载频量计算各个子干涉图间的载频误差，具体如下：步骤2.1、分别对四个子干涉图I0、I1、I2和I3，根据等效对比度计算公式计算对应的等效对比度γ0、γ1、γ2、γ3；其中干涉图的等效对比度计算公式为：式中γ表示干涉图的对比度，I(r)表示干涉图在空间位置矢量r下的灰度，表示干涉图的灰度均值，M
×
N表示参与对比度计算的干涉图尺寸大小，即像素数，Ω为干涉图有效区域；步骤2.2、对各个子干涉图进行条纹细化处理，并进行插值拟合相位，根据倾斜估计载频量，根据载频量计算各个子干涉图间的载频误差。4.根据权利要求3中所述的动态干涉仪同步移相干涉图参数失配的相位提取方法，其特征在于，步骤3所述的将计算的对比度、载频误差以及预设的移相间隔代入干涉图帧之间相关的方程组，计算干涉图表达式中的背景、调制度以及相位分布，具体如下：步骤3.1、预设移相间隔为(0
°
，90
°
，180
°
，270
°
)；
步骤3.2、将计算的对比度、载频误差以及预设的移相间隔代入四幅子干涉图模型方程，得到空间同步移相干涉图为：式中A(r)为干涉图背景，γ
n
为干涉图对比度，下标n为干涉图序号，表示待测相位，r为空间位置矢量，x空间位置坐标中的横坐标，y为纵坐标，a
n
和b
n
为相移项中两个正交方向上的倾斜误差因子，δ
n
为干涉图间的移相量；步骤3.3、定义Δ
n
(r，n)＝a
n
·
x+b
n
·
y+δ
n
，则式(2)转换为：I(r，n)＝A(r)+C(r)γ
n
cos(Δ
n
(r，n))+S(r)γ
n
sin(Δ
n
(r，n))
ꢀꢀꢀꢀ
(3)由最小二乘准则建立方程组：[A(r) C(r) S(r)]
T
＝Q
‑1B
ꢀꢀꢀꢀ
(4)式中：式中：其中I
mes
...

【专利技术属性】
技术研发人员：郑东晖，胡晨辉，陈磊，王云涛，张喆，马致遥，
申请(专利权)人：南京理工大学，
类型：发明
国别省市：