一种基于条纹图像频谱分析的干涉测量方法技术

本发明专利技术公开了一种基于条纹图像频谱分析的干涉测量方法。该方法首先获取待测元件经干涉仪生成的单幅加载频干涉条纹图像，然后对条纹图像进行延拓及边缘平滑处理，再对延拓条纹图像进行快速傅里叶变换(FFT)以获取其频谱。针对该频谱运用二维插值FFT算法估计频谱旁瓣峰值位置，并根据估计结果去除图像载频。最后通过计算得出待测波面的相位信息，实现待测光学元件的面形估计。本发明专利技术的方法基于经改进的条纹图像延拓方法与条纹载频估计方法进行条纹图像分析处理，具有更高的干涉测量精度。

A method of interferometry based on spectral analysis of stripe images

The invention discloses an interference measurement method based on spectral analysis of stripe images. First, we acquire the single loading frequency interference fringes generated by the interferometer, then extend the fringe pattern and smooth the edges. Then we use fast Fourier transform (FFT) to get the spectrum of the extended stripe images. The spectrum sidelobe peak position is estimated by using the two-dimensional interpolation FFT algorithm, and the image carrier frequency is removed according to the estimated results. At last, the phase information of the wave surface is calculated, and the surface shape of the optical element to be measured is realized. The method of the invention is based on the improved fringe image continuation method and fringe carrier frequency estimation method for fringe image analysis and processing, and has higher accuracy of interference measurement.

【技术实现步骤摘要】
一种基于条纹图像频谱分析的干涉测量方法
本专利技术属于光学测量领域，具体涉及到一种基于条纹图像频谱分析的干涉测量方法。
技术介绍
高精密度平面光学元件在光学工程领域有着极为广泛的应用。为了实现这些元件的面形检测、三维形貌获取及光学均匀性测试，通常采用基于光干涉原理的高精度干涉仪设备以及相应的干涉条纹测量方法。基于傅里叶变换处理干涉条纹图像是目前常用的一种干涉测量方法，其可以从单幅干涉条纹图像中获取待测光学元件的相关信息，降低了干涉系统的复杂程度，同时可以有效克服外界振动带来的干扰，适用于各类动态光学测量场合。然而在实际使用该方法进行干涉测量时，由于干涉条纹图像在进行傅里叶变换时会产生频谱能量泄露，并且基于频谱旁瓣峰值去除条纹图像载频的方法并不够准确，实际测算得出的波面相位值会存在较大的测量误差。近年来有多位学者针对图像频谱泄露问题提出了改进方案，但多以损失部分测量区域或增加大量运算处理时间为代价(楚兴春,吕海宝,赵尚弘.基于傅里叶变换的高准确度条纹细分方法；黄菁,朱日宏,陈磊.基于样本块匹配的干涉图延拓方法)。同时也有学者对条纹图像的载频准确估计方法进行了研究(粟银,范琦,王云飞.干涉条纹的高准确度傅里叶变换分析)，但其提出的基于旁瓣质心位置进行载频估计的方法存在物理意义不明确的不足，会对最终的相位值计算工作产生不利影响。
技术实现思路
本专利技术的目的在于提供一种基于条纹图像频谱分析的干涉测量方法。实现本专利技术目的的技术解决方案为：一种基于条纹图像频谱分析的干涉测量方法，具体步骤为：步骤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次迭代以得到对旁瓣谱峰位置本文档来自技高网...

【技术保护点】
一种基于条纹图像频谱分析的干涉测量方法，其特征在于，具体步骤为：步骤1、采集干涉条纹灰度图像；步骤2、对干涉条纹灰度图像进行延拓预处理得到干涉条纹延拓图像；步骤3、提取干涉条纹延拓图像延拓边缘区域的像素点，按条纹分布对其灰度值进行多项式线性回归处理，得到平滑处理后的干涉条纹延拓图像；步骤4、对平滑处理后的干涉条纹延拓图像进行频谱分析处理得到去载频图像频谱；步骤5、根据去载频图像频谱确定待测光学元件的面形分布。

【技术特征摘要】
1.一种基于条纹图像频谱分析的干涉测量方法，其特征在于，具体步骤为：步骤1、采集干涉条纹灰度图像；步骤2、对干涉条纹灰度图像进行延拓预处理得到干涉条纹延拓图像；步骤3、提取干涉条纹延拓图像延拓边缘区域的像素点，按条纹分布对其灰度值进行多项式线性回归处理，得到平滑处理后的干涉条纹延拓图像；步骤4、对平滑处理后的干涉条纹延拓图像进行频谱分析处理得到去载频图像频谱；步骤5、根据去载频图像频谱确定待测光学元件的面形分布。2.根据权利要求1所述的基于条纹图像频谱分析的干涉测量方法，其特征在于，步骤1中采集干涉条纹灰度图像，具体步骤为：步骤1.1、使用干涉仪装置针对待测光学元件进行测试，生成一幅加载频的干涉条纹图像；步骤1.2、对生成的干涉条纹图像进行采集并上传至计算机；步骤1.3、将干涉条纹图像转换为灰度图像，其灰度值分布表示如下：式中，a(x,y)和b(x,y)分别代表条纹图像的背景与调制度，fx0与fy0分别为x，y方向的空间载波频率，为待测相位。3.根据权利要求1所述的基于条纹图像频谱分析的干涉测量方法，其特征在于，步骤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，直至完成指定次数的迭代处理，得到干涉条纹延拓图像。4.根据权利要求1所述的基于条纹图像频谱分析的干涉测量方法，其特征在于，步骤3中提取干涉条纹延拓图像延拓边缘区域的像素点，按条纹分布对其灰度值进行多项式线性回归处理，得到边缘平滑处理后的干涉条纹延拓图像，具体步骤为：步骤3.1、查找干涉条纹图像中显示真实采集条纹的区域，并以此区域构建蒙版；步骤3.2、确定所有位于已知条纹区域边缘位置的w个像素点以这些边缘像素点为中心，向蒙版区域内外方向各提取n个像素点，其灰度值值记录为共计w组数据；步骤3.3、对每一组灰度值数据构建一元m阶多项式线性回归方程：Pik＝a0+a1i+a2i2+...+amim+εi(i＝-n,...,-1,0,1,...,n)各拟合点的关系用矩阵表示为：用矩阵符号表示为：按照最小二乘法求算多项式回归向量的系数估计值：分别求解w组数据中各像素点的灰度值估计值，具体计算公式为：步骤3.4、用各点的灰度值估计值替换原像素点灰度值数据，得到完成延拓及平滑预处理的干涉条纹图像。5.根据权利要求1所述的基于条纹图像频谱分析的干涉测量方法，其特征在于，步骤4中对平滑处理后的干涉条纹延拓图像进行频谱分析处理得到去载频图像频谱，具体步骤为：步骤4.1、对完成延拓及平滑预处理的干涉条纹图像进行二维快速傅里叶变换，得到如下式所示的频谱成分：

【专利技术属性】
技术研发人员：赵兆，阚凌志，
申请(专利权)人：南京理工大学，
类型：发明
国别省市：江苏,32