发明名称 基于横向剪切干涉结构的光学元件面形的测量方法
摘要 一种基于横向剪切干涉结构的光学元件面形的测量方法,测量系统包括光源、显微物镜、小孔光阑、准直物镜、剪切平板、光电探测器和计算机,测量时将待测光学元件置于准直物镜之后,并使光轴穿过待侧面的中心,测量步骤如下:①获取x方向和y方向的横向剪切干涉图;②在计算机中获取x方向和y方向的差分波前;③评价待测波前线性度,确定线性度最高的子矩形区域的初始点;④通过线性插值法计算重建的初始值;⑤通过区域法重建待测波前。本发明不增加测量次数,实现了高空间分辨率的光学元件面形的测量,并且提高了测量的精度和稳定性。
申请公布号 CN102878947B 申请公布日期 2014.12.17
申请号 CN201210356681.6 申请日期 2012.09.21
申请人 中国科学院上海光学精密机械研究所 发明人 戴凤钊;唐锋;王向朝;张敏
分类号 G01B11/24(2006.01)I 主分类号 G01B11/24(2006.01)I
代理机构 上海新天专利代理有限公司 31213 代理人 张泽纯
主权项 一种基于横向剪切干涉结构的光学元件面形测量方法,本方法采用的测量系统包括光源(1),沿该光源输出的光束方向依次为显微物镜(2)、小孔光阑(3)、准直物镜(4)、剪切平板(5)和准直物镜(6),光电探测器(8)置于剪切平板(5)下方,该光电探测器(8)的输出端与计算机(9)的输入端相连,测量时使剪切平板与光轴成45°角,将待测光学元件置于所述的准直物镜(6)之后,使光轴穿过待测面的中心,其特征在于该方法测量步骤如下:①获取x方向和y方向的横向剪切干涉图:将待测光学元件置于所述的准直物镜(6)之后,使光轴穿过待测面的中心,启动测量系统,在光轴方向上移动待测光学元件,使光束充满整个待测面,通过所述的光电探测器(8)采集x方向横向剪切干涉图,并输入至计算机(9)保存,然后将待测光学元件(7)绕光轴顺时针旋转90°,通过光电探测器采集y方向的横向剪切干涉图,并输入至计算机(9)保存;②在计算机中获取x方向和y方向的差分波前:在计算机中将所述的x方向的横向剪切干涉图通过傅里叶变换相位提取方法获取x方向的差分波前(11),设置一个N×(N‑S<sub>x</sub>)的全零矩阵,用x方向的差分波前的值替代该全零矩阵相应位置处的零值,获得x方向差分波前矩阵ΔW<sub>x</sub>(13),将该x方向的差分波前矩阵(13)的矩阵元逐行从左至右依次顺序编码,得到由N(N‑S<sub>x</sub>)个差分波前值构成的一维差分波前列向量ΔW<sub>x</sub>;在计算机中将所述的y方向的横向剪切干涉图通过傅里叶变换相位提取方法获取y方向的差分波前(12),设置一个(N‑S<sub>y</sub>)×N的全零矩阵,用y方向的差分波前的值替代该全零矩阵相应位置处的零值,获得y方向差分波前矩阵ΔW<sub>y</sub>(14),将所述的y方向的差分波前矩阵(14)的矩阵元逐行从左至右依次顺序编码,得到由N×(N‑S<sub>y</sub>)个差分波前值构成的一维差分波前列向量ΔW<sub>y</sub>;③评价待测波前线性度,确定线性度最高的子矩形区域的初始点:取x方向差分波前矩阵ΔW<sub>x</sub>的1到N‑S<sub>y</sub>行,1到N‑S<sub>x</sub>列的所有元素构成一个维度为(N‑S<sub>y</sub>)×(N‑S<sub>x</sub>)的差分波前子矩阵<img file="FDA0000571651750000011.GIF" wi="121" he="65" />即:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>&Delta;</mi><msubsup><mi>W</mi><mi>x</mi><mi>s</mi></msubsup><mo>[</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>]</mo><mo>=</mo><mi>&Delta;</mi><msub><mi>W</mi><mi>x</mi></msub><mo>[</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>]</mo></mrow>]]></math><img file="FDA0000571651750000012.GIF" wi="413" he="65" /></maths>其中:1≤m≤N‑S<sub>y</sub>,1≤n≤N‑S<sub>x</sub>,取y方向差分波前矩阵ΔW<sub>y</sub>的1到N‑S<sub>y</sub>行,1到N‑S<sub>x</sub>列的所有元素构成一个维度为(N‑S<sub>y</sub>)×(N‑S<sub>x</sub>)差分波前子矩阵<img file="FDA0000571651750000013.GIF" wi="119" he="71" />即:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>&Delta;</mi><msubsup><mi>W</mi><mi>y</mi><mi>s</mi></msubsup><mo>[</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>]</mo><mo>=</mo><mi>&Delta;</mi><msub><mi>W</mi><mi>y</mi></msub><mo>[</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>]</mo></mrow>]]></math><img file="FDA0000571651750000014.GIF" wi="415" he="71" /></maths>其中:1≤m≤N‑S<sub>y</sub>,1≤n≤N‑S<sub>x</sub>,由所述的差分波前子矩阵<img file="FDA0000571651750000021.GIF" wi="96" he="65" />和<img file="FDA0000571651750000022.GIF" wi="90" he="71" />计算平方和差分波前矩阵ΔW<sub>xy</sub>,ΔW<sub>xy</sub>的值由下式计算:ΔW<sub>xy</sub>[m,n]=(ΔW<sub>x</sub>[m,n])<sup>2</sup>+(ΔW<sub>y</sub>[m,n])<sup>2</sup>式中:1≤m≤N‑S<sub>y</sub>,1≤n≤N‑S<sub>x</sub>,搜索所述的平方和差分波前ΔW<sub>xy</sub>,找出除零值外最小值所在位置的坐标,标记该坐标为(r,c),该坐标表示的点即为本发明所确定的待测波前线性度最高的子矩形区域的初始点;④通过线性插值法计算重建的初始值:在计算机中设置一个S<sub>y</sub>×S<sub>x</sub>的初始值矩阵<img file="FDA0000571651750000023.GIF" wi="80" he="73" />该矩阵的矩阵元的值由差分波前矩阵ΔW<sub>x</sub>和ΔW<sub>y</sub>通过线性插值法计算,即:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mover><mi>W</mi><mo>~</mo></mover><mi>s</mi></msub><mo>[</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>]</mo><mo>=</mo><mfrac><mrow><mi>&Delta;</mi><msub><mi>W</mi><mi>y</mi></msub><mo>[</mo><mi>r</mi><mo>,</mo><mi>c</mi><mo>]</mo></mrow><msub><mi>S</mi><mi>y</mi></msub></mfrac><mo>&times;</mo><mrow><mo>(</mo><mi>i</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><mfrac><mrow><mi>&Delta;</mi><msub><mi>W</mi><mi>x</mi></msub><mo>[</mo><mi>r</mi><mo>,</mo><mi>c</mi><mo>]</mo></mrow><msub><mi>S</mi><mi>x</mi></msub></mfrac><mo>&times;</mo><mrow><mo>(</mo><mi>j</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000571651750000024.GIF" wi="795" he="130" /></maths>式中:i=1,2,...,S<sub>y</sub>,j=1,2,...,S<sub>x</sub>‑1,将所述的初始值矩阵<img file="FDA0000571651750000025.GIF" wi="58" he="71" />的矩阵元逐行从左至右依次编码,将初始值矩阵<img file="FDA0000571651750000026.GIF" wi="56" he="71" />表示为初始值列向量<img file="FDA0000571651750000027.GIF" wi="88" he="72" />⑤通过区域法重建待测波前:由正交方向上的两个差分波前列向量ΔW<sub>x</sub>和ΔW<sub>y</sub>与初始值列向量<img file="FDA0000571651750000028.GIF" wi="60" he="66" />计算待测波前在所有N<sup>2</sup>点上的值<img file="FDA0000571651750000029.GIF" wi="82" he="64" />计算方法如下:<img file="FDA00005716517500000210.GIF" wi="358" he="220" />式中:M<sup>x</sup>为N(N‑S<sub>x</sub>)×N<sup>2</sup>矩阵,其由下式给出:<img file="FDA00005716517500000211.GIF" wi="1260" he="263" />M<sup>y</sup>为N(N‑S<sub>y</sub>)×N<sup>2</sup>矩阵,其由下式给出:<img file="FDA00005716517500000212.GIF" wi="1010" he="236" />H为S<sub>x</sub>S<sub>y</sub>×N<sup>2</sup>矩阵,由下式给出:<img file="FDA00005716517500000213.GIF" wi="735" he="133" />式中:i=1,2,...,S<sub>x</sub>S<sub>y</sub>,k=1,2,...,S<sub>y</sub>,l=1,2,...,S<sub>x</sub>,将所述的N<sup>2</sup>×1的待测波前列向量<img file="FDA0000571651750000031.GIF" wi="52" he="64" />从第一个值按顺序分为N组,每组N个值表示为一个行向量,将所有N个行向量按顺序组成N×N待测波前<img file="FDA0000571651750000032.GIF" wi="72" he="72" /><img file="FDA0000571651750000033.GIF" wi="49" he="72" />即表示待测光学元件的面形。
地址 201800 上海市嘉定区800-211邮政信箱