基于二次图像补偿的相位差波前传感器 下载: 975次
1 引言
波前传感器(WFS)主要包括夏克-哈特曼传感器[1]、四棱锥传感器[2]、曲率传感器[3]、剪切干涉仪[4-5]和相位差波前传感器(PDWFS)[6]等,在自适应光学系统中起重要作用[7]。相比WFS,PDWFS只需要采集焦面和离焦面的图像信息,具有结构简单、易于实现的优点。且检测精度高、对环境要求低。随着计算机速度的发展,PDWFS在自适应光学系统中具有更大的应用潜力[8-9]。
1979年,Gonsalves和Chidlaw[10]提出了相位差(PD)算法,Dolne等[11]对PD算法进行了优化。Löfdahl等[12]成功使用PD算法重建了太阳表面图像。哈勃太空望远镜也采用了PD算法进行图像观测[13]。
PDWFS是电荷耦合器件(CCD),探测的光强信息中,会夹杂一些噪声数据[14-15],如背景噪声、信号光子噪声和CCD的读取噪声。其中读取噪声是一种高斯分布噪声,对图像的退化影响较大[16-18]。根据傅里叶变换,可将透镜的焦面图像看作是近场图像的频谱。由于焦面图像光强主要集中在图像中部的低空间频率区域,而周边的高频信息能量较弱,因此噪声对探测图像高频信息的影响较大。可以通过增加图像高频区域的信噪比(SNR),提高PDWFS的测量精度。Luo等[19]提出了在透镜前增加一个衍射光栅以提高图像信噪比的方法,根据CCD在接收面处不同衍射阶次的能量关系,计算得到焦距和离焦面上的光强分布,提高了PDWFS的探测精度。但其结构复杂、能量比例关系固定、光栅制作难度大以及对超大靶面CCD的要求限制了该方法在PDWFS中的实际应用。
当CCD相机接收的光强较弱时,采集到的光强分布受噪声影响严重。因此使用CCD相机时,要尽量增大曝光时间,以提高探测的信噪比。但曝光时间过大会造成CCD相机拍摄的图像饱和,导致无法探测。本文提出了一种基于图像补偿的相位差波前传感器,其光学结构比较简单。通过改变CCD相机的曝光时间,在焦面、离焦面上各采集两副图像,其中光强较大的图像,中间部分虽有饱和,但相比常规探测周围的高频信息信噪比较高,可准确探测。且该饱和区域可通过图像补偿算法由光强较小、未饱和的图像计算出来。将两幅图像融合成一幅高低频信息都非常准确的图像,可大大提高CCD相机的动态范围,有效抑制随机噪声。
2 基本原理
2.1 传统的相位反演算法
传统的PDWFS系统如
在忽略系数的情况下,可将焦面光场表示为[9]
式中,(x,y)为对应位置坐标,f为对应焦距,F为傅里叶变换,Enear为输入光场,φnear(x,y)为近场波前相位。近场光强分布可以通过焦面光场表示为
式中,F-为逆傅里叶变换,φf(x,y)为焦面相位。
可将离焦面的光场分布等效为在光瞳面引入一个附加相位差Δφd(x,y)
近场波前相位可以表示为
相位差波前复原算法流程如
当|Ed(x,y)|满足(5)式时,上述循环停止。
式中,XSSE为评价指标,σ为评价因子,通常设置为一个足够小的数值,实验中设为10-7,当达到判定条件,输出φnear(x,y)。
图 3. 带像差的远场图像。(a)焦面图像;(b)离焦面图像;(c)焦面饱和图像;(d)离焦面饱和图像;(e)焦面融合图像;(f)离焦面融合图像
Fig. 3. Far-field images with aberration. (a) Image of the focal plane; (b) image of the defocused plane; (c) image of the saturated focal plane; (d) image of the saturated defocused plane; (e) fusion image of the focal plane; (f) fusion image of defocused plane
受制作工艺、测试环境的影响,不可避免地会产生噪声。实际的焦面光强分布和离焦面光强分布可以表示为
式中,I'f(x,y)和I'd(x,y)分别为理想的焦面和离焦面光强分布,nf(x,y)和nd(x,y)为附加噪声。
2.2 基于图像补偿的相位差波前传感技术
通过基于图像补偿的相位差波前传感器(AET-PDWFS),在正常曝光时间下,即图像峰值光强接近相机饱和光强时,拍摄一幅图像;在饱和曝光时间下,即图像中间峰值部分饱和时,使图像对应点光强是原光强的数倍(一般取3~20之间),拍摄一幅图像,将两幅图像融合成一幅图像。
先模拟接近相机饱和情况下的正常焦面光强分布If0(x,y)和离焦面光强分布图像Id0(x,y),如
在焦面光强If0(x,y)和焦面饱和光强Ifs(x,y)中取对应的N个有效像素,可以得到饱和图像相对焦面图像的强度倍率
式中,N为两幅图像中未饱和区域均匀采集的对应点,一般取10,i∈[1,N]为序号。
Ids(x,y)和Id0(x,y)之间的比例系数可表示为
rf和rd取值范围为3~20。用If0(x,y)相应的像素光强乘rf代替Ifs(x,y)中饱和像素的光强得到Ifa(x,y),用Id0(x,y)相应的像素光强乘以rd代替Ids(x,y)中饱和像素的光强得到Ida(x,y)。新的光强分布Ifa(x,y)和Ida(x,y)分别如
使用计算得到的融合光强信息,可显著提高相位差法波前测量的精度。具体流程如
3 仿真分析
为验证AET-PDWFS的有效性和测量精度,进行了大量的数值模拟实验,并与传统PDWFS进行了对比。模拟参数:利用前65阶Zernike多项式生成了峰谷(PV)为0.5734λ和均方根(RMS)为0.0847λ的随机波前像差,如
采用PDWFS和AET-PDWFS重建
对随机波前在不同的信噪比下进行复原仿真,结果如
图 6. 光强分布。(a) PDWFS的焦面图像;(b) PDWFS的离焦面图像;(c) AET-PDWFS的焦面图像;(d) AET-PDWFS的离焦面图像
Fig. 6. Distribution of light intensity. (a) Focal plane image of PDWFS; (b) defocused plane image of PDWFS; (c) focal plane image of AET-PDWFS; (d) defocused plane image of AET-PDWFS
图 7. 波前复原结果。(a) PDWFS复原的波前;(b) AET-PDWFS复原的波前;(c) PDWFS复原的波前残差;(d) AET-PDWFS复原的波前残差;(e) PDWFS重构的Zernike系数;(f) AET-PDWFS重构的Zernike系数
Fig. 7. Results of wavefront restoration. (a) Wavefront of PDWFS restored; (b) wavefront of AET-PDWFS restored; (c) wavefront residual of PDWFS restored; (d) wavefront residual of AET-PDWFS restored; (e) Zernike coefficient reconstructed by PDWFS; (f) Zernike coefficient reconstructed by AET-PDWFS
图 8. 不同信噪比下两种方法的波面复原精度对比。(a) PV;(b) RMS
Fig. 8. Comparison of wavefront restoration accuracy of the two methods under different signal-to-noise ratios. (a) PV; (b) RMS
4 实验研究
为进一步验证AET-PDWFS的有效性,设计搭建了一套实验系统,如
采集到的图像如
由(8)式和(9)式,得到如
首先用干涉仪对像差板面形进行测量,实验使用的干涉仪为ZYGO公司生产的MARK IV型移相干涉仪,其标准镜头的PV为λ/20,RMS为λ/120。测量时温度范围为-30~60 ℃,相对湿度为95%,PV可重复性大于λ/100。MARK IV干涉仪测得的像差板面形如
分别通过PDWFS和AET-PDWFS进行波前复原,结果如
图 10. 实验采集的图像。(a)焦面常规图像;(b)离焦面常规图像;(c)焦面饱和图像;(d)离焦面饱和图像
Fig. 10. Images collected from experiment. (a) Regular image of the focal plane; (b) regular image of the defocused plane; (c) saturated image of the focal plane; (d) saturated image of the defocused plane
图 11. AET-PDWFS融合的图像。(a)焦面图像;(b)离焦面图像
Fig. 11. Images of the AET-PDWFS fusion. (a) Image of the focal plane; (b) image of the defocused plane
图 12. MARK IV 干涉仪测得的像差板面形
Fig. 12. Aberration plate shape measured by MARK IV interferometer
图 13. 波前复原结果。(a) PDWFS;(b) AET-PDWFS
Fig. 13. Results of wavefront reconstruction. (a) PDWFS; (b) AET-PDWFS
5 结论
阐述了AET-PDWFS的基本原理,通过数值仿真,证明了该传感器相比传统PDWFS在波前复原精度上有一定提高,特别在低信噪比的条件下,优势更加明显。设计搭建了实验系统,结果表明,该传感器操作简单,且无需改变光路、添加器件,只需要改变曝光时间即可实现较高波前的复原精度,具有较高的实用价值。
[1] Lane R G, Tallon M. Wave-front reconstruction using a Shack-Hartmann sensor[J]. Applied Optics, 1992, 31(32): 6902-6908.
[2] Chamot S R, Dainty C, Esposito S. Adaptive optics for ophthalmic applications using a pyramid wavefront sensor[J]. Optics Express, 2006, 14(2): 518-526.
[3] Roddier F. Curvature sensing and compensation: a new concept in adaptive optics[J]. Applied Optics, 1988, 27(7): 1223-1225.
[4] Wyant J C. Use of an ac heterodyne lateral shear interferometer with real-time wavefront correction systems[J]. Applied Optics, 1975, 14(11): 2622-2626.
[5] 纵榜铭, 栾嘉蕴, 蒋志龙, 等. 基于叠层扫描成像技术的成像透镜透射波前测量方法[J]. 中国激光, 2019, 46(8): 0804001.
[6] Gonsalves R A. Phase retrieval and diversity in adaptive optics[J]. Optical Engineering, 1982, 21(5): 829-832.
[7] Hardy J W. Adaptive optics: a progress review[J]. Proceedings of SPIE, 1991, 1542: 2-17.
[8] Zhang D, Zhang X B, Xu S Y, et al. Simplified phase diversity algorithm based on a first-order Taylor expansion[J]. Applied Optics, 2016, 55(28): 7872-7877.
[9] Jefferies S M, Lloyd-Hart M, Hege E K, et al. Sensing wave-front amplitude and phase with phase diversity[J]. Applied Optics, 2002, 41(11): 2095-2102.
[10] Gonsalves R A, Chidlaw R. Wavefront sensing by phase retrieval[J]. Proceedings of SPIE, 1979, 0207: 32-39.
[11] Dolne J J, Tansey R J, Black K A, et al. Practical issues in wave-front sensing by use of phase diversity[J]. Applied Optics, 2003, 42(26): 5284-5289.
[12] Löfdahl M G, Kendrick R L, Harwit A, et al. Phase diversity experiment to measure piston misalignment on the segmented primary mirror of the Keck II Telescope[J]. Proceedings of SPIE, 1998, 3356: 1190-1201.
[13] Fienup J R, Marron J C, Schulz T J, et al. Hubble space telescope characterized by using phase-retrieval algorithms[J]. Applied Optics, 1993, 32(10): 1747-1767.
[14] 李文力, 李凯, 彭迪, 等. AGRI遥感图像中条带噪声的分析与去除[J]. 光学学报, 2019, 39(12): 1228004.
[15] 郑成林, 何顶顶, 费庆国. 基于灰度梯度正则化去噪的改进数字图像相关法[J]. 光学学报, 2018, 38(8): 0812002.
[17] He Y X, Li X Y. Error analysis of laser beam quality measured with CCD sensor and choice of the optimal threshold[J]. Optics & Laser Technology, 2013, 45: 671-677.
[18] Snyder D L, Helstrom C W, Lanterman A D, et al. Compensation for readout noise in CCD images[J]. Journal of the Optical Society of America A, 1995, 12(2): 272-283.
[19] Luo Q, Huang L H, Gu N T, et al. A modified phase diversity wavefront sensor with a diffraction grating[J]. Chinese Physics B, 2012, 21(9): 094201.
Article Outline
孔庆峰, 王帅, 杨平, 林海奇, 刘永, 许冰. 基于二次图像补偿的相位差波前传感器[J]. 激光与光电子学进展, 2020, 57(12): 121103. Qingfeng Kong, Shuai Wang, Ping Yang, Haiqi Lin, Yong Liu, Bing Xu. Phase Diversity Wavefront Sensor Based on Secondary Image Compensation[J]. Laser & Optoelectronics Progress, 2020, 57(12): 121103.