激光与光电子学进展, 2020, 57 (12): 121103, 网络出版: 2020-06-03  

基于二次图像补偿的相位差波前传感器 下载: 975次

Phase Diversity Wavefront Sensor Based on Secondary Image Compensation
孔庆峰 1,2,3,4王帅 1,3,*杨平 1,3,**林海奇 1,3,4刘永 2许冰 1,3
作者单位
1 中国科学院自适应光学重点实验室, 四川 成都 610209
2 电子科技大学光电科学与工程学院, 四川 成都 610054
3 中国科学院光电技术研究所, 四川 成都 610209
4 中国科学院大学, 北京 100049
摘要
电荷耦合器件(CCD)的测量精度是影响波前复原精度的主要因素。受相机制造工艺的限制以及环境的影响,不可避免地会采集到噪声数据,从而影响低信噪比、高空间频率图像信息的采集。为了降低噪声对波前复原的影响,提出了一种基于图像补偿的改进相位差波前传感器。先在CCD相机的正常曝光时间和饱和曝光时间下分别采集两幅图像;然后融合成一幅包含准确高空间频率和低空间频率光强信息的图像。在相同信噪比下,对改进的相位差波前传感器和传统的相位差波前传感器进行了仿真。结果表明,改进的相位差波前传感器的波前反演精度有显著提高。最后通过实验进一步验证了该方法的有效性。
Abstract
The measurement precision of charge coupled device (CCD) is the main factor affecting the accuracy of wavefront reconstruction. Due to the limitations of camera manufacturing technology and environmental impact, it is inevitable that noise data will be collected, which will affect the image information collection with low signal-to-noise-ratio and high spatial frequency. An improved phase diversity wavefront sensor based on image compensation is proposed to reduce the influence of noise on wavefront reconstruction. First, two images are collected under the normal exposure time and saturated exposure time of CCD respectively. And then they are sewed into an image containing accurate high spatial frequency and low spatial frequency light intensity information. At the same signal-to-noise ratio, the improved phase diversity wavefront sensor and the traditional phase diversity wavefront sensor are simulated. The results show that the wavefront inversion accuracy of the improved phase difference wavefront sensor is significantly improved. Finally, the effectiveness of the method is further verified through experiments.

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系统如图1所示。待测光束经过透镜后被分光片(BS)分成两束,分别成像于相机焦面和相机离焦面靶面,其中d为离焦距离。

图 1. 光学系统模型

Fig. 1. Model of optical system

下载图片 查看所有图片

在忽略系数的情况下,可将焦面光场表示为[9]

Ef(x,y)=F{|Enear(x,y)|·exp[iφnear(x,y)]},(1)

式中,(x,y)为对应位置坐标,f为对应焦距,F为傅里叶变换,Enear为输入光场,φnear(x,y)为近场波前相位。近场光强分布可以通过焦面光场表示为

Enear(x,y)=F-|Ef(x,y)|·exp[iφf(x,y)]},(2)

式中,F-为逆傅里叶变换,φf(x,y)为焦面相位。

可将离焦面的光场分布等效为在光瞳面引入一个附加相位差Δφd(x,y)

Ed(x,y)=F{|Enear(x,y)|·exp[iφnear(x,y)+φd(x,y)]}(3)

近场波前相位可以表示为

Enear(x,y)=F-|Ed(x,y)|·exp[iφd(x,y)]}/exp[φd(x,y)](4)

相位差波前复原算法流程如图2所示。设初始φnear(x,y)为零,计算得到的|Enear(x,y)|、|Ef(x,y)|和|Ed(x,y)|分别用实际近场强度分布的平方根 Inear(x,y)、焦面强度分布的平方根 If(x,y)和离焦面强度分布的平方根 Id(x,y)代替。通常假定 Inear(x,y)在瞳孔内为1,其他区域为零。

当|Ed(x,y)|满足(5)式时,上述循环停止。

XSSE=|Ed(x,y)|-Id(x,y)]2/Id(x,y)<σ,(5)

式中,XSSE为评价指标,σ为评价因子,通常设置为一个足够小的数值,实验中设为10-7,当达到判定条件,输出φnear(x,y)。

图 2. PD算法流程图

Fig. 2. Flow chart of PD algorithm

下载图片 查看所有图片

图 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

下载图片 查看所有图片

受制作工艺、测试环境的影响,不可避免地会产生噪声。实际的焦面光强分布和离焦面光强分布可以表示为

If(x,y)=I'f(x,y)+nf(x,y),(6)

Id(x,y)=I'd(x,y)+nd(x,y),(7)

式中,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),如图3(a)和图3(b)所示,然后增加曝光时间至饱和,得到焦面饱和光强分布Ifs(x,y)和离焦面饱和光强分布图像Ids(x,y),如图3(c)和图3(d)所示。可以发现,增加曝光时间后,图像中间部分是饱和的,但其周围的高频光强信息却能被很好地探测到。

在焦面光强If0(x,y)和焦面饱和光强Ifs(x,y)中取对应的N个有效像素,可以得到饱和图像相对焦面图像的强度倍率

rf=i=1NIfs(xi,yi)/If0(xi,yi)N,(8)

式中,N为两幅图像中未饱和区域均匀采集的对应点,一般取10,i∈[1,N]为序号。

Ids(x,y)和Id0(x,y)之间的比例系数可表示为

rd=i=1NIds(xj,yj)/Id0(xj,yj)N(9)

rfrd取值范围为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(e)和图3(f)所示,既包含了图3(a)和图3(b)中光强峰值附近的光强信息,又融合了图3(c)和图3(d)中光强外围的高频光强信息。

使用计算得到的融合光强信息,可显著提高相位差法波前测量的精度。具体流程如图4所示。

图 4. AET-PDWFS流程图

Fig. 4. Flow chart of AET-PDWFS

下载图片 查看所有图片

3 仿真分析

为验证AET-PDWFS的有效性和测量精度,进行了大量的数值模拟实验,并与传统PDWFS进行了对比。模拟参数:利用前65阶Zernike多项式生成了峰谷(PV)为0.5734λ和均方根(RMS)为0.0847λ的随机波前像差,如图5所示。焦距为500 mm,离焦面在距离透镜505 mm处,通光孔径为4 mm,相机信噪比为26 dB,两种方法分别迭代50次。

图 5. 入射波前像差分布

Fig. 5. Aberration distribution of incident wavefront

下载图片 查看所有图片

图6为传统方法和改进方法的焦面采集图像,为方便观察,对光强值取对数。相比传统PDWFS的焦面和离焦面强度图像(图6(a)和图6(b)),基于AET-PDWFS的焦面和离焦面强度图像光强提高了10倍,通过图像融合得到新的焦面和离焦面强度图像,分别如图6(c)和图6(d)所示。可以发现,改进的AET-PDWFS光强高频细节比传统的PDWFS更丰富。

采用PDWFS和AET-PDWFS重建图5所示的随机入射波前,复原结果如图7所示。可以发现,传统PDWFS重构波前(PV为0.3258λ,RMS为0.0686λ)的细节与待测波前(PV为0.5734λ,RMS为0.0847λ)明显不同。对比图7(a),可以发现,图7(b)中AET-PDWFS的复原波前与待测波前非常接近。定义残余波前为复原波前与待测波前之差。图7(c)所示的传统PDWFS的恢复波前残余误差(PV为0.3134λ,RMS为0.0199λ)远大于图7(d)中AET-PDWFS的残余误差(PV为0.0609λ,RMS为0.003λ)。当复原波前用前65阶Zernike多项式表示时,可以发现,传统PDWFS的Zernike复原系数部分阶次的复原值与实际值差异较大,如图7(e)所示,黑色条柱为实际值,白色条柱为复原值。而图7(f)中AET-PDWFS的复原系数与实际值非常接近。特别是对于高阶Zernike多项式,AET-PDWFS的测量精度要比传统的PDWFS高很多,这表明在一定的噪声条件下,改进算法的波前复原精度有明显提高。

对随机波前在不同的信噪比下进行复原仿真,结果如图8所示。可以发现传统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的有效性,设计搭建了一套实验系统,如图9所示。系统中使用的光源是波长为632.8 nm的He-Ne激光器,光源发出的光经衰减片、空间滤波器滤波后入射到准直透镜,形成直径为70 mm的平行光。平行光经过光阑,形成口径为35 mm且能与像差板相匹配的光束。从像差板出射的光经缩束系统缩束后,进入成像透镜和远场相机。采用采样位数为12位的美国IMPERX TXG03型相机,通过前后移动CCD相机分别得到焦面和离焦面图像。在相同实验条件下,对PDWFS算法和AET-PDWFS算法进行对比。

图 9. 实验装置示意图

Fig. 9. Schematic diagram of the experimental device

下载图片 查看所有图片

采集到的图像如图10所示。图10(a)为焦面采集到的常规图像,图10(b)为离焦距离为2 mm时采集到的离焦面常规图像;图10(c)为提高CCD相机的曝光时间采集到的焦面饱和图像,对应光强为原光强的4倍;图10(d)为提高CCD相机曝光时间采集到的离焦面饱和图像,对应光强为原光强的8倍。

由(8)式和(9)式,得到如图11所示的新的焦面和离焦面的图像,可以明显发现,新图像相比常规采集的图像有更加丰富的高频信息。

首先用干涉仪对像差板面形进行测量,实验使用的干涉仪为ZYGO公司生产的MARK IV型移相干涉仪,其标准镜头的PV为λ/20,RMS为λ/120。测量时温度范围为-30~60 ℃,相对湿度为95%,PV可重复性大于λ/100。MARK IV干涉仪测得的像差板面形如图12所示,待测波前PV为1.287λ,RMS为0.258λ

分别通过PDWFS和AET-PDWFS进行波前复原,结果如图13所示。图13(a)为PDWFS算法的复原波前,PV为1.12λ,RMS为0.239λ。而AET-PDWFS算法的复原波前PV为1.291λ,RMS为0.248λ,与干涉仪测试结果非常接近。这表明AET-PDWFS相较于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.

    Zong B M, Luan J Y, Jiang Z L, et al. Method to measure transmitting wave-front of imaging lens based on extended ptychographical iterative engine[J]. Chinese Journal of Lasers, 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.

    Li W L, Li K, Peng D, et al. Analysis and removal of stripe noise in AGRI remote-sensing images[J]. Acta Optica Sinica, 2019, 39(12): 1228004.

[15] 郑成林, 何顶顶, 费庆国. 基于灰度梯度正则化去噪的改进数字图像相关法[J]. 光学学报, 2018, 38(8): 0812002.

    Zheng C L, He D D, Fei Q G. Improved digital image correlation method based on gray gradient denoised by regularization method[J]. Acta Optica Sinica, 2018, 38(8): 0812002.

[16] Laligant O, Truchetet F, Fauvet E. Noise estimation from digital step-model signal[J]. IEEE Transactions on Image Processing, 2013, 22(12): 5158-5167.

[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.

孔庆峰, 王帅, 杨平, 林海奇, 刘永, 许冰. 基于二次图像补偿的相位差波前传感器[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.

引用该论文: TXT   |   EndNote

相关论文

加载中...

关于本站 Cookie 的使用提示

中国光学期刊网使用基于 cookie 的技术来更好地为您提供各项服务,点击此处了解我们的隐私策略。 如您需继续使用本网站,请您授权我们使用本地 cookie 来保存部分信息。
全站搜索
您最值得信赖的光电行业旗舰网络服务平台!