简单的个人主页网站制作,做服装网站的意义是什么,公司图案设计,wordpress如何上传Zixin Zhao1#xff0c;Hong Zhao1、Lu Zhang 1#xff0c;Fen Gao2#xff0c;Yuwei Qin3#xff0c;Hubing Du
摘要:
我们提出了一种适用于一般干涉测量应用的相位展开方法。所提出的方法依赖于导数泽尼克多项式拟合#xff08;DZPF#xff09;技术#xff0c;其中相…Zixin Zhao1Hong Zhao1、Lu Zhang 1Fen Gao2Yuwei Qin3Hubing Du
摘要:
我们提出了一种适用于一般干涉测量应用的相位展开方法。所提出的方法依赖于导数泽尼克多项式拟合DZPF技术其中相位近似为泽尼克多项式的组合。然后使用最小二乘法来估计拟合系数。因此相位展开问题被简化为这些系数的计算。由于采用了全场运算因此该方法快速有效。此外该方法直接提供所需相位而无需进一步计算未对准误差。该方法结合了相位展开和波前拟合过程。仿真和实验结果验证了该方法的潜力。
关键词相位展开算法导数Zernike多项式拟合干涉测量应用 1.引言
光学干涉仪可用于测量广泛的物理量。干涉仪提供的有趣数据中有条纹图案它是由被测量的波前像差调制的相位的余弦函数。在分析这些干涉图的过程中相位解调和相位展开都是必要的阶段。第一阶段通过相位测量算法[1]从原始干涉图中提取相位。通常第一阶段导出的相位是包裹的通常限制在-π到π的范围内。应使用相位展开器去除包裹相位中隐藏的不连续性。相位展开是相位提取过程中最具挑战性的最后一步。许多研究人员对此进行了探索并提出了几种方法如路径跟随方法[2-5]、最小二乘法[6-10]、区域法[11-13]或Regula。
在这里我们主要讨论由Ghiglia等人[6]首次引入的最小二乘技术将相位展开问题重新表述为离散泊松方程的解。已经提出了几种方法[910]来求解该方程。然后梁[7]介绍了一种基于模型的相位展开方法。该方法通过截断泰勒级数和残差函数来表示解包裹的相位函数。即使它是为展开MR图像而设计的它也应该同样适用于可以用低阶多项式近似的任何相位表面。之后Arines[8]提出了一种使用正交多项式如Zernike多项式进行Shack–Hartmann检验的迭代最小二乘模态估计方法。这更适合于一般的干涉测量应用但不能求解拟合过程中的第一个常数项。因此引入了一个迭代过程。
在综合研究上述方法的基础上本文提出了一种基于导数Zernike多项式拟合DZPF技术的高效非迭代相位展开算法。它本质上是一种最小二乘法但工作速度更快因为在展开过程中不需要迭代。该方法是在圆形条纹图案的背景下提出的尽管它也适用于其他类型的条纹图案。本文的第二节介绍了该算法的原理。第3节和第4节分别介绍并讨论了将所提出的算法应用于展开模拟和真实相位图的结果
2.拟用方法说明
图1显示了新的相位展开方法的流程图。下面解释详细的代数公式。通常相位展开问题可以建模如下 其中x y是真相位φx y是包裹相位k x y是要求解的整数。
对于一般的干涉测量测试测量的相位总是可以分解为在单位圆中定义的一系列正交Zernike多项式[15]的总和 其中Z x yi是以极坐标书写的方程3中定义的第i个Zernike多项式ci是相应的系数M是要使用的多项式的数量并且通常设置为等于36。 其中i 是模式排序编号而 其中指数n和m分别是径向度数和方位角频率并且满足m≤ n n− |m|even偶数。M与n的关系为M(n1)/(n2)12/2。通过计算真实相位和包裹相位的一阶差相位展开问题导致计算方程2中定义的系数c i(i1,2,3,...m)。因此我们正在寻找的阶段将是使以下成本函数最小化的阶段为了符号的简单性省略了xy依赖性 其中C(c1,c2,c3...cM)是包含拟合系数的向量W·是满足方程6的包裹算子。 计算成本函数的偏导数并将其设为零就能得到如图所示的正态方程。将其设为零我们就得到了如下所示的正则方程 式中: 图2:a 模拟条纹图案。b 以弧度表示的包裹相位。c 以弧度表示的展开相位。d 所提出方法的弧度残差。
其中A是对方程2中定义的相位场进行一阶导数运算的矩阵Δφ是包含x和y方向上的包裹相位差的向量T是矩阵换位。显式表达式如下: 其中Z_(ix)和Z_(iy)是第i个Zernike多项式在x和y方向上的导数W·是方程6中定义的包裹算子。k1,2,....m -1,l1,2,...n-1 其中 m n 分别是x方向和y方向上的像素数。
在最小二乘意义上求解正规方程我们可以得到方程10左侧定义的估计相位_cx y。
