CN105527027B - 亚微米孔径光纤点衍射波前测量的结构误差校正方法 - Google Patents
亚微米孔径光纤点衍射波前测量的结构误差校正方法 Download PDFInfo
- Publication number
- CN105527027B CN105527027B CN201610008014.7A CN201610008014A CN105527027B CN 105527027 B CN105527027 B CN 105527027B CN 201610008014 A CN201610008014 A CN 201610008014A CN 105527027 B CN105527027 B CN 105527027B
- Authority
- CN
- China
- Prior art keywords
- diffraction
- point
- sub
- optical fiber
- pore size
- 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.)
- Expired - Fee Related
Links
- 239000013307 optical fiber Substances 0.000 title claims abstract description 79
- 239000011148 porous material Substances 0.000 title claims abstract description 61
- 238000005259 measurement Methods 0.000 title claims abstract description 19
- 238000010008 shearing Methods 0.000 claims abstract description 55
- 239000000523 sample Substances 0.000 claims abstract description 37
- 238000000034 method Methods 0.000 claims abstract description 25
- 238000003384 imaging method Methods 0.000 claims abstract description 5
- 238000012937 correction Methods 0.000 claims description 21
- 238000006073 displacement reaction Methods 0.000 claims description 9
- 230000003287 optical effect Effects 0.000 claims description 8
- 238000001514 detection method Methods 0.000 claims description 4
- 238000005457 optimization Methods 0.000 claims description 4
- 230000017105 transposition Effects 0.000 claims description 4
- 230000001427 coherent effect Effects 0.000 claims description 3
- 238000000605 extraction Methods 0.000 claims description 3
- 239000000835 fiber Substances 0.000 description 9
- 238000012360 testing method Methods 0.000 description 5
- 239000000919 ceramic Substances 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 230000009897 systematic effect Effects 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 241001269238 Data Species 0.000 description 2
- 230000004075 alteration Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 241001614291 Anoplistes Species 0.000 description 1
- 241000208340 Araliaceae Species 0.000 description 1
- 206010073261 Ovarian theca cell tumour Diseases 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000005305 interferometry Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 238000004064 recycling Methods 0.000 description 1
- 208000001644 thecoma Diseases 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J9/00—Measuring optical phase difference; Determining degree of coherence; Measuring optical wavelength
- G01J9/02—Measuring optical phase difference; Determining degree of coherence; Measuring optical wavelength by interferometric methods
- G01J9/0215—Measuring optical phase difference; Determining degree of coherence; Measuring optical wavelength by interferometric methods by shearing interferometric methods
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J9/00—Measuring optical phase difference; Determining degree of coherence; Measuring optical wavelength
- G01J9/02—Measuring optical phase difference; Determining degree of coherence; Measuring optical wavelength by interferometric methods
- G01J2009/0223—Common path interferometry; Point diffraction interferometry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J9/00—Measuring optical phase difference; Determining degree of coherence; Measuring optical wavelength
- G01J9/02—Measuring optical phase difference; Determining degree of coherence; Measuring optical wavelength by interferometric methods
- G01J2009/0242—Compensator
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Physics & Mathematics (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明涉及一种亚微米孔径光纤点衍射波前测量的结构误差校正方法,通过无成像镜头的CCD探测器获取原始剪切波前数据,利用三维坐标重构对剪切波面数据进行预校正,再获取180度旋转测量探头后的预校正数据,将二者叠加,消除两点衍射源偏移引入的结构误差,得到真实剪切波面数据,再根据差分泽尼克多项式拟合方法拟合得到待测的亚微米孔径光纤点衍射波前,实现亚微米孔径光纤点衍射波前的高精度测量。
Description
技术领域
本发明属于光学测量技术领域,尤其涉及一种亚微米孔径光纤点衍射波前测量的结构误差校正方法。
背景技术
随着光学制造技术的发展,干涉测量技术已得到广泛应用。传统的干涉仪,如斐索干涉仪和泰曼-格林干涉仪,一般采用标准透镜来产生参考波前,其可达到的测量精度受标准透镜精度限制。点衍射干涉仪利用一点衍射波前作为理想的球面参考波前,克服了由于标准元件的精度限制,在高精度球面和平面等面形检测以及三维坐标绝对测量中得到了广泛应用。相比于传统的针孔点衍射干涉仪以及单模光纤点衍射干涉仪,亚微米孔径光纤点衍射干涉仪由于可同时获得大数值孔径波前和高点衍射光能量,极大的扩展了点衍射干涉仪的应用范围。
由于实现点衍射干涉仪测试的精度主要取决于衍射波前的球形,球形的衍射波前分析已成为评估点衍射干涉仪的性能的根本途径。各种实验测试提出了如何测量点衍射波前的球形,其中大部分都是基于剪切法和混合法。通常情况下,混合方法需要几个测量与测试光学旋转和位移,它对环境干扰很敏感,不能完全独立的系统误差。剪切法是一个自我参照的方法,它被广泛地应用于重建点衍射波前校正点衍射干涉仪。然而,它需要系统结构参数等先验知识,以便于去除高阶像差,特别在高数值孔径的情况下,点衍射源偏移量不可忽略,会造成较大的测量误差。
发明内容
本发明为了解决上述技术问题,提供一种亚微米孔径光纤点衍射波前测量的结构误差校正方法,其针对两点衍射源间距结构参量与结构误差存在奇函数关系,提出了基于三维坐标重构以及对称位置补偿的结构误差校正方法,解决了点衍射源偏移引入的结构误差问题,并使用差分泽尼克多项式拟合方法,进而实现亚微米孔径光纤点衍射波前的高精度测量。
本发明的上述技术问题主要是通过下述技术方案得以解决的:本发明亚微米孔径光纤点衍射波前测量的结构误差校正方法,包括如下步骤:
(1-1)将两根具有相同结构尺寸的亚微米孔径光纤安装于亚波长孔径光纤点衍射干涉仪的测量探头中,所述的亚微米孔径光纤的前端头呈锥形,构成锥形探针出光端面,两个锥形探针出光端面并排且共面,水平放置所述的测量探头,并且使两个所述的锥形探针出光端面的连线呈水平方向;
(1-2)将所述的亚波长孔径光纤点衍射干涉仪的参考路和检测路的两路相干光分别导入两根亚微米孔径光纤,采用不带成像镜头的CCD探测器采集来自于两根亚微米孔径光纤的点衍射干涉图,利用多步移相算法解调出对应的点衍射剪切波面数据ΔW(s),其中s为两亚微米孔径光纤点衍射源的横向偏移量;
(1-3)在点衍射剪切波面数据ΔW(s)中选取6个以上像素点处的相位值,组成非线性超定方程组,运用迭代优化算法重构出两亚微米孔径光纤点衍射源的三维坐标为(x1,y1,z1)和(x2,y2,z2),已知所述的CCD探测器上每像素点坐标为(x,y,z),根据光程差公式:
生成理想波面数据ΔW(s)′;
(1-4)将所述的测量探头绕其轴线旋转180度,使两亚微米孔径光纤点衍射源的横向偏移量变为-s,保持CCD探测器位置不变,重复步骤(1-2),得到另一组点衍射剪切波面数据ΔW(-s),再用点衍射剪切波面数据ΔW(-s)重复步骤(1-3),得到对应的另一组理想波面数据ΔW(-s)′;
(1-5)利用基于三坐标重构得到的两组理想波面数据(ΔW(s)′,ΔW(-s)′),对测量得到的原始点衍射剪切波面数据(ΔW(s),ΔW(-s))进行第一步校正:
得到初步校正后的点衍射剪切波面数据
(1-6)将初步校正后的点衍射剪切波面数据进行叠加,实现点衍射剪切波面数据的第二步校正,得到不存在结构误差的真实点衍射剪切波面ΔW:
作为优选,所述的亚微米孔径光纤点衍射波前测量的结构误差校正方法,包括亚波长孔径光纤点衍射波前重构方法:
(2-1)利用步骤(1-3)所重构得到的两亚微米孔径光纤点衍射源的三维坐标(x1,y1,z1)和(x2,y2,z2),计算得到两亚微米孔径光纤点衍射源的横向偏移量s,
并据此得到相应的差分泽尼克多项式ΔZ={ΔZi},
ΔZi(x,y)=Zi(x,y)-Zi(x+s,y)
其中,下标i表示多项式的项数,{Zi}表示泽尼克多项式,坐标(x,y)为波面上各点坐标数据;
(2-2)利用最小二乘法对步骤(1-6)所获得的真实点衍射剪切波面ΔW进行差分泽尼克多项式拟合,得到差分泽尼克多项式系数a={ai},
a=(ΔZTΔZ)-1ΔZTΔW,
其中T表示转置运算;
(2-3)利用差分泽尼克多项式系数a={ai},即可重构得到实际待测的亚微米孔径光纤点衍射波面数据W0,
本发明通过CCD探测器获取原始剪切波前数据,利用三维坐标对剪切波前数据进行了预校正,再获取旋转测量探头后的预校正数据,将二者叠加,即可消除点衍射源偏移引入的高阶像差,得到真实的点衍射剪切波前数据,并根据差分泽尼克多项式拟合方法得到亚微米孔径光纤点衍射波前的高精度测量。
本发明的有益效果是:操作简单、高效,可通用化,可实现超高精度的校正,适用于任何数值孔径的点衍射探头的测量,在无需了解关于干涉仪特性、点衍射源数值孔径等先验信息的情况下,即可实现结构误差的高精度校正,实现亚微米孔径光纤点衍射波前的高精度测量,可在实际点衍射干涉仪测量中得到很好的应用。
附图说明
图1是本发明中测量探头的一种立体结构示意图。
图2是本发明中亚波长孔径光纤点衍射干涉仪的一种连接结构示意图。
图3是点衍射源偏移引入结构误差的几何模型示意图。
图4是点衍射源偏移量为s的原始的点衍射剪切波面数据图。
图5是点衍射源偏移量为-s的原始的点衍射剪切波面数据图。
图6是点衍射源偏移量为s预校正剪切波面数据。
图7是点源偏移量为-s预校正剪切波面数据。
图8是本发明实施例经过两步校正后的实验结果。
图9是本发明实施例经过亚波长孔径光纤点衍射波前重构方法的实验结果。
图中1.激光器,2.第一半波片,3.第二半波片,4.分光棱镜,5.第一光纤耦合器,6.第二光纤耦合器,7.测量探头,8.1/4波片,9.反射镜,10.压电陶瓷,11.电脑,12.亚微米孔径光纤,13.CCD探测器。
具体实施方式
下面通过实施例,并结合附图,对本发明的技术方案作进一步具体的说明。
实施例:本实施例的亚微米孔径光纤点衍射波前测量的结构误差校正方法,采用亚波长孔径光纤点衍射干涉仪,如图2所示,亚波长孔径光纤点衍射干涉仪包括激光器1、第一半波片2、分光棱镜4、第二半波片3、第一光纤耦合器5、第二光纤耦合器6、测量探头7和1/4波片8、反射镜9、压电陶瓷10,激光器1发出的光经第一半波片2到分光棱镜4,经分光棱镜4分光,一路光路经第二半波片3到第一光纤耦合器5,另一路光路到第二光纤耦合器6,还有一路经1/4波片8到反射镜9,反射镜9放置在压电陶瓷10上,压电陶瓷10和电脑11相连。如图1所示,测量探头7中安装有两根具有相同结构尺寸的亚微米孔径光纤12,亚微米孔径光纤12的前端头呈锥形,构成锥形探针出光端面,两个锥形探针出光端面并排且共面,两根亚微米孔径光纤的后端头分别和第一光纤耦合器5、第二光纤耦合器6相连。不带成像镜头的CCD探测器13放置在和测量探头7的出光端面相对的位置处,CCD探测器13的数据线和电脑11相连。
亚微米孔径光纤点衍射波前测量的结构误差校正方法,包括如下步骤:
(1-1)如图1所示,将两根具有相同结构尺寸的亚微米孔径光纤安装于亚波长孔径光纤点衍射干涉仪的测量探头中,亚微米孔径光纤的前端头呈锥形,构成锥形探针出光端面,两个锥形探针出光端面并排且共面,水平放置测量探头,并且使两个锥形探针出光端面的连线呈水平方向;
(1-2)将亚波长孔径光纤点衍射干涉仪的参考路和检测路的两路相干光分别导入两根亚微米孔径光纤,采用不带成像镜头的CCD探测器采集来自于两根亚微米孔径光纤的点衍射干涉图,利用多步移相算法解调出对应的点衍射剪切波面数据ΔW(s),其中s为两亚微米孔径光纤点衍射源的横向偏移量;
(1-3)在点衍射剪切波面数据ΔW(s)中选取6个以上像素点处的相位值,组成非线性超定方程组,运用迭代优化算法重构出两亚微米孔径光纤点衍射源的三维坐标为(x1,y1,z1)和(x2,y2,z2),已知CCD探测器上每像素点坐标为(x,y,z),根据光程差公式:
生成理想波面数据ΔW(s)′;
(1-4)将测量探头绕其轴线旋转180度,使两亚微米孔径光纤点衍射源的横向偏移量变为-s,保持CCD探测器位置不变,重复步骤(1-2),得到另一组点衍射剪切波面数据ΔW(-s),再利用该点衍射剪切波面数据ΔW(-s)重复步骤(1-3),得到对应的另一组理想波面数据ΔW(-s)′;
(1-5)利用基于三坐标重构得到的两组理想波面数据(ΔW(s)′,ΔW(-s)′),对测量得到的原始点衍射剪切波面数据(ΔW(s),ΔW(-s))进行第一步校正:
得到初步校正后的点衍射剪切波面数据
如图3所示,是两根亚微米孔径光纤的几何误差分析图,S1、S2代表两根亚微米孔径光纤,定义测量探头出光端面和CCD探测器的探测面之间的距离为D,S1是坐标系统的原点。S1与S2之间的距离(点源间距)为s,S1与S2到CCD探测器上一个随机点P的距离分别为R1与R2,可以得到S1、S2和P点坐标分别为(0,0,0)、(s,0,0)和(x,y,D),R1与R2光程差ΔR:
为了简化分析,定义r为极半径,θ为极角,将直角坐标系下光程差变为极坐标系光程差得:
定义Rm为CCD探测器的半径,NA为测量光束的最大数值孔径,假设:
其中,ρ代表归一化半径,采用泽尼克多项式表示光程差ΔR得:
ΔR=a2Z2+a9Z9+a19Z19+a33Z33
其中Z2为x方向的倾斜项,Z9、Z19和Z33是泽尼克多项式的一级、二级和x方向的彗差项,a2,a9,a19和a33分别为对应的系数,并得到:
由上式可知,泽尼克系数均为奇函数,定义ΔR(s)为点源偏移量为s所引入的剪切波前像差,ΔR(-s)为点源偏移量为-s所引入的剪切波前像差,得ΔR(s)和ΔR(-s)是点源偏移量s的奇函数,进而得到:
ΔR(s)+ΔR(-s)≌0
定义偏移量s和-s的对应初步校正的剪切波前数据为和待测剪切波前为ΔW可以得到:
(1-6)将步骤(1-5)处理得到初步校正后的点衍射剪切波面数据进行叠加,实现点衍射剪切波面数据的第二步校正,得到不含结构误差的真实点衍射剪切波面ΔW:
亚波长孔径光纤点衍射波前重构方法为:
(2-1)利用上述步骤(1-3)所重构得到的两亚微米孔径光纤点衍射源的三维坐标(x1,y1,z1)和(x2,y2,z2),计算得到两亚微米孔径光纤点衍射源的横向偏移量s,
并据此得到相应的差分泽尼克多项式ΔZ={ΔZi},
ΔZi(x,y)=Zi(x,y)-Zi(x+s,y)
其中,下标i表示多项式的项数,{Zi}表示泽尼克多项式,坐标(x,y)为波面上各点坐标数据。
(2-2)利用最小二乘法对上述步骤(1-6)所获得的真实点衍射剪切波面ΔW进行差分泽尼克多项式拟合,得到差分泽尼克多项式系数a={ai},
a=(ΔZTΔZ)-1ΔZTΔW
其中T表示转置运算;
(2-3)利用步骤(2-2)得到的差分泽尼克多项式系数a={ai},即可高精度重构得到实际待测的亚微米孔径光纤点衍射波面数据W0,
本发明由点源间距的结构参量与结构误差的奇函数关系,提出了对称的系统误差消除方法,解决了点衍射源偏移引入的结构误差问题。在点衍射源偏移量超过50μm,数值孔径NA为0.1的前提下,用传统泽尼克拟合的方法,残余误差均方根大于2.0×10-3λ,并随着数值孔径的不断增大,传统方法的残余误差均方根将显著增大。在点源偏移量超过300μm,数值孔径NA为0.6情况下,用本发明提出的方法进行点衍射波前测量的残余误差均方根小于1.0×10-4λ。
再以一个具体的测量实例作进一步说明。
采用本发明的方法对点源偏移量是250μm、数值孔径NA为0.60的亚波长孔径光纤点衍射干涉仪的测量探头进行测量校正,包括如下步骤:
(1-1)如图1所示,将两根具有相同结构尺寸的亚微米孔径光纤安装于亚波长孔径光纤点衍射干涉仪的测量探头中,亚微米孔径光纤的前端头呈锥形,构成锥形探针出光端面,两个锥形探针出光端面并排且共面安装于一个测量探头中,水平放置测量探头,并且使两个锥形探针出光端面的连线呈水平方向;
(1-2)如图2所示,将亚波长孔径光纤点衍射干涉仪的参考路和检测路的两路相干光分别导入两根亚微米孔径光纤,并利用丹麦JAI RMC-4200GE CCD探测器采集来自于两亚微米孔径光纤的点衍射干涉图,利用多步移相算法解调出对应的点衍射剪切波面数据ΔW(s),如图4所示,两亚微米孔径光纤点衍射源的横向偏移量s=250μm;
(1-3)在点衍射剪切波面数据ΔW(s)中选取6个以上像素点处的相位值,组成非线性超定方程组,运用迭代优化算法重构出两亚微米孔径光纤点衍射源的三维坐标为(0.125,0,300)和(-0.125,0,300),已知CCD探测器上每像素点坐标为(x,y,z),根据光程差公式:
生成理想波面数据ΔW(s)′;
(1-4)将测量探头绕其轴线旋转180度,使两亚微米孔径光纤点衍射源的横向偏移量变为-s,保持CCD探测器位置不变,重复步骤(1-2),得到另一组点衍射剪切波面数据ΔW(-s),如图5所示,再利用该波面数据ΔW(-s)重复步骤(1-3),得到对应的另一组理想波面数据ΔW(-s)′;
(1-5)利用基于三坐标重构得到的两组理想波面数据(ΔW(s)′,ΔW(-s)′),对测量得到的原始点衍射剪切波面数据(ΔW(s),ΔW(-s))进行第一步校正:
得到如图6与图7所示的初步校正后点衍射剪切波面数据
(1-6)将步骤(1-5)处理得到初步校正后的点衍射剪切波面数据进行叠加,实现点衍射剪切波面数据的第二步校正,得到不含结构误差的真实点衍射剪切波面ΔW,如图8所示,
亚波长孔径光纤点衍射波前重构方法为:
(2-1)利用步骤(1-3)所重构得到的两亚微米孔径光纤点衍射源的三维坐标(0.125,0,300)和(-0.125,0,300),计算得到两亚微米孔径光纤点衍射源的横向偏移量s=250μm,并据此得到相应的差分泽尼克多项式ΔZ={ΔZi},
ΔZi(x,y)=Zi(x,y)-Zi(x+s,y)
其中,下标i表示多项式的项数,{Zi}表示泽尼克多项式,坐标(x,y)为波面上各点坐标数据;
(2-2)利用最小二乘法对所获得的真实点衍射剪切波面ΔW进行差分泽尼克多项式拟合,得到差分泽尼克多项式系数a={ai},
a=(ΔZTΔZ)-1ΔZTΔW,
其中T表示转置运算;
(2-3)利用步骤(2-2)得到的差分泽尼克多项式系数a={ai},即可高精度重构得到实际待测的亚微米孔径光纤点衍射波面数据W0,如图9所示,
通过实验验证,在所提出的方法中的预校正波前剪切数据相对于传统的方法,可以显著减少系统误差,均方根值为0.0077λ。为了最大限度地减少残余的系统误差,实现对点衍射波前的高精度测量,以对称点源偏移补偿为基础,再对预校正数据进行叠加,均方根值可以减少到0.8253×10-4λ。这表明采用本发明的方法进行测量,消除了结构误差,实现了点衍射波前的高精度测量。
Claims (1)
1.一种亚微米孔径光纤点衍射波前测量的结构误差校正方法,其特征在于包括如下步骤:
(1-1)将两根具有相同结构尺寸的亚微米孔径光纤安装于亚波长孔径光纤点衍射干涉仪的测量探头中,所述的亚微米孔径光纤的前端头呈锥形,构成锥形探针出光端面,两个锥形探针出光端面并排且共面,水平放置所述的测量探头,并且使两个所述的锥形探针出光端面的连线呈水平方向;
(1-2)将所述的亚波长孔径光纤点衍射干涉仪的参考路和检测路的两路相干光分别导入两根亚微米孔径光纤,采用不带成像镜头的CCD探测器采集来自于两根亚微米孔径光纤的点衍射干涉图,利用多步移相算法解调出对应的点衍射剪切波面数据ΔW(s),其中s为两亚微米孔径光纤点衍射源的横向偏移量;
(1-3)在点衍射剪切波面数据ΔW(s)中选取6个以上像素点处的相位值,组成非线性超定方程组,运用迭代优化算法重构出两亚微米孔径光纤点衍射源的三维坐标为(x1,y1,z1)和(x2,y2,z2),已知所述的CCD探测器上每像素点坐标为(x,y,z),根据光程差公式:
生成理想波面数据ΔW(s)′;
(1-4)将所述的测量探头绕其轴线旋转180度,使两亚微米孔径光纤点衍射源的横向偏移量变为-s,保持CCD探测器位置不变,重复步骤(1-2),得到另一组点衍射剪切波面数据ΔW(-s),再用点衍射剪切波面数据ΔW(-s)重复步骤(1-3),得到对应的另一组理想波面数据ΔW(-s)′;
(1-5)利用基于三坐标重构得到的两组理想波面数据(ΔW(s)′,ΔW(-s)′),对测量得到的原始点衍射剪切波面数据(ΔW(s),ΔW(-s))进行第一步校正:
ΔW1 (s)=ΔW(s)-ΔW(s)′
ΔW1 (-s)=ΔW(-s)-ΔW(-s)′
得到初步校正后的点衍射剪切波面数据(ΔW1 (s),ΔW1 (-s));
(1-6)将初步校正后的点衍射剪切波面数据(ΔW1 (s),ΔW1 (-s))进行叠加,实现点衍射剪切波面数据的第二步校正,得到不存在结构误差的真实点衍射剪切波面ΔW:
包括亚波长孔径光纤点衍射波前重构方法:
(2-1)利用步骤(1-3)所重构得到的两亚微米孔径光纤点衍射源的三维坐标(x1,y1,z1)和(x2,y2,z2),计算得到两亚微米孔径光纤点衍射源的横向偏移量s,
并据此得到相应的差分泽尼克多项式ΔZ={ΔZi},
ΔZi(x,y)=Zi(x,y)-Zi(x+s,y)
其中,下标i表示多项式的项数,{Zi}表示泽尼克多项式,坐标(x,y)为波面上各点坐标数据;
(2-2)利用最小二乘法对步骤(1-6)所获得的真实点衍射剪切波面ΔW进行差分泽尼克多项式拟合,得到差分泽尼克多项式系数a={ai},
a=(ΔZTΔZ)-1ΔZTΔW,
其中T表示转置运算;
(2-3)利用差分泽尼克多项式系数a={ai},即可重构得到实际待测的亚微米孔径光纤点衍射波面数据W0,
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610008014.7A CN105527027B (zh) | 2016-01-05 | 2016-01-05 | 亚微米孔径光纤点衍射波前测量的结构误差校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610008014.7A CN105527027B (zh) | 2016-01-05 | 2016-01-05 | 亚微米孔径光纤点衍射波前测量的结构误差校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105527027A CN105527027A (zh) | 2016-04-27 |
CN105527027B true CN105527027B (zh) | 2018-11-02 |
Family
ID=55769396
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610008014.7A Expired - Fee Related CN105527027B (zh) | 2016-01-05 | 2016-01-05 | 亚微米孔径光纤点衍射波前测量的结构误差校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105527027B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105933060B (zh) * | 2016-06-24 | 2018-02-13 | 温州大学 | 一种基于动力学反馈神经网络的波前重构方法 |
CN109341523B (zh) * | 2018-10-08 | 2020-08-14 | 中国计量大学 | 一种点衍射干涉三坐标测量装置探头端面调平方法 |
CN109528161B (zh) * | 2018-11-06 | 2021-01-12 | 浙江大学 | 一种基于机器学习的高速高分辨扫描显微成像系统与方法 |
CN113218519B (zh) * | 2021-06-04 | 2023-03-28 | 中国工程物理研究院应用电子学研究所 | 一种基于双层亚波长孔结构的径向剪切波前测量系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103217104A (zh) * | 2013-03-15 | 2013-07-24 | 中国计量学院 | 一种点衍射三维绝对位移无导轨测量装置及方法 |
CN104330039A (zh) * | 2014-09-28 | 2015-02-04 | 中国计量学院 | 一种用于三坐标测量的大数值孔径光纤点衍射干涉装置及方法 |
-
2016
- 2016-01-05 CN CN201610008014.7A patent/CN105527027B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103217104A (zh) * | 2013-03-15 | 2013-07-24 | 中国计量学院 | 一种点衍射三维绝对位移无导轨测量装置及方法 |
CN104330039A (zh) * | 2014-09-28 | 2015-02-04 | 中国计量学院 | 一种用于三坐标测量的大数值孔径光纤点衍射干涉装置及方法 |
Non-Patent Citations (6)
Title |
---|
一种新的相移点衍射干涉仪系统误差标定方法;刘克 等;《光学学报》;20101031;第30卷(第10期);全文 * |
亚波长孔径光纤点衍射波前质量分析;陈茜茜 等;《光学学报》;20150930;第35卷(第9期);全文 * |
基于Levenbery-Marquardt算法的点衍射三维坐标测量方法;王道档 等;《光学学报》;20140831;第34卷(第8期);全文 * |
点衍射干涉仪基准波前质量测评研究;邵晶;《中国博士学位论文全文数据库 基础科学辑》;20131031;第5.2节 * |
点衍射干涉波前参考源误差及公差分析;代晓珂 等;《中国光学》;20141031;第7卷(第5期);第3节,图4-5 * |
点衍射球面检测中的斜反射波前像差校正;王道档 等;《光学学报》;20110630;第31卷(第6期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN105527027A (zh) | 2016-04-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106052583B (zh) | 基于可变补偿透镜的非球面面形干涉测量方法及装置 | |
CN105527027B (zh) | 亚微米孔径光纤点衍射波前测量的结构误差校正方法 | |
CN101876540B (zh) | 基于多波前透镜补偿器的非球面绝对测量系统 | |
CN106840027A (zh) | 光学自由曲面的像散补偿型干涉检测装置与检测方法 | |
CN101865670B (zh) | 光纤点衍射移相干涉仪的平面面形测量方法 | |
CN107421436B (zh) | 基于空间光调制器参考面的非球面干涉测量系统及方法 | |
CN101949691A (zh) | 非零位补偿浅度光学非球面面形检测方法 | |
CN106595529B (zh) | 基于虚拟牛顿环的大曲率半径非零位干涉测量方法及装置 | |
CN106517086A (zh) | 一种大面积高分辨率宽视场在线测量装置及其测量方法 | |
US11635289B2 (en) | Surface shape measurement device and surface shape measurement method | |
CN108061639A (zh) | 一种结合自适应光学技术的大动态范围、高精度相位差法波前测量仪 | |
CN104655053B (zh) | 基于针孔式点衍射干涉仪球面镜曲率半径测量方法 | |
CN111751012A (zh) | 动态高分辨光学波前相位测量装置及测量方法 | |
WO2018000943A1 (zh) | 一种凹柱面及柱面发散镜的检测方法及装置 | |
CN105784129A (zh) | 一种用于激光波前检测的低频外差干涉仪 | |
Wu et al. | Precise 3-D microscopic profilometry using diffractive image microscopy and artificial neural network in single-exposure manner | |
CN105823563A (zh) | 动态高分辨率波前测量装置及方法 | |
CN111397634B (zh) | 星敏感器固定端面热变形的高分辨干涉检测装置及方法 | |
JP2005524065A (ja) | 光学検査方法および装置 | |
CN102589472A (zh) | 在球面面形干涉检测中高精度消除调整误差的方法 | |
CN102494631B (zh) | 一种绕光轴旋转对准误差分析方法 | |
CN102927930B (zh) | 采用平行光管拼接检测超大口径反射镜面形误差的方法 | |
KR20110065365A (ko) | 비구면체 측정 방법 및 장치 | |
CN112611548B (zh) | 一种基于数字全息的透镜焦距测量装置及方法 | |
CN205642308U (zh) | 基于双波长干涉的微表面形貌测量装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB02 | Change of applicant information | ||
CB02 | Change of applicant information |
Address after: The town of Anping Road Wenling city Taizhou city Zhejiang province 317523 No. 120 Applicant after: CHINA JILIANG UNIVERSITY Address before: Hangzhou City, Zhejiang province 310018 Xiasha Higher Education Park source Street No. 258 Applicant before: China Jiliang University |
|
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20181102 |