CN113295287B - 一种针对光瞳动态强度分布的哈特曼子孔径减阈值方法 - Google Patents
一种针对光瞳动态强度分布的哈特曼子孔径减阈值方法 Download PDFInfo
- Publication number
- CN113295287B CN113295287B CN202110579014.3A CN202110579014A CN113295287B CN 113295287 B CN113295287 B CN 113295287B CN 202110579014 A CN202110579014 A CN 202110579014A CN 113295287 B CN113295287 B CN 113295287B
- Authority
- CN
- China
- Prior art keywords
- sub
- aperture
- threshold value
- hartmann
- weight
- 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 43
- 210000001747 pupil Anatomy 0.000 title claims abstract description 26
- 238000009826 distribution Methods 0.000 title claims abstract description 18
- 230000009467 reduction Effects 0.000 title claims description 6
- 238000001514 detection method Methods 0.000 claims abstract description 13
- 238000004364 calculation method Methods 0.000 claims abstract description 10
- 230000008569 process Effects 0.000 claims description 10
- 238000013528 artificial neural network Methods 0.000 claims description 2
- 238000013135 deep learning Methods 0.000 claims description 2
- 238000005192 partition Methods 0.000 claims description 2
- 239000011148 porous material Substances 0.000 claims description 2
- 238000003672 processing method Methods 0.000 claims description 2
- 230000004048 modification Effects 0.000 claims 1
- 238000012986 modification Methods 0.000 claims 1
- 230000008859 change Effects 0.000 abstract description 9
- 238000000605 extraction Methods 0.000 abstract description 5
- 238000009828 non-uniform distribution Methods 0.000 abstract description 2
- 230000009286 beneficial effect Effects 0.000 abstract 1
- 230000000694 effects Effects 0.000 abstract 1
- 230000003044 adaptive effect Effects 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 238000011084 recovery Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000007418 data mining Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 229910044991 metal oxide Inorganic materials 0.000 description 1
- 150000004706 metal oxides Chemical class 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000005693 optoelectronics Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
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
- 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)
- Spectroscopy & Molecular Physics (AREA)
- General Physics & Mathematics (AREA)
- Studio Devices (AREA)
Abstract
本发明提供一种针对光瞳动态强度分布的哈特曼子孔径减阈值方法，该方法结合不同时刻的子孔径光强信息动态设置阈值权重，降低光强时域动态变化导致的子孔径质心提取误差波动，可望改善哈特曼子孔径斜率计算准确性和稳定性。这种减阈值方法可降低光瞳强度分布时域非均匀性对子孔径质心计算的影响，同时也有助于光瞳强度空间非均匀分布下的子孔径质心精确计算，对于光瞳闪烁等情况下哈特曼波前探测的准确性和稳定性有重要作用。
Description
技术领域
本发明属于光学信息测量技术领域，具体涉及一种针对光瞳动态强度分布的哈特曼子孔径减阈值方法。
背景技术
哈特曼波前传感器是一种兼具高精度和高灵敏度的波前测量手段，已成功应用于自适应光学、光电系统集成、光学检测等领域。哈特曼波前传感器以子孔径阵列的形式，对光学波面进行微分，由各子孔径斜率解算完整波面的相位信息(即，波前信息)。因此，子孔径斜率误差对波前复原精度有直接影响。如何减小子孔径斜率计算误差，以确保波前复原精度，始终是这一领域所关注的重点问题。
由于哈特曼波前传感器采用的探测成像器件多采用CCD或CMOS相机，工作过程中不可避免地携带热噪声、散粒噪声及产生-复合噪声等多种误差因素的影响，因此，哈特曼子孔径质心提取及斜率计算必然涉及背景噪声的有效处理。目前，从算法和硬件层面已经发展了多种方法，如“哈特曼波前传感器质心测量精度优化方法”(专利号CN101055223)及“自适应系统中哈特曼波前传感器光斑质心最佳标定位置”(马晓燠,郑瀚清等,光电工程,2009,36(4))等。当待测波前光瞳强度分布出现时间和空间非均匀性时，同一时刻，各子孔径内的信噪比可能出现差异，而同一子孔径在不同时刻的信噪比也可能出现差异。这将给子孔径质心提取带来新的技术挑战。近年来先后发展了“一种采用分时曝光的哈特曼波前传感器”(专利公布号CN102607718A)和“一种适应背景变化点源目标波前探测的夏克-哈特曼波前传感器”(专利公布号CN1971222)等方法。但这些方法均使系统或算法复杂度有所增加。
本发明涉及方法中，针对传统的背景噪声处理方式可能导致斜率计算误差，提出具备动态适应能力的减阈值方法。该方法过程更为简单，可有效降低光瞳强度分布时域非均匀性对子孔径质心计算的影响，同时也有助于光瞳强度空间非均匀分布下的子孔径质心精确计算，可望在光瞳闪烁等情况下波前探测中实现更高精度和稳定性。
发明内容
本发明要解决的技术问题是：动态强度变化光瞳的时间非均匀性和空间非均匀性将导致哈特曼波前传感器在探测过程中各子孔径信噪比动态变化。针对这一问题，提出一种具备动态适应能力的减阈值方法，确保动态强度分布光瞳的高精度波前探测或复原控制。
本发明要解决其技术问题所采用的技术方案是：一种针对光瞳动态强度分布的哈特曼子孔径减阈值方法，本方法以光瞳动态强度分布下的高精度质心提取为目标，提出动态阈值权重技术，通过对光瞳时间/空间非均匀性的数理分析，确定阈值权重变更条件和变更方法，减少子孔径信噪比动态变化导致的误差起伏，提高哈特曼波前探测或复原控制的精度。
具体实现步骤如下：
步骤(1)、采集哈特曼传感器子图像前，先设定各子孔径的阈值权重初始值N(j,I,0)；
其中，N为阈值权重，j为子孔径序号，I为子孔径光强峰值，t为采集样本数。
子孔径阈值权重初始值N(j,I,0)，可有如下几种设定方式：
1)统一设定为某一特定值，即，N(j,I,0)＝N0；
2)分子孔径或分区域设定为固定值N(j,I,0)＝f(j)；
3)分子孔径或分区域设定为子孔径光强的统一权值N(j,I,0)＝f(I)；
4)分子孔径或分区域设定为子孔径光强的不同权值N(j,I,0)＝f(j,I)。
步骤(2)、开始采集哈特曼子图像，按照预定的阈值动态更新规则，逐个子孔径进行光强判断。符合阈值权重更新条件则相应更改阈值权重N(j,I,t)，否则阈值权重N(j,I,t)保持不变；
权重阈值更新条件与采集样本数t直接相关。该条件可以只与t相关，也可以与子孔径序号、光强中的一种或两种相关。
权重阈值更新条件，可以结合系统特性人工设定，也可以采用深度学习、模型辨识、神经网络等多种智能化处理方法设定。
步骤(3)、各子孔径根据相应的阈值权重，减阈值后计算质心和子孔径斜率，完成波前探测。直至采集流程结束。
子孔径质心和斜率计算过程中，可并行质心修正、斜率加权等操作。
本发明的原理是：结合不同时刻的子孔径光强信息动态设置各子孔径阈值权重，各子孔径可望近似等效于固定信噪比，以此提升子孔径质心提取的精度和稳定性。
本发明与现有技术相比具有如下优点：本方法综合考虑了光瞳分布的时间/空间非均匀特性，可望提升哈特曼子孔径质心提取精度，并确保斜率计算的准确性和稳定性。本方法可提高哈特曼波前传感器的实用能力，使异形、闪烁等光瞳条件下的波前准确探测或复原控制成为可能。
附图说明
图1为本发明涉及的某一子孔径光强变化曲线图；
图2为本发明涉及的某一子孔径在光强动态变化时本方法相对于传统方法的改善效果；
图3为本发明的一种针对光瞳动态强度分布的哈特曼子孔径减阈值方法的流程图。
具体实施方式
下面结合附图以及具体实施方式进一步说明本发明。
如图3所示，本发明的一种针对光瞳动态强度分布的哈特曼子孔径减阈值方法的具体实施方法如下：
步骤(1)、采集哈特曼传感器子图像前，先设定各子孔径的阈值权重初始值N(j,I,0)＝0.3Ij，即，各子孔径阈值为各自光强峰值的30％；
步骤(2)、开始采集哈特曼子图像，按照预定的阈值动态更新规则，逐个子孔径进行光强判断。符合阈值权重更新条件则相应更改阈值权重N(j,I,t)，否则阈值权重N(j,I,t)保持不变；
采集过程中，某子孔径光强变化规律如图1所示，该子孔径光强随t有波动趋势。通过大量数据分析，设定如下的权值阈值变更规则：当子孔径光强峰值≥1500，权重设为0.3，当峰值＜1500，权重改为0.1。
步骤(3)、各子孔径根据相应的阈值权重，减阈值后计算质心和子孔径斜率，完成波前探测。直至采集流程结束。
本方法以光瞳动态强度分布下的高精度质心提取为目标，提出动态阈值权重技术，通过对光瞳时间/空间非均匀性的数理分析，结合数据挖掘、数据融合及实验测试等过程，确定阈值权重变更条件和变更方法，减少子孔径信噪比动态变化导致的误差起伏，提高哈特曼波前探测或复原控制的精度。
本实施例中，本发明涉及方法和传统方法的对比如图2所示。本发明涉及方法在子孔径光强动态变化时，计算的斜率更为平稳，计算斜率均方根值为0.120；而传统方法在光强起伏时斜率也相应地产生起伏，与光强变化有负相关的关系，计算斜率均方根值为0.236。本发明涉及方法计算的斜率稳定性有明显提升。
本发明说明书中未作详细描述的内容属于本领域专业技术人员公知的现有技术。
Claims (4)
1.一种针对光瞳动态强度分布的哈特曼子孔径减阈值方法，其特征在于，实现步骤如下：
步骤(1)、采集哈特曼传感器子图像前，先设定各子孔径的阈值权重初始值N(j,I,t)，t＝0，以下记为N(j,I,0)；N为阈值权重，j为子孔径序号，I为子孔径光强峰值，t为采集样本数；
步骤(2)、开始采集哈特曼子图像，按照预定的阈值动态更新规则，逐个子孔径进行光强判断，符合阈值权重更新条件则相应更改阈值权重N(j,I,t)，否则阈值权重N(j,I,t)保持不变；
步骤(3)、各子孔径根据相应的阈值权重，减阈值后计算质心和子孔径斜率，完成波前探测，直至采集流程结束；
其中，步骤(1)所述的子孔径阈值权重初始值N(j,I,0)，可有如下几种设定方式：
1)统一设定为某一特定值，即，N(j,I,0)＝N0；
2)分子孔径或分区域设定为固定值N(j,I,0)＝f(j)；
3)分子孔径或分区域设定为子孔径光强的统一权值N(j,I,0)＝f(I)；
4)分子孔径或分区域设定为子孔径光强的不同权值N(j,I,0)＝f(j,I)；
该方法结合不同时刻的子孔径光强信息动态设置各子孔径阈值权重，各子孔径可望近似等效于固定信噪比，以此提升子孔径质心提取的精度和稳定性。
2.根据权利要求1所述的一种针对光瞳动态强度分布的哈特曼子孔径减阈值方法，其特征在于：步骤(2)所述的阈值权重 更新条件与采集样本数t直接相关，该条件可以只与t相关，也可以与子孔径序号、光强中的一种或两种相关。
3.根据权利要求1所述的一种针对光瞳动态强度分布的哈特曼子孔径减阈值方法，其特征在于：步骤(2)所述的阈值权重 更新条件，可以结合系统特性人工设定，也可以采用深度学习、模型辨识、神经网络智能化处理方法设定。
4.根据权利要求1所述的一种针对光瞳动态强度分布的哈特曼子孔径减阈值方法，其特征在于：步骤(3)所述的子孔径质心和斜率计算过程中，可并行进一步操作：质心修正和/或斜率加权。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110579014.3A CN113295287B (zh) | 2021-05-26 | 2021-05-26 | 一种针对光瞳动态强度分布的哈特曼子孔径减阈值方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110579014.3A CN113295287B (zh) | 2021-05-26 | 2021-05-26 | 一种针对光瞳动态强度分布的哈特曼子孔径减阈值方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113295287A CN113295287A (zh) | 2021-08-24 |
| CN113295287B true CN113295287B (zh) | 2022-11-11 |
Family
ID=77325179
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110579014.3A Active CN113295287B (zh) | 2021-05-26 | 2021-05-26 | 一种针对光瞳动态强度分布的哈特曼子孔径减阈值方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113295287B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114252163B (zh) * | 2021-12-21 | 2023-12-22 | 中国科学院光电技术研究所 | 一种基于图像噪声去除的低信噪比子光斑波前复原方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5287165A (en) * | 1991-09-30 | 1994-02-15 | Kaman Aerospace Corporation | High sensitivity-wide dynamic range optical tilt sensor |
| CN101055223A (zh) * | 2007-04-26 | 2007-10-17 | 中国科学院光电技术研究所 | 哈特曼波前传感器质心测量精度优化方法 |
| CN104239740A (zh) * | 2014-09-26 | 2014-12-24 | 中国科学院光电技术研究所 | 一种基于哈特曼波前传感器的模式波前复原方法 |
| CN104596650A (zh) * | 2015-01-27 | 2015-05-06 | 中国科学院光电技术研究所 | 一种哈特曼波前传感器超分辨波前复原方法 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE10146643C1 (de) * | 2001-09-21 | 2003-08-14 | Eads Deutschland Gmbh | Verfahren zur Kalibrierung der Radarsignale an den Subaperturen der Antenne eines zweikanaligen SAR/MTI Radarsystems |
| US6924899B2 (en) * | 2002-05-31 | 2005-08-02 | Optical Physics Company | System for measuring wavefront tilt in optical systems and method of calibrating wavefront sensors |
| US20040227932A1 (en) * | 2003-02-13 | 2004-11-18 | Geunyoung Yoon | Large dynamic range shack-hartmann wavefront sensor |
| US7457545B2 (en) * | 2004-02-12 | 2008-11-25 | Northrop Grumman Corporation | Process for controlling a Hartmann wavefront sensor (WFS) in an adaptive optic (AO) system |
| CN106408019A (zh) * | 2016-09-14 | 2017-02-15 | 南京信息工程大学 | 强天光背景下的自适应光学天体目标的检测方法 |
| CN112782844B (zh) * | 2021-01-12 | 2023-03-31 | 中国科学院光电技术研究所 | 一种自适应光学系统稳定闭环控制方法 |
-
2021
- 2021-05-26 CN CN202110579014.3A patent/CN113295287B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5287165A (en) * | 1991-09-30 | 1994-02-15 | Kaman Aerospace Corporation | High sensitivity-wide dynamic range optical tilt sensor |
| CN101055223A (zh) * | 2007-04-26 | 2007-10-17 | 中国科学院光电技术研究所 | 哈特曼波前传感器质心测量精度优化方法 |
| CN104239740A (zh) * | 2014-09-26 | 2014-12-24 | 中国科学院光电技术研究所 | 一种基于哈特曼波前传感器的模式波前复原方法 |
| CN104596650A (zh) * | 2015-01-27 | 2015-05-06 | 中国科学院光电技术研究所 | 一种哈特曼波前传感器超分辨波前复原方法 |
Non-Patent Citations (1)
| Title |
|---|
| Scene-based nonuniformity correction algorithm based on optical flow;Chen Peng 等;《Display Technologies III》;20130301;第8643卷;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113295287A (zh) | 2021-08-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113295287B (zh) | 一种针对光瞳动态强度分布的哈特曼子孔径减阈值方法 | |
| KR100780464B1 (ko) | 별센서에서의 별 중심찾기 방법 | |
| NO316849B1 (no) | Adaptiv ikke-uniform kompensasjonsalgoritme | |
| CN104101297B (zh) | 一种基于光电观测的空间物体尺度获取方法 | |
| CN111985143B (zh) | 一种基于Zernike多项式分解的全视场望远镜主动准直方法 | |
| CN103916599B (zh) | 一种针对遥感相机外景成像的快速调焦方法 | |
| CN105203213B (zh) | 一种计算复合波前传感自适应光学系统复原电压的方法 | |
| CN106989689A (zh) | 大口径平面光学元件面形的子孔径拼接检测技术及装置 | |
| Plazas et al. | Transverse electric fields' effects in the Dark Energy Camera CCDs | |
| NL2010457A (en) | Hartmann wavefront measuring instrument adapted for non-uniform light illumination. | |
| CN102761360A (zh) | 一种光信号处理方法及装置 | |
| CN114186664B (zh) | 一种基于神经网络的模式波前复原方法 | |
| CN102928196A (zh) | 自由曲面镜片的检测方法和装置 | |
| CN109872315B (zh) | 一种光学天文望远镜杂散光均匀性实时检测方法 | |
| CN113188671B (zh) | 一种基于交叉迭代自动位置矫正的波前检测方法 | |
| CN101581556A (zh) | 基于复合图形靶板的瞄具零位走动数字化测量装置 | |
| CN107421717B (zh) | 一种红外成像仪最小可探测温差自动测试方法及装置 | |
| CN109709460A (zh) | 一种日盲型紫外成像系统光子计数方法 | |
| Gudehus et al. | The design and construction of a charge-coupled device imaging system | |
| Chulani et al. | Preliminary performance results of the weighted Fourier phase slope centroiding method for Shack–Hartmann wavefront sensors obtained with the OOMAO simulator | |
| CN115356777A (zh) | 一种搜索天体测量型微引力透镜事件最大观测信号及星对最接近时刻的方法 | |
| CN114964522A (zh) | 一种基于光瞳映射模型的哈特曼波前复原方法 | |
| CN104007619B (zh) | 一种光瞳面光强分布的测量方法 | |
| Basden | Shack–Hartmann Wavefront Sensors | |
| CN113379617B (zh) | 一种流道内流体温度场图像的后处理方法 |
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 |