CN114322829A - 一种基于双波长涡旋光束相移干涉的拼接主镜共相误差检测系统及其工作方法 - Google Patents

一种基于双波长涡旋光束相移干涉的拼接主镜共相误差检测系统及其工作方法 Download PDF

Info

Publication number
CN114322829A
CN114322829A CN202111545618.2A CN202111545618A CN114322829A CN 114322829 A CN114322829 A CN 114322829A CN 202111545618 A CN202111545618 A CN 202111545618A CN 114322829 A CN114322829 A CN 114322829A
Authority
CN
China
Prior art keywords
mirror
splicing
error
light
sub
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.)
Granted
Application number
CN202111545618.2A
Other languages
English (en)
Other versions
CN114322829B (zh
Inventor
杨忠明
杨丽丽
刘兆军
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shandong University
Original Assignee
Shandong University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shandong University filed Critical Shandong University
Priority to CN202111545618.2A priority Critical patent/CN114322829B/zh
Publication of CN114322829A publication Critical patent/CN114322829A/zh
Application granted granted Critical
Publication of CN114322829B publication Critical patent/CN114322829B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Instruments For Measurement Of Length By Optical Means (AREA)

Abstract

本发明涉及一种基于双波长涡旋光束相移干涉的拼接主镜共相误差检测系统及其工作方法,包括:第一激光器、第一扩束器、第二激光器、第二扩束器、第一分束器、第二分束器、四分之一波片、零级涡旋波片、第三分束器、平面反射镜、PZT位移平台、聚焦透镜、第四分束器、调制透镜、第五分束器、CCD相机、拼接镜、PC终端;本发明将涡旋光束相移干涉与双波长算法结合用于拼接镜的共相误差检测中,本发明避免了一些加工难度大的器件的使用,并且测量精度高、测量范围大、检测效率高、实时性好、原理简单、不受2π模糊影响等,为拼接镜共相误差的测量提供了一种新颖有效的思路。

Description

一种基于双波长涡旋光束相移干涉的拼接主镜共相误差检测 系统及其工作方法
技术领域
本发明涉及光学共相检测领域,特别是涉及一种基于双波长涡旋光束相移干涉的拼接主镜共相 误差检测系统及其工作方法。
背景技术
随着科学技术的不断发展,大口径天文望远镜在重大前沿科学问题研究、宇宙天体研究等领域 展现出越来越重要的作用。增大望远镜的主镜口径是提高观测分辨率和集光能力的有效手段,然而 由于加工运输、镜面检测、制造成本以及后期维护等因素的制约,研究人员普遍认为单体式望远镜 的口径难以超过10米。
为了解决大型望远镜口径技术难点的制约,研究人员提出了拼接镜技术。然而,拼接镜技术带 来了在装调、检测等方面的新挑战,只有实现所有子镜的光学上共相位,才能使拼接式望远镜和单 体式望远镜具有相同的成像质量。共相误差包括平移误差和倾斜误差的探测与校正是实现各个子镜 共相、保障拼接式望远镜成像质量的关键技术。
现有的拼接镜共相检测技术有窄带/宽带哈特曼传感技术、色散条纹传感技术、四棱锥传感法、 相位差法等多种探测方法,但在探测范围、测量精度或计算速度等方面存在各种问题,如窄带/宽 带哈特曼传感技术需要加入精密的光瞳掩膜,且存在掩模与子镜光瞳边界对准难等问题;色散条纹 法探测范围广,但测量精度低;四棱锥传感法中的四棱锥存在顶点对准难、加工难度大等问题;相 位差法精度高,但运算量大,存在难以实现实时检测等问题。
发明内容
针对现有技术的不足,本发明提供一种基于双波长涡旋光束相移干涉的拼接主镜共相误差检测 系统及其工作方法,实现测量精度高、范围大的拼接镜共相误差实时测量。
本发明首次将涡旋光束用于拼接镜的共相误差检测,将涡旋光束干涉测量与相移法及双波长算 法结合测量拼接镜的共相误差,具有精度高、实时性好、原理简单、测量范围大、不受2π模糊影 响等优点。
术语解释:
CCD,是电荷耦合器件(charge coupled device)的简称,也可以称为CCD图像传感器。它使 用一种高感光度的半导体材料制成,能把光线转变成电荷,然后通过模数转换器芯片将电信号转换 成数字信号,数字信号经过压缩处理经USB接口传到电脑上就形成所采集的图像。即CCD相机用 于图像采集,捕获干涉条纹图像。
PZT位移平台,PZT是锆钛酸铅压电陶瓷的简称,是一种高精度微位移器件,利用压电陶瓷 逆压电效应,能实现亚微米级的微位移。当在电介质的极化方向上施加电场,这些电介质也会发生 变形,电场去掉后,电介质的变形随之消失,这种现象称为逆压电效应。
本发明的技术方案为:
一种基于双波长涡旋光束相移干涉的拼接主镜共相误差检测系统,包括:第一激光器、第一扩 束器、第二激光器、第二扩束器、第一分束器、第二分束器、四分之一波片、零级涡旋波片、第三 分束器、平面反射镜、PZT位移平台、聚焦透镜、第四分束器、调制透镜、第五分束器、CCD相 机、拼接镜、PC终端;
由第一激光器、第二激光器发出的不同波长的单色偏振光分别经过第一扩束器、第二扩束器扩 束成平行光,第一激光器、第二激光器不同时使用,第一分束器将经过第一扩束器、第二扩束器扩 束成的平行光垂直入射至第二分束器,分成透射光T1与反射光F1;
反射光F1经过四分之一波片变成圆偏振光,再经过零级涡旋波片转化为涡旋光束,涡旋光束 入射到第三分束器上,并由第三分束器反射至平面反射镜上后经过第三分束器原路返回作为参考光 入射至第五分束器上;
透射光T1经聚焦透镜聚焦后入射至第四分束器,并透射到拼接镜上,拼接镜原路返回的携带 拼接镜共相误差信息的球面波作为测试光由第四分束器反射,后再经由调制透镜调制后入射到第五 分束器上;
参考光即涡旋光束与测试光即球面光束在第五分束器处合束并发生干涉,CCD相机捕获干涉 条纹;PC终端先将干涉条纹处理得到包裹相位图,然后提取包裹相位图的梯度线拟合为费马螺线, 通过费马螺线的旋转角与位置特征反解出各个拼接子镜的平移误差与倾斜误差,再由PC终端将共 相误差反馈给各拼接子镜的促动器实现共相误差的校正,拼接镜包括多个拼接子镜。
根据本发明优选的,平面反射镜设置于PZT位移平台上用于移相。
根据本发明优选的,第一分束器、第二分束器、第三分束器、第四分束器、第五分束器的反射 面分别与所在光路的方向呈顺时针或逆时针45°夹角。
根据本发明优选的,所述聚焦透镜的焦点与拼接镜的曲率中心重合。
上述基于双波长涡旋光束相移干涉的拼接主镜共相误差检测系统的工作方法,包括步骤如下:
(1)第一激光器、第二激光器发出的不同波长的单色偏振光分别经过第一扩束器、第二扩束 器扩束成平行光,第一分束器将经过第一扩束器、第二扩束器扩束成的平行光垂直入射至第二分束 器,分成透射光T1与反射光F1;
反射光F1经过四分之一波片变成圆偏振光,再经过零级涡旋波片转化为涡旋光束,涡旋光束 入射到第三分束器上,并由第三分束器反射至平面反射镜上后经过第三分束器原路返回作为参考光 入射至第五分束器上;
透射光T1经聚焦透镜聚焦后入射至第四分束器,并透射到拼接镜上,拼接镜原路返回的携带 拼接镜共相误差信息的球面波作为测试光由第四分束器反射,再经由调制透镜调制后入射到第五分 束器上;
(2)参考光即涡旋光束与测试光即球面光束在第五分束器处合束并发生干涉,CCD相机捕获 干涉条纹;
(3)PC终端将干涉条纹处理得到包裹相位图;
(4)提取包裹相位图的梯度线拟合为费马螺线,通过费马螺线的旋转角与位置特征反解出各 个拼接子镜的平移误差与倾斜误差;
(5)PC终端将共相误差反馈给各拼接子镜的促动器实现共相误差的校正。
根据本发明优选的,涡旋光束作为参考光束,其光场表达式如式(Ⅰ)所示:
Figure BDA0003415730910000031
测试光即球面光束,其光场表达式如式(Ⅱ)所示:
Figure BDA0003415730910000032
其中,Av、As分别是涡旋光与球面波的振幅,i是指表示第i个拼接子镜,m是拓扑电荷数,θ 是涡旋光束的方位角,k是波数,z是传输距离,r是极坐标系下的径向分量,
Figure BDA0003415730910000033
是由拼接镜以及调 制透镜等引入的球面波前相位,常数
Figure BDA0003415730910000034
Figure BDA0003415730910000035
分别为对应于涡旋光与球面波的初相。
根据本发明优选的,当拼接镜存在共相误差,且取涡旋光的拓扑荷数为1时,干涉光强表达 式如式(Ⅲ)所示:
Figure BDA0003415730910000036
式(Ⅲ)中,
Figure BDA0003415730910000037
k是波数,x,y表示干涉图像中每个像素点对应的 坐标值,i表示第i个拼接子镜,n为拼接子镜的个数,pi是拼接子镜的平移误差,η*i是x轴倾斜误 差,ηyi是y轴倾斜误差,2(pi+xsinη*i+ysinηyi)是由拼接镜共相误差而带来的光路差。
根据本发明优选的,步骤(3)中,PC终端通过三步相移法将干涉条纹处理得到包裹相位图, 相移量分别为-2π/3,0,2π/3,三步相移法的表示式如式(Ⅳ)所示:
Figure BDA0003415730910000041
式(Ⅳ)中,
Figure BDA0003415730910000042
I1、I2、I3分别对应于三步相移过程中的干涉光强;
通过式(Ⅲ)与式(Ⅳ),得到相位表达式如式(Ⅴ)所示:
Figure BDA0003415730910000043
式(Ⅴ)中,
Figure BDA0003415730910000044
k是波数,x,y表示包裹相位图中每个像素点对应的坐标值,i 表示第i个拼接子镜,n为拼接子镜的个数,pi是拼接子镜的平移误差,η*i是x轴倾斜误差,ηyi是 y轴倾斜误差,2(pi+xsinη*i+ysinηyi)是由拼接镜共相误差而带来的光路差,I1、I2、I3分别对 应于三步相移过程中的干涉光强。
根据本发明优选的,步骤(4)中,相位被包裹在-π到π区间,包裹相位图的相位分布有一个 断层式的变化,即产生了一个梯度或边缘,这个梯度或边缘满足式(Ⅵ):
Figure BDA0003415730910000045
式(Ⅵ)中,
Figure BDA0003415730910000046
k是波数,x,y表示包裹相位图中每个像素点对应的坐标值,i 表示第i个拼接子镜,n为拼接子镜的个数,pi是拼接子镜的平移误差,η*i是x轴倾斜误差,ηyi是 y轴倾斜误差,2(pi+xsinη*i+ysinηyi)是由拼接镜共相误差而带来的光路差,j为整数。
当拼接镜中仅存在平移误差时,边缘提取的梯度曲线符合标准的费马螺旋线,其旋转角度θi随 平移误差pi的变化而变化,拼接镜拟合的费马螺线的旋转角度与拼接镜的平移误差的关系如(Ⅶ) 所示:
pi=Δθi/(2k) (Ⅶ)
式(Ⅶ)中,pi是第i个拼接子镜的平移误差,Δθi是由第i个拼接子镜与参考子镜拟合的费马 螺线的旋转角的差值,k是波数。
费马螺线的旋转角与位置特征通过曲线拟合得到;
合成波长的表达式如式(Ⅷ)所示:
Figure BDA0003415730910000047
式(Ⅷ)中,λ1、λ2分别是指第一激光器、第二激光器的波长;
假如拼接子镜的实际平移误差为p,利用式(Ⅶ)得到的两波长的平移误差分别为p1和p2, 则平移误差p的求取公式如式(Ⅸ)所示:
Figure BDA0003415730910000051
式(Ⅸ)中,u、v是待求解的常数,对两波长上求得的平移误差p进行比较,当两波长上p 最匹配时,此时对应的u与v的值就为求解值;
当拼接镜中仅存在倾斜误差时,将式(Ⅵ)整理,得到边缘提取的梯度曲线即变形的费马螺线 如式(Ⅹ)所示:
Figure BDA0003415730910000052
式(Ⅹ)中,
Figure BDA0003415730910000053
k是波数,x,y表示包裹相位图中每个像素点对应的坐标值,i 表示第i个拼接子镜,n为拼接子镜的个数,η*i是x轴倾斜误差,ηyi是y轴倾斜误差,j为整数;
将式(Ⅹ)称为变形的费马螺线,并且设定(2R sinη*i,2R sinηyi)是变形的费马螺线的中 心,通过曲线拟合得到的变形的费马螺线中心位置从原点移动了x0和y0的距离,则倾斜误差 如式(Ⅺ)所示:
Figure BDA0003415730910000054
式(Ⅺ)中,η*是x轴倾斜误差,ηy是y轴倾斜误差,x0和y0是通过曲线拟合得到的变形的 费马螺线的中心位置。
根据本发明优选的,步骤(5)之前先对倾斜误差进行校正,具体是指:
选择一个校正过倾斜误差的拼接子镜为参考子镜,即该拼接子镜拟合的费马曲线为标准的费马 螺线,标准的费马螺线的中心位置为原点;
将该拼接子镜拟合的费马螺线的旋转角与中心位置特征作为参考值,将其它拼接子镜拟合的费 马螺线的旋转角、中心位置特征与参考值作比较,先通过中心位置特征反解出对应拼接子镜的倾斜 误差,并将其校正,倾斜误差校正后再通过曲线拟合得到费马螺线的旋转角,并反解出对应拼接子 镜的平移误差。将其余拼接子镜拟合得到的变形的费马螺线的中心位置特征代入式(Ⅺ),即将x0和y0代入式(Ⅺ)得到对应拼接子镜的倾斜误差;在倾斜误差校正后,将其余拼接子镜与参考子 镜拟合得到的费马螺线的旋转角特征相减后代入式(Ⅶ),即将Δθi代入(Ⅶ)得到对应拼接子镜的 平移误差pi
将使用双波长测量计算的共相误差反馈给促动器来实时校正各拼接子镜的共相误差,具体是指: 将计算的共相误差作为校正值,引导促动器对各拼接子镜的位姿进行调整,促动器控制各拼接子镜 的左右旋转、前后俯仰与沿着光轴的前后移动。
本发明的有益效果为:
1、本发明想法新颖,首次将涡旋光束相移干涉用于拼接镜的共相误差检测中,相比传统干涉 图像的竖直条纹或牛顿环,涡旋光束干涉图像中的螺旋条纹的信息量更多,干涉图中费马螺线一直 存在,没有传统干涉方法如牛顿环中圆心出现暗条纹可能存在信息难以确定的情况。
2、本发明测量精度高、测量范围大、检测效率高、可实时得到子镜之间的倾斜误差和平移误 差,避免了一些加工难度大的器件的使用,方法简便,原理简单。
附图说明
图1为本发明基于双波长涡旋光束相移干涉的拼接主镜共相误差检测系统的结构示意图;
图2为本发明基于双波长涡旋光束相移干涉的拼接主镜共相误差检测系统的工作方法的流程 示意图;
图3(a)为右侧拼接子镜不存在共相误差下的曲线拟合示例示意图;
图3(b)为右侧拼接子镜仅存在平移误差下的曲线拟合示例示意图;
图3(c)为右侧拼接子镜仅存在x轴倾斜误差下的曲线拟合示例示意图;
图3(d)为右侧拼接子镜仅存在y轴倾斜误差下的曲线拟合示例示意图。
1、第一激光器,2、第一扩束器,3、第二激光器,4、第二扩束器,5、第一分束器,6、第二分束器,7、四分之一波片,8、零级涡旋波片,9、第三分束器,10、平面反射镜,11、PZT位 移平台,12、聚焦透镜,13、第四分束器,14、调制透镜,15、第五分束器,16、CCD相机,17、 拼接镜,18、PC终端。
具体实施方式
下面结合说明书附图和实施例对本发明作进一步限定,但不限于此。
实施例1
一种基于双波长涡旋光束相移干涉的拼接主镜共相误差检测系统,如图1所示,包括:第一激 光器1、第一扩束器2、第二激光器3、第二扩束器4、第一分束器5、第二分束器6、四分之一波 片7、零级涡旋波片8、第三分束器9、平面反射镜10、PZT位移平台11、聚焦透镜12、第四分束 器13、调制透镜14、第五分束器15、CCD相机16、拼接镜17、PC终端18;
由第一激光器1、第二激光器3发出的不同波长的单色偏振光分别经过第一扩束器2、第二扩 束器4扩束成平行光,第一激光器1、第二激光器3不同时使用,第一分束器5将经过第一扩束器 2、第二扩束器4扩束成的平行光垂直入射至第二分束器6,分成透射光T1与反射光F1;
反射光F1经过四分之一波片7变成圆偏振光,再经过零级涡旋波片8转化为涡旋光束,涡旋 光束入射到第三分束器9上,并由第三分束器9反射至平面反射镜10上后经过第三分束器9原路 返回作为参考光入射至第五分束器15上;
透射光T1经聚焦透镜12聚焦后入射至第四分束器13,并透射到拼接镜17上,拼接镜17原 路返回的携带拼接镜17共相误差信息的球面波作为测试光由第四分束器13反射,后再经由调制透 镜14调制后入射到第五分束器15上;在干涉测量领域,一般需要干涉条纹不要太密集,以方便数 据处理。调制透镜14用于调制条纹密度。
参考光即涡旋光束与测试光即球面光束在第五分束器15处合束并发生干涉,CCD相机16捕 获干涉条纹;PC终端18先将干涉条纹处理得到包裹相位图,然后提取包裹相位图的梯度线拟合为 费马螺线,通过费马螺线的旋转角与位置特征反解出各个拼接子镜的平移误差与倾斜误差,再由 PC终端18将共相误差反馈给各拼接子镜的促动器实现共相误差的校正,拼接镜17包括多个拼接 子镜。拼接子镜由促动器控制调节平移误差与倾斜误差。
参考光即涡旋光束与测试光即球面光束在第五分束器15处合束并发生干涉,干涉光束形成干 涉条纹被CCD相机16捕获。干涉条纹呈费马螺旋的形式分布。当拼接子镜存在平移或倾斜误差 时,拼接子镜拟合的费马螺线的旋转角、中心位置特征不相同,通过检测螺旋条纹的旋转角度和位 置特征可以反馈拼接镜17的共相误差。
平面反射镜10设置于PZT位移平台11上用于移相。
第一分束器5、第二分束器6、第三分束器9、第四分束器13、第五分束器15的反射面分别与 所在光路的方向呈顺时针或逆时针45°夹角。
聚焦透镜12的焦点与拼接镜17的曲率中心重合。
实施例2
实施例1所述的双波长涡旋光束相移干涉的拼接主镜共相误差检测系统的工作方法,如图2 所示,包括步骤如下:
(1)第一激光器1、第二激光器3发出的不同波长的单色偏振光分别经过第一扩束器2、第二 扩束器4扩束成平行光,第一分束器5将经过第一扩束器2、第二扩束器4扩束成的平行光垂直入 射至第二分束器6,分成透射光T1与反射光F1;
反射光F1经过四分之一波片变7成圆偏振光,再经过零级涡旋波片8转化为涡旋光束,涡旋 光束入射到第三分束器9上,并由第三分束器9反射至平面反射镜10上后经过第三分束器9原路 返回作为参考光入射至第五分束器15上;
透射光T1经聚焦透镜12聚焦后入射至第四分束器13,并透射到拼接镜17上,拼接镜17原 路返回的携带拼接镜17共相误差信息的球面波作为测试光由第四分束器13反射,再经由调制透镜 14调制后入射到第五分束器15上;
涡旋光束作为参考光束,其光场表达式如式(Ⅰ)所示:
Figure BDA0003415730910000081
测试光即球面光束,其光场表达式如式(Ⅱ)所示:
Figure BDA0003415730910000082
其中,Av、As分别是涡旋光与球面波的振幅,i是指表示第i个拼接子镜,m是拓扑电荷数,θ 是涡旋光束的方位角,k是波数,z是传输距离,r是极坐标系下的径向分量,
Figure BDA0003415730910000083
是由拼接镜17以 及调制透镜14等引入的球面波前相位,常数
Figure BDA0003415730910000084
Figure BDA0003415730910000085
分别为对应于涡旋光与球面波的初相。
(2)参考光即涡旋光束与测试光即球面光束在第五分束器15处合束并发生干涉,CCD相机 16捕获干涉条纹;CCD相机16上捕获的干涉光强表达式为:
Figure BDA0003415730910000086
其中,
Figure BDA0003415730910000087
当拼接镜17存在共相误差,且取涡旋光的拓扑荷数为1时,CCD相机16上捕获的干涉光强 表达式如式(Ⅲ)所示:
Figure BDA0003415730910000088
式(Ⅲ)中,
Figure BDA0003415730910000089
k是波数,x,y表示干涉图像中每个像素点对应的 坐标值,i表示第i个拼接子镜,n为拼接子镜的个数,pi是拼接子镜的平移误差,η*i是x轴倾斜误 差,ηyi是y轴倾斜误差,2(pi+xsinη*i+ysinηyi)是由拼接镜17共相误差而带来的光路差。
(3)PC终端18将干涉条纹处理得到包裹相位图;
由式(Ⅲ)可以看出,干涉条纹呈螺旋状分布,通过检测螺旋条纹的旋转角度和位置特征,可 以反馈出分段镜的共相误差。但是,由于干涉条纹具有一定的宽度,因此难以直接精确提取干涉条 纹的极大值或条纹中心,此外,条纹的极大值及条纹中心的确定受多种因素的影响,如测试光本身 的强度分布、噪声等。因此,难以拟合螺旋条纹的旋转角度与位置特征,即难以满足高精度测量的 要求。
因此,本发明引入了三步相移干涉法来提高相位和空间分辨率,其中相移量被选为2π/3,因 为它对非线性误差的敏感性较低。步骤(3)中,PC终端18通过三步相移法将干涉条纹处理得到 包裹相位图,相移量分别为-2π/3,0,2π/3,因为等间隔相移对非线性误差的敏感性较低。三步相 移法的表示式如式(Ⅳ)所示:
Figure BDA0003415730910000091
式(Ⅳ)中,
Figure BDA0003415730910000092
I1、I2、I3分别对应于三步相移过程中的干涉光强;
相移法将相位限制在[-π,π),因此,在π+2jπ,j为整数处,包裹相位图将出现断层式变化,即存 在梯度,该梯度拟合曲线呈费马螺线的形式分布。相移法具有抗干扰能力强,测量精度较高,对条 纹噪声敏感性较低等优点,因此,本发明选择通过相移来提取螺旋条纹的旋转角度和位置特征。
通过式(Ⅲ)与式(Ⅳ),得到相位表达式如式(Ⅴ)所示:
Figure BDA0003415730910000093
式(Ⅴ)中,
Figure BDA0003415730910000094
k是波数,x,y表示包裹相位图中每个像素点对应的坐标值,i 表示第i个拼接子镜,n为拼接子镜的个数,pi是拼接子镜的平移误差,η*i是x轴倾斜误差,ηyi是 y轴倾斜误差,2(pi+xsinη*i+ysinηyi)是由拼接镜17共相误差而带来的光路差,I1、I2、I3分 别对应于三步相移过程中的干涉光强。
(4)提取包裹相位图的梯度线拟合为费马螺线,通过费马螺线的旋转角与位置特征反解出各 个拼接子镜的平移误差与倾斜误差;
由式(Ⅴ)可以看出,由于求解过程中反正切函数的区间是从-π到π,因此,相位被包裹在-π 到π区间,包裹相位图的相位分布有一个断层式的变化,即产生了一个梯度或边缘,这个梯度或边 缘满足式(Ⅵ):
Figure BDA0003415730910000095
式(Ⅵ)中,
Figure BDA0003415730910000096
k是波数,x,y表示包裹相位图中每个像素点对应的坐标值,i 表示第i个拼接子镜,n为拼接子镜的个数,pi是拼接子镜的平移误差,η*i是x轴倾斜误差,ηyi是 y轴倾斜误差,2(pi+xsinη*i+ysinηyi)是由拼接镜17共相误差而带来的光路差,j为整数。
当拼接镜17中仅存在平移误差时,边缘提取的梯度曲线符合标准的费马螺旋线(r2=a(θ+b),r为极坐标下的径向分量,θ为极角。),其旋转角度θi随平移误差pi的变化而变化,拼接镜17拟合 的费马螺线的旋转角度与拼接镜17的平移误差的关系如(Ⅶ)所示:
pi=Δθi/(2k) (Ⅶ)
式(Ⅶ)中,pi是第i个拼接子镜的平移误差,Δθi是由第i个拼接子镜与参考子镜拟合的费马 螺线的旋转角的差值,k是波数。
费马螺线的旋转角与位置特征通过曲线拟合得到;通过提取的包裹相位图梯度线的数据点的曲 线拟合,得到费马螺线的表达式,表达式中包含了费马螺线的旋转角与位置特征。
干涉条纹和费马螺旋的分布是周期性的。当用单色波长测量时,相差λ/2的平移误差将产生相 同的条纹分布,即存在2π模糊。因此,单色光检测时,平移误差的检测量程仅为λ/2。因此,采 用双波长光源切换测量即双波长算法来解决2π模糊问题,唯一的确定拼接镜17的平移误差,并使 得测量范围扩展到合成波长的一半即几微米。合成波长的表达式如式(Ⅷ)所示:
Figure BDA0003415730910000101
式(Ⅷ)中,λ1、λ2分别是指第一激光器1、第二激光器3的波长;
假如拼接子镜的实际平移误差为p,利用式(Ⅶ)得到的两波长的平移误差分别为p1和p2, 则平移误差p的求取公式如式(Ⅸ)所示:
Figure BDA0003415730910000102
式(Ⅸ)中,u、v是待求解的常数,对两波长上求得的平移误差p进行比较,当两波长上p 最匹配时,此时对应的u与v的值就为求解值;
通过求解u与v就可以唯一的确定实际平移误差p,即解决2π模糊的限制问题,将平移误差 的测量量程由λ1/2或λ2/2扩展到λΛ/2。
当拼接镜17中仅存在倾斜误差时,将式(Ⅵ)整理,得到边缘提取的梯度曲线即变形的费马 螺线如式(Ⅹ)所示:
Figure BDA0003415730910000103
式(Ⅹ)中,
Figure BDA0003415730910000104
k是波数,x,y表示包裹相位图中每个像素点对应的坐标值,i 表示第i个拼接子镜,n为拼接子镜的个数,η*i是x轴倾斜误差,ηyi是y轴倾斜误差,j为整数;
该梯度曲线不再符合标准的费马螺旋,为了便于描述,将式(Ⅹ)称为变形的费马螺线, 并且设定(2R sinη*i,2R sinηyi)是变形的费马螺线的中心,通过曲线拟合得到的变形的费马螺 线中心位置从原点移动了x0和y0的距离,则倾斜误差如式(Ⅺ)所示:
Figure BDA0003415730910000111
式(Ⅺ)中,η*是x轴倾斜误差,ηy是y轴倾斜误差,x0和y0是通过曲线拟合得到的变形的 费马螺线的中心位置。
为直观展示,以两块半圆形子镜为例,说明了平移误差与倾斜误差对拟合的螺线旋转角及位 置特征的影响。
图3(a)为右侧拼接子镜不存在共相误差下的曲线拟合示例;图3(b)为右侧拼接子镜仅存在平移 误差下的曲线拟合示例;图3(c)为右侧拼接子镜仅存在x轴倾斜误差下的曲线拟合示例;图3(d) 为右侧拼接子镜仅存在y轴倾斜误差下的曲线拟合示例。
当子镜不存在共相误差时,它的包裹相位图梯度拟合的曲线为标准的费马螺线,当子镜存在 平移误差时,它的包裹相位图梯度拟合的曲线同样为标准的费马螺线,费马螺线的旋转角随平移误 差的变化而变化,当子镜存在倾斜误差时,它的包裹相位图梯度拟合的曲线不再是标准的费马曲线, 为了方面描述,称它为变形的费马螺旋,变形费马螺线的中心位置随倾斜误差的变化而变化。
当拼接镜17共相时,由其余子镜包裹相位图的梯度拟合的费马螺线与参考子镜拟合的费马螺 线重合,即旋转角与中心位置特征都相同。当其余拼接子镜与参考子镜仅存在平移误差时,由其余 拼接子镜包裹相位图的梯度拟合的费马螺线与参考子镜拟合的费马螺线中心位置相同,而旋转角不 同,将旋转角代入公式(Ⅶ)就反解得到其余子镜的平移误差。当其余拼接子镜与参考子镜仅存在 倾斜误差时,由其余拼接子镜包裹相位图的梯度拟合的费马螺线与参考子镜拟合的费马螺线中心位 置不同,旋转角也不同,将拟合得到的中心位置代入公式(Ⅺ)就可以反解得到其余拼接子镜的倾 斜误差。由于倾斜误差也会使螺旋的旋转角发生改变,因此在检测平移误差前需要先将倾斜误差校 正。
先对倾斜误差进行校正,具体是指:
选择一个校正过倾斜误差的拼接子镜为参考子镜,即该拼接子镜拟合的费马曲线为标准的费马 螺线,标准的费马螺线的中心位置为原点;
将该拼接子镜拟合的费马螺线的旋转角与中心位置特征(即x=0,y=0)作为参考值,将其它拼 接子镜拟合的费马螺线的旋转角、中心位置特征与参考值作比较,先通过中心位置特征反解出对应 拼接子镜的倾斜误差,并将其校正,倾斜误差校正后再通过曲线拟合得到费马螺线的旋转角,并反 解出对应拼接子镜的平移误差,将其余拼接子镜拟合得到的变形的费马螺线的中心位置特征代入式 (Ⅺ),即将x0和y0代入式(Ⅺ)得到对应拼接子镜的倾斜误差;在倾斜误差校正后,将其余拼接 子镜与参考子镜拟合得到的费马螺线的旋转角特征相减后代入式(Ⅶ),即将Δθi代入(Ⅶ)得到对 应拼接子镜的平移误差pi。曲线拟合过程、旋转角与位置特征提取过程可以通过计算机很容易的实 时实现。
干涉条纹和费马螺旋的分布是周期性的。当用单色波长测量时,相差λ/2的平移误差将产生相 同的条纹分布,即存在2π模糊。因此,单色光检测时,平移误差的检测量程仅为λ/2。因此,本 发明采用双波长光源切换测量即双波长算法来解决2π模糊问题,唯一的确定拼接镜17的平移误差, 并使得测量范围扩展到合成波长的一半即几微米。
(5)PC终端18将共相误差反馈给各拼接子镜的促动器实现共相误差的校正。具体是指:将 计算的共相误差作为校正值,引导促动器对各拼接子镜的位姿进行调整,促动器控制各拼接子镜的 左右旋转、前后俯仰与沿着光轴的前后移动。

Claims (10)

1.一种基于双波长涡旋光束相移干涉的拼接主镜共相误差检测系统,其特征在于,包括:第一激光器、第一扩束器、第二激光器、第二扩束器、第一分束器、第二分束器、四分之一波片、零级涡旋波片、第三分束器、平面反射镜、PZT位移平台、聚焦透镜、第四分束器、调制透镜、第五分束器、CCD相机、拼接镜、PC终端;
由第一激光器、第二激光器发出的不同波长的单色偏振光分别经过第一扩束器、第二扩束器扩束成平行光,第一激光器、第二激光器不同时使用,第一分束器将经过第一扩束器、第二扩束器扩束成的平行光垂直入射至第二分束器,分成透射光T1与反射光F1;
反射光F1经过四分之一波片变成圆偏振光,再经过零级涡旋波片转化为涡旋光束,涡旋光束入射到第三分束器上,并由第三分束器反射至平面反射镜上后经过第三分束器原路返回作为参考光入射至第五分束器上;
透射光T1经聚焦透镜聚焦后入射至第四分束器,并透射到拼接镜上,拼接镜原路返回的携带拼接镜共相误差信息的球面波作为测试光由第四分束器反射,后再经由调制透镜调制后入射到第五分束器上;
参考光即涡旋光束与测试光即球面光束在第五分束器处合束并发生干涉,CCD相机捕获干涉条纹;PC终端先将干涉条纹处理得到包裹相位图,然后提取包裹相位图的梯度线拟合为费马螺线,通过费马螺线的旋转角与位置特征反解出各个拼接子镜的平移误差与倾斜误差,再由PC终端将共相误差反馈给各拼接子镜的促动器实现共相误差的校正,拼接镜包括多个拼接子镜。
2.根据权利要求1所述的一种基于双波长涡旋光束相移干涉的拼接主镜共相误差检测系统,其特征在于,平面反射镜设置于PZT位移平台上用于移相。
3.根据权利要求1所述的一种基于双波长涡旋光束相移干涉的拼接主镜共相误差检测系统,其特征在于,第一分束器、第二分束器、第三分束器、第四分束器、第五分束器的反射面分别与所在光路的方向呈顺时针或逆时针45°夹角。
4.根据权利要求1所述的一种基于双波长涡旋光束相移干涉的拼接主镜共相误差检测系统,其特征在于,所述聚焦透镜的焦点与拼接镜的曲率中心重合。
5.权利要求1-4任一所述的基于双波长涡旋光束相移干涉的拼接主镜共相误差检测系统的工作方法,其特征在于,包括步骤如下:
(1)第一激光器、第二激光器发出的不同波长的单色偏振光分别经过第一扩束器、第二扩束器扩束成平行光,第一分束器将经过第一扩束器、第二扩束器扩束成的平行光垂直入射至第二分束器,分成透射光T1与反射光F1;
反射光F1经过四分之一波片变成圆偏振光,再经过零级涡旋波片转化为涡旋光束,涡旋光束入射到第三分束器上,并由第三分束器反射至平面反射镜上后经过第三分束器原路返回作为参考光入射至第五分束器上;
透射光T1经聚焦透镜聚焦后入射至第四分束器,并透射到拼接镜上,拼接镜原路返回的携带拼接镜共相误差信息的球面波作为测试光由第四分束器反射,再经由调制透镜调制后入射到第五分束器上;
(2)参考光即涡旋光束与测试光即球面光束在第五分束器处合束并发生干涉,CCD相机捕获干涉条纹;
(3)PC终端将干涉条纹处理得到包裹相位图;
(4)提取包裹相位图的梯度线拟合为费马螺线,通过费马螺线的旋转角与位置特征反解出各个拼接子镜的平移误差与倾斜误差;
(5)PC终端将共相误差反馈给各拼接子镜的促动器实现共相误差的校正。
6.根据权利要求5所述的基于双波长涡旋光束相移干涉的拼接主镜共相误差检测系统的工作方法,其特征在于,涡旋光束作为参考光束,其光场表达式如式(Ⅰ)所示:
Figure FDA0003415730900000021
测试光即球面光束,其光场表达式如式(Ⅱ)所示:
Figure FDA0003415730900000022
其中,Av、As分别是涡旋光与球面波的振幅,i是指表示第i个拼接子镜,m是拓扑电荷数,θ是涡旋光束的方位角,k是波数,z是传输距离,r是极坐标系下的径向分量,
Figure FDA0003415730900000023
是由拼接镜以及调制透镜引入的球面波前相位,常数
Figure FDA0003415730900000024
Figure FDA0003415730900000025
分别为对应于涡旋光与球面波的初相。
7.根据权利要求5所述的基于双波长涡旋光束相移干涉的拼接主镜共相误差检测系统的工作方法,其特征在于,当拼接镜存在共相误差,且取涡旋光的拓扑荷数为1时,干涉光强表达式如式(Ⅲ)所示:
Figure FDA0003415730900000026
式(Ⅲ)中,
Figure FDA0003415730900000027
k是波数,x,y表示干涉图像中每个像素点对应的坐标值,i表示第i个拼接子镜,n为拼接子镜的个数,pi是拼接子镜的平移误差,η*i是x轴倾斜误差,ηyi是y轴倾斜误差,2(pi+x sinη*i+y sinηyi)是由拼接镜共相误差而带来的光路差。
8.根据权利要求5所述的基于双波长涡旋光束相移干涉的拼接主镜共相误差检测系统的工作方法,其特征在于,步骤(3)中,PC终端通过三步相移法将干涉条纹处理得到包裹相位图,相移量分别为-2π/3,0,2π/3,三步相移法的表示式如式(Ⅳ)所示:
Figure FDA0003415730900000031
式(Ⅳ)中,
Figure FDA0003415730900000032
I1、I2、I3分别对应于三步相移过程中的干涉光强;
通过式(Ⅲ)与式(Ⅳ),得到相位表达式如式(Ⅴ)所示:
Figure FDA0003415730900000033
式(Ⅴ)中,
Figure FDA0003415730900000034
k是波数,x,y表示包裹相位图中每个像素点对应的坐标值,i表示第i个拼接子镜,n为拼接子镜的个数,pi是拼接子镜的平移误差,η*i是x轴倾斜误差,ηyi是y轴倾斜误差,2(pi+x sinη*i+y sinηyi)是由拼接镜共相误差而带来的光路差,I1、I2、I3分别对应于三步相移过程中的干涉光强。
9.根据权利要求5所述的基于双波长涡旋光束相移干涉的拼接主镜共相误差检测系统的工作方法,其特征在于,步骤(4)中,相位被包裹在-π到π区间,包裹相位图的相位分布有一个断层式的变化,即产生了一个梯度或边缘,这个梯度或边缘满足式(Ⅵ):
Figure FDA0003415730900000035
式(Ⅵ)中,
Figure FDA0003415730900000036
k是波数,x,y表示包裹相位图中每个像素点对应的坐标值,i表示第i个拼接子镜,n为拼接子镜的个数,pi是拼接子镜的平移误差,η*i是x轴倾斜误差,ηyi是y轴倾斜误差,2(pi+x sinη*i+y sinηyi)是由拼接镜共相误差而带来的光路差,j为整数;
当拼接镜中仅存在平移误差时,边缘提取的梯度曲线符合标准的费马螺旋线,其旋转角度θi随平移误差pi的变化而变化,拼接镜拟合的费马螺线的旋转角度与拼接镜的平移误差的关系如(Ⅶ)所示:
pi=Δθi/(2k) (Ⅶ)
式(Ⅶ)中,pi是第i个拼接子镜的平移误差,Δθi是由第i个拼接子镜与参考子镜拟合的费马螺线的旋转角的差值,k是波数;
费马螺线的旋转角与位置特征通过曲线拟合得到;
合成波长的表达式如式(Ⅷ)所示:
Figure FDA0003415730900000041
式(Ⅷ)中,λ1、λ2分别是指第一激光器、第二激光器的波长;
假如拼接子镜的实际平移误差为p,利用式(Ⅶ)得到的两波长的平移误差分别为p1和p2,则平移误差p的求取公式如式(Ⅸ)所示:
Figure FDA0003415730900000042
式(Ⅸ)中,u、v是待求解的常数,对两波长上求得的平移误差p进行比较,当两波长上p最匹配时,此时对应的u与v的值就为求解值;
当拼接镜中仅存在倾斜误差时,将式(Ⅵ)整理,得到边缘提取的梯度曲线即变形的费马螺线如式(Ⅹ)所示:
Figure FDA0003415730900000043
式(Ⅹ)中,
Figure FDA0003415730900000044
k是波数,x,y表示包裹相位图中每个像素点对应的坐标值,i表示第i个拼接子镜,n为拼接子镜的个数,η*i是x轴倾斜误差,ηyi是y轴倾斜误差,j为整数;
将式(Ⅹ)称为变形的费马螺线,并且设定(2R sinη*i,2R sinηyi)是变形的费马螺线的中心,通过曲线拟合得到的变形的费马螺线中心位置从原点移动了x0和y0的距离,则倾斜误差如式(Ⅺ)所示:
Figure FDA0003415730900000045
式(Ⅺ)中,η*是x轴倾斜误差,ηy是y轴倾斜误差,x0和y0是通过曲线拟合得到的变形的费马螺线的中心位置。
10.根据权利要求5所述的基于双波长涡旋光束相移干涉的拼接主镜共相误差检测系统的工作方法,其特征在于,步骤(5)之前先对倾斜误差进行校正,具体是指:
选择一个校正过倾斜误差的拼接子镜为参考子镜,即该拼接子镜拟合的费马曲线为标准的费马螺线,标准的费马螺线的中心位置为原点;
将该拼接子镜拟合的费马螺线的旋转角与中心位置特征作为参考值,将其它拼接子镜拟合的费马螺线的旋转角、中心位置特征与参考值作比较,先通过中心位置特征反解出对应拼接子镜的倾斜误差,并将其校正,倾斜误差校正后再通过曲线拟合得到费马螺线的旋转角,并反解出对应拼接子镜的平移误差,将其余拼接子镜拟合得到的变形的费马螺线的中心位置特征代入式(Ⅺ),即将x0和y0代入式(Ⅺ)得到对应拼接子镜的倾斜误差;在倾斜误差校正后,将其余拼接子镜与参考子镜拟合得到的费马螺线的旋转角特征相减后代入式(Ⅶ),即将Δθi代入(Ⅶ)得到对应拼接子镜的平移误差pi
将使用双波长测量计算的共相误差反馈给促动器来实时校正各拼接子镜的共相误差,具体是指:将计算的共相误差作为校正值,引导促动器对各拼接子镜的位姿进行调整,促动器控制各拼接子镜的左右旋转、前后俯仰与沿着光轴的前后移动。
CN202111545618.2A 2021-12-16 2021-12-16 一种基于双波长涡旋光束相移干涉的拼接主镜共相误差检测系统及其工作方法 Active CN114322829B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111545618.2A CN114322829B (zh) 2021-12-16 2021-12-16 一种基于双波长涡旋光束相移干涉的拼接主镜共相误差检测系统及其工作方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111545618.2A CN114322829B (zh) 2021-12-16 2021-12-16 一种基于双波长涡旋光束相移干涉的拼接主镜共相误差检测系统及其工作方法

Publications (2)

Publication Number Publication Date
CN114322829A true CN114322829A (zh) 2022-04-12
CN114322829B CN114322829B (zh) 2022-12-20

Family

ID=81052194

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111545618.2A Active CN114322829B (zh) 2021-12-16 2021-12-16 一种基于双波长涡旋光束相移干涉的拼接主镜共相误差检测系统及其工作方法

Country Status (1)

Country Link
CN (1) CN114322829B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115046476A (zh) * 2022-06-08 2022-09-13 中国科学院长春光学精密机械与物理研究所 一种稀疏孔径望远镜边缘传感方法
CN116385547A (zh) * 2023-06-02 2023-07-04 中国科学院云南天文台 一种天文望远镜主镜面共相检测方法、系统及计算机设备

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101276056A (zh) * 2007-12-04 2008-10-01 中国科学院光电技术研究所 高精度拼接子镜相对倾斜误差光电检测系统
TW201033583A (en) * 2009-03-10 2010-09-16 Univ Nat Cheng Kung Interferometer and interference method for generating stable vortex beam
CN105700128A (zh) * 2016-05-03 2016-06-22 中国科学院上海天文台 一种拼接望远镜共相位控制装置及控制方法
CN108007572A (zh) * 2017-11-23 2018-05-08 南京大学 一种基于涡旋光束和萨格拉克干涉仪的旋转扰动测量系统
CN109297930A (zh) * 2018-11-14 2019-02-01 深圳大学 一种基于涡旋光束共轭干涉的三阶非线性测量装置及方法
CN110487212A (zh) * 2019-08-02 2019-11-22 中北大学 一种基于涡旋光螺旋相位相移干涉检测物体面型装置

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101276056A (zh) * 2007-12-04 2008-10-01 中国科学院光电技术研究所 高精度拼接子镜相对倾斜误差光电检测系统
TW201033583A (en) * 2009-03-10 2010-09-16 Univ Nat Cheng Kung Interferometer and interference method for generating stable vortex beam
CN105700128A (zh) * 2016-05-03 2016-06-22 中国科学院上海天文台 一种拼接望远镜共相位控制装置及控制方法
CN108007572A (zh) * 2017-11-23 2018-05-08 南京大学 一种基于涡旋光束和萨格拉克干涉仪的旋转扰动测量系统
CN109297930A (zh) * 2018-11-14 2019-02-01 深圳大学 一种基于涡旋光束共轭干涉的三阶非线性测量装置及方法
CN110487212A (zh) * 2019-08-02 2019-11-22 中北大学 一种基于涡旋光螺旋相位相移干涉检测物体面型装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
LILI YANG等: "Co-phase state detection for segmented mirrors", 《OPTICS EXPRESS》 *
张禹: "大口径多光谱通道波前测量系统设计", 《红外与激光工程》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115046476A (zh) * 2022-06-08 2022-09-13 中国科学院长春光学精密机械与物理研究所 一种稀疏孔径望远镜边缘传感方法
CN116385547A (zh) * 2023-06-02 2023-07-04 中国科学院云南天文台 一种天文望远镜主镜面共相检测方法、系统及计算机设备
CN116385547B (zh) * 2023-06-02 2023-08-01 中国科学院云南天文台 一种天文望远镜主镜面共相检测方法、系统及计算机设备

Also Published As

Publication number Publication date
CN114322829B (zh) 2022-12-20

Similar Documents

Publication Publication Date Title
CN102589416B (zh) 用于非球面测量的波长扫描干涉仪及方法
CN114322829B (zh) 一种基于双波长涡旋光束相移干涉的拼接主镜共相误差检测系统及其工作方法
US6690474B1 (en) Apparatus and methods for surface contour measurement
CN102289152B (zh) 光学系统波像差检测装置
CA2397095C (en) Apparatus and methods for surface contour measurement
CN102589440B (zh) 连续变角度数字全息干涉测量的方法
WO2005052502A2 (en) Pixelated phase-mask interferometer
WO1999063304A1 (en) Apparatus and methods for surface contour measurement
CN108844492A (zh) 一种基于光谱调制度深度编码的微结构形貌测量方法及其装置
CN101324421A (zh) 同步移相菲索干涉仪
CN113196003B (zh) 分别用于确定输入光场的输入相位和/或输入振幅的方法,干涉仪和信号处理设备
CN108895986B (zh) 基于条纹成像投影的显微三维形貌测量装置
CN105300273B (zh) 条纹对比度可调的动态点衍射干涉仪
CN208635741U (zh) 一种基于光谱调制度深度编码的微结构形貌测量装置
CN108680108B (zh) 线激光移相干涉三角微位移测量装置及方法
CN104296676A (zh) 基于低频差声光移频器移相的外差点衍射干涉仪
CN102401630B (zh) 空间移相菲索球面干涉仪
KR100916593B1 (ko) 실시간 3차원 형상 측정 시스템
CN103712569A (zh) 基于偏转角度的单幅图像快速相移系统及相位检测方法
CN114354141A (zh) 一种基于频域测量超表面相位的方法及系统
CN101907445A (zh) 大口径单层薄膜厚度全场检测装置
CN105674875B (zh) 一种全视场低频外差点衍射干涉仪
CN106644076A (zh) 移相式干涉光谱成像系统和成像方法
CN201251428Y (zh) 同步移相菲索干涉仪
CN108362222A (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