CN112197877B - 大视场巡天望远镜波前曲率传感方法、装置、设备及介质 - Google Patents

大视场巡天望远镜波前曲率传感方法、装置、设备及介质 Download PDF

Info

Publication number
CN112197877B
CN112197877B CN202011246688.3A CN202011246688A CN112197877B CN 112197877 B CN112197877 B CN 112197877B CN 202011246688 A CN202011246688 A CN 202011246688A CN 112197877 B CN112197877 B CN 112197877B
Authority
CN
China
Prior art keywords
curvature
wavefront
boundary
telescope
slope
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
CN202011246688.3A
Other languages
English (en)
Other versions
CN112197877A (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.)
Changchun Institute of Optics Fine Mechanics and Physics of CAS
Original Assignee
Changchun Institute of Optics Fine Mechanics and Physics 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 Changchun Institute of Optics Fine Mechanics and Physics of CAS filed Critical Changchun Institute of Optics Fine Mechanics and Physics of CAS
Priority to CN202011246688.3A priority Critical patent/CN112197877B/zh
Publication of CN112197877A publication Critical patent/CN112197877A/zh
Application granted granted Critical
Publication of CN112197877B publication Critical patent/CN112197877B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • G—PHYSICS
    • G02—OPTICS
    • G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B23/00—Telescopes, e.g. binoculars; Periscopes; Instruments for viewing the inside of hollow bodies; Viewfinders; Optical aiming or sighting devices
    • 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
    • G01J2009/002—Wavefront phase distribution

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Optics & Photonics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Astronomy & Astrophysics (AREA)
  • Telescopes (AREA)

Abstract

本申请公开了一种大视场巡天望远镜波前曲率传感方法、装置、设备及介质，该方法包括：包括：对大口径主焦点组件进行粗对准后，通过位于焦面两侧的错位型曲率传感器成像，得到离焦星点图；获取离焦星点图的边界，将离焦星点图分割成多个子孔径区域；子孔径区域被分成未含边界的纯内部区域和含边界的区域；在纯内部区域内，求解出各纯内部区域对应的波前曲率；在含边界的区域内，求解出含边界的区域对应的波前曲率和波前斜率，并利用极大似然估计约束波前斜率的参数；将求解出的波前曲率和约束后的波前斜率进行综合整理分析，获得望远镜的波前相位信息。这样可以使最终求解的波前传感的结果更精确，降低校正时间，提高校正精度。

Description

大视场巡天望远镜波前曲率传感方法、装置、设备及介质
技术领域
本发明涉及望远镜技术领域，特别是涉及一种大视场巡天望远镜波前曲率传感方法、装置、设备及介质。
背景技术
望远镜口径的增加不仅可以有效地提高对临近目标的分辨能力，同时还以平方规律提升望远镜的集光能力，可有效提升暗弱目标成像信噪比、拓展极限探测能力，最终实现对更加深远的宇宙的探索。因此，大口径大视场望远镜是未来验证宇宙学最新理论、增加时域天文等领域学术话语权的关键。
大口径大视场望远镜在近二十年来发展获得了飞速发展，为了获得更高的巡天效率与集光能力，其口径与视场都在不断扩大。主动光学作为大口径大视场望远镜的关键技术，已经获得了广泛的应用。国外已经研制并成功运行多台大口径大视场望远镜，八米级的LSST已经投入建设，而国内尚未开展两米以上的大视场望远镜研究。不论是在占领“太空高地”保障国土安全方面，还是在探测存在撞击威胁的小行星等天文学邻域，均存在较大的差距。为了进一步发挥大口径大视场望远镜的探测能力，通过主动光学对望远镜中的各个主要部件进行独立、实时的面形校正与姿态控制，不仅可以降低对光学加工、系统装配精度的要求，还可以有效地放宽对大型跟踪架刚度的要求，降低系统运动惯量。相对于高分辨成像望远镜，大口径大视场望远镜观测任务更加紧张，更长的观测时间会直接影响结果为所面临的外部观测环境更加的恶劣。
目前，一般利用相位差异技术求解大视场巡天望远镜系统的波前相位，该技术是基于波前经过存在差异的光学系统(如引入少量的离焦)后产生的变化，利用优化迭代的方法进行波前相位的求解。但是，对于大视场巡天系统而言，相位差异技术单一像差量程小于一个波长的特点，导致其无法有效地针对巡天望远镜像差较大的边缘视场；同时，天光背景与渐晕造成的光瞳非均匀照明，对相位恢复的精度也会造成影响。因此，相位差异技术在巡天主动光学望远镜波前传感系统的应用，目前还存在着诸多限制因素。
因此，如何精确地求解波前相位，以提高整个主动光学系统的校正效果，是本领域技术人员亟待解决的技术问题。
发明内容
有鉴于此，本发明的目的在于提供一种大视场巡天望远镜波前曲率传感方法、装置、设备及介质，可以更精确地求解波前相位信息，降低校正时间，提高校正精度。其具体方案如下：
一种大视场巡天望远镜波前曲率传感方法，包括：
对大口径主焦点组件进行粗对准后，通过位于焦面两侧的错位型曲率传感器成像，得到离焦星点图；
获取所述离焦星点图的边界，将所述离焦星点图分割成多个子孔径区域；所述子孔径区域被分成未含边界的纯内部区域和含边界的区域；
在所述纯内部区域内，求解出所述纯内部区域对应的波前曲率；
在所述含边界的区域内，求解出所述含边界的区域对应的波前曲率和波前斜率，并利用极大似然估计约束波前斜率的参数；
将求解出的波前曲率和约束后的波前斜率进行综合整理分析，获得所述望远镜的波前相位信息。
优选地，在本发明实施例提供的上述大视场巡天望远镜波前曲率传感方法中，获取所述离焦星点图的边界，包括：
利用膨胀运算将所述离焦星点图的孔洞进行平滑；
利用腐蚀运算将所述离焦星点图的边缘进行还原；所述边缘对应的像素同时含有光瞳边界斜率信息以及光瞳内部曲率信息；
利用膨胀运算将所述离焦星点图的边界逐层取出。
优选地，在本发明实施例提供的上述大视场巡天望远镜波前曲率传感方法中，所述子孔径区域是按照像素阵列排布的方向分割而成。
优选地，在本发明实施例提供的上述大视场巡天望远镜波前曲率传感方法中，波前斜率的参数包括斜率方向和斜率大小；所述斜率方向垂直于波前曲率的方向。
优选地，在本发明实施例提供的上述大视场巡天望远镜波前曲率传感方法中，所述主焦点组件包括科学终端；所述错位型曲率传感器设置在所述科学终端的内部。
本发明实施例还提供了一种大视场巡天望远镜波前曲率传感装置，包括：
曲率传感成像模块，用于对大口径主焦点组件进行粗对准后，通过位于焦面两侧的错位型曲率传感器成像，得到离焦星点图；
子孔径分割模块，用于获取所述离焦星点图的边界，将所述离焦星点图分割成多个子孔径区域；所述子孔径区域被分成未含边界的纯内部区域和含边界的区域；
纯内部区域计算模块，用于在所述纯内部区域内，求解出所述纯内部区域对应的波前曲率；
边界区域计算模块，用于在所述含边界的区域内，求解出所述含边界的区域对应的波前曲率和波前斜率，并利用极大似然估计约束波前斜率的参数；
信息整理分析模块，用于将求解出的波前曲率和约束后的波前斜率进行综合整理分析，获得所述望远镜的波前相位信息。
本发明实施例还提供了一种大视场巡天望远镜波前曲率传感设备，包括处理器和存储器，其中，所述处理器执行所述存储器中保存的计算机程序时实现如本发明实施例提供的上述大视场巡天望远镜波前曲率传感方法。
本发明实施例还提供了一种计算机可读存储介质，用于存储计算机程序，其中，所述计算机程序被处理器执行时实现如本发明实施例提供的上述大视场巡天望远镜波前曲率传感方法。
从上述技术方案可以看出，本发明所提供的一种大视场巡天望远镜波前曲率传感方法，包括：对大口径主焦点组件进行粗对准后，通过位于焦面两侧的错位型曲率传感器成像，得到离焦星点图；获取离焦星点图的边界，将离焦星点图分割成多个子孔径区域；子孔径区域被分成未含边界的纯内部区域和含边界的区域；在纯内部区域内，求解出各纯内部区域对应的波前曲率；在含边界的区域内，求解出含边界的区域对应的波前曲率和波前斜率，并利用极大似然估计约束波前斜率的参数；将求解出的波前曲率和约束后的波前斜率进行综合整理分析，获得望远镜的波前相位信息。
本发明在系统完成粗对准后，通过曲率传感的方式得到离焦星点图，再通过分区域的方式获取不同的子孔径区域，分别对仅包含曲率的区域和同时包含曲率和斜率的区域进行求解，最后综合整理分析，可以使最终求解的波前传感的结果更精确，降低了整个主动光学系统的校正时间，提高了整个主动光学系统的校正精度，进而提高望远镜对深空领域的观测精度，满足实际需求。此外，本发明还针对大视场巡天望远镜波前曲率传感方法提供了相应的装置、设备及计算机可读存储介质，进一步使得上述方法更具有实用性，该装置、设备及计算机可读存储介质具有相应的优点。
附图说明
为了更清楚地说明本发明实施例或相关技术中的技术方案，下面将对实施例或相关技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据提供的附图获得其他的附图。
图1为本发明实施例提供的大视场巡天望远镜波前曲率传感方法的流程图；
图2为本发明实施例提供的望远镜系统的结构示意图；
图3为本发明实施例提供的子孔径区域的平面示意图；
图4为本发明实施例提供的子孔径区域的立体结构示意图；
图5为本发明实施例提供的曲率传感的原理示意图；
图6为本发明实施例提供的大视场巡天望远镜波前曲率传感装置的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明提供一种大视场巡天望远镜波前曲率传感方法，如图1所示，包括以下步骤：
S101、对大口径主焦点组件进行粗对准后，通过位于焦面两侧的错位型曲率传感器成像，得到离焦星点图；
需要说明的是，提升大视场巡天望远镜口径，可以有效地提高暗弱目标探测信噪比、实现对更加深远宇宙的探索；针对大视场巡天望远镜越来越大的口径，为了充分发挥望远镜的极限探测能力、保证光学成像质量，本发明通过波前传感系统实现主要光学元件的面形校正与姿态控制，可以克服大视场巡天望远镜光学系统自身特点(如成像质量对光学元件失调误差敏感)与望远镜结构刚度不足的耦合影响；
另外，如图2所示，以主焦点式的大口径大视场巡天望远镜为例，利用粗对准环节对大口径主焦点组件进行对准后，波前传感器实现成像，得到的离焦星点图与望远镜的主镜面是一一对应的。这里的波前传感器选用错位型曲率传感器，均位于主焦点组件对应的焦面两侧，主焦点组件包括科学终端，该错位型曲率传感器可以设置在科学终端的内部，即在终端中内置错位型曲率传感器，设光源照明强度为无饱和，主焦点组件进行离焦运动时，错位型传感器基于两侧离焦光强分布的差分，可以对波前曲率进行估计；
S102、获取离焦星点图的边界，将离焦星点图分割成多个子孔径区域；子孔径区域被分成未含边界的纯内部区域和含边界的区域；
需要说明的是，曲率传感是通过离焦的星点像，对曲率进行估计，实现重建波前。它的边缘处跟波前的斜率相关，理论上可以裁掉边缘，但是在实际检测过程中，由于用像素取样，在边缘处就会产生曲率跟斜率都有的情况下，因此，为了使波前传感更加高精度，本发明将星点像分割成多个区域，对不同区域进行分类计算。如图3所示，虚线表示边界，虚线最下方的区域全部为纯内部区域，该纯内部区域仅包括曲率，虚线所在的区域全部为含边界的区域，该区域为曲率斜率混合的区域；如图4所示，圆环状正下方对应的像素为纯内部区域的像素(Pixel in Donut)，在圆环状的边缘线正下方对应的像素为含边界的区域的像素(Pixel on Boundary)，圆环状外边区域对应的像素为外部区域的像素(Pixel out ofDonut)，该外部区域与波前不产生任何的关系，可不做考虑；
S103、在纯内部区域内，求解出纯内部区域对应的波前曲率；
S104、在含边界的区域内，求解出含边界的区域对应的波前曲率和波前斜率，并利用极大似然估计约束波前斜率的参数；
可以理解的是，只将纯内部区域进行波前信息计算，方法较为简单，但如果单纯只进行曲率计算的话，它对一阶导数的这个边界条件是没有约束的，所以本发明为了进一步使波前传感更加高精度，同时考虑了要将斜率影响的区域的影响降到最低，在计算完纯内部区域的波前曲率后，还需同时兼顾含边界的区域，利用极大似然估计约束波前斜率的参数，这样得到的结果和实际光强分布会尽量地吻合；
S105、将求解出的波前曲率和约束后的波前斜率进行综合整理分析，获得望远镜的波前相位信息。
在本发明实施例提供的上述大视场巡天望远镜波前曲率传感方法中，在系统完成粗对准后，通过曲率传感的方式得到离焦星点图，再通过分区域的方式获取不同的子孔径区域，分别对仅包含曲率的区域和同时包含曲率和斜率的区域进行求解，最后综合整理分析，可以使最终求解的波前传感的结果更精确，降低了整个主动光学系统的校正时间，提高了整个主动光学系统的校正精度，进而提高望远镜对深空领域的观测精度，满足实际需求。
需要说明的是，曲率传感器的基本原理为光瞳处波前局部的曲率变化，所对应的焦内像与焦外像的光强分布会发生对应的变化。
具体地，如图5所示，P0为微透镜阵列的所在平面，P1为焦前离焦平面，P2为焦后离焦平面，F点为微透镜阵列的焦点，z为光轴方向。根据近场电磁波的传输方程，可以解算出波前信息，传输方程如下：
Figure BDA0002770262240000061
其中，ρ为微透镜阵列的光瞳内M(x,y)点的向量坐标，I(ρ)为M(x,y)的光强度，Φ(ρ)为M(x,y)点波前相位，
Figure BDA0002770262240000062
为梯度算子，
Figure BDA0002770262240000063
为波前斜率，
Figure BDA0002770262240000064
为拉普拉斯算子，
Figure BDA0002770262240000065
为波前曲率。
这里可以看出，曲率传感器得到的焦内像、焦外像的光强分布，与波前相位的二阶导数(即曲率)和一阶导数(即斜率)都有关系，主要区域还是和二阶导数(即曲率)有关。但是图像的边缘处和整个波前相位的一阶导数(即斜率)是息息相关的，为了解决这个一阶导数(即斜率)对这个内部区域的影像，所以本发明用分割区域的方式来进行计算。
对于自适应光学系统，一般离焦量(即离焦平面和焦平面的距离)仅为几个焦深，在离焦平面上形成的离焦星点像十分接近光瞳形状，由此可以认为：
Figure BDA0002770262240000071
其中，δ为光强变化曲率，r为光斑的半径；
故近场电磁波的传输方程通过近似可得：
Figure BDA0002770262240000072
其中，Δz为焦前离焦平面P1和焦后离焦平面P2共轭位置相对入瞳的距离，因此，Δz＝f(f-l)/l；其中，f为微透镜阵列的焦距，l为光波在离焦平面和焦平面之间的光程。
又因为在离焦平面上形成的离焦星点像十分接近光瞳形状，因此可认为：2I0＝I1(ρ)+I2(ρ)；
Figure BDA0002770262240000073
设
Figure BDA0002770262240000074
可以得到
Figure BDA0002770262240000075
对上式两端同时进行傅里叶变换，可得：
Figure BDA0002770262240000076
由此，对上式进行傅里叶变换逆运算，即可得：
Figure BDA0002770262240000077
其中，
Figure BDA0002770262240000078
为空间向量坐标，m和n为空间频率，Δz为离焦平面和焦平面在光轴方向上的距离，I1为焦前光强分布，I2为焦后光强分布，FFT为傅里叶变换，IFFT为逆傅里叶变换。
根据离焦光强分布公式，可获得离焦平面上随各个成像位置点(x,y)的变化的离焦光强分布I(x,y)，离焦光强分布公式如下：
Figure BDA0002770262240000081
其中，I0为焦面光强分布，(u,v)为微透镜焦平面上的成像位置点坐标，α为第一间距和第二间距的比值，第一间距为主镜头所在平面和微透镜阵列所在平面之间的间距，第二间距为主镜头所在平面和微透镜阵列的焦平面之间的间距。
在具体实施时，在本发明实施例提供的上述大视场巡天望远镜波前曲率传感方法中，在进行形态学分析前，需要针对离焦星点像一些探测器坏点、饱和点、以及散斑噪声造成的图形缺陷，因此，步骤S102获取离焦星点图的边界，具体可以包括：首先利用“膨胀”运算将离焦星点图的孔洞进行平滑；然后利用“腐蚀”运算将离焦星点图的真正的边缘进行还原；边缘对应的像素同时含有光瞳边界斜率信息以及光瞳内部曲率信息；最后利用“膨胀运算”将离焦星点图的边界逐层取出。
在具体实施时，在本发明实施例提供的上述大视场巡天望远镜波前曲率传感方法中，如图3所示，子孔径区域可以按照像素阵列排布的方向分割而成。可以理解的是，错位型传感器是基于两侧离焦光强分布的差分，对波前曲率进行估计，在建立波前传感过程中可以对差分区域进行分割，其分割方法与像素阵列方向相同。
在具体实施时，在本发明实施例提供的上述大视场巡天望远镜波前曲率传感方法中，波前斜率的参数包括斜率方向和斜率大小，也就是说，用一个向量对波前斜率进行描述，这样可以确定边界的法向在哪个方向；这里的斜率方向垂直于波前曲率的方向。
具体地，在每个子孔径之内，采用傅里叶变换方法进行求解。将子孔径分为两类，分别为：1、纯内部区域Si,j，由于曲率传感器获得的为波前曲率估计，因此，可以化简为齐次方程与二维一次函数fi,j(x,y)，fi,j(x,y)中共包括两个方向斜率与活塞误差三个参数记为ηi,j(αx,αy，βz)。相邻子孔径需要共用一条边界；2、含边界的区域，曲率传感所获得的光强图像的边缘像素中，同时含有光瞳边界斜率信息以及光瞳内部曲率信息。首先利用“膨胀”这一图形学常用的运算，将边界逐层取出，利用两层的差分和内部的层数进行比较，当二者接近时，即证明取得了斜率部分与曲率部分的真实边界，根据曲率传感的基本原理，其曲率垂直于边界，因此需要在估计法向方向的基础上，将斜率分解到正交方向上，Bηi,j(θi,j，βz)。
接下来，利用极大似然估计进行子孔径综合，首先在利用边界条件获得参数Bηi,j(θi,j，βz)i,j∈(1,N),(1,M)的估计，进一步获得ηi,j(αx,αy，βz)i,j∈(1,N),(1,M)，并在内部区域求解泊松方程；
设观测值y是一个随机变量，其概率密度p依赖于未知参数Bηi,j。为了由y估计Bηi,j，选取使似然函数L(D/Bηi,j)＝p(D/Bηi,j)极大化的Bηi,j值。即存在
Figure BDA0002770262240000091
是L(D/Bηi,j)的最大值，则
Figure BDA0002770262240000092
是准确数值的可能性最大，此时称
Figure BDA0002770262240000093
为Bηi,j的极大似然估计。
设观测值y∈Rm与随机噪声ε∈Rm，当{ε(k)}具有同协方差Σ的m维独立高斯分布序列时，则似然函数为：
Figure BDA0002770262240000094
当Σ未知时，首先利用最小二乘估计得到
Figure BDA0002770262240000095
进而得到
Figure BDA0002770262240000096
求
Figure BDA0002770262240000097
可得极大似然估计。
需要说明的是，本发明也可以利用最小二乘估计法来约束波前斜率的参数。但是，利用极大似然估计法，不论残差偏离正态分布是否为白噪声时，都可达到相同的渐进方差，因此，极大似然法的鲁棒性要远远好于最小二乘估计法。
在具体实施时，在执行步骤S105获得望远镜的波前相位信息之后，还可以包括：建立变换矩阵，将波前相位信息采用Zernike系数表达与执行元件驱动量之间的关系。具体地，首先获取主焦点组件相对主镜的偏离信息，然后以获得的波前相位信息(即各个视场所测得的低阶像差)为输入，以获取的主焦点组件相对主镜的偏离信息为输出，建立并训练神经网络，之后通过训练好的神经网络解算各个视场所测得的低阶像差，获得主焦点组件对应的偏离量。这样可以使神经网络输出与实际位置误差的二次方和最小。利用该神经网络可以得到改正镜组姿态误差，包括倾斜、平移、面内旋转。
基于同一发明构思，本发明实施例还提供了一种大视场巡天望远镜波前曲率传感装置，由于该装置解决问题的原理与前述一种大视场巡天望远镜波前曲率传感方法相似，因此该装置的实施可以参见大视场巡天望远镜波前曲率传感方法的实施，重复之处不再赘述。
在具体实施时，本发明实施例提供的大视场巡天望远镜波前曲率传感装置，如图6所示，包括：
曲率传感成像模块11，用于对大口径主焦点组件进行粗对准后，通过位于焦面两侧的错位型曲率传感器成像，得到离焦星点图；
子孔径分割模块12，用于获取离焦星点图的边界，将离焦星点图分割成多个子孔径区域；子孔径区域被分成未含边界的纯内部区域和含边界的区域；
纯内部区域计算模块13，用于在纯内部区域内，求解出各纯内部区域对应的波前曲率；
边界区域计算模块14，用于在含边界的区域内，求解出含边界的区域对应的波前曲率和波前斜率，并利用极大似然估计约束波前斜率的参数；
信息整理分析模块15，用于将求解出的波前曲率和约束后的波前斜率进行综合整理分析，获得望远镜的波前相位信息。
在本发明实施例提供的上述大视场巡天望远镜波前曲率传感装置中，可以通过上述五个模块的相互作用，可以使最终求解的波前传感的结果更精确，降低了整个主动光学系统的校正时间，提高了整个主动光学系统的校正精度，进而提高望远镜对深空领域的观测精度。
关于上述各个模块更加具体的工作过程可以参考前述实施例公开的相应内容，在此不再进行赘述。
相应的，本发明实施例还公开了一种大视场巡天望远镜波前曲率传感设备，包括处理器和存储器；其中，处理器执行存储器中保存的计算机程序时实现前述实施例公开的大视场巡天望远镜波前曲率传感方法。
关于上述方法更加具体的过程可以参考前述实施例中公开的相应内容，在此不再进行赘述。
进一步的，本发明还公开了一种计算机可读存储介质，用于存储计算机程序；计算机程序被处理器执行时实现前述公开的大视场巡天望远镜波前曲率传感方法。
关于上述方法更加具体的过程可以参考前述实施例中公开的相应内容，在此不再进行赘述。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其它实施例的不同之处，各个实施例之间相同或相似部分互相参见即可。对于实施例公开的装置、设备、存储介质而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
专业人员还可以进一步意识到，结合本文中所公开的实施例描述的各示例的单元及算法步骤，能够以电子硬件、计算机软件或者二者的结合来实现，为了清楚地说明硬件和软件的可互换性，在上述说明中已经按照功能一般性地描述了各示例的组成及步骤。这些功能究竟以硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本申请的范围。
结合本文中所公开的实施例描述的方法或算法的步骤可以直接用硬件、处理器执行的软件模块，或者二者的结合来实施。软件模块可以置于随机存储器(RAM)、内存、只读存储器(ROM)、电可编程ROM、电可擦除可编程ROM、寄存器、硬盘、可移动磁盘、CD-ROM、或技术领域内所公知的任意其它形式的存储介质中。
本发明实施例提供的一种大视场巡天望远镜波前曲率传感方法，包括：对大口径主焦点组件进行粗对准后，通过位于焦面两侧的错位型曲率传感器成像，得到离焦星点图；获取离焦星点图的边界，将离焦星点图分割成多个子孔径区域；子孔径区域被分成未含边界的纯内部区域和含边界的区域；在纯内部区域内，求解出各纯内部区域对应的波前曲率；在含边界的区域内，求解出含边界的区域对应的波前曲率和波前斜率，并利用极大似然估计约束波前斜率的参数；将求解出的波前曲率和约束后的波前斜率进行综合整理分析，获得望远镜的波前相位信息。这样可以使最终的波前传感的结果求解的更精确，降低了整个主动光学系统的校正时间，提高了整个主动光学系统的校正精度，进而提高望远镜对深空领域的观测精度，满足实际需求。此外，本发明还针对大视场巡天望远镜波前曲率传感方法提供了相应的装置、设备及计算机可读存储介质，进一步使得上述方法更具有实用性，该装置、设备及计算机可读存储介质具有相应的优点。
最后，还需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
以上对本发明所提供的大视场巡天望远镜波前曲率传感方法、装置、设备及介质进行了详细介绍，本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。

Claims (7)

1.一种大视场巡天望远镜波前曲率传感方法，其特征在于，包括：
对大口径主焦点组件进行粗对准后，通过位于焦面两侧的错位型曲率传感器成像，得到离焦星点图；
利用膨胀运算将所述离焦星点图的孔洞进行平滑；
利用腐蚀运算将所述离焦星点图的边缘进行还原；所述边缘对应的像素同时含有光瞳边界斜率信息以及光瞳内部曲率信息；
利用膨胀运算将所述离焦星点图的边界逐层取出；
将所述离焦星点图分割成多个子孔径区域；所述子孔径区域被分成未含边界的纯内部区域和含边界的区域；
在所述纯内部区域内，求解出所述纯内部区域对应的波前曲率；
在所述含边界的区域内，求解出所述含边界的区域对应的波前曲率和波前斜率，并利用极大似然估计约束波前斜率的参数；
将求解出的波前曲率和约束后的波前斜率进行综合整理分析，获得所述望远镜的波前相位信息。
2.根据权利要求1所述的大视场巡天望远镜波前曲率传感方法，其特征在于，所述子孔径区域是按照像素阵列排布的方向分割而成。
3.根据权利要求2所述的大视场巡天望远镜波前曲率传感方法，其特征在于，波前斜率的参数包括斜率方向和斜率大小；所述斜率方向垂直于波前曲率的方向。
4.根据权利要求3所述的大视场巡天望远镜波前曲率传感方法，其特征在于，所述主焦点组件包括科学终端；所述错位型曲率传感器设置在所述科学终端的内部。
5.一种大视场巡天望远镜波前曲率传感装置，其特征在于，包括：
曲率传感成像模块，用于对大口径主焦点组件进行粗对准后，通过位于焦面两侧的错位型曲率传感器成像，得到离焦星点图；
子孔径分割模块，用于利用膨胀运算将所述离焦星点图的孔洞进行平滑；利用腐蚀运算将所述离焦星点图的边缘进行还原；所述边缘对应的像素同时含有光瞳边界斜率信息以及光瞳内部曲率信息；利用膨胀运算将所述离焦星点图的边界逐层取出；将所述离焦星点图分割成多个子孔径区域；所述子孔径区域被分成未含边界的纯内部区域和含边界的区域；
纯内部区域计算模块，用于在所述纯内部区域内，求解出所述纯内部区域对应的波前曲率；
边界区域计算模块，用于在所述含边界的区域内，求解出所述含边界的区域对应的波前曲率和波前斜率，并利用极大似然估计约束波前斜率的参数；
信息整理分析模块，用于将求解出的波前曲率和约束后的波前斜率进行综合整理分析，获得所述望远镜的波前相位信息。
6.一种大视场巡天望远镜波前曲率传感设备，其特征在于，包括处理器和存储器，其中，所述处理器执行所述存储器中保存的计算机程序时实现如权利要求1至4任一项所述的大视场巡天望远镜波前曲率传感方法。
7.一种计算机可读存储介质，其特征在于，用于存储计算机程序，其中，所述计算机程序被处理器执行时实现如权利要求1至4任一项所述的大视场巡天望远镜波前曲率传感方法。
CN202011246688.3A 2020-11-10 2020-11-10 大视场巡天望远镜波前曲率传感方法、装置、设备及介质 Active CN112197877B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011246688.3A CN112197877B (zh) 2020-11-10 2020-11-10 大视场巡天望远镜波前曲率传感方法、装置、设备及介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011246688.3A CN112197877B (zh) 2020-11-10 2020-11-10 大视场巡天望远镜波前曲率传感方法、装置、设备及介质

Publications (2)

Publication Number Publication Date
CN112197877A CN112197877A (zh) 2021-01-08
CN112197877B true CN112197877B (zh) 2021-10-08

Family

ID=74034449

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011246688.3A Active CN112197877B (zh) 2020-11-10 2020-11-10 大视场巡天望远镜波前曲率传感方法、装置、设备及介质

Country Status (1)

Country Link
CN (1) CN112197877B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113280932B (zh) * 2021-04-19 2022-04-12 西安理工大学 相干光通信波前校正中去除与补偿子孔径光斑坏点的方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7619191B1 (en) * 2006-04-28 2009-11-17 University Of Central Florida Research Foundation, Inc. Increase spatial sampling for wave front mid-spatial frequency error recovery
CN102721478A (zh) * 2012-07-10 2012-10-10 中国科学院光电技术研究所 一种应用于曲率波前传感器的波前复原方法
US20160209646A1 (en) * 2013-09-20 2016-07-21 Inter-University Research Institute Corporation National Institutes Of Natural Sciences Adaptive optics system and optical device
CN106969844A (zh) * 2017-01-06 2017-07-21 北京理工大学 基于子孔径波前振幅调制的斜率和曲率信号提取方法
CN109343212A (zh) * 2018-11-02 2019-02-15 中国科学院紫金山天文台 大视场巡天望远镜主动光学系统实现方法
CN109407311A (zh) * 2018-12-28 2019-03-01 哈尔滨工业大学 一种基于相位差异法的衍射光学成像系统波前反演算法
CN111751013A (zh) * 2020-07-07 2020-10-09 中国科学院长春光学精密机械与物理研究所 光学成像的像差测量方法及光学成像方法
CN111811785A (zh) * 2020-07-21 2020-10-23 中国科学院长春光学精密机械与物理研究所 大动态范围望远镜像差的检测方法、装置、设备及介质

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3635079B2 (ja) * 2002-08-20 2005-03-30 規 有賀 撮像装置とその光学系

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7619191B1 (en) * 2006-04-28 2009-11-17 University Of Central Florida Research Foundation, Inc. Increase spatial sampling for wave front mid-spatial frequency error recovery
CN102721478A (zh) * 2012-07-10 2012-10-10 中国科学院光电技术研究所 一种应用于曲率波前传感器的波前复原方法
US20160209646A1 (en) * 2013-09-20 2016-07-21 Inter-University Research Institute Corporation National Institutes Of Natural Sciences Adaptive optics system and optical device
CN106969844A (zh) * 2017-01-06 2017-07-21 北京理工大学 基于子孔径波前振幅调制的斜率和曲率信号提取方法
CN109343212A (zh) * 2018-11-02 2019-02-15 中国科学院紫金山天文台 大视场巡天望远镜主动光学系统实现方法
CN109407311A (zh) * 2018-12-28 2019-03-01 哈尔滨工业大学 一种基于相位差异法的衍射光学成像系统波前反演算法
CN111751013A (zh) * 2020-07-07 2020-10-09 中国科学院长春光学精密机械与物理研究所 光学成像的像差测量方法及光学成像方法
CN111811785A (zh) * 2020-07-21 2020-10-23 中国科学院长春光学精密机械与物理研究所 大动态范围望远镜像差的检测方法、装置、设备及介质

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于结构函数的子孔径拼接算法研究;安其昌等;《红外与激光工程》;20150331;第44卷(第3期);第929-933页 *

Also Published As

Publication number Publication date
CN112197877A (zh) 2021-01-08

Similar Documents

Publication Publication Date Title
US8305485B2 (en) Digital camera with coded aperture rangefinder
US8436912B2 (en) Range measurement using multiple coded apertures
US8432479B2 (en) Range measurement using a zoom camera
Mahmood et al. Optimal depth estimation by combining focus measures using genetic programming
US10317597B2 (en) Light-field microscopy with phase masking
US8330852B2 (en) Range measurement using symmetric coded apertures
CN106845024B (zh) 一种基于波前反演的光学卫星在轨成像仿真方法
US20110267485A1 (en) Range measurement using a coded aperture
US20130216124A1 (en) Spatial Reconstruction of Plenoptic Images
EP2702556A1 (en) A method and apparatus for providing image data for constructing an image of a region of a target object
CN111811785B (zh) 大动态范围望远镜像差的检测方法、装置、设备及介质
US7113268B2 (en) Scintillation tolerant optical field sensing system and associated method
CN111812836A (zh) 大口径大视场相机的标校装调方法、装置、设备及介质
Li et al. Efficient assessment method of on-board modulation transfer function of optical remote sensing sensors
US10230940B2 (en) Method for determining the complex amplitude of the electromagnetic field associated with a scene
CN112197877B (zh) 大视场巡天望远镜波前曲率传感方法、装置、设备及介质
CN115523839A (zh) 一种基于观测目标的拼接子镜边缘误差高精度实时探测方法
CN112113514B (zh) 一种波前信息的检测方法、装置以及系统
KR20190138632A (ko) 전자기장의 파면의 단층 영상 분포를 획득하는 방법 및 시스템
CN112525496B (zh) 一种巡天望远镜波前曲率传感方法、装置、设备及介质
JP2021060353A (ja) 光学系の内部誤差の推定方法および計測装置
Wang et al. Parameterized modeling of spatially varying PSF for lens aberration and defocus
CN115683360A (zh) 一种大口径望远镜波前曲率传感方法、电子设备及计算机可读存储介质
Salvaggio Image Quality Modeling and Optimization for Non-Conventional Aperture Imaging Systems
Bonse et al. Use the 4S (Signal-Safe Speckle Subtraction): Explainable Machine Learning reveals the Giant Exoplanet AF Lep b in High-Contrast Imaging Data from 2011

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