两种提取Fabry-Perot干涉条纹圆心新方法.docVIP

两种提取Fabry-Perot干涉条纹圆心新方法 　　摘要本文介绍两种提取Fabry-Perot干涉条纹圆心点坐标的新方法。首先对干涉图像依次进行二值化处理，对所得到的条纹强度曲线进行均平滤波和自适应滤波，根据条纹灰度值强度余弦函数分布的特点，对条纹灰度值数据进行最小二乘法拟合，获得条纹强度峰值坐标，然后通过本文介绍的算法，进而获得Fabry-Perot干涉条纹圆心点的坐标。 　　【关键词】Fabry-Perot干涉条纹数字图像处理 最小二乘法 　　1 引言 　　在现代光学检测中经常要对干涉条纹图进行分析和研究。由于所使用的干涉测量原理、装置和对象的不同，条纹图像的差异很大，这给精确判定干涉条纹中心位置和级次带来很大困难，而且会严重影响测量精度。通过参考已有文献：“一种基于曲线拟合提取干涉条纹中心点的新方法” 和“一种确定干涉条纹中心的新方法” ，本文提出了一种迭代算法和一种均分干涉圆算法获取Fabry-Perot干涉条纹的圆心坐标，并将两种方法结果进行比较。实验结果表明，迭代算法具有高精度和运算速度快的特点。 　　2 基本原理 　　2.1 Fabry-Perot干涉仪成像原理 　　 法布里-帕罗多光束干涉仪是由平行放置的两块平面板组成的，在两板相对的平面上镀薄银膜或其它有较高反射系数的薄膜，要求镀膜面的平面度优于λ/20，而为消除两平板相背的平面上的反射光的干涉与我们所研究的干涉的重叠，故每块板都不是平行平面板，板的两面成一很小的夹角。 　　图1示出了Fabry-Perot干涉仪的工作原理，S为一扩展光源，L为准直透镜，G1和G2是两块平行平面板（板的两面有一很小的夹角），在两板相对的平面上镀有高反射率的薄银膜，G1和G2构成法-帕（F-P）腔，L1为聚焦透镜。S 上任一点发出的光线在腔内形成多次反射情况，结果构成多个平行的透射光，如图1所示。 　　在透射的诸光束中，相邻光束的光程差均为： 式中n为折射率，d为两块平行平面板间隔，i为入射光与板平面法线的夹角，由#8710;引起的位相差δ则为：#8710;=2ndcosi。 　　当光程差为波长的整数倍时，则透射光线的叠加形成干涉极大。由于两反射面是平行的，且S是扩展面光源，所以产生等倾干涉，即所有相同光程差且满足：mλ=2ndcosi 的入射光线在垂直于观察方向的平面上的轨迹是一组同心圆，因而在透镜L1焦面上得到的是一组同心环形条纹的图样（如图1中F’干涉环）。 　　2.2 确定干涉条纹圆心 　　确定干涉条纹圆心所使用的原理为Fabry-Perot干涉仪成像原理，以上已经讲过，在此不做多余阐述，使用的图像与数据处理软件为matlab7.0，在对图像进行预处理后，方可用一下方法确定圆心： 　　2.2.1 方法1 　　首先在干涉图上大体确定中心圆的坐标范围，然后在范围内引一条水平直线x=xi，与干涉条纹最内环的2个交点附近的强度峰值位置坐标为B11，B12，然后计算圆心的纵坐标为：。接着依次引水平直线 x=xi+nC（n=1,2,3……，C为一常量，受n的取值影响。本次实验中取n=10，计算后得C=30，即将圆的水平范围平均分成10份，n取的越大，实验结果越精确），用上述方法可求得n组圆心的纵坐标，求和并取平均值得。用同样的方法即可求得圆心横坐标，多次计算求平均，即可得到圆心坐标（）即为干涉条纹圆心点坐标。 　　 2.2.2 方法2 　　由于入射光来自空中同一区域，原则上Fabry-Perot干涉条纹是一系列同心圆（如图2），其中心位置的又一确定方法如下：首先同样对Fabry-Perot干涉图像进行二值化处理，对所得到的条纹强度曲线进行均平滤波和自适应滤波，在干涉图上选择一个近似的中心点，用一条过该点的水平直线将图像分为两半。对直线与干涉条纹交点附近的光强分布进行余弦函数拟合得到各个峰值的横坐标，通过非线性最小二乘法： 　　A，B，C和D的数值可以由拟合得到。假设水平直线与干涉条纹（最内环）的2个交点附近的强度峰值位置坐标为A11，A12，则一次迭代后圆心的横坐标为： 。再过横坐标Cx做一条垂直线，由上面相同的方法，我们可以得垂直直线与干涉条纹（最内环）的2个交点附近的强度峰值位置坐标为B11，B12，则可以的到圆心纵坐标为： 。再过纵坐标 Cy做一水平直线，用同样的方式重新确定Cx，重复这个过程，直至每相邻两次Cx，Cy的变化在设置的条件范围内，迭代可结束，（Cx，Cy）即为干涉条纹圆心点坐标。 　　3 实验与数据处理结果 　　对于方法1： 　　如图3为取x=430时的在X方向的条纹强度分布，可看到波峰最大的两个值所对应的坐标即可看做水平直线与灰度二维图最中间圆的两个交点，取交点的平均值即为该水平直线下的圆心纵坐标。实验数据最终如表1、2所示。 　　如表1、2，表格左边第一列