CN106324854B - 一种基于二元方形衍射元件的相位反演方法 - Google Patents

一种基于二元方形衍射元件的相位反演方法 Download PDF

Info

Publication number
CN106324854B
CN106324854B CN201610952805.5A CN201610952805A CN106324854B CN 106324854 B CN106324854 B CN 106324854B CN 201610952805 A CN201610952805 A CN 201610952805A CN 106324854 B CN106324854 B CN 106324854B
Authority
CN
China
Prior art keywords
far
phase
field
binary
diffraction element
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
Application number
CN201610952805.5A
Other languages
English (en)
Other versions
CN106324854A (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.)
Institute of Optics and Electronics of CAS
Original Assignee
Institute of Optics and Electronics of CAS
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 Institute of Optics and Electronics of CAS filed Critical Institute of Optics and Electronics of CAS
Priority to CN201610952805.5A priority Critical patent/CN106324854B/zh
Publication of CN106324854A publication Critical patent/CN106324854A/zh
Application granted granted Critical
Publication of CN106324854B publication Critical patent/CN106324854B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/42—Diffraction optics, i.e. systems including a diffractive element being designed for providing a diffractive effect
    • G02B27/4205—Diffraction optics, i.e. systems including a diffractive element being designed for providing a diffractive effect having a diffractive optical element [DOE] contributing to image formation, e.g. whereby modulation transfer function MTF or optical aberrations are relevant
    • 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/0012—Optical design, e.g. procedures, algorithms, optimisation routines

Abstract

本发明提供了一种基于二元方形衍射元件的相位反演方法，在复原算法中，描述光强分布的常规方法是x、y方向的二维矩阵，所以二元圆形衍射元件建模时其圆周边沿有很多方形锯齿，很难建立准确的数学模型，从而影响了其复原波前的精度，而二元方形衍射元件则很好地解决了上述问题，且二元方形衍射元件本身对输入波面进行了调制，较传统单帧波前反演的基于玻璃透镜的GS算法，波前复原精度有较大提高。本发明元件简单，轻便，适合大孔径及超大孔径光波相位探测，可实现基于单帧远场图像的准确近场相位反演。

Description

一种基于二元方形衍射元件的相位反演方法
技术领域
本发明涉及一种基于二元方形衍射元件的相位反演方法，属于波前传感技术领域。
背景技术
相位反演技术是波前传感技术的重要分支，特点是直接利用采集到的光束远场光斑图像信息，通过衍射光学理论，反演推算出光束的近场波前相位分布信息。
基于单帧图像的相位反演光学实现结构非常简单，但其主要缺欠是由于点扩散函数不唯一，波前复原存在多解性，因而基于更多远场信息的相位反演方法应运而生。基于多帧图像的相位反演技术通过获取更多的远场图像作为已知信息，增加相位反演限定条件，从而解决解的唯一性问题。然而，同时获取多帧远场图像意味着光学实现结构复杂度的增加和可靠性的降低。因此，若能通过新的技术手段和数学方法解决多解问题，利用单帧远场图像实现准确的相位反演，那么光学实现结构能够十分简洁，并且无需分光元件，光能利用率和信号信噪比高，是较为理想的相位反演实施方式。
发明内容
(一)要解决的技术问题
鉴于上述技术问题，本发明提供了一种基于二元方形衍射元件的相位反演方法，在复原算法中，描述光强分布的常规方法是x、y方向的二维矩阵，所以二元圆形衍射元件建模时其圆周边沿有很多方形锯齿，很难建立准确的数学模型，从而影响了其复原波前的精度，而二元方形衍射元件则很好地解决了上述问题，且二元方形衍射元件本身对输入波面进行了调制，较传统单帧波前反演的GS算法波前复原精度有较大提高。本发明元件简单，轻便，适合大孔径及超大孔径波前探测，可实现基于单帧远场图像的准确近场相位反演。
(二)技术方案
根据本发明的一个方面，提供了一种基于二元方形衍射元件的相位反演方法，由于二元衍射元件相较传统的玻璃透镜更容易特制，任何图形布局都可以在二元衍射元件上得到实现，可以在远场形成了特殊的光斑形状，使原光波在远近场有唯一的点扩散函数。光学实现结构包括二元方形衍射元件1和光电探测器2，所述的二元方形衍射元件1由0、pi的相位间隔构成，根据设计不同，可以形成独特的远场光强分布，光电探测器用于探测调制后光波所成的远场光斑图像，进而通过以下步骤处理远场光斑图像，实现对入射波前相位复原：
步骤1：已知入射光束近场强度分布Inear和对应的远场光斑图像强度分布Ifar，并设定相位反演方法中近场波前相位分布的初始值为0；
步骤2：计算二元方形衍射元件的远场复振幅分布：
式中Afar为计算远场光波振幅分布，为计算的远场光波相位分布；
步骤3：对比计算的远场光波振幅|Efar|和实际的远场光波振幅分布计算表征两者之间差异的评价指标：
若SSE小于设定的判定标准，表明本次计算所用近场光波复振幅与实际入射光波拥有一致的远场光强分布，则当前的近场波前即为实际的光束近场相位分布，作为相位反演结果输出，基于二元方形衍射元件的相位反演方法结束；若SSE大于设定的判定标准，则方法继续；
步骤4：将远场实际光强分布Ifar平方根作为远场光波振幅，变换后远场光波复振幅为：
利用变化后远场光波复振幅E′far，计算逆向衍射后对应二元方形衍射元件的近场光波复振幅：
式中Anear为计算近场光波振幅分布，为计算的近场调制后光波相位分布；
步骤5：将Inear的平方根作为近场光波振幅代替Anear，从而构成新的近场调制后光波复振幅，重新进入复原方法步骤2，开始新一轮的迭代计算，直至某次迭代复原运算的步骤3满足判定标准，则基于二元方形衍射元件的相位反演方法结束，输出反演的光束近场相位分布结果。
(三)有益效果
本发明一种基于二元方形衍射元件的相位反演方法。简单易于实现，轻薄可靠，在复原算法中，描述光强分布的常规方法是x、v方向的二维矩阵，所以二元圆形衍射元件建模时其圆周边沿有很多方形锯齿，很难建立准确的数学模型，从而影响了其复原波前的精度，而二元方形衍射元件则很好地解决了上述问题，且二元方形衍射元件本身对输入波面进行了调制，较传统单帧波前反演的GS算法波前复原精度有较大提高。对大孔径波前探测及单帧相位反演探测拥有重要意义。
附图说明
图1为本发明方法原理流程图；
图2为本发明方法光路图；
图3为传统GS算法示意图，其中，图3(a)为初始波前面型，图3(b)为算法复原面型，图3(c)为波前残余误差，图3(d)为算法收敛曲线；
图4为二元方形衍射元件方形结构示意图；
图5为实施例一中本发明方法，其中，图5(a)为初始波前面型，图5(b)为算法复原面型，图5(c)为波前残余误差，图5(d)为算法收敛曲线；
图6为二元方形衍射元件四区域边长不同的方形结构示意图；
图7为实施例二中本发明方法，其中，图7(a)为初始波前面型，图7(b)为算法复原面型，图7(c)为波前残余误差，图7(d)为算法收敛曲线；
图中附图标记含义为：1为二元方形衍射元件，2为光电探测器件。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白，以下结合具体实施例，并参照附图，对本发明进一步详细说明。
本发明是基于二元方形衍射元件的相位反演方法，核心思想是通过二元方形衍射元件对波前进行调制，在远场产生特殊的光场分布，增加了算法的约束条件，元件简单，轻便，适合大孔径及超大孔径光波探测，可实现基于单帧远场图像的准确近场相位反演。
实施例一
在本发明的示例性实施例中，提供了一种基于方形结构二元衍射元件的相位反演方法。
请参照图1～5，本实施例方法包括：
图2是系统光路图，实施例中待测的入射光波波前相位是由65阶Zernike多项式组成的随机波前(PV＝3.7415rad，rms＝0.8082rad)，图4为方形结构二元方形衍射元件结构图。
光学实现结构包括二元方形衍射元件1和光电探测器2，二元方形衍射元件1，由0、pi的相位间隔构成，根据设计不同，可以形成独特的远场光强分布，光电探测器2用于探测调制后光波所成的远场光斑图像，进而通过以下步骤处理远场光斑图像，实现对入射波前相位复原：
步骤1：已知入射光束近场强度分布Inear和对应的远场光斑图像强度分布Ifar，并设定相位反演方法中近场波前相位分布的初始值为0。
步骤2：计算二元方形衍射元件的远场复振幅分布：
式中Afar为计算远场光波振幅分布，为计算的远场光波相位分布；
步骤3：对比计算的远场光波振幅|Efar|和实际的远场光波振幅分布计算表征两者之间差异的评价指标：
若SSE小于设定的判定标准，表明本次计算所用近场光波复振幅与实际入射光波拥有一致的远场光强分布，则当前的近场波前即为实际的光束近场相位分布，作为相位反演结果输出，基于二元方形衍射元件的相位反演方法结束；若SSE大于设定的判定标准，则方法继续；
步骤4：将远场实际光强分布Ifar平方根作为远场光波振幅，变换后远场光波复振幅为：
利用变化后远场光波复振幅E′far，计算逆向衍射后对应二元方形衍射元件的近场光波复振幅：
式中Anear为计算近场光波振幅分布，为计算的近场调制后光波相位分布；
步骤5：将Inear的平方根作为近场光波振幅代替Anear，从而构成新的近场调制后光波复振幅，重新进入复原方法步骤2，开始新一轮的迭代计算，直至50次迭代复原运算后满足步骤4判定标准，则基于二元方形衍射元件的相位反演方法结束，输出反演的光束近场相位分布结果。
将传统单帧GS算法迭代200次后仍远达不到步骤4判定标准，所以同样迭代50次，与本发明方法进行对比。图3为传统GS算法复原结果，图3(a)为初始波前面型(PV＝3.7415rad，rms＝0.8082rad)，图3(b)为算法复原面型(PV＝3.721rad，rms＝0.8062rad)，图3(c)为波前残余误差(PV＝1.2482rad，rms＝0.1267rad)，图3(d)为算法收敛曲线；图5为实施例一中本发明方法复原结果，图5(a)为初始波前面型(PV＝3.7415rad，rms＝0.8082rad)，图5(b)为算法复原面型(PV＝3.7785rad，rms＝0.7974rad)，图5(c)为波前残余误差(PV＝0.2332rad，rms＝0.0291rad)，图5(d)为算法收敛曲线。可以看出本发明方法波前复原细节、精度都较传统方法有较大提高。收敛更快，且不容易停滞，相同迭代次数下，本发明方法相较传统方法精度提高近5倍。
本实施例中，基于二元方形衍射元件的相位反演方法，核心思想是通过二元方形衍射元件对波前进行调制，在远场产生特殊的光场分布，能对单帧远场图像进行较准确的相位反演。
实施例二
在本发明的示例性实施例中，提供了一种基于四区域边长不同的方形结构二元方形衍射元件的相位反演方法。
请参照图1，图2，图3，图6，图7，本实施例方法包括：
图2是系统光路图，实施例中待测的入射光波波前相位是由65阶Zernike多项式组成的随机波前(PV＝3.7415rad，rms＝0.8082rad)，图6为四区域边长不同的方形结构二元方形衍射元件结构图。
步骤1：已知入射光束近场强度分布Inear和对应的远场光斑图像强度分布Ifar，并设定相位反演方法中近场波前相位分布的初始值为0。
步骤2：计算二元方形衍射元件在焦面附近的远场复振幅分布：
式中Afar为计算远场光波振幅分布，为计算的远场光波相位分布；
步骤3：对比计算的远场光波振幅|Efar|和远场实际光强分布Ifar计算表征两者之间差异的评价指标：
若SSE小于设定的判定标准，表明本次计算所用近场光波复振幅与实际入射光波拥有一致的远场光强分布，则当前的近场波前即为实际的光束近场相位分布，作为相位反演结果输出，基于二元方形衍射元件的相位反演方法结束；若SSE大于设定的判定标准，则方法继续；
步骤4：将远场实际光强分布Ifar平方根作为远场光波振幅，变换后远场光波复振幅为：
利用变化后远场光波复振幅E′far，计算逆向衍射后对应二元方形衍射元件的近场光波复振幅：
式中Anear为计算近场光波振幅分布，为计算的近场调制后光波相位分布；
步骤5：将Inear的平方根作为近场光波振幅代替Anear，从而构成新的近场调制后光波复振幅，重新进入复原方法步骤2，开始新一轮的迭代计算，直至50次迭代复原运算后满足步骤3判定标准，则基于二元方形衍射元件的相位反演方法结束，输出反演的光束近场相位分布结果。
将传统单帧GS算法迭代200次后仍远达不到步骤4判定标准，所以同样迭代50次，与本发明方法进行对比。图3为传统GS算法复原结果，图3(a)为初始波前面型(PV＝3.7415rad，rms＝0.8082rad)，图3(b)为算法复原面型(PV＝3.721rad，rms＝0.8062rad)，图3(c)为波前残余误差(PV＝1.2482rad，rms＝0.1267rad)，图3(d)为算法收敛曲线；图7为实施例二中本发明方法复原结果，图7(a)为初始波前面型(PV＝3.7415rad，rms＝0.8082rad)，图7(b)为算法复原面型(PV＝3.7519rad，rms＝0.8038rad)，图7(c)为波前残余误差(PV＝0.0934rad，rms＝0.0102rad)，图7(d)为算法收敛曲线。可以看出本发明方法波前复原细节、精度都较传统方法有较大提高。收敛更快，且不容易停滞，相同迭代次数下，本发明方法相较传统方法精度提高十几倍。
本实施例中，基于四区域边长不同的方形结构二元方形衍射元件的相位反演方法，核心思想是通过不同区域的衍射元件排布不同对波前进一步调制，在远场产生更加特殊的光场分布，能对单帧远场图像进行较准确的相位反演。
至此，已经结合附图对本发明实施例进行了详细描述。依据以上描述，本领域技术人员应当对本发明基于二元方形衍射元件的相位反演方法有了清楚的认识。
以上所述，仅为本发明中的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉该技术的人在本发明所揭露的技术范围内，可理解想到的变换或替换，都应涵盖在本发明的包含范围之中。

Claims (1)

1.一种基于四区域边长不同的方形结构二元方形衍射元件的相位反演方法，其特征在于：该方法包括：待测的入射光波波前相位是由65阶Zernike多项式组成的随机波前，PV＝3.7415rad，rms＝0.8082rad，该方法利用四区域边长不同的方形结构二元方形衍射元件，
步骤1：已知入射光束近场强度分布Inear和对应的远场光斑图像强度分布Ifar，并设定相位反演方法中近场波前相位分布的初始值为0；
步骤2：计算二元方形衍射元件在焦面附近的远场复振幅分布：
式中Afar为计算远场光波振幅分布，为计算的远场光波相位分布；
步骤3：对比计算的远场光波振幅|Efar|和远场实际光强分布Ifar计算表征两者之间差异的评价指标：
若SSE小于设定的判定标准，表明本次计算所用近场光波复振幅与实际入射光波拥有一致的远场光强分布，则当前的近场波前即为实际的光束近场相位分布，作为相位反演结果输出，基于二元方形衍射元件的相位反演方法结束；若SSE大于设定的判定标准，则方法继续；
步骤4：将远场实际光强分布Ifar平方根作为远场光波振幅，变换后远场光波复振幅为：
利用变化后远场光波复振幅E′far，计算逆向衍射后对应二元方形衍射元件的近场光波复振幅：
式中Anear为计算近场光波振幅分布，为计算的近场调制后光波相位分布；
步骤5：将Inear的平方根作为近场光波振幅代替Anear，从而构成新的近场调制后光波复振幅，重新进入复原方法步骤2，开始新一轮的迭代计算，直至50次迭代复原运算后满足步骤3判定标准，则基于二元方形衍射元件的相位反演方法结束，输出反演的光束近场相位分布结果；
将传统单帧GS算法迭代200次后仍远达不到步骤4判定标准，所以同样迭代50次，可以看出该方法波前复原细节、精度都较传统方法有较大提高，收敛更快，且不容易停滞，相同迭代次数下，该方法相较传统方法精度提高十几倍。
CN201610952805.5A 2016-10-27 2016-10-27 一种基于二元方形衍射元件的相位反演方法 Active CN106324854B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610952805.5A CN106324854B (zh) 2016-10-27 2016-10-27 一种基于二元方形衍射元件的相位反演方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610952805.5A CN106324854B (zh) 2016-10-27 2016-10-27 一种基于二元方形衍射元件的相位反演方法

Publications (2)

Publication Number Publication Date
CN106324854A CN106324854A (zh) 2017-01-11
CN106324854B true CN106324854B (zh) 2019-06-25

Family

ID=57818947

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610952805.5A Active CN106324854B (zh) 2016-10-27 2016-10-27 一种基于二元方形衍射元件的相位反演方法

Country Status (1)

Country Link
CN (1) CN106324854B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111240009B (zh) * 2019-12-31 2020-12-29 嘉兴驭光光电科技有限公司 可用于投射斜线的衍射光学元件、投射装置及其设计方法
US20210232093A1 (en) * 2020-01-27 2021-07-29 Texas Instruments Incorporated Projector with phase hologram modulator
CN113625379B (zh) * 2020-05-06 2023-09-26 中国科学院微电子研究所 一种菲涅尔波带片的设计方法、制作方法和设计装置
CN111854981B (zh) * 2020-07-10 2022-09-20 中国科学院光电技术研究所 一种基于单帧焦面光强图像的深度学习波前复原方法

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3854791A (en) * 1971-02-25 1974-12-17 Hitachi Ltd Holographic memory with random phase illumination hologram
US20040165268A1 (en) * 2001-07-16 2004-08-26 Jari Turunen Diffractive shaping of the intensity distribution of a spatially partially coherent light beam
JP2008034080A (ja) * 2006-06-29 2008-02-14 Sharp Corp 光集積ユニットおよび光ピックアップ装置
CN103676498B (zh) * 2013-11-18 2017-04-05 中国科学院上海光学精密机械研究所 光刻机光瞳整形单元结构及其衍射光学元件设计方法
CN103777472B (zh) * 2014-01-23 2015-07-29 中国科学院上海光学精密机械研究所 用于准分子激光光束整形的衍射光学元件设计方法
CN104215339B (zh) * 2014-09-14 2017-02-15 中国科学院光电技术研究所 一种基于连续远场的波前复原系统及方法
CN104280880A (zh) * 2014-09-27 2015-01-14 郑敏 基于干涉法的衍射光学元件设计方法
CN104596650B (zh) * 2015-01-27 2017-07-11 中国科学院光电技术研究所 一种哈特曼波前传感器超分辨波前复原方法
CN104596651B (zh) * 2015-02-02 2017-07-07 中国科学院光电技术研究所 一种基于四象限二元相位调制的相位反演方法

Also Published As

Publication number Publication date
CN106324854A (zh) 2017-01-11

Similar Documents

Publication Publication Date Title
CN106324854B (zh) 一种基于二元方形衍射元件的相位反演方法
CN106338823B (zh) 一种基于混合焦距菲涅耳波带片的相位反演方法
CN103323396B (zh) 一种基于符合测量的二维压缩鬼成像系统及方法
CN102183214B (zh) 一种大口径非球面镜结构光检测方法
CN102288305B (zh) 一种自适应光学系统波前传感器及其探测方法
Zammit et al. Extended depth-of-field imaging and ranging in a snapshot
CN102967380B (zh) 一种基于单元光敏探测器阵列的哈特曼波前传感器
Trouvé et al. Passive depth estimation using chromatic aberration and a depth from defocus approach
Starikov et al. Wavefront reconstruction of an optical vortex by a Hartmann-Shack sensor
CN106526839B (zh) 一种基于模式的同步无波前自适应光学系统
CN103308189B (zh) 一种基于双压缩符合测量的纠缠成像系统及方法
CN104596650B (zh) 一种哈特曼波前传感器超分辨波前复原方法
Yuan et al. High-accuracy phase demodulation method compatible to closed fringes in a single-frame interferogram based on deep learning
CN104483752A (zh) 一种反射式数字成像系统的设计方法
Yue et al. Co-phasing of the segmented mirror and image retrieval based on phase diversity using a modified algorithm
Védrenne et al. Laser beam complex amplitude measurement by phase diversity
Cheffot et al. Measuring the cophasing state of a segmented mirror with a wavelength sweep and a Zernike phase contrast sensor
CN108303038A (zh) 基于二维光学点阵的反射型面形测量方法和装置
CN104596651B (zh) 一种基于四象限二元相位调制的相位反演方法
Jingshan et al. Sparse ACEKF for phase reconstruction
Lawson Group-delay tracking in optical stellar interferometry with the fast Fourier transform
CN115523839A (zh) 一种基于观测目标的拼接子镜边缘误差高精度实时探测方法
Chen et al. Distortion spot correction and center location base on deep neural network and MBAS in measuring large curvature aspheric optical element
Shi et al. Rapid all-in-focus imaging via physical neural network optical encoding
CN111652372B (zh) 一种基于衍射光学神经网络的波前复原方法及系统

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant