一种基于条纹图像频谱分析的干涉测量方法与流程

本发明属于光学测量领域，具体涉及到一种基于条纹图像频谱分析的干涉测量方法。

背景技术：

高精密度平面光学元件在光学工程领域有着极为广泛的应用。为了实现这些元件的面形检测、三维形貌获取及光学均匀性测试，通常采用基于光干涉原理的高精度干涉仪设备以及相应的干涉条纹测量方法。基于傅里叶变换处理干涉条纹图像是目前常用的一种干涉测量方法，其可以从单幅干涉条纹图像中获取待测光学元件的相关信息，降低了干涉系统的复杂程度，同时可以有效克服外界振动带来的干扰，适用于各类动态光学测量场合。

然而在实际使用该方法进行干涉测量时，由于干涉条纹图像在进行傅里叶变换时会产生频谱能量泄露，并且基于频谱旁瓣峰值去除条纹图像载频的方法并不够准确，实际测算得出的波面相位值会存在较大的测量误差。近年来有多位学者针对图像频谱泄露问题提出了改进方案，但多以损失部分测量区域或增加大量运算处理时间为代价(楚兴春,吕海宝,赵尚弘.基于傅里叶变换的高准确度条纹细分方法；黄菁,朱日宏,陈磊.基于样本块匹配的干涉图延拓方法)。同时也有学者对条纹图像的载频准确估计方法进行了研究(粟银,范琦,王云飞.干涉条纹的高准确度傅里叶变换分析)，但其提出的基于旁瓣质心位置进行载频估计的方法存在物理意义不明确的不足，会对最终的相位值计算工作产生不利影响。

技术实现要素：

本发明的目的在于提供一种基于条纹图像频谱分析的干涉测量方法。

实现本发明目的的技术解决方案为：一种基于条纹图像频谱分析的干涉测量方法，具体步骤为：

步骤1、采集干涉条纹灰度图像；

步骤2、对干涉条纹灰度图像进行延拓预处理得到干涉条纹延拓图像；

步骤3、提取干涉条纹延拓图像延拓边缘区域的像素点，按条纹分布对其灰度值进行多项式线性回归处理，得到平滑处理后的干涉条纹延拓图像；

步骤4、对平滑处理后的干涉条纹延拓图像进行频谱分析处理得到去载频图像频谱；

步骤5、根据去载频图像频谱计算出待测光学元件的面形分布。

本发明的方法与现有技术相比，其优点为：1)本发明的方法基于经改进的条纹图像延拓方法与条纹载频估计方法进行条纹图像分析处理，具有更高的干涉测量精度。2)本发明的方法处理流程简单，算法耗时短，有利于开发高效的软件产品。

下面结合附图对本发明作进一步详细描述。

附图说明

图1是本发明使用的基于条纹图像频谱分析的干涉测量方法流程图。

图2是本发明使用的干涉条纹图像延拓平滑预处理方法流程图。

图3是本发明使用的干涉条纹图像载频估计方法流程图。

图4是本发明实施例1采集的原始干涉条纹图像。

图5是本发明实施例1原始干涉条纹二维傅里叶变换频谱图。

图6为本发明实施例1中滤波后的频谱作逆变换得到的新图像。

图7是本发明实施例1原图像中非拓展条纹区域的像素值替换新图像中的相同位置得到的新图。

图8是本发明实施例1迭代50次后的干涉条纹延拓图像。

图9为本发明实施例1其中一组像素点分别在平滑前后的灰度值分布图，图9(a)为平滑前的灰度值分布图，图9(b)为平滑后的灰度值分布图。

图10是本发明实施例1完成延拓及平滑处理后的最终效果图。

图11为每一次迭代运算后的新偏移量δ随迭代次数q的变化情况。

图12是本发明实施例1最终得到的面形分布结果图及波差图，其中图12(a)为最终得到的面形分布结果图，图12(b)为最终得到的波差图。

下面结合具体实施例对本发明作进一步详细描述

具体实施方式

本发明提出了一种基于条纹图像频谱分析的干涉测量方法，具体步骤为：

步骤1、采集干涉条纹灰度图像，具体步骤为：

步骤1.1、使用干涉仪装置针对待测光学元件进行测试，生成一幅加载频的干涉条纹图像；

步骤1.2、对生成的干涉条纹图像进行采集并上传至计算机；

步骤1.3、将干涉条纹图像转换为灰度图像，其灰度值分布表示如下：

式中，a(x,y)和b(x,y)分别代表条纹图像的背景与调制度，fx0与fy0分别为x，y方向的空间载波频率，为待测相位。

步骤2、对干涉条纹灰度图像进行延拓预处理得到干涉条纹延拓图像，具体步骤为：

步骤2.1、对干涉条纹灰度图像进行二维快速傅里叶变换，干涉条纹图像灰度值分布的复数表达形式为：

i(x,y)＝a(x,y)+c(x,y)exp(j2πfx0x+j2πfy0y)+c^*(x,y)exp(-j2πfx0x-j2πfy0y)

其中，

对其进行二维傅里叶变换得到干涉条纹灰度图像的频谱成分：

i(fx,fy)＝a(fx,fy)+c(fx-fx0,fy-fy0)+c^*(fx+fx0,fy+fy0)

步骤2.2、从干涉条纹灰度图像的频谱成分中提取代表条纹信息的两个旁瓣成分c，c^*，并对其余频谱成分予以滤除，然后对滤除过后的频谱作逆变换得到一幅新的图像；

步骤2.3、用原条纹图像中的已知条纹区域像素灰度值替换新图像中相同位置的像素灰度值，得到替换处理后的图像；

步骤2.4、对替换处理后的图像重复进行步骤2.1、2.2、2.3，直至完成指定次数的迭代处理，得到干涉条纹延拓图像。

步骤3、提取干涉条纹延拓图像延拓边缘区域的像素点，按条纹分布对其灰度值进行多项式线性回归处理，得到平滑处理后的干涉条纹延拓图像，具体步骤为：

步骤3.1、查找干涉条纹图像中显示真实采集条纹的区域，并以此区域构建蒙版；

步骤3.2、确定所有位于已知条纹区域边缘位置的w个像素点以这些边缘点为中心，向蒙版区域内外方向各提取n个像素点，其灰度值值记录为共计w组数据；

步骤3.3、对每一组灰度值数据构建一元m阶多项式线性回归方程：

各拟合点的关系用矩阵表示为：

用矩阵符号表示为：

按照最小二乘法求算多项式回归向量的系数估计值：

分别求解w组数据中各像素点的灰度值估计值：

步骤3.4、用各点的灰度值估计值替换原像素点灰度值数据，得到完成延拓及平滑预处理的干涉条纹图像。

步骤4、对平滑处理后的干涉条纹延拓图像进行频谱分析处理得到去载频图像频谱；具体步骤为：

步骤4.1、对完成延拓及平滑预处理的干涉条纹图像进行二维快速傅里叶变换，得到如下式所示的频谱成分：

步骤4.2、滤除代表图像背景的基频成分a0(fx,fy)；

步骤4.3、基于二维插值fft算法对频谱进行精确的旁瓣峰值估计；具体步骤为：

步骤4.3.1、对干涉条纹图像二维频谱进行谱峰搜索，得到其正一级旁瓣谱峰坐标点i(x0,y0)；

步骤4.3.2、计算真实频率峰值在x、y方向上偏离i(x0,y0)的偏差值，得到对旁瓣谱峰位置的最终估计结果，具体方法为：

步骤4.3.2.1、在(x0-0.5,y0)与(x0+0.5,y0)、(x0,y0-0.5)与(x0,y0+0.5)、(x0-0.5,y0+0.5)与(x0+0.5,y0-0.5)、(x0-0.5,y0-0.5)与(x0+0.5,y0+0.5)添加4组共8个插值点；

步骤4.3.2.2、设定各组插值点的偏差值分量初值为

步骤4.3.2.3、根据各差值点的偏差值分量计算各组插值点的dft值，具体计算公式为：

步骤4.3.2.4、根据上述dft值计算各组插值点下的偏差修正值分量，具体计算公式为：

其中与分别表示每组数据中位于坐标轴正(负)方向上的插值点dft值；

各组插值点x、y方向上的偏差修正值，具体为：

得到各组插值点x、y方向上的偏差值分量，具体为：

步骤4.3.2.5、重复步骤4.3.2.3、4.3.2.4，直至完成预设q次迭代以得到对旁瓣谱峰位置的最终估计结果。

步骤4.4、依据步骤4.3中得到的旁瓣谱峰坐标对图像频谱旁瓣进行频移至零频位置，得到后续计算所需的去载频图像频谱。

步骤5、根据去载频图像频谱计算出待测光学元件的面形分布，具体步骤为：

步骤5.1、对去载频图像频谱中的条纹频谱成分进行傅里叶逆变换得到c0(x,y)；

步骤5.2、计算待测波面的相位分布情况，具体计算公式为：

其中re[c0(x,y)]和im[c0(x,y)]分别为c0(x,y)的实部和虚部；

步骤5.3、对步骤5.2计算得到的相位分布进行解包裹操作，即可得到代表待测光学元件面形分布的真实波面相位值分布情况。

下面结合实施例对本发明做进一步详细描述。

实施例1

本实施例中使用英特飞光电技术有限公司(nanjinginterferoopto-electronicstenhnologyco.,ltd)自主研发的25.4mm口径菲索型干涉仪搭建测试光路，对一组光学平面镜进行了干涉测量实验。该干涉仪采用635nm波长半导体激光器作为光源，同时搭载了高精度条纹图像分析装置。干涉仪实物如图1所示。具体的测量步骤如下所示：

步骤1、采集到干涉条纹图像，本实施例中采集的干涉条纹图像如图4所示。

步骤2、对干涉条纹灰度图像进行延拓预处理得到干涉条纹延拓图像，具体步骤为：

步骤2.1、对采集到的干涉条纹图像进行一次二维傅里叶变换处理，得到的频谱如图5所示。

步骤2.2、选定代表条纹信息的两个旁瓣区域，除去选定区域的其余所有频谱成分滤除；然后对滤波后的频谱作逆变换得到新图像，如图6所示。

步骤2.3、用原图像中非拓展条纹区域的像素值替换新图像中的相同位置，如图7所示。

步骤2.4、对替换处理后的图像重复进行步骤2.1、2.2、2.3，直至完成50次的迭代处理，得到干涉条纹延拓图像图，图8所示为迭代50次后的结果。

步骤3、对延拓图像进行平滑处理。本实施例中设置单边取样宽度n为20，并选取了所有边缘位置像素点灰度值数据组进行平滑。图9为其中一组像素点分别在平滑前后的灰度值分布图，通过对比可以看出平滑处理有效缓解了延拓图像边缘区域灰度值突变问题。完成延拓及平滑处理后的最终效果图如图10所示。

步骤4、对完成延拓及平滑处理的干涉条纹图像再次进行二维傅里叶变换，去除图像频谱零频背景后，对其旁瓣谱峰进行插值fft估计。图11所示为每一次迭代运算后的新偏移量δ随迭代次数q的变化情况，可以看出，x与y方向上的新偏移量在迭代次数超过15次后均趋向于0，此时真实谱峰位置的估计值趋于稳定，表明插值fft估计算法具有较好的收敛特性。最后按照插值fft估计的结果移除条纹载频。

步骤5、根据步骤4中得到的去载频图像频谱进行相位计算，得到的面形分布结果及波差图如图12所示。经计算得知，该光学平面镜的pv值为0.0632λ，rms值为0.0110λ。