CN112099229B - 基于远场的高速自适应光学闭环控制方法 - Google Patents
基于远场的高速自适应光学闭环控制方法 Download PDFInfo
- Publication number
- CN112099229B CN112099229B CN202011190608.7A CN202011190608A CN112099229B CN 112099229 B CN112099229 B CN 112099229B CN 202011190608 A CN202011190608 A CN 202011190608A CN 112099229 B CN112099229 B CN 112099229B
- Authority
- CN
- China
- Prior art keywords
- wavefront
- distribution
- far
- far field
- pupil
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 33
- 230000003287 optical effect Effects 0.000 title claims abstract description 13
- 238000009826 distribution Methods 0.000 claims abstract description 86
- 230000004075 alteration Effects 0.000 claims abstract description 20
- 230000008859 change Effects 0.000 claims abstract description 13
- 230000008569 process Effects 0.000 claims abstract description 8
- 210000001747 pupil Anatomy 0.000 claims description 36
- 230000003044 adaptive effect Effects 0.000 claims description 15
- 230000000295 complement effect Effects 0.000 claims description 5
- 238000011084 recovery Methods 0.000 claims description 3
- 238000009827 uniform distribution Methods 0.000 claims description 3
- 238000001514 detection method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G02—OPTICS
- G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
- G02B27/00—Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00
- G02B27/0025—Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00 for optical correction, e.g. distorsion, aberration
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Physics & Mathematics (AREA)
- Optics & Photonics (AREA)
- Eye Examination Apparatus (AREA)
Abstract
本发明公开了一种基于远场的高速自适应光学闭环控制方法,该方法利用的装置由远场采集模块,远场到波前复原模块,高压放大器和波前校正器模块组成,远场采集模块实时记录远场图像变化过程,波前复原模块根据当前帧远场光强分布和前一帧的远场光强分布复原波前,根据复原获得的波前信息驱动高压放大器控制波前校正器校正光路中的相差。采用单远场光强分布进行闭环控制,无需增加额外的远场信息,如离焦面上的远场光斑分布,可以简化光路结构,同时提高光能利用率;可以根据前后两帧远场确定入射波前的位相,特别是经过次数较少的自适应光学校正或输入像差较小的情况下,可以准确的获得像差数据,从而实现快速的对像差进行准确校正,实时性更强。
Description
技术领域
本发明涉及一种波前闭环控制方法,特别是一种基于远场的高速自适应光学闭环控制方法。
背景技术
自适应光学波前控制方法是自适应光学技术的一个重要组成部分。目前的自适应光学波前控制主要采用有:基于哈特曼波前传感器的控制方法、基于相位差(PD)远场的波前控制方法,以及基于四棱锥和曲率传感器等波前控制方法。
与其他的波前测量方法相比,相位差控制方法具有更高的光能利用率,在天文、生物医学等需要较高光能量利用率的领域具有比较大的优势。基于单远场的自适应光学波前控制方法(无波前传感器自适应光学控制算法)在近些年有很大的发展,如随机梯度下降控制算法,遗传算法,蚁群算法等在自适应光学的应用,进一步拓宽了自适应光学在弱光领域的应用。
然而,作为光能利用率最高的无波前传感器自适应光学控制算法在进行闭环控制过程中存在的主要困难是:闭环收敛速度慢,虽然近些年来提出的各种优化的无波前控制算法已经在不同程度的提高了闭环的收敛速度,如黄林海等在文献《Wavefrontsensorless adaptive optics:a general model-based approach》提出了一种高效的基于远场的波前控制算法,但是,与传统的基于哈特曼波前探测器的闭环速度相比仍然较低。对于一个由前十项Zernike组成的像差,无波前传感器自适应光学控制算法最少需要探测11次才能收敛(而哈特曼理论上只需1次),运算用时长,因而无法运用到诸如天文等快速动态变化的领域中。
发明内容
本发明要解决的技术问题是:针对无波前传感器自适应光学控制算法中收敛速度慢的问题,提供一种快速的闭环收敛算法。本方法直接利用成像远场数据进行闭环控制,无需增加额外的远场信息。
本发明采用的技术方案为:基于远场的高速自适应光学闭环控制方法,该方法利用的装置由远场图像采集模块、波前复原模块、高压放大器和波前校正器模块组成。该方法工作过程是:远场图像采集模块采集第1帧远场光强变化图像,并将采集图像保存在计算机中;高压放大器驱动波前校正器产生一个特定的波前分布,使远场光强分布产生变化,远场图像采集模块采集第2帧远场光强变化图像,根据上述采集的两帧远场光强变化图像,波前复原模块利用相位差波前复原算法迭代产出校正波前Φ1,并驱动高压放大器使波前校正器模块形成与校正波前Φ1互补的面型分布,最终得到经过第1次校正后的远场光强分布,远场图像采集模块记录第3帧远场光强分布,利用当前帧远场光强分布、上一帧的远场光强分布以及互补的波前Φ1,结合相位差波前复原算法迭代产出校正波前Φ2,并驱动高压放大器使波前校正器模块形成与校正波前Φ2互补的面型分布,得到经过第2次校正后的远场光强分布,重复上述校正步骤,持续对动态变化的像差进行实时的校正。
进一步地,所述的一个特定的波前分布指利用变形镜的单个或多个影响函数组成的波前分布。
进一步地,所述的相位差波前复原算法,具体过程可以描述如下:
第一步,首先构造入射光瞳处的广义光瞳函数Pk(r)=p(r)*exp[i(Ψk(r)+Φk(r))],其中,r表示瞳面坐标,Ψk(r)表示光瞳处待测畸变波前的初始相位分布(一般假定为0),Φk(r)表示光瞳处波前校正器产生的第k帧相位分布,p(r)为入射光瞳的孔径函数,孔径内值为1,其余为0;
第二步,对光瞳的复振幅分布进行傅里叶变换计算出焦面光场的复振幅分布Uk=|Uk(w)|*exp[i*Ωk],w为像面坐标,|Uk(w)|表示第k次迭代中当前帧的光场振幅分布,Ωk为焦面光场相位分布;
第三步,根据远场图像采集模块(1)记录到的当前帧光强分布Ik(w),计算其平方根值sqrt[Ik(w)],取代第二步傅里叶变换得到的振幅分布,组成新的焦面光场复振幅Ek(k)=sqrt[Ik(w)]*exp[i*Ωk];
第四步,对Ek(k)进行傅里叶逆变换,可计算得到光瞳处复振幅分布Qk=qk*exp[i*(Θk+Φk)],qk表示光瞳函数的振幅分布,Θk表示光瞳处波前畸变的相位分布;
第五步,将上一步计算得到的光瞳处复振幅的振幅分布修改为均匀分布的孔径函数p(r),同时将位相中波前校正器产生的像差Φk移除,可组成一个新的光瞳函数Pk N(r)=p(r)*exp[i*(Θk)];
第六步,对新的光瞳函数Pk N(r)进行傅里叶变换,计算出新的光场振幅分布Uk N(w)=|Uk N(w)|*exp[i*Ωk N];
第七步,将上一帧采集到的光强分布I k-1(w)的平方根值取代光场振幅分布|Uk N(w)|,组成新的复振幅Ek N(k)=sqrt[Ik-1(w)]*exp[i*Ωk N];
第八步,对新的复振幅Ek N(k)进行傅里叶变换,可计算得到光瞳处复振幅分布Qk N=qk N*exp[i*Θk N];按照上述步骤完成一次迭代后,比较前后两次迭代位相Θk N的差异的均方根值,如果两次位相差异小于目标函数E则迭代结束,并输出最终位相Θk -Φk(r),其中,E为预先给定的大于零实数,E一般小于0.01。
本发明的原理是:
a、本发明采用了与相位差相似的波前探测方法,波前探测过程中均采用了两个远场光强分布数据,不同的是,相位差方法需要同时采集两个远场,本发明采用的远场为前后两帧的远场,并且相位差方法中使用的两个远场具有固定的相位差,而本发明使用的前后两帧远场不具有固定的相位差,相位差随闭环过程变化;
b、本发明工作的前提是前后两帧动态像差之间存在较大相关性,即在高速采样条件下,两帧像差之间是连续变化的,并且这两帧像差之间的差异较波前校正器产生的像差小;
因此,在上述两个前提情况下,通过前后两帧的信息连续闭环控制逐步消除光路中的像差,包括***静态像差和动态像差。
本发明与现有技术相比有如下优点:
(1)本发明采用单远场光强分布进行闭环控制,无需增加额外的远场信息,如离焦面上的远场光斑分布,可以简化光路结构,同时提高光能利用率;
(2)本发明可以根据前后两帧远场确定入射波前的位相,特别是经过次数较少的自适应光学校正或输入像差较小的情况下,可以准确的获得像差数据,从而实现快速的对像差进行准确校正,因此,本发明较其他迭代控制算法具有更快的收敛速度,实时性更强。
附图说明
图1为基于远场的高速自适应光学闭环控制方法原理示意图;
图2为开环图样;
图3为经过10次闭环控制结果。
具体实施方式
下面结合附图以及具体实施方式进一步说明本发明。
如图1所示,本发明基于远场的高速自适应光学闭环控制方法,包括:该方法利用的装置由远场图像采集模块1、波前复原模块2、高压放大器3和波前校正器模块4组成。其中,远场图像采集模块1一般采用CCD等图像记录设备,波前复原模块2指高性能计算机或者信号处理电路,高压放大器3和波前校正器模块4指用于改变光路中波前分布规律的变形镜及其响应的驱动装置。***准备完成后,远场图像采集模块1采集第1帧远场光强变化图像,并将采集图像保存在计算机中;高压放大器3驱动波前校正器模块4产生一个近似离焦的波前分布ΔΦ(不需要精确的生成离焦面型),使远场光强分布产生变化,远场图像采集模块1采集第2帧远场光强变化图像,根据上述采集的两帧远场光强变化图像和波前分布ΔΦ,波前复原模块2利用相位差波前复原算法迭代产出校正波前Φ1,并驱动高压放大器3使波前校正器模块4形成与校正波前Φ1互相补偿的面型分布,最终得到经过第1次校正后的远场光强分布,远场图像采集模块1记录第3帧远场光强分布,利用当前帧远场光强分布、上一帧的远场光强分布以及互补的波前Φ1,结合相位差波前复原算法迭代产出校正波前Φ2,并驱动高压放大器3使波前校正器模块4形成与校正波前Φ2互相补偿的面型分布,得到经过第2次校正后的远场光强分布,重复上述校正步骤,持续对动态变化的像差进行实时的校正。
所述的相位差波前复原算法,具体过程可以描述如下:
第一步,首先构造入射光瞳处的广义光瞳函数Pk(r)=p(r)*exp[i(Ψk(r)+Φk(r))],其中,r表示瞳面坐标,Ψk(r)表示光瞳处待测畸变波前的初始相位分布(一般假定为0),Φk(r)表示光瞳处波前校正器产生的第k帧相位分布,p(r)为入射光瞳的孔径函数,孔径内值为1,其余为0;
第二步,对光瞳的复振幅分布进行傅里叶变换计算出焦面光场的复振幅分布Uk=|Uk(w)|*exp[i*Ωk],w为像面坐标,|Uk(w)|表示第k次迭代中当前帧的光场振幅分布,Ωk为焦面光场相位分布;
第三步,根据远场图像采集模块(1)记录到的当前帧光强分布Ik(w),计算其平方根值sqrt[Ik(w)],取代第二步傅里叶变换得到的振幅分布,组成新的焦面光场复振幅Ek(k)=sqrt[Ik(w)]*exp[i*Ωk];
第四步,对Ek(k)进行傅里叶逆变换,可计算得到光瞳处复振幅分布Qk=qk*exp[i*(Θk+Φk)],qk表示光瞳函数的振幅分布,Θk表示光瞳处波前畸变的相位分布;
第五步,将上一步计算得到的光瞳处复振幅的振幅分布修改为均匀分布的孔径函数p(r),同时将位相中波前校正器产生的像差Φk移除,可组成一个新的光瞳函数Pk N(r)=p(r)*exp[i*(Θk)];
第六步,对新的光瞳函数Pk N(r)进行傅里叶变换,计算出新的光场振幅分布Uk N(w)=|Uk N(w)|*exp[i*Ωk N];
第七步,将上一帧采集到的光强分布I k-1(w)的平方根值取代光场振幅分布|Uk N(w)|,组成新的复振幅Ek N(k)=sqrt[Ik-1(w)]*exp[i*Ωk N];
第八步,对新的复振幅Ek N(k)进行傅里叶变换,可计算得到光瞳处复振幅分布Qk N=qk N*exp[i*Θk N];按照上述步骤完成一次迭代后,比较前后两次迭代位相Θk N的差异的均方根值,如果两次位相差异小于目标函数E则迭代结束,并输出最终位相Θk -Φk(r),其中,E为预先给定的大于零实数,E一般小于0.01。
Claims (2)
1.基于远场的高速自适应光学闭环控制方法,其特征在于:该方法利用的装置由远场图像采集模块(1)、波前复原模块(2),高压放大器(3)和波前校正器模块(4)组成,该方法工作过程是:远场图像采集模块(1)采集第1帧远场光强变化图像,并将采集图像保存在计算机中;高压放大器(3)驱动波前校正器产生利用变形镜的单个或多个影响函数组成的波前分布,使远场光强分布产生变化,远场图像采集模块(1)采集第2帧远场光强变化图像,根据上述采集的两帧远场光强变化图像,波前复原模块(2)利用相位差波前复原算法迭代生成校正波前Φ1,并驱动高压放大器(3)使波前校正器模块(4)形成与校正波前Φ1互补的面型分布,最终得到经过第1次校正后的远场光强分布,远场图像采集模块(1)记录第3帧远场光强分布,利用当前帧远场光强分布和上一帧的远场光强分布以及互补的波前Φ1,结合相位差波前复原算法迭代产出校正波前Φ2,并驱动高压放大器(3)使波前校正器模块(4)形成与校正波前Φ2互补的面型分布,得到经过第2次校正后的远场光强分布,重复上述步骤,持续对动态变化的像差进行实时的校正。
2.根据权利要求1所述的基于远场的高速自适应光学闭环控制方法,其特征在于:所述的相位差波前复原算法,具体过程可以描述如下:
第一步,首先构造入射光瞳处的广义光瞳函数Pk(r)=p(r)*exp[i(Ψk(r)+Φk(r))],其中,r表示瞳面坐标,Ψk(r)表示光瞳处待测畸变波前的初始相位分布,Φk(r)表示光瞳处波前校正器产生的第k帧相位分布,p(r)为入射光瞳的孔径函数,孔径内值为1,其余为0;
第二步,对光瞳的复振幅分布进行傅里叶变换计算出焦面光场的复振幅分布Uk=|Uk(w)|*exp[i*Ωk],w为像面坐标,|Uk(w)|表示第k次迭代中当前帧的光场振幅分布,Ωk为焦面光场相位分布;
第三步,根据远场图像采集模块(1)记录到的当前帧光强分布Ik(w),计算其平方根值sqrt[Ik(w)],取代第二步傅里叶变换得到的振幅分布,组成新的焦面光场复振幅Ek(k)=sqrt[Ik(w)]*exp[i*Ωk];
第四步,对Ek(k)进行傅里叶逆变换,可计算得到光瞳处复振幅分布Qk=qk*exp[i*(Θk+Φk)],qk表示光瞳函数的振幅分布,Θk表示光瞳处波前畸变的相位分布;
第五步,将上一步计算得到的光瞳处复振幅的振幅分布修改为均匀分布的孔径函数p(r),同时将位相中波前校正器产生的像差Φk移除,可组成一个新的光瞳函数Pk N(r)=p(r)*exp[i*(Θk)];
第六步,对新的光瞳函数Pk N(r)进行傅里叶变换,计算出新的光场振幅分布Uk N(w)=|Uk N(w)|*exp[i*Ωk N];
第七步,将上一帧采集到的光强分布Ik-1(w)的平方根值取代光场振幅分布|Uk N(w)|,组成新的复振幅Ek N(k)=sqrt[Ik-1(w)]*exp[i*Ωk N];
第八步,对新的复振幅Ek N(k)进行傅里叶变换,可计算得到光瞳处复振幅分布Qk N=qk N*exp[i*Θk N];按照上述步骤完成一次迭代后,比较前后两次迭代位相Θk N的差异的均方根值,如果两次位相差异小于目标函数E则迭代结束,并输出最终位相Θk -Φk(r),其中,E为预先给定的大于零实数,E小于0.01。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011190608.7A CN112099229B (zh) | 2020-10-30 | 2020-10-30 | 基于远场的高速自适应光学闭环控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011190608.7A CN112099229B (zh) | 2020-10-30 | 2020-10-30 | 基于远场的高速自适应光学闭环控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112099229A CN112099229A (zh) | 2020-12-18 |
CN112099229B true CN112099229B (zh) | 2022-06-28 |
Family
ID=73785788
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011190608.7A Active CN112099229B (zh) | 2020-10-30 | 2020-10-30 | 基于远场的高速自适应光学闭环控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112099229B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114200664B (zh) * | 2021-11-11 | 2023-09-05 | 常州北邮新一代信息技术研究院有限公司 | 基于改进相位差算法的自适应光学*** |
CN114137814B (zh) * | 2021-11-30 | 2022-08-16 | 武汉大学 | 一种实现独立全息图像复用的超表面器件及其构建方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104215339A (zh) * | 2014-09-14 | 2014-12-17 | 中国科学院光电技术研究所 | 一种基于连续远场的波前复原***及方法 |
CN108646407A (zh) * | 2018-05-08 | 2018-10-12 | 宜宾学院 | 一种无波前感知自适应光学成像***及其成像方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6721510B2 (en) * | 2001-06-26 | 2004-04-13 | Aoptix Technologies, Inc. | Atmospheric optical data transmission system |
US7406263B2 (en) * | 2002-10-17 | 2008-07-29 | Aoptix Technologies | Combined wavefront sensor and data detector for a free space optical communications system with adaptive optics |
CN102095504B (zh) * | 2010-12-07 | 2012-07-11 | 四川大学 | 基于空间相位调制的环形共光路点衍射干涉仪 |
CN102889935A (zh) * | 2012-09-14 | 2013-01-23 | 中国科学院光电技术研究所 | 基于相位差法的自适应光学***近场波前传感器标定装置及标定方法 |
-
2020
- 2020-10-30 CN CN202011190608.7A patent/CN112099229B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104215339A (zh) * | 2014-09-14 | 2014-12-17 | 中国科学院光电技术研究所 | 一种基于连续远场的波前复原***及方法 |
CN108646407A (zh) * | 2018-05-08 | 2018-10-12 | 宜宾学院 | 一种无波前感知自适应光学成像***及其成像方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112099229A (zh) | 2020-12-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109031654B (zh) | 一种基于卷积神经网络的自适应光学校正方法及*** | |
CN107843982B (zh) | 基于实时相位差异技术的无波前探测自适应光学*** | |
CN112099229B (zh) | 基于远场的高速自适应光学闭环控制方法 | |
US7274442B2 (en) | Closed-loop wavefront sensor using field programmable gate array | |
CN110207835A (zh) | 一种基于离焦图像训练的波前校正方法 | |
US11368608B2 (en) | Compressed sensing based object imaging system and imaging method therefor | |
CN106338823B (zh) | 一种基于混合焦距菲涅耳波带片的相位反演方法 | |
CN106526839B (zh) | 一种基于模式的同步无波前自适应光学*** | |
CN110675451B (zh) | 基于相空间光学的数字自适应校正方法及*** | |
CN105933060B (zh) | 一种基于动力学反馈神经网络的波前重构方法 | |
CN116046183A (zh) | 基于图像增强的夏克-哈特曼波前传感器波前复原方法 | |
CN113837947B (zh) | 一种获得光学相干层析大焦深图像的处理方法 | |
CN109581849B (zh) | 一种同轴全息重建方法及*** | |
CN101162294B (zh) | 用于惯性约束聚变装置全光路像差的测量和校正方法 | |
CN108646406B (zh) | 一种基于棱锥波前传感器本征模式控制的自适应光学装置 | |
CN105044906B (zh) | 一种基于图像信息的快速扩展目标成像校正方法 | |
CN112484866B (zh) | 一种基于夏克-哈特曼波前传感器的波前复原方法 | |
CN104596651B (zh) | 一种基于四象限二元相位调制的相位反演方法 | |
CN111695676B (zh) | 一种基于生成对抗网络的波前复原方法及*** | |
CN116400495B (zh) | 一种基于run优化算法的波前校正*** | |
CN207882564U (zh) | 基于实时相位差异技术的无波前探测自适应光学*** | |
CN104215339B (zh) | 一种基于连续远场的波前复原***及方法 | |
CN102621688B (zh) | 基于差分传感器的自适应光学*** | |
US12025793B1 (en) | Wavefront correction system based on RUN optimization algorithm | |
CN112179503A (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 |