CN105716725B - 一种基于叠层扫描的相位差波前探测和图像复原方法 - Google Patents

一种基于叠层扫描的相位差波前探测和图像复原方法 Download PDF

Info

Publication number
CN105716725B
CN105716725B CN201610136602.9A CN201610136602A CN105716725B CN 105716725 B CN105716725 B CN 105716725B CN 201610136602 A CN201610136602 A CN 201610136602A CN 105716725 B CN105716725 B CN 105716725B
Authority
CN
China
Prior art keywords
phase difference
lamination
scanning
wavefront
subgraph
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
CN201610136602.9A
Other languages
English (en)
Other versions
CN105716725A (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 CN201610136602.9A priority Critical patent/CN105716725B/zh
Publication of CN105716725A publication Critical patent/CN105716725A/zh
Application granted granted Critical
Publication of CN105716725B publication Critical patent/CN105716725B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J9/00Measuring optical phase difference; Determining degree of coherence; Measuring optical wavelength

Landscapes

  • Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Testing Of Optical Devices Or Fibers (AREA)
  • Eye Examination Apparatus (AREA)
  • Lenses (AREA)

Abstract

本发明涉及一种基于叠层扫描的相位差波前探测和图像复原方法,可用于恢复波前畸变,并且能够有效复原受像差影响的模糊图像。本发明将一个小孔径的光阑在成像系统的光瞳平面内以叠层扫描的方式进行移动,并用图像传感器记录下相应的受不同像差影响的子图像。通过使用基于叠层扫描的相位差波前探测和图像复原算法对记录的系列子图像进行处理,可以探测出系统的波前畸变并复原图像。本发明采用小孔径的光阑进行空间扫描以产生含有相位差异的子图像,相对于目前用于相位差法的各种技术,更容易满足奈斯奎特采样频率,像差探测范围大,并且无需离焦光路和各种衍射器件,具有系统紧凑,使用方便等优点。

Description

一种基于叠层扫描的相位差波前探测和图像复原方法
技术领域
本发明涉及波前探测和图像复原的技术领域,特别涉及一种基于叠层扫描的相位差波前探测和图像复原方法。
背景技术
相位差法首先由Gonsalves R A提出,通过在待测波前中添加已知的离焦像差,根据得到的焦面光强分布和多幅位于不同离焦面的光强分布,基于最小二乘,采用GS算法可有效探测波前畸变,同时恢复受像差污染的模糊图像。该技术因结构简单、测量精度高以及对光源无特殊要求,在遥感和波前探测领域得到了广泛的应用,受到了国内外科研工作者的广泛关注和研究。例如美国亚利桑那光学中心、洛克希德-马丁公司下属实验室等国外科研机构相继搭建了实验平台对相位差技术进行研究,力图完善其算法性能并拓展其应用;在国内,中科院光电所自适应光学重点实验室也开展了相关研究,提出了多种改良的相位差波前探测器。
但是,目前相位差法仍然存在多种局限:相位差法本质上是一种基于图像的反演技术,需要根据采集的光强信息来提取波前信息和复原图像,作为先验信息而人为引入的像差势必会使图像丢失部分信息,从而影响波前探测的精度以及复原图像的质量,相位差法对波前畸变的探测能力因此受到制约;在对于分辨率有很高要求的应用领域,例如天文观测,成像系统的口径需要做得很大,常见CCD的像元尺寸很难满足奈斯奎特采样定律,在后端算法处理中往往需要做插值,从而影响相位差算法的复原精度;目前广泛应用的相位差波前探测和成像系统都需要额外的光学元件,常见的如分光镜和反射镜,特殊的如发明专利CN102331303使用的光栅,以及发明专利CN102564612B提出的组合棱镜,这些光学元件的使用可以保证图像采集的实时性、增强对高频像差的探测能力,却增加了系统的体积、复杂度和成本,并引入额外的系统误差以及未知像差,严重影响算法精度。
本发明提出一种基于叠层扫描的相位差新方法及装置,采用小孔径光阑进行空间调制来产生相位差异图像,没有引入额外像差,使原图像得以保留更多细节信息,小孔径光阑的使用,也使一般的CCD像素尺寸都可以满足光学系统要求的奈斯奎特采样频率,从而提高了波前探测的精度和图像复原的质量。另外,本发明无需离焦光路和额外的光学元件,具有光路紧凑、使用便捷、成本低廉等优点。
发明内容
为了克服现有技术存在的问题和实现的复杂程度,本发明提出了一种基于叠层扫描的相位差波前探测和图像复原方法。
本发明采用的技术方案是:一种基于叠层扫描的相位差波前探测和图像复原方法,该方法包含以下步骤:
第一步,将小孔径的光阑在成像系统的光瞳平面内以叠层扫描的方式进行平移,并用图像传感器依次记录下与各扫描位置所对应的受不同像差影响的子图像;
第二步,使用基于叠层扫描的相位差波前探测和图像复原算法对序列子图像进行处理,得到待测的畸变波前和复原的图像,处理算法流程如下:
1)根据极大似然估计或最小二乘优化理论建立目标函数E:
其中,Dk(fx,fy)是调制光阑处在第k个扫描位置时图像传感器采集的子图像的傅里叶变换,Hk(fx,fy)是调制光阑处在第k个扫描位置时对应的光学通道的光学传递函数,是Hk(fx,fy)的共轭K为扫描的次数,fx,fy为像平面频域的坐标。
2)根据傅里叶光学,估计出第k个光学通道的点扩散函数,系统的波前相位分布可以由泽尼克多项式来表征:
其中Zn(u,v)表示第n阶泽尼克多项式,αn表示第n阶泽尼克多项式的系数;
3)目标函数E演变为关于变量αn的函数E(α),采用诸如SPGD,遗传算法,模拟退火,神经网络等最优化算法寻找使目标函数最小的解,即可求得待测的波前畸变
4)根据求得的波前畸变,利用下式可以复原图像R:
其中,该方法利用的装置包括成像系统,空间调制单元,图像传感器和计算机;空间调制单元放置于成像系统的出瞳面上,图像传感器置于成像系统的像平面上,并输入到计算机中;成像目标发出的光波进入到成像系统中,受到空间调制单元的叠层调制后,在图像传感器上成像,该图像存储在计算机中,最后利用算法对所有子图像进行处理得到探测的波前以及恢复的图像。
其中,该方法需要对成像系统的光瞳平面进行叠层扫描的空间调制,用于调制的小孔径光阑在光瞳二维平面内依次平移一定的距离完成扫描,并保证相邻平移位置间光瞳平面上被光阑扫描的区域内存在一定的交叠,交叠部分对应的像差使采集的子图像含有部分相同的相位信息,而非交叠部分对应的像差则在采集的子图像中引入了相位差异。
其中,该方法子图像之间的相位差异隐藏在已知的扫描位置信息里,而不是引入的额外像差,因此,采集的子图像得以保留更多的细节信息,从而提高波前探测的范围和精度。
其中,该方法叠层扫描不仅可以由机械移动平台完成,也可以由数字微镜器件(Digital Micromirror Device,DMD)或者空间光调制器(Spatial Light Modulator,SLM)等快速光学器件来完成,其相应速度也足够保证本发明对准静态像差的探测。
其中,该方法既适用于点光源也适用于扩展目标;既可以用于相干光成像也可以用于非相干光成像。
其中,该方法既可以用来探测波前畸变,也可以用来恢复图像。
本发明相对于现有方法,具有如下优点:
(1)相对于传统相位差法通过引入离焦量来产生不同的光学通道,本发明采用空间调制的方式得到含有相位差异的子图像,使图像得以保留更多细节信息。由于相位差法本质上是一种基于光强图像复原待测波前的间接波前探测技术,因此本发明能够对更大范围的波前畸变进行更为准确的测量,提高了相位差法的波前探测能力。
(2)本发明采用小孔径光阑对光瞳平面进行空间调制,使得一般的CCD像素尺寸都可以满足光学系统要求的奈斯奎特采样频率,保证了相位差算法的有效性,扩大了相位差法的应用领域。
(3)相对于传统相位差波前传感器的光路结构,本发明采用空间调制和相位差法相结合,以较小的时间代价消除了系统对其他光学元件的需求,避免引入额外的相位误差,使整套系统也更紧凑,成本更低廉,使用更方便。
附图说明
图1为基于叠层扫描的相位差波前探测和图像复原新方法的装置和原理示意图,其中,1为成像系统,2为空间调制单元,3为图像传感器,4为计算机;
图2为叠层扫描示意图。
具体实施方式
下面结合附图以及具体实施方式进一步说明本发明。
本发明一种基于叠层扫描的相位差波前探测和图像复原方法,该方法的步骤如下:
(1)第一步,将小孔径的光阑在成像系统的光瞳平面内以叠层扫描的方式进行平移,并用图像传感器依次记录下与各扫描位置所对应的受不同像差影响的子图像。
所述叠层扫描是这样一种调制方式,用于扫描的小孔径光阑在光瞳二维平面内依次平移一定的距离完成调制,并保证相邻平移位置间光瞳平面上被光阑扫描的区域内存在一定的交叠。
叠层扫描不仅可以由机械移动平台完成,在对实时性有一定要求的应用领域,也可以由数字微镜器件(Digital Micromirror Device,DMD)或者空间光调制器(SpatialLight Modulator,SLM)等快速调制光学器件来完成,其相应速度也足够保证本发明对准静态像差的探测。
子图像之间的相位差异隐藏在已知的调制信息里,而不是引入的额外像差:交叠部分对应的像差使采集的子图像含有部分相同的相位信息,而非交叠部分对应的像差则在采集的子图像中引入了相位差异。
(2)第二步,使用基于叠层扫描的相位差波前探测和图像复原算法对子图像进行处理,得到待测的畸变波前和复原的图像。处理算法流程如下:
1)根据极大似然估计或最小二乘优化理论建立目标函数E:
其中,Dk(fx,fy)是调制光阑处在第k个扫描位置时图像传感器采集的子图像的傅里叶变换Hk(fx,fy)是调制光阑处在第k个扫描位置时对应的光学通道的光学传递函数,是Hk(fx,fy)的共轭,K为扫描的次数,fx,fy为像平面频域的坐标。
2)根据傅里叶光学,估计出第k个光学通道的点扩散函数,系统的波前相位分布可以由泽尼克多项式来表征:
其中Zn(u,v)表示第n阶泽尼克多项式,αn表示第n阶泽尼克多项式的系数。
3)因此,目标函数E演变为关于变量αn的函数E(α),采用诸如SPGD,遗传算法,模拟退火,神经网络等最优化算法寻找使目标函数最小的解,即可求得待测的波前畸变4)根据求得的波前畸变,利用下式可以复原图像R:
实施例:
基于本发明波前探测和图像复原新方法的光学系统由成像系统1、空间调制单元2、图像传感器3和计算机4组成,整个系统的工作原理如图1所示:
其具体的工作过程为:
(1)将小孔径的光阑在成像系统1的光瞳平面内以叠层扫描的方式进行平移。叠层扫描如图2所示,用于扫描的小孔径光阑在光瞳二维平面内依次平移一定的距离完成调制,并保证相邻平移位置间光瞳平面上被光阑扫描的区域内存在一定的交叠。叠层扫描由空间调制单元2来完成,空间调制单元2可以是机械移动平台,也可以是数字微镜器件(DigitalMicromirror Device,DMD)或者空间光调制器(Spatial Light Modulator,SLM)等快速调制器件。
(2)图像传感器3依次记录下与各扫描位置所对应的受不同像差影响的子图像。
(3)在计算机4中,使用基于叠层扫描的相位差波前探测和图像复原算法对子图像进行处理,得到待测的畸变波前和复原的图像。
以上所述,仅为本发明中的具体实施方式,但本发明的保护范围并不局限于此。只要是通过叠层扫描空间调制来产生不同光学通道并利用相位差算法进行处理的成像方法、装置和系统,均属于本发明的保护范围。

Claims (4)

1.一种基于叠层扫描的相位差波前探测和图像复原方法,其特征在于该方法包含以下步骤:
第一步,将小孔径的光阑在成像系统的光瞳平面内以叠层扫描的方式进行平移,并用图像传感器依次记录下与各扫描位置所对应的受不同像差影响的子图像;
第二步,使用基于叠层扫描的相位差波前探测和图像复原算法对序列子图像进行处理,得到待测的畸变波前和复原的图像,处理算法流程如下:
1)根据极大似然估计或最小二乘优化理论建立目标函数E:
其中,Dk(fx,fy)是调制光阑处在第k个扫描位置时图像传感器采集的子图像的傅里叶变换,Hk(fx,fy)是调制光阑处在第k个扫描位置时对应的光学通道的光学传递函数,是Hk(fx,fy)的共轭,K为扫描的次数,fx,fy为像平面频域的坐标;
2)根据傅里叶光学,估计出第k个光学通道的点扩散函数,系统的波前相位分布由泽尼克多项式来表征:
其中Zn(u,v)表示第n阶泽尼克多项式,αn表示第n阶泽尼克多项式的系数;
3)目标函数E演变为关于变量αn的函数E(α),采用最优化算法寻找使目标函数最小的解,求得待测的波前畸变
4)根据求得的波前畸变,利用下式复原图像R:
该方法利用的装置包括成像系统(1),空间调制单元(2),图像传感器(3)和计算机(4);其中空间调制单元(2)放置于成像系统(1)的出瞳面上,图像传感器置于成像系统(1)的像平面上,并输入到计算机(4)中;成像目标发出的光波进入到成像系统(1)中,受到空间调制单元(2)的叠层调制后,在图像传感器(3)上成像,该图像存储在计算机(4)中,最后利用算法对所有子图像进行处理得到探测的波前以及恢复的图像;
该方法需要对成像系统的光瞳平面进行叠层扫描的空间调制,用于调制的小孔径光阑在光瞳二维平面内依次平移一定的距离完成扫描,并保证相邻平移位置间光瞳平面上被光阑扫描的区域内存在一定的交叠,交叠部分对应的像差使采集的子图像含有部分相同的相位信息,而非交叠部分对应的像差则在采集的子图像中引入了相位差异;
该方法子图像之间的相位差异隐藏在已知的扫描位置信息里,而不是引入的额外像差,因此,采集的子图像得以保留更多的细节信息,从而提高波前探测的范围和精度;
该方法叠层扫描不仅可以由机械移动平台完成,也可以由数字微镜器件(DigitalMicromirror Device,DMD)或者空间光调制器(Spatial Light Modulator,SLM)快速光学器件来完成,其相应速度也足够保证本发明对准静态像差的探测;
该方法通过引入离焦量来产生不同的光学通道,采用空间调制的方式得到含有相位差异的子图像,使图像得以保留更多细节信息,能够对更大范围的波前畸变进行更为准确的测量,提高了相位差法的波前探测能力;
该方法采用小孔径光阑对光瞳平面进行空间调制,使得一般的CCD像素尺寸都满足光学系统要求的奈斯奎特采样频率,保证了相位差算法的有效性,扩大了相位差法的应用领域;
该方法采用空间调制和相位差法相结合,以较小的时间代价消除了系统对其他光学元件的需求,避免引入额外的相位误差,使整套系统也更紧凑,成本更低廉,使用更方便。
2.根据权利要求1所述的基于叠层扫描的相位差波前探测和图像复原方法,其特征在于,该方法既适用于点光源也适用于扩展目标。
3.根据权利要求1所述的基于叠层扫描的相位差波前探测和图像复原方法,其特征在于,该方法既可以用于相干光成像也可以用于非相干光成像。
4.根据权利要求1所述的基于叠层扫描的相位差波前探测和图像复原方法,其特征在于,该方法既可以用来探测波前畸变,也可以用来恢复图像。
CN201610136602.9A 2016-03-10 2016-03-10 一种基于叠层扫描的相位差波前探测和图像复原方法 Active CN105716725B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610136602.9A CN105716725B (zh) 2016-03-10 2016-03-10 一种基于叠层扫描的相位差波前探测和图像复原方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610136602.9A CN105716725B (zh) 2016-03-10 2016-03-10 一种基于叠层扫描的相位差波前探测和图像复原方法

Publications (2)

Publication Number Publication Date
CN105716725A CN105716725A (zh) 2016-06-29
CN105716725B true CN105716725B (zh) 2019-02-12

Family

ID=56157616

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610136602.9A Active CN105716725B (zh) 2016-03-10 2016-03-10 一种基于叠层扫描的相位差波前探测和图像复原方法

Country Status (1)

Country Link
CN (1) CN105716725B (zh)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3270404A1 (en) * 2016-07-13 2018-01-17 FEI Company Method of imaging a specimen using ptychography
CN109407311B (zh) * 2018-12-28 2021-06-08 哈尔滨工业大学 一种基于相位差异法的衍射光学成像系统波前反演算法
CN109785266A (zh) * 2019-01-17 2019-05-21 哈尔滨工业大学 基于空间光调制器的single-shot叠层相位恢复技术
CN109782432A (zh) * 2019-01-17 2019-05-21 哈尔滨工业大学 基于空间光调制器的汇聚透镜single-shot叠层相位恢复技术
CN110044498B (zh) * 2019-04-18 2021-02-19 中国科学院光电技术研究所 一种基于深度学习的哈特曼波前传感器模式波前复原方法
CN110458901A (zh) * 2019-06-26 2019-11-15 西安电子科技大学 一种基于计算成像的光电成像系统的全局性优化设计方法
CN111121969B (zh) * 2019-12-31 2021-12-10 杭州电子科技大学 一种基于孔径扫描的多光谱叠层成像方法
CN111179371B (zh) * 2019-12-31 2022-05-17 南京大学 一种扫描电子显微成像方法
CN117419815B (zh) * 2023-12-19 2024-03-08 苏州吉天星舟空间技术有限公司 一种太空望远镜波前探测装置及其波前探测方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102636270A (zh) * 2012-04-09 2012-08-15 中国人民解放军63655部队 一种大气相干长度的光学测量仪器及方法
CN102866133A (zh) * 2012-09-11 2013-01-09 中国科学院上海光学精密机械研究所 透射型样品振幅和相位成像装置和方法
CN204788651U (zh) * 2015-05-22 2015-11-18 中国科学院光电研究院 一种微孔衍射波前质量的测量装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
HU170623B (zh) * 1974-06-25 1977-07-28

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102636270A (zh) * 2012-04-09 2012-08-15 中国人民解放军63655部队 一种大气相干长度的光学测量仪器及方法
CN102866133A (zh) * 2012-09-11 2013-01-09 中国科学院上海光学精密机械研究所 透射型样品振幅和相位成像装置和方法
CN204788651U (zh) * 2015-05-22 2015-11-18 中国科学院光电研究院 一种微孔衍射波前质量的测量装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Synthetic aperture imaging by using spatial modulation diversity technology with stochastic parallel gradient descent algorithm;Haotong Ma等;《Optics Express》;20150528;第23卷(第11期);正文第1页第1段-第6页第2段,第12页第1段,第14页第1-2段及图1 *
小孔扫描傅里叶叠层成像的关键参量研究;谢宗良等;《光学学报》;20151031;第35卷(第10期);正文第1011002-1页第1段-第1011002-3页第1段及图1 *

Also Published As

Publication number Publication date
CN105716725A (zh) 2016-06-29

Similar Documents

Publication Publication Date Title
CN105716725B (zh) 一种基于叠层扫描的相位差波前探测和图像复原方法
US8305485B2 (en) Digital camera with coded aperture rangefinder
Mahalanobis et al. Recent results of medium wave infrared compressive sensing
CN106803892B (zh) 一种基于光场测量的光场高清晰成像方法
US8587686B1 (en) Hybrid differential optical sensing imager
CN105589210B (zh) 一种基于光瞳调制的数字化合成孔径成像方法
McCormac et al. DONUTS: A science frame autoguiding algorithm with sub-pixel precision, capable of guiding on defocused stars
US20110267485A1 (en) Range measurement using a coded aperture
US20070278386A1 (en) Measurement-diverse imaging and wavefront sensing with amplitude and phase estimation
US20110267477A1 (en) Range measurement using symmetric coded apertures
CN110675451B (zh) 基于相空间光学的数字自适应校正方法及系统
Vidal et al. Analysis of on-sky MOAO performance of CANARY using natural guide stars
WO2010018515A1 (en) Measuring and correcting lens distortion in a multispot scanning device.
Li et al. Self-measurements of point-spread function for remote sensing optical imaging instruments
US8749640B2 (en) Blur-calibration system for electro-optical sensors and method using a moving multi-focal multi-target constellation
JP2020506390A (ja) 電磁場の波頭の断層撮影分布を取得する方法及び光学システム
Li et al. Modulation transfer function measurements using a learning approach from multiple diffractive grids for optical cameras
CN104142131B (zh) 一种相位成像系统
Seshadri et al. Initial results from a laboratory emulation of weak gravitational lensing measurements
Du Bosq et al. An overview of joint activities on computational imaging and compressive sensing systems by NATO SET-232
Neuner III et al. Digital adaptive optics with interferometric homodyne encoding for mitigating atmospheric turbulence
US20220018718A1 (en) Wavefront curvature sensor involving temporal sampling of the image intensity distribution
Driggers et al. Target identification performance of infrared imager models as a function of blur and sampling
Monaldi et al. Rolling shutter effect aberration compensation in digital holographic microscopy
Xie et al. Restoration of degraded images using pupil-size diversity technology with stochastic parallel gradient descent algorithm

Legal Events

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