CN108872153B - 基于非均匀傅里叶变换的平行平板光学均匀性的测量方法 - Google Patents
基于非均匀傅里叶变换的平行平板光学均匀性的测量方法 Download PDFInfo
- Publication number
- CN108872153B CN108872153B CN201810950390.7A CN201810950390A CN108872153B CN 108872153 B CN108872153 B CN 108872153B CN 201810950390 A CN201810950390 A CN 201810950390A CN 108872153 B CN108872153 B CN 108872153B
- Authority
- CN
- China
- Prior art keywords
- interference
- uniform
- phase
- fourier transform
- optical uniformity
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/41—Refractivity; Phase-affecting properties, e.g. optical path length
- G01N21/45—Refractivity; Phase-affecting properties, e.g. optical path length using interferometric methods; using Schlieren methods
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Instruments For Measurement Of Length By Optical Means (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公开了一种基于非均匀傅里叶变换的平行平板光学均匀性的测量方法,首先通过在菲索型波长移相干涉仪中,分别移相采样放置待测平行平板和空腔下的干涉图;其次,对干涉图各点的干涉光强数据,采用基于快速高斯网格法(FGG)的非均匀傅里叶变换,转化至频域进行频谱分析;最后,提取不同干涉腔长下对应的频率分量进行逆傅里叶变换获得干涉图的相位信息,将其恢复为波面信息后,计算得到待测平行平板的光学均匀性。本发明采用非均匀傅里叶变换频谱分析法,降低了波长移相时的非线性误差引起的测量误差,操作简单,步骤少,可用于大口径光学材料光学均匀性的高精度测量。
Description
技术领域
本发明涉及光干涉计量测试领域,特别是一种基于非均匀傅里叶变换的平行平板光学均匀性的测量方法。
背景技术
光学透射材料作为光学***中最常用的光学元件,其光学性能对整个光学***有着重大影响。光学均匀性是评价光学透射材料光学性能的一项重要指标,它反映了同一块光学材料内部折射率的不一致性。光学透射材料的光学均匀性会直接导致透射波前的变化,改变光学***的波像差,进而使光学***的输出产生偏差。因此,高精度的光学均匀性测量具有重要的理论和实际意义。
由于光学均匀性会引起透射光波的变形,因此可以通过干涉测量的方法,测量由于光学均匀性引起的波前改变量,从而得到材料的光学均匀性。干涉测量法包括直接测量法、样品翻转法和绝对测量法等。其中绝对测量法由于将待测材料的前后表面面形误差、干涉仪的***误差对结果的影响考虑在内,具有很高的检测精度。
基于波长调谐的时域傅里叶变换算法,通过波长移相干涉仪,对透射平行平板干涉腔及空腔进行两次测量,实现待测平板和参考面之间的干涉条纹的分离,从而获得平板前后表面信息和材料均匀性信息,是一种行之有效的绝对测量法。由于快速傅里叶变换(FFT)算法对时域的采样具有严格的要求,其数据分布局限于均匀分布的标准网格上,因此在测量中应保持移相量等于预设值,即在时域的采样点为等间隔均匀变化的,否则将导致波面测量误差。然而,由于波长调谐移相时存在一定的非线性误差,导致移相量非均匀变化,直接使用快速傅里叶变换算法会影响移相精度,需要额外的光学方法对波长改变量和电压之间的关系进行标定,其过程过于复杂、繁琐,使得测量效率下降。
非均匀快速傅里叶变换(NUFFT)是一种针对数据采样分布在非均匀网格上的情形,实现时域傅里叶变换算法的途径。NUFFT概念的首次出现是由Dutt等在《Fast Fouriertransforms for nonequispaced data》(SIAM J Sci Comput,1993,14(6):1368-1393)中提出的,他们对非均匀傅里叶变换的五种不同情况进行了研究与分析,并使用高斯窗函数插值解算其中的一种最常见的情况,给出了计算结果和误差分析,但由于插值过程窗函数插值过程重复性过高,因此计算量大,耗时长。Fessler等在《Nonuniform fast Fouriertransforms using min-max interpolation》(Signal Process,IEEE Trans.2003,51(2):560-574)中提出了基于min-max插值的NUFFT算法,该算法的主要是在最大误差最小的框架下来计算插值因子,使得NUFFT插值误差更小。针对NUFFT算法过采样过程需要大运算量和大物理储存的情况,Kunis在《Nonequispaced fast Fourier transforms withoutoversampling》(Proceedings in Applied Mathematics and Mechanics,2008)中提出了使用多维Taylor展开方法的不需要过采样的快速算法,节省了过采样方法需要消耗的大量存储空间,但仅仅限于某些点性能较好,其原因在于泰勒展开是自变量的局部领域内展开,只在较小的领域上保持较高的近似精度。
发明内容
本发明的目的在于提供一种基于非均匀傅里叶变换的平行平板光学均匀性的测量方法,用以降低波长移相过程中的非线性误差引起的测量误差,避免繁杂的标定过程,使得测量结果不受参考面面形精度的影响,精度高,易于实现。
实现本发明目的的技术解决方案为:一种基于非均匀傅里叶变换的平行平板光学均匀性的测量方法,包括以下步骤:
步骤1,在菲索波长移相干涉仪的干涉腔内,于干涉腔透射参考平面T和反射参考平面R之间放置待测平行平板;设置采集干涉图数目N及移相步长后进行移相采样,得到N幅干涉图的干涉光强数据,N=2n,n为满足测试精度所需的干涉图数目所取的正整数,所述每幅干涉图的干涉光强数据包括非多表面干涉叠加区域和多表面干涉叠加区域,转入步骤2;。
步骤2,将上述每幅干涉图的非多表面干涉叠加区域,通过相位恢复算法计算其中的移相量,所得结果记为原始非均匀采样序列xj,转入步骤3;
步骤3,将原始非均匀采样序列xj进行过采样,得到均匀采样序列xi,以此预计算高斯核函数用于卷积时的各指数部分值,并使用高斯核函数与步骤1中每幅干涉图中的多表面干涉区域的单个像素点的非均匀干涉数据fj卷积,得到高斯网格下的均匀离散值fτ,转入步骤4;
步骤4,对fτ做一维快速傅里叶变换和退卷积处理,得到频谱数据Fdeconv,完成对该像素点的非均匀傅里叶变换,转入步骤5;
步骤5,返回步骤3,直至遍历多表面干涉区域所有像素点,加窗提取频谱数据Fdeconv中的峰值并对其做逆傅里叶变换,得到各组干涉条纹对应的相位信息,将其转化为波面信息W1~W6,转入步骤6;
步骤6,空腔测量,撤去待测平行平板,保持干涉腔透射参考平面T和反射参考平面R的位置固定,进行波长调谐移相测量操作,得到波面信息W7,转入步骤7;
步骤7,综合步骤5和步骤6的测量结果,计算得到待测平行平板光学均匀性信息Δn。
本发明与现有测量技术相比,其显著优点:(1)采用两步绝对测量法,消除了参考面面形精度的影响,测量精度高;(2)测量过程操作简单,只需平板测量和空腔测量两个步骤;(3)使用快速高斯网格法的非均匀快速傅里叶变换处理移相采样干涉数据,抑制了波长调谐移相时的非线性误差对测量结果的影响,避免了人为标定移相量所带来的额外工作量。
附图说明
图1为本发明的基于非均匀傅里叶变换的平行平板光学均匀性的测量方法流程示意图。
图2为本发明的波长移相干涉原理示意图。
图3为本发明的波长移相采样所得干涉图。
图4为本发明的波长移相采样过程中存在非线性误差的移相量示意图。
图5为本发明的干涉信号经非均匀傅里叶变换后的频谱图。
具体实施方式
下面结合附图和具体实施方式对本发明做进一步说明。
结合图1,本发明的一种基于非均匀傅里叶变换的平行平板光学均匀性的测量方法,包括以下步骤:
步骤1,在菲索波长移相干涉仪的干涉腔内,于干涉腔透射参考平面T和反射参考平面R之间放置待测平行平板;设置采集干涉图数目N及移相步长后进行移相采样,得到N幅干涉图的干涉光强数据,N=2n,n为满足测试精度所需的干涉图数目所取的正整数,所述每幅干涉图的干涉光强数据包括非多表面干涉叠加区域和多表面干涉叠加区域。固定波长移相的移相步长进行N步移相采样,不考虑多次反射,则所述每幅干涉图中间部分的干涉光强数据由多表面干涉叠加形成,即为多表面干涉叠加区域,其余部分为非多表面干涉叠加区域。如图2,TF为透射参考平面,RF为反射参考平面,用T、A、B、R分别表示TF的反射面、待测平行平板的前、后表面、RF的反射面的面形,则共产生T和A、A和B、T和B、T和R、A和R与B和R六组干涉条纹,每组干涉光强可写成:
式中,hi为干涉腔长,λ0为初始波长,Δλ为每次移相时的波长变化量,t为移相次数。
步骤2,将上述每幅干涉图的非多表面干涉叠加区域,通过相位恢复算法计算其中的移相量,所得结果记为原始非均匀采样序列xj。如图3,干涉腔截面面积大于待测平行平板截面面积,仅有待测平行平板截面区域为多表面干涉叠加区域,使用掩膜将其提取出来,则剩余部分为简单的非多表面(二表面,T与R)干涉叠加区域,即单组干涉条纹数据。使用傅里叶变换相位测量法对该部分干涉数据进行相位重建。傅里叶变换相位测量法又称作附加空间线性载频相位测量法,使同为平行光的两束光波(分别为T与R的反射光)之间产生一个固定的夹角,在条纹图中引入线性变化的相位,为了简化计算,引入的相位一般只沿着x方向变化,因此,第i(i=1~N)幅干涉图中的单组干涉条纹光强数据为:
Ii(x)=ai+bi cos(φi+2πfxx) (4)
式中,φi为初始相位与移相量之和。fx表示单组干涉条纹光强数据的空间频率。运用欧拉公式对(4)式进行变形得:
对(5)式两边同时进行傅里叶变换,得到:
式中,aiδ(0)为直流项,干涉图的相位φi被锁定在幅度谱的峰值位置f=±fx中。选取合适的滤波器,取出频谱分布的正一级(或负一级)谱,令其为C(f),即:
对(7)式做逆傅里叶变换,得到:
对(8)式取实部与虚部进行运算,有:
对N幅干涉图进行同样的操作,可以将单组干涉条纹光强数据转化为相位数据φi,由此可以得到每幅干涉图之间的移相量为:
得到的原始非均匀采样序列如图4。
步骤3,将原始非均匀采样序列xj进行过采样,得到均匀采样序列xi,以此预计算高斯核函数用于卷积时的各指数部分值,并使用高斯核函数与步骤1中每幅干涉图中的多表面干涉区域的单个像素点的非均匀干涉数据fj卷积,得到高斯网格下的均匀离散值fτ。设置过采样率R、过采样点数Mτ、过采样间隔等参数,建立高斯均匀网格点。对原始非均匀采样序列xj中的各点,离xj最近的高斯均匀网格点(xj≤xi)的集合即为所求均匀采样序列xi。设置高斯核函数对非均匀干涉数据fj进行卷积平滑处理,高斯核函数可表示为:
gτ(x)=exp(-x2/4τ) (12)
式中,τ为高斯核函数参数,它决定了核函数的指数衰减速率。τ由下式决定:
式中,M为非均匀干涉数据长度,其值与N相同,MSP为高斯核函数单边可覆盖的均匀网格点数,其取值由所需的计算精度确定。利用高斯核函数衰减速度快的特点,计算高斯均匀采样点数据值时只需考虑该点附近2MSP个点的干涉强度的影响,因此需要预计算高斯核函数在卷积过程中的指数部分用以下一步骤的计算,公式如下:
E1j=exp[-(xj-xi)2/4τ] (14)
E2j=exp[(xj-xi)π/Mττ] (15)
E3m′=exp[-(πm′/Mτ)2/τ],-MSP<m′≤MSP (16)
卷积过程公式如下:
其中,i=1~Mτ。遍历整个xj序列后将fτ各点的值累加,获得高斯网格下的均匀离散值fτ。
步骤4,对高斯网格值fτ做一维快速傅里叶变换和退卷积处理,得到频谱数据Fdeconv,完成对该点的非均匀傅里叶变换。高斯网格下的均匀离散值fτ,对fτ进行均匀快速傅里叶变换(FFT)得到Fτ,再对Fτ做退卷积处理消除高斯核函数的平滑影响,公式如下:
式中,k为退卷积区间内的坐标值。
步骤5,返回步骤3,直至遍历多表面干涉区域所有像素点,加窗提取频谱数据Fdeconv中的峰值并对其做逆傅里叶变换,得到各组干涉条纹对应的相位信息,将其转化为波面信息W1~W6。如图5,每个像素点的频谱数据Fdeconv中的峰值对应各组干涉条纹的腔长值,不同的频率峰值加窗提取的复值包含干涉图的相位信息,通过逆傅里叶变换从中解调,对所有像素点的Fdeconv中的每个频率峰值进行上述操作,得到整个多表面干涉区域分别在六个不同条纹频率下的相位信息。再经过相位解包、消倾斜处理,获得对应的波面信息W1~W6。W1~W6与干涉腔长的关系如下表所示:
表1 各组条纹的光程差以及干涉腔长
表中,n0为待测平行平板折射率,t0为待测平行平板厚度。
步骤6,空腔测量,撤去待测平行平板,保持干涉腔透射参考平面T和反射参考平面R的位置固定,进行波长调谐移相测量操作,得到波面信息W7。易知此时干涉表面为T与R,干涉腔长为h6,则:
W7=2R-2T (19)
h6=2(L1+L2+t) (20)
步骤7,综合步骤5、6所得的测量结果,计算得到待测平行平板光学均匀性信息Δn。由表1与(19)式可得:
Claims (5)
1.一种基于非均匀傅里叶变换的平行平板光学均匀性的测量方法,其特征在于,包括以下步骤:
步骤1,在菲索波长移相干涉仪的干涉腔内,于干涉腔透射参考平面T和反射参考平面R之间放置待测平行平板;设置采集干涉图数目N及移相步长后进行移相采样,得到N幅干涉图的干涉光强数据,N=2n,n为满足测试精度所需的干涉图数目所取的正整数,所述每幅干涉图的干涉光强数据包括非多表面干涉叠加区域和多表面干涉叠加区域,转入步骤2;
步骤2,将上述每幅干涉图的非多表面干涉叠加区域,通过相位恢复算法计算其中的移相量,所得结果记为原始非均匀采样序列xj,转入步骤3;
步骤3,将原始非均匀采样序列xj进行过采样,得到均匀采样序列xi,以此预计算高斯核函数用于卷积时的各指数部分值,并使用高斯核函数与步骤1中每幅干涉图中的多表面干涉区域的单个像素点的非均匀干涉数据fj卷积,得到高斯网格下的均匀离散值fτ,公式如下:
式中,i=1~Mτ,Mτ为过采样点数,MSP为高斯核函数单边可覆盖的均匀网格点数,E1、E2、E3为高斯核函数用于卷积时预计算的指数部分值;
转入步骤4;
步骤4,对fτ做一维快速傅里叶变换和退卷积处理,得到频谱数据Fdeconv,完成对该像素点的非均匀傅里叶变换,转入步骤5;
步骤5,返回步骤3,直至遍历多表面干涉区域所有像素点,加窗提取频谱数据Fdeconv中的峰值并对其做逆傅里叶变换,得到各组干涉条纹对应的相位信息,将其转化为波面信息W1~W6,转入步骤6;
步骤6,空腔测量,撤去待测平行平板,保持干涉腔透射参考平面T和反射参考平面R的位置固定,进行波长调谐移相测量操作,得到波面信息W7,转入步骤7;
步骤7,综合步骤5和步骤6的测量结果,计算得到待测平行平板光学均匀性信息Δn。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810950390.7A CN108872153B (zh) | 2018-08-20 | 2018-08-20 | 基于非均匀傅里叶变换的平行平板光学均匀性的测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810950390.7A CN108872153B (zh) | 2018-08-20 | 2018-08-20 | 基于非均匀傅里叶变换的平行平板光学均匀性的测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108872153A CN108872153A (zh) | 2018-11-23 |
CN108872153B true CN108872153B (zh) | 2020-12-11 |
Family
ID=64320966
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810950390.7A Active CN108872153B (zh) | 2018-08-20 | 2018-08-20 | 基于非均匀傅里叶变换的平行平板光学均匀性的测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108872153B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113446930A (zh) * | 2020-03-26 | 2021-09-28 | 山东大学 | 一种基于希尔伯特变换的白光干涉信号非均匀采样的校正方法 |
CN111562088A (zh) * | 2020-04-30 | 2020-08-21 | 南京理工大学 | 基于采样函数的平行平板光学参数的测量方法 |
CN111664800B (zh) * | 2020-05-19 | 2022-01-07 | 上海大学 | 一种平行平板多表面检测方法及夹具 |
CN111964876B (zh) * | 2020-07-29 | 2021-11-23 | 南京理工大学 | 基于lrte-nufft的平行平板光学均匀性的测量方法 |
CN112097678B (zh) * | 2020-08-31 | 2022-07-08 | 上海大学 | 基于频率盲估计的多表面面形测量方法 |
CN112880981B (zh) * | 2021-01-25 | 2022-12-23 | 上海大学 | 一种自适应调节腔长的多表面干涉测量方法 |
CN112880569B (zh) * | 2021-01-25 | 2022-11-08 | 上海大学 | 一种基于腔长校正的多表面测量方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101000464A (zh) * | 2006-01-10 | 2007-07-18 | 比斯泰克半导体***股份有限公司 | 用于照亮并检查表面的装置 |
EP1817542A1 (en) * | 2004-10-06 | 2007-08-15 | Rudolph Technologies, Inc. | Measuring elastic moduli of dielectric thin films using an optical metrology system |
JP6035794B2 (ja) * | 2012-03-14 | 2016-11-30 | 大日本印刷株式会社 | ムラ情報生成装置、ムラ情報生成プログラム及びムラ情報生成方法 |
CN105092530B (zh) * | 2015-05-21 | 2018-01-05 | 南京理工大学 | 平行平晶光学非均匀性的绝对测量方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5974160A (en) * | 1993-10-26 | 1999-10-26 | Asahi Kasei Kogyo Kabushiki Kaisha | Measuring method and apparatus of gloss irregularity and printing unevenness |
FR2953029B1 (fr) * | 2009-11-25 | 2011-11-18 | Draka Comteq France | Fibre optique multimode a tres large bande passante avec une interface coeur-gaine optimisee |
CN103335982B (zh) * | 2013-06-21 | 2016-05-11 | 中国科学院上海光学精密机械研究所 | 利用波长调谐移相干涉仪测量平行平板光学均匀性的方法 |
CN103454249B (zh) * | 2013-09-16 | 2016-03-02 | 南京理工大学 | 基于白光干涉的光学玻璃均匀性检测方法及装置 |
CN106501216B (zh) * | 2016-12-27 | 2019-04-16 | 南京理工大学 | 一种光学平行平板玻璃的均匀性绝对测量方法 |
-
2018
- 2018-08-20 CN CN201810950390.7A patent/CN108872153B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1817542A1 (en) * | 2004-10-06 | 2007-08-15 | Rudolph Technologies, Inc. | Measuring elastic moduli of dielectric thin films using an optical metrology system |
CN101000464A (zh) * | 2006-01-10 | 2007-07-18 | 比斯泰克半导体***股份有限公司 | 用于照亮并检查表面的装置 |
JP6035794B2 (ja) * | 2012-03-14 | 2016-11-30 | 大日本印刷株式会社 | ムラ情報生成装置、ムラ情報生成プログラム及びムラ情報生成方法 |
CN105092530B (zh) * | 2015-05-21 | 2018-01-05 | 南京理工大学 | 平行平晶光学非均匀性的绝对测量方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108872153A (zh) | 2018-11-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108872153B (zh) | 基于非均匀傅里叶变换的平行平板光学均匀性的测量方法 | |
Surrel | Fringe analysis | |
Kujawinska et al. | High accuracy Fourier transform fringe pattern analysis | |
CN107917676B (zh) | 一种基于条纹图像频谱分析的干涉测量方法 | |
AU659271B2 (en) | Spatial wavefront evaluation by intensity relationship | |
CN112097678B (zh) | 基于频率盲估计的多表面面形测量方法 | |
JP6207250B2 (ja) | サンプルの性状を測定する方法及び装置 | |
CN111879730A (zh) | 基于矩形窗函数优化的光学相干层析成像信号处理方法 | |
CN111562088A (zh) | 基于采样函数的平行平板光学参数的测量方法 | |
Schmit et al. | Spatial and temporal phase-measurement techniques: a comparison of major error sources in one dimension | |
CN108170636B (zh) | 一种从线性载频移相干涉图中恢复出真实相位的方法 | |
CN111811398B (zh) | 一种基于相移特征多项式高精度拟合的多表面测量方法 | |
CN104316204B (zh) | 一种高精度的变频干涉的四步相移标定方法 | |
CN111964876B (zh) | 基于lrte-nufft的平行平板光学均匀性的测量方法 | |
CN101968383B (zh) | 一种抗扰动的时频域波前检测方法 | |
CN114964033B (zh) | 一种使用波长移相法测量超表面形貌分布的方法 | |
Chang et al. | Wavelength-tuning phase-shifting interferometry of transparent plates using sub-signal frequency correction | |
CN110500968B (zh) | 基于稀疏傅里叶变换的数字莫尔干涉相位实时测量方法 | |
Barajas et al. | Towards an on-chip signal processing solution for the online calibration of SS-OCT systems | |
Yang et al. | A radial-shearing interference system of testing laser-pulse wavefront distortion and the original wavefront reconstructing | |
CN118089533A (zh) | 基于频偏修正的多干涉谐波叠加信号频率及相位参数重建 | |
CN116586640B (zh) | 球面测试板及其制作方法以及干涉仪传递函数的标定方法 | |
Cui et al. | Real time phase unwrapping implementation based on Unscented Kalman filter for interferometric fiber sensor applications | |
CN115307746B (zh) | 一种单色光的干涉测量相位计算方法 | |
CN111397642B (zh) | 基于相位保留滤波和峰位置追踪的低相干干涉解调方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |