[go: up one dir, main page]

CN115183886A - Wavefront sensor based on defocused grating array - Google Patents

Wavefront sensor based on defocused grating array Download PDF

Info

Publication number
CN115183886A
CN115183886A CN202210809078.2A CN202210809078A CN115183886A CN 115183886 A CN115183886 A CN 115183886A CN 202210809078 A CN202210809078 A CN 202210809078A CN 115183886 A CN115183886 A CN 115183886A
Authority
CN
China
Prior art keywords
wavefront
array
aperture
sub
grating array
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.)
Granted
Application number
CN202210809078.2A
Other languages
Chinese (zh)
Other versions
CN115183886B (en
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 CN202210809078.2A priority Critical patent/CN115183886B/en
Publication of CN115183886A publication Critical patent/CN115183886A/en
Application granted granted Critical
Publication of CN115183886B publication Critical patent/CN115183886B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • GPHYSICS
    • G02OPTICS
    • G02BOPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B27/00Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00
    • G02B27/42Diffraction optics, i.e. systems including a diffractive element being designed for providing a diffractive effect
    • G02B27/44Grating systems; Zone plate systems
    • 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
    • G01J2009/002Wavefront phase distribution

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Optics & Photonics (AREA)
  • Testing Of Optical Devices Or Fibers (AREA)

Abstract

本发明公开了一种基于离焦光栅阵列的波前传感器,利用离焦光栅阵列,同时实现对待测光束的孔径分割和子孔径内±1级离焦光斑的生成,利用每个子孔径中的±1级离焦光斑图像可提取子孔径中波前斜率信息和高阶像差模式系数,实现对子孔径内波前的高分辨测量。与传统夏克‑哈特曼波前传感器相比,本发明可以提取子孔径内更多的细节分布信息,因而可以在同等子孔径条件下以更高精度复原波前,或在相同复原精度下以更少的子孔径数目完成波前复原,从而为弱光、高精度波前探测等领域提供了一种新的技术途径。

Figure 202210809078

The invention discloses a wavefront sensor based on a defocus grating array, which utilizes the defocus grating array to simultaneously realize the aperture division of the beam to be measured and the generation of ±1-order defocused light spots in the sub-apertures. The first-order defocus spot image can extract the wavefront slope information and high-order aberration mode coefficients in the sub-aperture, and realize high-resolution measurement of the wavefront in the sub-aperture. Compared with the traditional Shack-Hartmann wavefront sensor, the present invention can extract more detailed distribution information in the sub-aperture, so the wavefront can be recovered with higher accuracy under the same sub-aperture condition, or under the same recovery accuracy. The wavefront restoration is completed with a smaller number of sub-apertures, thus providing a new technical approach for the fields of low-light and high-precision wavefront detection.

Figure 202210809078

Description

一种基于离焦光栅阵列的波前传感器A wavefront sensor based on defocus grating array

技术领域technical field

本发明属于波前探测技术领域,尤其涉及一种基于离焦光栅阵列的波前传感器。The invention belongs to the technical field of wavefront detection, and in particular relates to a wavefront sensor based on a defocus grating array.

背景技术Background technique

波前探测一直是光学界研究热点之一,在天文观测、光学检测、医学成像、自适应光学等诸多领域都涉及到相位测量的问题。现有主流相位测量方法主要分为干涉法、直接测量法和基于强度分布的间接测量法三大类,每类方法都有其独特的优势,分别被应用于不同的场合。Wavefront detection has always been one of the research hotspots in the optical field, and phase measurement is involved in many fields such as astronomical observation, optical detection, medical imaging, and adaptive optics. The existing mainstream phase measurement methods are mainly divided into three categories: interferometry, direct measurement and indirect measurement based on intensity distribution. Each type of method has its own unique advantages and is used in different occasions.

离焦光栅实质上是一个离轴的菲涅尔波带片,具有普通光栅的棱镜作用,将入射波前在光栅的不同的衍射级上分束,同时具有菲涅尔波带片的透镜作用,在不同的衍射级上引入不同的透镜效应。最终使得输入波前有对称分布的±1级衍射光轴,且在±1级上具有大小相等,一正一负的焦距。The defocus grating is essentially an off-axis Fresnel zone plate, which has the prism function of an ordinary grating, splits the incident wavefront on different diffraction orders of the grating, and has the lens function of the Fresnel zone plate at the same time. , introducing different lens effects at different diffraction orders. Finally, the input wavefront has a symmetrical distribution of ±1st-order diffraction optical axes, and has equal focal lengths, one positive and one negative, on the ±1st order.

夏克-哈特曼波前传感器主要是通过微透镜阵列对波前进行分割采样,子波前经过微透镜阵列后被聚焦于光电探测器上,形成光斑阵列图,其测量精度与空间采样率存在固有矛盾,限制了夏克-哈特曼波前传感器的测量性能。曲率波前传感器是通过计算聚焦透镜前后光强的归一化之差,在光瞳内部与波前曲率成正比,但曲率波前传感器无法对高阶信息进行测量,限制了其发展。因此,目前亟需寻找一种探测精度高、鲁棒性好的低空间采样波前探测技术路线。The Shack-Hartmann wavefront sensor mainly divides and samples the wavefront through the microlens array. After passing through the microlens array, the sub-wavefront is focused on the photodetector to form a light spot array diagram. The measurement accuracy is related to the spatial sampling rate. There are inherent contradictions that limit the measurement performance of the Shack-Hartmann wavefront sensor. The curvature wavefront sensor calculates the normalized difference of the light intensity before and after the focusing lens, which is proportional to the wavefront curvature inside the pupil, but the curvature wavefront sensor cannot measure high-order information, which limits its development. Therefore, it is urgent to find a low-space sampling wavefront detection technology route with high detection accuracy and robustness.

结合离焦光栅的特性、曲率波前传感器的原理及哈特曼波前传感器的结构,提出一种基于离焦光栅阵列波前传感器,在单个子孔径内引入离焦光栅对入射波前进行相位调制,在其焦面上获得±1级离焦光斑,结合曲率波前传感器的原理,得到单个波前的信息,最终根据复原矩阵和单个子孔径内的波前信息复原波前。Combining the characteristics of the defocus grating, the principle of the curvature wavefront sensor and the structure of the Hartmann wavefront sensor, a defocus grating array wavefront sensor is proposed. The defocus grating is introduced into a single sub-aperture to phase the incident wavefront. Modulation to obtain a ±1-order defocused spot on its focal plane, combined with the principle of the curvature wavefront sensor, to obtain the information of a single wavefront, and finally restore the wavefront according to the restoration matrix and the wavefront information in a single sub-aperture.

发明内容SUMMARY OF THE INVENTION

本发明要解决的技术问题是:夏克-哈特曼波前传感器测量精度和空间采样率的固有矛盾以及曲率波前传感器无法对波前的高阶信息进行探测的问题。在低空间采样条件下完成高精度波前复原,在同等子孔径条件下以更高精度复原更高阶的波前像差信息。The technical problems to be solved by the present invention are: the inherent contradiction between the measurement accuracy of the Shack-Hartmann wavefront sensor and the spatial sampling rate, and the problem that the curvature wavefront sensor cannot detect the high-order information of the wavefront. Complete high-precision wavefront restoration under the condition of low spatial sampling, and restore higher-order wavefront aberration information with higher precision under the same sub-aperture condition.

本发明解决上述技术问题采用的技术方案是:一种基于离焦光栅阵列的波前传感器,所述波前传感器采用离焦光栅阵列对待测光束进行孔径分割和相位调制形成±1级衍射阵列,再利用微透镜阵列对±1级衍射阵列进行聚焦,形成±1级离焦光斑阵列;所述波前传感器主要包括离焦光栅阵列、微透镜阵列和光电探测器三部分;The technical solution adopted by the present invention to solve the above technical problems is: a wavefront sensor based on a defocus grating array, the wavefront sensor adopts the defocus grating array to perform aperture division and phase modulation of the beam to be measured to form a ±1 order diffraction array, Then use the microlens array to focus the ±1-order diffraction array to form a ±1-order defocused spot array; the wavefront sensor mainly includes three parts: a defocus grating array, a microlens array and a photodetector;

所述离焦光栅阵列中的光栅为相位型离焦光栅,用于对待测光束进行孔径分割以形成若干子孔径并对子孔径内波前进行相位调制;设单个子孔径调制前的波前为φj(x,y),j为子孔径编号,对子孔径内波前调制的相移函数为φm(x,y),则对第j个子孔径调制后的波前φ'j(x,y)为:The gratings in the defocus grating array are phase-type defocus gratings, which are used for aperture division of the beam to be measured to form a number of sub-apertures and phase modulation of the wavefront in the sub-aperture; let the wavefront before modulation of a single sub-aperture be φ j (x, y), j is the sub-aperture number, and the phase shift function of the wavefront modulation in the sub-aperture is φ m (x, y), then for the jth sub-aperture modulated wavefront φ' j (x ,y) is:

φ'j(x,y)=φj(x,y)+φm(x,y) φ'j (x,y)= φj (x,y)+ φm (x,y)

其中,x,y为离焦光栅的坐标;Among them, x, y are the coordinates of the defocus grating;

所述微透镜阵列位于所述离焦光栅阵列之后,且与离焦光栅阵列一一对应,所述微透镜阵列对通过离焦光栅阵列调制后的光束进行聚焦,在其焦面上形成±1级离焦光斑阵列;The microlens array is located behind the defocusing grating array and corresponds to the defocusing grating array one-to-one. The microlens array focuses the light beam modulated by the defocusing grating array, forming ±1 on the focal plane. stage defocused spot array;

所述光电探测器位于微透镜阵列的焦面,探测每个子孔径内的±1级离焦光斑图像数据,所述离焦光斑图像数据用于利用曲率波前复原方法计算子孔径中的高阶像差信息和重建全口径的高分辨波前。The photodetector is located at the focal plane of the microlens array, and detects ±1-order defocus spot image data within each sub-aperture, and the de-focus spot image data is used to calculate higher-order orders in the sub-aperture using the curvature wavefront restoration method Aberration information and reconstruction of full aperture high-resolution wavefronts.

进一步的,所述离焦光栅阵列大于等于2×2。Further, the defocus grating array is greater than or equal to 2×2.

进一步的,所述的曲率波前复原方法是基于两幅±1级离焦光斑图像的本征函数波前复原方法,包括所有基于±1级离焦光斑提取波前信息的方法。Further, the curvature wavefront restoration method is an eigenfunction wavefront restoration method based on two ±1-order defocused light spot images, including all methods for extracting wavefront information based on ±1-order defocused light spots.

进一步的,所述光电探测器可以是CCD、CMOS,或其他阵列型探测器。Further, the photodetector may be a CCD, CMOS, or other array type detectors.

本发明与现有技术相比有以下优点:Compared with the prior art, the present invention has the following advantages:

本发明利用离焦光栅对输入波前进行调制,调制后的波前在微透镜阵列的焦面可以得到正负焦面的信息,结合曲率波前传感器的波前复原算法,可以获得单个子孔径对应波前的高阶信息,以提高波前复原精度;本发明与传统夏克-哈特曼波前传感器复原方法相比,可以在相同子孔径条件下以更高精度复原波前,稀疏子孔径条件下可以以更高精度复原更高阶畸变波前,有望用于对弱光、高精度等领域的波前探测。In the present invention, the input wavefront is modulated by the defocus grating, and the modulated wavefront can obtain the information of the positive and negative focal planes on the focal plane of the microlens array. Combined with the wavefront recovery algorithm of the curvature wavefront sensor, a single sub-aperture can be obtained. Corresponding to the high-order information of the wavefront to improve the wavefront restoration accuracy; compared with the traditional Shack-Hartmann wavefront sensor restoration method, the present invention can restore the wavefront with higher accuracy under the same sub-aperture condition, and the sparse sub-aperture Under the condition of aperture, higher-order distorted wavefronts can be recovered with higher precision, and it is expected to be used for wavefront detection in the fields of weak light and high precision.

附图说明Description of drawings

图1为本发明实施例一2×2离焦光栅阵列波前传感器的结构示意图;FIG. 1 is a schematic structural diagram of a 2×2 defocus grating array wavefront sensor according to an embodiment of the present invention;

图2为本发明实施例一2×2离焦光栅阵列的相位图;FIG. 2 is a phase diagram of a 2×2 defocus grating array according to an embodiment of the present invention;

图3为本发明实施例一中的待测波前和焦面光强图,其中,(a)为待测波前图,(b)为焦面的光强图;Fig. 3 is the light intensity diagram of the wavefront to be measured and the focal plane in Embodiment 1 of the present invention, wherein (a) is the wavefront diagram to be measured, and (b) is the light intensity diagram of the focal plane;

图4为本发明实施例一波前复原结果,其中,(a)为本发明实施例一复原的波前图,(b)为本发明实施例一波前复原残差图;4 is a wavefront restoration result according to Embodiment 1 of the present invention, wherein (a) is a restored wavefront graph according to Embodiment 1 of the present invention, and (b) is a residual graph of wavefront restoration according to Embodiment 1 of the present invention;

图5为传统夏克-哈特曼波前传感器复原算法的复原结果,其中,(a)为模式法复原的波前图,(b)为模式法波前复原残差图。Fig. 5 is the restoration result of the traditional Shack-Hartmann wavefront sensor restoration algorithm, wherein (a) is the wavefront image restored by the mode method, and (b) is the residual image of the mode method wavefront restoration.

具体实施方式Detailed ways

为使本发明的目的和技术方案更加清楚明白,以下结合具体实施例一,并参照附图,对本发明进一步详细说明。In order to make the purpose and technical solution of the present invention clearer, the present invention will be further described in detail below with reference to the specific embodiment 1 and the accompanying drawings.

本发明为一种基于离焦光栅阵列波前传感器,实施例一中采用2×2离焦光栅阵列波前传感器,其光学结构如图1所示,离焦光栅对输入波前进行相位调制,离焦光栅阵列1位于微透镜阵列2前,CCD 3位于微透镜阵列2的焦平面处。其中,离焦光栅阵列、微透镜阵列均采用2×2排布,微透镜阵列焦距为50mm,单个子孔径尺寸为1200μm,光波波长为632.8nm。The present invention is a defocus grating array wavefront sensor. In the first embodiment, a 2×2 defocus grating array wavefront sensor is used. Its optical structure is shown in FIG. 1. The defocus grating performs phase modulation on the input wavefront. The defocus grating array 1 is located in front of the microlens array 2 , and the CCD 3 is located at the focal plane of the microlens array 2 . Among them, the defocus grating array and the microlens array are arranged in 2×2, the focal length of the microlens array is 50mm, the size of a single sub-aperture is 1200μm, and the wavelength of the light wave is 632.8nm.

具体实施分为两个部分,如下:The specific implementation is divided into two parts, as follows:

1.离焦光栅对入射波前的调制作用表示为:1. The modulation effect of the defocus grating on the incident wavefront is expressed as:

离焦光栅为二元相位调制板,相移函数φ'(x,y)满足如下关系式:The defocus grating is a binary phase modulation plate, and the phase shift function φ'(x, y) satisfies the following relationship:

Figure BDA0003739754690000031
Figure BDA0003739754690000031

其中,x,y为离焦光栅的坐标,x0为离焦光栅的偏移量。Among them, x, y are the coordinates of the defocused grating, and x 0 is the offset of the defocused grating.

实施例一中,待测波前包含前16阶的本征函数像差,光瞳为长度为a宽度为b的矩形,且矩形区域上的本征函数的解析表达式为:In the first embodiment, the wavefront to be measured includes the first 16 eigenfunction aberrations, the pupil is a rectangle with a length a and a width b, and the analytical expression of the eigenfunction on the rectangular area is:

Figure BDA0003739754690000032
Figure BDA0003739754690000032

对应的本征值为

Figure BDA0003739754690000033
其中m,n=0、1、2、3…The corresponding eigenvalues are
Figure BDA0003739754690000033
where m,n=0, 1, 2, 3...

则输入波前可以由本征函数表示为:Then the input wavefront can be represented by the eigenfunction as:

Figure BDA0003739754690000034
Figure BDA0003739754690000034

其中,Wi表示第i阶的本征函数的相差,ai为第i阶的本征函数的系数,波前总共由L阶的本征函数表示。Among them, Wi represents the phase difference of the eigenfunction of the ith order, a i is the coefficient of the eigenfunction of the ith order, and the wavefront is represented by the eigenfunction of the L order in total.

待测波前通过离焦光栅对相位进行调制,并经过微透镜阵列后在其焦面的CCD上形成正负离焦的光强图像,离焦光栅阵列的相位图2所示,输入的待测波前如图3中(a)所示,焦面的CCD上形成正负离焦的光强图像如图3中(b)所示。The phase of the wavefront to be measured is modulated by the defocus grating, and after passing through the microlens array, positive and negative defocus light intensity images are formed on the CCD of its focal plane. The phase of the defocus grating array is shown in Figure 2. The measured wavefront is shown in Fig. 3(a), and the positive and negative defocused light intensity images formed on the CCD on the focal plane are shown in Fig. 3(b).

设第j个子孔径调制前的波前为φj(x,y),对子孔径内波前调制的相移函数为φm(x,y),则对第j个子孔径调制后的波前φ'j(x,y)为:Let the wavefront before the modulation of the jth sub-aperture be φ j (x, y), and the phase shift function of the modulation of the wave front in the sub-aperture is φ m (x, y), then the wavefront after modulation of the jth sub-aperture is φ'j (x,y) is:

φ'j(x,y)=φj(x,y)+φm(x,y) φ'j (x,y)= φj (x,y)+ φm (x,y)

2.基于离焦光栅阵列的波前探测器波前复原过程如下:2. The wavefront recovery process of the wavefront detector based on the defocus grating array is as follows:

步骤1:求解子波前信息Step 1: Solve for the wavefront information

输入的子波前可以由本征函数表示为:The input wavefront can be represented by the eigenfunction as:

Figure BDA0003739754690000041
Figure BDA0003739754690000041

其中,Wi表示第i阶的本征函数的相差,aij为第i阶的本征函数的系数,子波前总共由k阶的本征函数表示。Wherein, Wi represents the phase difference of the eigenfunction of the ith order, a ij is the coefficient of the eigenfunction of the ith order, and the wavefront is represented by the eigenfunction of the kth order in total.

经过调制后的波前,在单个子孔径的焦面上可以得到传输距离为δz的正负离焦面的光强信息I+,I-。定义S(x,y)为:After the modulated wavefront, the light intensity information I + , I- of the positive and negative defocal planes with a transmission distance of δz can be obtained on the focal plane of a single sub-aperture. Define S(x,y) as:

Figure BDA0003739754690000042
Figure BDA0003739754690000042

其中,k为波矢,I0为近场光强信息。Among them, k is the wave vector, and I 0 is the near-field light intensity information.

则根据曲率波前传感器的波前复原算法,即可复原单个子孔径内的波前信息,具体的方法为:Then, according to the wavefront recovery algorithm of the curvature wavefront sensor, the wavefront information in a single sub-aperture can be recovered. The specific method is as follows:

Figure BDA0003739754690000043
Figure BDA0003739754690000043

其中,Wi为第i阶本征函数基,ci=∫σWi(x,y)dσ为第i阶本征函数基的归一化系数,λi对应于本征函数基的本征值。Among them, Wi is the eigenfunction basis of the ith order, c i =∫ σ Wi ( x,y)dσ is the normalization coefficient of the eigenfunction basis of the ith order, and λ i corresponds to the basis of the eigenfunction basis eigenvalue.

由此便可得到每一阶本征函数的系数,进而得到单个子孔径内的波前信息。From this, the coefficient of each order eigenfunction can be obtained, and then the wavefront information in a single sub-aperture can be obtained.

步骤2:求解子孔径的系数向量Step 2: Solve the coefficient vector for the subaperture

根据步骤1中的方法,可以求解每个子孔径内的波前信息,将单个子孔径内的波前信息表示为ai,j,i表示第几阶信息,j表示子孔径编号,则第j个子孔径中的波前信息表示为

Figure BDA0003739754690000044
将各个子孔径的系数按列排列为一个列向量A,则A矩阵的转置矩阵
Figure BDA0003739754690000051
计算得到的子孔径的系数向量A将用于步骤4中的全口径波前重构。According to the method in step 1, the wavefront information in each sub-aperture can be solved, and the wavefront information in a single sub-aperture can be expressed as a i,j , where i is the order information, j is the sub-aperture number, then the jth The wavefront information in the sub-apertures is expressed as
Figure BDA0003739754690000044
Arrange the coefficients of each sub-aperture into a column vector A, then the transpose matrix of the A matrix
Figure BDA0003739754690000051
The calculated coefficient vector A of the sub-aperture will be used for the full-aperture wavefront reconstruction in step 4.

步骤3:求解全口径的复原矩阵Step 3: Solve the full aperture recovery matrix

求解全口径的本征函数的复原矩阵,有效子孔径的数目为j,则j=1、2、3…n。同上,将每个子孔径的波前信息表示为Aj。全口径的波前总共用m阶的本征函数的来表示,将不同阶数m的本征函数对应的每个子波前的信息Ai表示为Ai,j,i=1、2、3…m,则由Ai,j构成的全口径的复原矩阵D为:To solve the restoration matrix of the eigenfunction of the full aperture, the number of effective sub-apertures is j, then j=1, 2, 3...n. As above, denote the wavefront information of each sub-aperture as A j . The full aperture wavefront is represented by the eigenfunctions of order m in total, and the information A i of each wavefront corresponding to the eigenfunctions of different orders m is represented as A i,j , i=1, 2, 3 ...m, then the full aperture restoration matrix D composed of A i, j is:

Figure BDA0003739754690000052
Figure BDA0003739754690000052

此处,j对应离焦光栅阵列的波前探测器中的有效子孔径序数j=1、2、3…n,i对应本征函数的阶序数i=1、2、3…m。每一个矩阵元素都包含有一个子波前的本征函数的系数列向量Aj,D是一个6×i行j列的矩阵。Here, j corresponds to the effective sub-aperture number j=1, 2, 3...n in the wavefront detector of the defocus grating array, and i corresponds to the order number i=1, 2, 3...m of the eigenfunction. Each matrix element contains a column vector A j of coefficients of the eigenfunction of the wavefront, and D is a matrix of 6×i rows and j columns.

将全口径的复原矩阵D,换算为伪逆矩阵D+,并将伪逆矩阵D+的数据存入波前复原的计算机系统之中。将用于步骤4中的全口径波前重构。Convert the full aperture restoration matrix D into a pseudo-inverse matrix D + , and store the data of the pseudo-inverse matrix D + into the computer system for wavefront restoration. Will be used for full aperture wavefront reconstruction in step 4.

步骤4:全口径波前重构Step 4: Full Aperture Wavefront Reconstruction

将全口径的波前φ0(x,y)表示为由全口径复原矩阵计算得到的各阶本征函数的系数的线性叠加模式:The full-aperture wavefront φ 0 (x, y) is expressed as a linear superposition pattern of coefficients of various order eigenfunctions calculated from the full-aperture restoration matrix:

Figure BDA0003739754690000053
Figure BDA0003739754690000053

其中,a'i为计算得到的每一项本征函数的系数,i=1、2、3…m。Among them, a' i is the coefficient of each eigenfunction obtained by calculation, i=1, 2, 3...m.

离焦光栅阵列的波前探测器探测全口径的波前φ0(x,y)时,根据步骤2中利用曲率波前传感器的复原算法计算得到的系数向量和步骤3中计算得到的复原矩阵D的伪逆矩阵D+满足如下关系式:When the wavefront detector of the defocus grating array detects the full-aperture wavefront φ 0 (x, y), according to the coefficient vector calculated by the restoration algorithm of the curvature wavefront sensor in step 2 and the restoration matrix calculated in step 3 The pseudo-inverse matrix D + of D satisfies the following relation:

a'i=D+·Aa' i =D + ·A

波前复原的计算根据上述的数学关系式便可以计算重构出全口径的波前信息φ0(x,y)。The calculation of the wavefront restoration can calculate and reconstruct the full-aperture wavefront information φ 0 (x, y) according to the above mathematical relationship.

最终输出复原的波前如图4中(a)(RMS=2.6344rad,PV=10.4121rad)所示,波前复原残差如图4中(b)(RMS=0.0308rad,PV=0.1161rad)所示,其RMS值和PV值分别是输入波前的1.17%、1.12%,可以较好的复原波前,为了突出本发明的优势,实施例一中在同等条件下,利用模式法对该输入波前进行复原,复原结果如图5所示,其中(a)为模式法复原的波前(RMS=0.4434rad,PV=1.9403rad),(b)为模式法波前复原残差(RMS=2.5154rad,PV=8.5746rad),RMS值和PV值分别是输入波前的567.3%、441.92%,已无法有效复原该畸变波前。以上结果充分证明:本发明在2×2子孔径条件下可以高精度复原波前,波前复原精度几乎不受子孔径数目影响。The final output restored wavefront is shown in Figure 4 (a) (RMS=2.6344rad, PV=10.4121rad), and the wavefront restoration residual is shown in Figure 4 (b) (RMS=0.0308rad, PV=0.1161rad) As shown, the RMS value and PV value are respectively 1.17% and 1.12% of the input wavefront, which can restore the wavefront well. In order to highlight the advantages of the present invention, in the first embodiment, under the same conditions, the model method is used The input wavefront is restored, and the restoration result is shown in Figure 5, where (a) is the wavefront restored by the mode method (RMS=0.4434rad, PV=1.9403rad), (b) is the residual error of the mode method restoration (RMS =2.5154rad, PV=8.5746rad), the RMS value and the PV value are respectively 567.3% and 441.92% of the input wavefront, and the distorted wavefront cannot be effectively restored. The above results fully prove that the present invention can restore the wavefront with high precision under the condition of 2×2 sub-apertures, and the wavefront restoration accuracy is hardly affected by the number of sub-apertures.

以上所述,仅为本发明中的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉该技术的人在本发明所揭露的技术范围内,可理解想到的变换或替换,都应涵盖在本发明的包含范围之内。The above is only a specific embodiment of the present invention, but the protection scope of the present invention is not limited to this, any person familiar with the technology can understand the transformation or replacement that comes to mind within the technical scope disclosed by the present invention, All should be included within the scope of the present invention.

Claims (4)

1. A wave front sensor based on an out-of-focus grating array is characterized in that: the wavefront sensor adopts an out-of-focus grating array to perform aperture segmentation and phase modulation on a light beam to be detected to form a +/-1-level diffraction array, and then utilizes a micro-lens array to focus the +/-1-level diffraction array to form a +/-1-level out-of-focus light spot array; the wavefront sensor mainly comprises three parts, namely a defocusing grating array, a micro-lens array and a photoelectric detector;
the gratings in the defocusing grating array are phase-type defocusing gratings and are used for carrying out aperture segmentation on the light beam to be detected so as to form a plurality of sub apertures and carrying out phase modulation on wavefront in the sub apertures; let the wavefront before single subaperture modulation be phi j (x, y), j is the sub-aperture number, and the phase shift function for the wavefront modulation within the sub-aperture is φ m (x, y), wavefront φ 'after modulating the jth sub-aperture' j (x, y) is:
φ′ j (x,y)=φ j (x,y)+φ m (x,y),
wherein x and y are coordinates of the defocused grating;
the micro lens array is positioned behind the defocusing grating array and corresponds to the defocusing grating array one by one, and focuses the light beams modulated by the defocusing grating array to form a +/-1-level defocusing light spot array on a focal plane of the micro lens array;
the photoelectric detector is positioned on the focal plane of the micro-lens array and is used for detecting +/-1-level defocusing spot image data in each sub-aperture, and the defocusing spot image data is used for calculating high-order aberration information in the sub-apertures and reconstructing full-aperture high-resolution wavefront by using a curvature wavefront restoration method.
2. The wavefront sensor based on an out-of-focus grating array of claim 1, wherein: the defocused grating array is larger than or equal to 2 x 2.
3. The wavefront sensor based on an out-of-focus grating array of claim 1, wherein: the curvature wave front restoration method is an eigen function wave front restoration method based on two +/-1-level defocused light spot images, and comprises all methods for extracting wave front information based on +/-1-level defocused light spots.
4. The wavefront sensor based on an out-of-focus grating array of claim 1, wherein: the photodetector may be a CCD, CMOS, or other array type detector.
CN202210809078.2A 2022-07-11 2022-07-11 A wavefront sensor based on defocused grating array Active CN115183886B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210809078.2A CN115183886B (en) 2022-07-11 2022-07-11 A wavefront sensor based on defocused grating array

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210809078.2A CN115183886B (en) 2022-07-11 2022-07-11 A wavefront sensor based on defocused grating array

Publications (2)

Publication Number Publication Date
CN115183886A true CN115183886A (en) 2022-10-14
CN115183886B CN115183886B (en) 2024-10-18

Family

ID=83517125

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210809078.2A Active CN115183886B (en) 2022-07-11 2022-07-11 A wavefront sensor based on defocused grating array

Country Status (1)

Country Link
CN (1) CN115183886B (en)

Citations (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2665955A1 (en) * 1985-11-20 1992-02-21 Onera (Off Nat Aerospatiale) Opto-electronic analyser for wave surfaces with microlens mosaics
US5233174A (en) * 1992-03-11 1993-08-03 Hughes Danbury Optical Systems, Inc. Wavefront sensor having a lenslet array as a null corrector
US6057913A (en) * 1997-02-05 2000-05-02 Mems Optical Inc. Compact shearing wavefront sensor and method
CN1607379A (en) * 2003-10-15 2005-04-20 中国科学院光电技术研究所 Hartmann wavefront sensor for pulse light beam quality detection based on microprism array
US7232999B1 (en) * 2003-07-30 2007-06-19 Kestrel Corporation Laser wavefront characterization
CN101278874A (en) * 2007-12-28 2008-10-08 中国科学院光电技术研究所 A Transmissive Hartmann Measuring Instrument for Intraocular Lens Aberration
US20080265130A1 (en) * 2005-02-23 2008-10-30 Tristan Colomb Wave Front Sensing Method and Apparatus
US20090180115A1 (en) * 2008-01-11 2009-07-16 Wilson Daniel W Single-lens computed tomography imaging spectrometer and method of capturing spatial and spectral information
CN101701847A (en) * 2009-11-27 2010-05-05 中国科学院光电技术研究所 Wide Dynamic Range Imaging System Based on Grating and CCD Imaging Detector
US20100141959A1 (en) * 2008-12-09 2010-06-10 Kuechel Michael Two grating lateral shearing wavefront sensor
US20110134436A1 (en) * 2009-04-29 2011-06-09 Adrian Podoleanu Method for depth resolved wavefront sensing, depth resolved wavefront sensors and method and apparatus for optical imaging
CN102331303A (en) * 2011-08-05 2012-01-25 中国科学院光电技术研究所 A grating-based phase difference wavefront sensor
CN102419213A (en) * 2011-12-20 2012-04-18 四川大学 Hartmann wavefront sensor based on diffraction grating arrays
CN102507020A (en) * 2011-11-01 2012-06-20 南京理工大学 Microlens array-based synchronized phase-shifting interference test method and test device
CN102607820A (en) * 2012-04-05 2012-07-25 中国科学院光电技术研究所 Method for measuring focal length of micro-lens array
CN202420688U (en) * 2011-12-20 2012-09-05 四川大学 Hartmann wavefront sensor based on diffraction grating arrays
CN102914373A (en) * 2012-11-20 2013-02-06 天津理工大学 Hartmann wave-front sensor based on micro-cylindrical lens array
CN104198164A (en) * 2014-09-19 2014-12-10 中国科学院光电技术研究所 Focus detection method based on Hartmann wavefront detection principle
US20190361226A1 (en) * 2018-05-25 2019-11-28 Canon Kabushiki Kaisha Wavefront sensor, wavefront measurement apparatus, method of manufacturing optical element, and method of manufacturing optical system
CN111561864A (en) * 2020-04-29 2020-08-21 西安电子科技大学 A point diffraction digital holographic microscope device and method based on polarization grating
CN111998962A (en) * 2020-08-19 2020-11-27 中国科学院光电技术研究所 Hartmann wavefront sensor based on array type binary phase modulation
CN112484866A (en) * 2020-11-24 2021-03-12 中国科学院光电技术研究所 Wavefront restoration method based on shack-Hartmann wavefront sensor
CN112484865A (en) * 2020-11-20 2021-03-12 中国科学院光电技术研究所 Real-time polarization modulation Hartmann-shack wavefront detection device
CN112629680A (en) * 2020-12-07 2021-04-09 中国科学院长春光学精密机械与物理研究所 Aviation camera focus detection device and method based on shack-Hartmann wavefront sensing

Patent Citations (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2665955A1 (en) * 1985-11-20 1992-02-21 Onera (Off Nat Aerospatiale) Opto-electronic analyser for wave surfaces with microlens mosaics
US5233174A (en) * 1992-03-11 1993-08-03 Hughes Danbury Optical Systems, Inc. Wavefront sensor having a lenslet array as a null corrector
US6057913A (en) * 1997-02-05 2000-05-02 Mems Optical Inc. Compact shearing wavefront sensor and method
US7232999B1 (en) * 2003-07-30 2007-06-19 Kestrel Corporation Laser wavefront characterization
CN1607379A (en) * 2003-10-15 2005-04-20 中国科学院光电技术研究所 Hartmann wavefront sensor for pulse light beam quality detection based on microprism array
US20080265130A1 (en) * 2005-02-23 2008-10-30 Tristan Colomb Wave Front Sensing Method and Apparatus
CN101278874A (en) * 2007-12-28 2008-10-08 中国科学院光电技术研究所 A Transmissive Hartmann Measuring Instrument for Intraocular Lens Aberration
US20090180115A1 (en) * 2008-01-11 2009-07-16 Wilson Daniel W Single-lens computed tomography imaging spectrometer and method of capturing spatial and spectral information
US20100141959A1 (en) * 2008-12-09 2010-06-10 Kuechel Michael Two grating lateral shearing wavefront sensor
US20110134436A1 (en) * 2009-04-29 2011-06-09 Adrian Podoleanu Method for depth resolved wavefront sensing, depth resolved wavefront sensors and method and apparatus for optical imaging
CN101701847A (en) * 2009-11-27 2010-05-05 中国科学院光电技术研究所 Wide Dynamic Range Imaging System Based on Grating and CCD Imaging Detector
CN102331303A (en) * 2011-08-05 2012-01-25 中国科学院光电技术研究所 A grating-based phase difference wavefront sensor
CN102507020A (en) * 2011-11-01 2012-06-20 南京理工大学 Microlens array-based synchronized phase-shifting interference test method and test device
CN102419213A (en) * 2011-12-20 2012-04-18 四川大学 Hartmann wavefront sensor based on diffraction grating arrays
CN202420688U (en) * 2011-12-20 2012-09-05 四川大学 Hartmann wavefront sensor based on diffraction grating arrays
CN102607820A (en) * 2012-04-05 2012-07-25 中国科学院光电技术研究所 Method for measuring focal length of micro-lens array
CN102914373A (en) * 2012-11-20 2013-02-06 天津理工大学 Hartmann wave-front sensor based on micro-cylindrical lens array
CN104198164A (en) * 2014-09-19 2014-12-10 中国科学院光电技术研究所 Focus detection method based on Hartmann wavefront detection principle
US20190361226A1 (en) * 2018-05-25 2019-11-28 Canon Kabushiki Kaisha Wavefront sensor, wavefront measurement apparatus, method of manufacturing optical element, and method of manufacturing optical system
CN111561864A (en) * 2020-04-29 2020-08-21 西安电子科技大学 A point diffraction digital holographic microscope device and method based on polarization grating
CN111998962A (en) * 2020-08-19 2020-11-27 中国科学院光电技术研究所 Hartmann wavefront sensor based on array type binary phase modulation
CN112484865A (en) * 2020-11-20 2021-03-12 中国科学院光电技术研究所 Real-time polarization modulation Hartmann-shack wavefront detection device
CN112484866A (en) * 2020-11-24 2021-03-12 中国科学院光电技术研究所 Wavefront restoration method based on shack-Hartmann wavefront sensor
CN112629680A (en) * 2020-12-07 2021-04-09 中国科学院长春光学精密机械与物理研究所 Aviation camera focus detection device and method based on shack-Hartmann wavefront sensing

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
PAUL M. BLANCHARD: "Simultaneous multiplane imaging with a distorted diffraction grating", 《APPLIED OPTICS》, vol. 38, no. 32, 10 November 1999 (1999-11-10), pages 6992 - 6999, XP000893705, DOI: 10.1364/AO.38.006692 *
张锐;杨金生;田雨;饶学军;: "焦面哈特曼传感器波前相位复原", 光电工程, no. 02, 15 February 2013 (2013-02-15), pages 36 - 43 *
李斐: "基于离焦光栅的最小二乘波前复原算法", 《中国激光》, vol. 37, no. 11, 30 November 2010 (2010-11-30), pages 2813 - 2820 *
邱学晶: "基于四象限离散相位调制的单帧图像深度学习 相位反演算法研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》, 15 August 2021 (2021-08-15), pages 138 - 419 *

Also Published As

Publication number Publication date
CN115183886B (en) 2024-10-18

Similar Documents

Publication Publication Date Title
CN102331303B (en) Phase difference wavefront sensor based on grating
CN106950195B (en) Programmable optical element and light field control system and method based on scattering medium
CN100580501C (en) Adaptive Optics Stellar Target Imaging System Based on Image Clearing Principle
CN102967380B (en) Hartmann wavefront sensor based on unit photosensitive detector array
CN102288305B (en) Adaptive optical system wavefront sensor and detection method thereof
CN111220361B (en) Method for measuring focal length of micro-lens array
CN104596650B (en) Hartmann wavefront sensor super-resolution wavefront restoration method
CN106338343B (en) Wavefront detection method based on Fresnel zone plate
CN106338823B (en) Phase inversion method based on mixed focal length Fresnel zone plate
US20060082882A1 (en) Achromatic imaging lens with extended depth of focus
CN112484866B (en) A Wavefront Restoration Method Based on Shack-Hartmann Wavefront Sensor
CN101285712B (en) Linear phase inversion wavefront sensor based on discrete light intensity measuring device
CN111829671A (en) High-resolution wavefront detection device and wavefront restoration method
CN106525255A (en) Phase difference wavefront sensor based on image compensation
CN113433688A (en) Interference imaging method and system based on micro-lens array and photonic integrated chip
CN114186664B (en) A method of pattern wavefront restoration based on neural network
CN113805335B (en) A method to reduce the influence of the chromatic aberration of the beam splitter on the wavefront detection accuracy of the phase-difference method
CN111998962B (en) A Hartmann Wavefront Sensor Based on Array Binary Phase Modulation
CN113670456B (en) A Wavefront Restoration Method Using a Hartmann Wavefront Sensor with Adjustable Spatial Resolution
Ge et al. Target-independent dynamic wavefront sensing method based on distorted grating and deep learning
CN115183886B (en) A wavefront sensor based on defocused grating array
CN103105235B (en) Method for measuring oblique aberration in sub-aperture of Hartmann wavefront sensor
CN101285711A (en) Linear phase inversion wavefront sensor based on area array CCD
CN113432731B (en) A Compensation Method in the Reconstruction Process of Grating Transverse Shear Interference Wavefront
CN112013973B (en) Fibonacci photon sieve based variable shear ratio four-wave shearing interferometer

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