CN107402074B - 一种基于光学稀疏子孔径区域内标准正交多项式计算进行重构波前的方法 - Google Patents
一种基于光学稀疏子孔径区域内标准正交多项式计算进行重构波前的方法 Download PDFInfo
- Publication number
- CN107402074B CN107402074B CN201710563807.XA CN201710563807A CN107402074B CN 107402074 B CN107402074 B CN 107402074B CN 201710563807 A CN201710563807 A CN 201710563807A CN 107402074 B CN107402074 B CN 107402074B
- Authority
- CN
- China
- Prior art keywords
- aperture region
- sparse sub
- orthogonal
- sparse
- orthogonal polynomial
- 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 24
- 230000003287 optical effect Effects 0.000 title abstract description 5
- 238000004364 calculation method Methods 0.000 title abstract description 4
- 238000005259 measurement Methods 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 238000001514 detection method Methods 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 210000001747 pupil Anatomy 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/17—Function evaluation by approximation methods, e.g. inter- or extrapolation, smoothing, least mean square method
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- Algebra (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Databases & Information Systems (AREA)
- Exposure And Positioning Against Photoresist Photosensitive Materials (AREA)
- Testing Of Optical Devices Or Fibers (AREA)
Abstract
本发明公开了一种基于光学稀疏子孔径区域内标准正交多项式计算进行重构波前的方法。首先把第一个Zernike多项式作为稀疏子孔径区域内正交多项式的第一个,将第i个Zernike多项式去除其在稀疏子孔径区域内正交多项式的前i‑1上的投影,得到的残差作为稀疏子孔径区域内正交多项式的第i个,并对其进行归一化处理直到得到N个稀疏子孔径区域内标准正交基,然后,利用N个稀疏子孔径区域内标准正交基对波前相位数据进行重构。本发明解决了圆Zernike在稀疏子孔径区域的非正交性问题。
Description
技术领域
本发明属于波前检测领域,具体涉及一种基于光学稀疏子孔径区域内标准正交多项式计算进行重构波前的方法。
背景技术
在光学元件的检测技术中,对离散化的波前数据进行拟合分析从而得到所关注的信息是数据处理阶段的重要部分,其本质是利用合适的数学波前函数来拟合离散的采样数据,拟合的波前与实际波前残差尽可能小。模式法是基于待求波面相位建立正交函数系的解析模型,用最小二乘法求出模型的各项系数。大多数情况下,检测的光学元件具有圆形光瞳或通光孔,反射或透射的波前总是趋于光滑且连续的,所以,常采用模式法进行检测,在圆域内选用Zernike多项式作为基底函数进行波前数据拟合,因为它在单位圆上加权正交,函数正交的性质使其系数相互独立,有利于消除偶然误差的影响,并且Zernike系数能与Seidel像差建立联系。
Zernike多项式是一组定义在单位圆连续区域内的完备正交基。在圆域进行采样时,常采用Zernike多项式进行波前重构。本发明公开了当数据的采样区域不再是圆域,而是稀疏子孔径阵列时,在圆域Zernike的基础上构造在稀疏子孔径区域内的正交的基底函数。
发明内容
本发明要解决在稀疏子孔径区域采样情况下,采用模式法重构波前选用的正交基底函数问题。本发明是为了克服上述圆Zernike多项式的非正交性,提供了一种稀疏子孔径区域内标准正交多项式计算方法。
本发明采用的技术方案为:一种基于光学稀疏子孔径区域内标准正交多项式计算进行重构波前的方法,按以下步骤实现:
步骤一,确定所要计算的稀疏子孔径区域内正交多项式的个数N,即波前重构的阶数;
步骤二,同样取N个单位圆内归一化的Zernike多项式,表示为{Zi(x,y)};
步骤三,将第一个Zernike多项式作为稀疏子孔径区域内正交多项式的第一个;
步骤四,将第二个Zernike多项式去除其在稀疏子孔径区域内正交多项式的第一个上的投影,得到的残差作为稀疏子孔径区域内正交多项式的第二个;
步骤五,将第三个Zernike多项式去除其在稀疏子孔径区域内正交多项式的第一个和第二个上的投影,得到的残差作为稀疏子孔径区域内正交多项式的第三个;
步骤六,以此类推,将N个Zernike多项式正交得到N个稀疏子孔径区域内正交多项式,并对其进行归一化处理得到N个稀疏子孔径区域内标准正交基,利用N个稀疏子孔径区域内标准正交基对波前相位数据进行重构。
进一步地,步骤一中所述所要计算的稀疏子孔径区域内正交多项式的个数N为任意值,都可计算。
进一步地,步骤四具体过程为:
取第二个Zernike多项式Z2,去除其在稀疏子孔径区域内正交多项式的第一个上的投影,得到的残差作为稀疏子孔径区域内正交多项式的第二个S2,即其中,<>表示在稀疏子孔径区域的内积。
进一步地,步骤五具体过程为:
取第三个Zernike多项式Z3,减去其在稀疏子孔径区域内正交多项式的第一个和第二个上的投影,得到的残差作为稀疏子孔径区域内正交多项式的第三个S3,即其中,<>表示在稀疏子孔径区域的内积。
进一步地,步骤六具体过程为:
取第i个Zernike多项式Zi,减去其在稀疏子孔径区域内正交多项式的第1个、第2个…第i-1个上的投影,得到的残差作为稀疏子孔径区域内正交多项式的第i个Si,即以上步骤计算出的{Si(x,y)},做归一化处理,即其中,<>表示在稀疏子孔径区域的内积。
进一步地,可以计算出任意稀疏子孔径阵列的基底函数表达式,得到一组新的正交多项式多项式,记为{Si(x,y)},利用N个稀疏子孔径区域内标准正交基对波前相位数据进行重构。
本发明与现有技术相比的优点在于:
(1)Zernike多项式仅能作为圆域内的正交基。对于稀疏子孔径区域的光学***,Zernike多项式则不具有正交性,本发明解决了圆Zernike在稀疏子孔径区域的非正交性问题。
(2)圆Zernike在稀疏子孔径区域的非正交性会造成模式耦合和模式混淆误差,本发明有效地消除了模式重构法中的这两种误差。
附图说明
图1为本发明一种基于光学稀疏子孔径区域内标准正交多项式计算进行重构波前的方法的流程图;
图2为选用的稀疏子孔径阵列。
具体实施方式
下面结合附图以及具体实施方式进一步说明本发明。
本实施方式的一种稀疏子孔径区域内标准正交多项式计算方法按以下步骤实现:
步骤一,确定所要计算的稀疏子孔径区域内正交多项式的个数N,即波前重构的阶数;
步骤二,同样取N个单位圆内归一化的Zernike多项式,表示为{Zi(x,y)};
步骤三,将第一个Zernike多项式作为稀疏子孔径区域内正交多项式的第一个;
步骤四,将第二个Zernike多项式去除其在稀疏子孔径区域内正交多项式的第一个上的投影,得到的残差作为稀疏子孔径区域内正交多项式的第二个;
步骤五,将第三个Zernike多项式去除其在稀疏子孔径区域内正交多项式的第一个和第二个上的投影,得到的残差作为稀疏子孔径区域内正交多项式的第三个;
步骤六,以此类推,将N个Zernike多项式正交得到N个稀疏子孔径区域内正交多项式,并对其进行归一化处理得到N个稀疏子孔径区域内标准正交基,记为{Si(x,y)},利用N个稀疏子孔径区域内标准正交基对波前相位数据进行重构。
实施例:
1、波前重构通常分为区域法和模式法,区域法通过相位斜率方向的测量来拟合特定子孔径上的波前。模式法中将全孔径内的波前相位展开成不同的模式,即正交多项式组,然后用全孔径内的测量数据去求解各模式的系数,得到完整的波前展开式。Zernike多项式是一组定义在单位圆连续区域内的完备正交基,记为{Zi(x,y)},在圆域进行采样时,模式中常采用Zernike多项式对圆域内的波前数据进行拟合。但是,当在稀疏子孔径区域进行采样时,就需要采用稀疏子孔径区域内标准正交多项式对采样区域内的波前数据进行拟合。
本实施例对所选用的稀疏子孔径阵列,如图2所示,利用本发明所计算出的9项标准正交多项式如表1所示。
表1选用的稀疏子孔径阵列所计算出的9项标准正交多项式表
本发明中涉及到的本领域公知技术未详细阐述。
Claims (1)
1.一种基于光学稀疏子孔径区域内标准正交多项式计算进行重构波前的方法,其特征在于:它按以下步骤实现:
步骤一,确定所要计算的稀疏子孔径区域内正交多项式的个数N,即波前重构的阶数;
步骤二,同样取N个单位圆内归一化的Zernike多项式,表示为{Zi(x,y)};
步骤三,将第一个Zernike多项式作为稀疏子孔径区域内正交多项式的第一个;
步骤四,将第二个Zernike多项式去除其在稀疏子孔径区域内正交多项式的第一个上的投影,得到的残差作为稀疏子孔径区域内正交多项式的第二个;
步骤五,将第三个Zernike多项式去除其在稀疏子孔径区域内正交多项式的第一个和第二个上的投影,得到的残差作为稀疏子孔径区域内正交多项式的第三个;
步骤六,将N个Zernike多项式正交得到N个稀疏子孔径区域内正交多项式,并对其进行归一化处理得到N个稀疏子孔径区域内标准正交基,利用N个稀疏子孔径区域内标准正交基对波前相位数据进行重构;
步骤一中所述所要计算的稀疏子孔径区域内正交多项式的个数N为任意值,都可计算;
步骤四具体过程为:
稀疏子孔径区域内正交多项式表示为{Si(x,y)};
取第二个Zernike多项式Z2,去除其在稀疏子孔径区域内正交多项式的第一个上的投影,得到的残差作为稀疏子孔径区域内正交多项式的第二个S2,即其中,<>表示在稀疏子孔径区域的内积;
步骤五具体过程为:
取第三个Zernike多项式Z3,减去其在稀疏子孔径区域内正交多项式的第一个和第二个上的投影,得到的残差作为稀疏子孔径区域内正交多项式的第三个S3,即:其中,<>表示在稀疏子孔径区域的内积;
步骤六具体过程为:
取第i个Zernike多项式Zi,减去其在稀疏子孔径区域内正交多项式的第1个、第2个…第i-1个上的投影,得到的残差作为稀疏子孔径区域内正交多项式的第i个Si,即:以上步骤计算出的{Si(x,y)},做归一化处理,即其中,<>表示在稀疏子孔径区域的内积;
可以计算出任意稀疏子孔径阵列的基底函数表达式,得到一组新的正交多项式,利用N个稀疏子孔径区域内标准正交基对波前相位数据进行重构。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710563807.XA CN107402074B (zh) | 2017-07-12 | 2017-07-12 | 一种基于光学稀疏子孔径区域内标准正交多项式计算进行重构波前的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710563807.XA CN107402074B (zh) | 2017-07-12 | 2017-07-12 | 一种基于光学稀疏子孔径区域内标准正交多项式计算进行重构波前的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107402074A CN107402074A (zh) | 2017-11-28 |
CN107402074B true CN107402074B (zh) | 2019-08-02 |
Family
ID=60405367
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710563807.XA Active CN107402074B (zh) | 2017-07-12 | 2017-07-12 | 一种基于光学稀疏子孔径区域内标准正交多项式计算进行重构波前的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107402074B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109613346B (zh) * | 2018-12-28 | 2020-11-13 | 中国科学院新疆天文台 | 反射面天线远场方向图的重构方法和*** |
CN111486974A (zh) * | 2020-04-22 | 2020-08-04 | 中国科学院上海光学精密机械研究所 | 任意孔径形状高阶自由曲面波前的重构方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104089583A (zh) * | 2014-06-24 | 2014-10-08 | 北京空间机电研究所 | 一种光学***波前的子孔径反演方法 |
-
2017
- 2017-07-12 CN CN201710563807.XA patent/CN107402074B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104089583A (zh) * | 2014-06-24 | 2014-10-08 | 北京空间机电研究所 | 一种光学***波前的子孔径反演方法 |
Non-Patent Citations (2)
Title |
---|
基于稀疏孔径的波前重构算法;王棪 等;《航天返回与遥感》;20151031;第36卷(第5期);全文 * |
波面重构中非圆域Zernike正交基底构造方法;赵齐 等;《光学技术》;20170531;第43卷(第3期);第0-1节,图1 * |
Also Published As
Publication number | Publication date |
---|---|
CN107402074A (zh) | 2017-11-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Koleva et al. | ULySS: a full spectrum fitting package | |
CN107402074B (zh) | 一种基于光学稀疏子孔径区域内标准正交多项式计算进行重构波前的方法 | |
CN107490397B (zh) | 高精度自适应滤波fbg光谱快速寻峰方法 | |
CN111221123B (zh) | 一种基于模型的无波前传感器自适应光学校正方法 | |
CN110146258A (zh) | 一种Poisson noise模型下对扩展目标成像时的相位恢复方法 | |
Laidlaw et al. | Optimizing the accuracy and efficiency of optical turbulence profiling using adaptive optics telemetry for extremely large telescopes | |
CN110307962A (zh) | 检测光学***任意波长透射波前的方法 | |
CN105266784A (zh) | 一种基于示波法测量血压的方法和装置 | |
Martin et al. | William Herschel Telescope site characterization using the MOAO pathfinder CANARY on-sky data | |
JP2014055882A (ja) | 収差推定方法、プログラムおよび撮像装置 | |
CN1793805B (zh) | 一种消除Shack-Hartmann波前传感器模型误差的方法 | |
Bardou et al. | ELT-scale elongated LGS wavefront sensing: on-sky results | |
Martin et al. | Temporal convergence of phase spatial covariance matrix measurements in tomographic adaptive optics | |
JP5554234B2 (ja) | 光学系の波面の少なくとも一つの変形または光学系によって観察される対象を推定する方法および関連する装置 | |
CN106530351A (zh) | 一种基于图像传感器像素内量子效率获取的质心定位方法 | |
Zielinski et al. | Determination of the sampling factor in a phase-diverse phase retrieval algorithm | |
Le Bigot et al. | Branch point reconstructors for discontinuous light phase functions | |
Klee et al. | Digital holographic tomography for path-resolved turbulence measurements | |
CN103258129B (zh) | 准分子激光器干涉条纹的数据处理方法 | |
Chulani | Weighted fourier phase slope as a centroiding method in a shack-hartmann wavefront sensor for adaptive optics in astronomy | |
CN102818640B (zh) | 用包含4f***误差的图像频谱确定哈特曼阵列数的方法 | |
CN107977680A (zh) | 一种信号稀疏重构的方法 | |
Dai | Wavefront reconstruction from slope data within pupils of arbitrary shapes using iterative Fourier transform | |
Wang et al. | A transformation approach for aberration-mode coefficients of Walsh functions and Zernike polynomials | |
Ellerbroek | Adaptive optics without borders: performance evaluation in the infinite aperture limit |
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 |