CN109115350A - 一种基于压缩感知的波前探测系统 - Google Patents
一种基于压缩感知的波前探测系统 Download PDFInfo
- Publication number
- CN109115350A CN109115350A CN201811029190.4A CN201811029190A CN109115350A CN 109115350 A CN109115350 A CN 109115350A CN 201811029190 A CN201811029190 A CN 201811029190A CN 109115350 A CN109115350 A CN 109115350A
- Authority
- CN
- China
- Prior art keywords
- wavefront
- slope signal
- signals
- information
- matrix
- 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
Links
- 238000005070 sampling Methods 0.000 claims abstract description 22
- 230000006835 compression Effects 0.000 claims abstract description 6
- 238000007906 compression Methods 0.000 claims abstract description 6
- 239000011159 matrix material Substances 0.000 claims description 22
- 238000000034 method Methods 0.000 claims description 18
- 230000003044 adaptive effect Effects 0.000 abstract description 7
- 238000005259 measurement Methods 0.000 abstract description 6
- 230000000694 effects Effects 0.000 abstract description 3
- 238000001514 detection method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 238000011084 recovery Methods 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
- 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)
- Mechanical Light Control Or Optical Switches (AREA)
- Testing Of Optical Devices Or Fibers (AREA)
Abstract
本发明涉及自适应光学系统中波前信号的测量。一种基于压缩感知的波前探测器，包括变形镜、波前传感器和计算机，变形镜对波前信息进行下采样然后通过波前传感器得到压缩采样的波前斜率信息传输到计算机，计算机通过重构算法得到完整的波前信息。本发明降低了对硬件的高要求，解决了大量数据存储和传输的问题，提供一种低成本，效果优良的波前探测系统。
Description
技术领域
本发明涉及自适应光学系统中波前信号的测量。
背景技术
在自适应光学系统，为有效地克服湍流波前的干扰，提高图像的分辨率，在自适应光学系统中，波前探测器承担着对湍流动态变化造成的波前扰动实时测量的任务。随着观测图像的分辨率提高，波前探测器必须具备高帧频和微观探测能力，同时自适应光学系统的实时控制依托于波前测量的后期处理技术，这要求波前探测数据量尽可能的小以便传输和存储。在自适应光学系统中，如果变形镜的驱动单元增加，要求高空间分辨率的波前探测器与之匹配，波前探测器的分辨率的提高会使其波前探测能力大幅度下降，这对波前探测器有更加高要求的硬件需求。如果探测器的数量增多，波前探测的数据量增大，也带来大量的数据传输和存储的压力。这都给自适应光学系统实时控制及基于波前信息的数据后期处理提出了苛刻的要求。
传统的全采样方法，不但会产生大量的波前探测数据，且对数据的传输和存储有巨大压力，会降低信息的实时性；压缩感知理论为解决上述问题提供了新的途径。压缩感知理论充分利利用波前信息的可压缩性实现探测数据的压缩采集。在波前信息满足可以稀疏表示的理论前提下，将信号的稀疏性和采样相结合，可在低数据量的前提下，实现波前信息的精确重建。将压缩感知理论，对波前湍流的采样方法进行突破，在信号采集的过程中打破传统的Nyquist采样定理，在信号采集时直接对实时信号压缩采集，只将少量的有用信息进行传输。在压缩重建模块，利用采集到的少量的波前测量信息精确地重建出精确的波前相位。这样不但可以降低数据的传输与存储，同时提高了波前探测器的微光探测能力，降低对光电信号采集系统像素数目的要求，可以提高波前探测的帧频。因此，压缩感知技术在波前探测规程中可以发挥很重要的作用。
对基于压缩感知的波前探测系统，最主要的部分包括：1)设计能够很好稀疏表示波前斜率的稀疏字典。2)在变形镜加载采样矩阵。3)对斜率信号采样得到测量信号，4)测量信号进行重构，得到波前斜率信号，5)通过区域法复原波前信号，完成波前信号的探测。
发明内容
本发明所要解决的技术问题是：如何克服变形镜的驱动单元增加带来的波前探测器的高的硬件需求，以及大量数据传输和存储的问题，提供一种基于压缩感知技术的波前探测系统，降低对硬件的高要求和解决大量数据存储和传输的问题。本发明所采用的技术方案是：一种基于压缩感知的波前探测器，包括变形镜、波前传感器和计算机，变形镜对波前信息进行下采样然后通过波前传感器得到压缩采样的波前斜率信息传输到计算机，计算机通过重构算法得到完整的波前信息，完整的波前信息计算步骤如下：
步骤一、建立直角坐标系，以平行于波前传感器中心轴线的方向为Z轴方向，以波前传感器的微透镜阵列平面内的任意一条横向直线为X轴方向，以波前传感器的微透镜阵列平面内的任意一条纵向直线为Y轴方向，射到波前传感器的微透镜阵列平面内的波与Z轴的偏移量即与Z轴的夹角能通过该波与X方向夹角和该波与Y方向夹角确定，波前信号φ对应的X轴方向的斜率信号Gx，波前信号φ对应的Y轴方向的斜率信号Gy，波前传感器的微透镜阵列平面上有N×N个子孔径，通过波前信息的特征建立起能够稀疏表示出波前斜率信号的稀疏字典Ψ，使得波前斜率信号等于极少量的系数乘以稀疏字典Ψ，即Gx＝Ψα,Gy＝Ψβ，α,β表示对应X轴方向和Y轴方向的斜率信号的系数，Gx、Gy、α、β、Ψ都为一个N×N的矩阵，
斜率信号Gx表示为
斜率信号Gy表示为
其中α、β中不为零的数据有k个，k<<N×N，
波前信号φ具有低频多高频少的特点，利用这一特点，通过黄金分割法在频率域对频率分布进行分割，得到符合黄金分割法(见李娟娟等."大气湍流畸变波前斜率的稀疏分解."光电工程2(2018))的频率序列F＝[f1,f2,f3,...,fn]，将F代入
得到大小为N×N的矩阵，将该矩阵分割为的小块，得到N个块组成稀疏字典Ψ，其中m、n都为大于0小于等于N的正整数，u、v都为F中的某一个频率，j为虚数单位；
步骤二、变形镜上有与波前传感器子孔径一一对应的变形镜子孔径，通过变形镜的上加载的采样矩阵Φ对波前信息进行下采样，采样矩阵Φ为一个M×N的矩阵，矩阵中每个元素的取值为0或者1，其中M<N，其中1代表通，0代表断，波前传感器得到波前斜率信号的测量信号，X轴方向的测量信号CSx＝ΦGx，Y轴方向的测量信号CSy＝ΦGy，CSx和CSy都为一个M×N矩阵；
步骤三、由公式CSx＝ΦΨα＝Aα和CSy＝ΦΨβ＝Aβ，通过迭代权重最小二乘法(见Chartrand,Rick,and W.Yin."Iteratively reweighted algorithms for compressivesensing."(2008):3869-3872)，求解得到α、β，通过Gx＝Ψα,Gy＝Ψβ获得波前信号φ对应的X轴方向的斜率信号Gx和波前信号φ对应的Y轴方向的斜率信号Gy；
步骤四、通过对波前信号φ对应的X轴方向的斜率信号Gx和波前信号φ对应的Y轴方向的斜率信号Gy使用区域法复原出波前信号φ(区域法复原见张强,姜文汉,and许冰."用于Hartmann-Shack波前探测器的区域法算法研究."强激光与粒子束10.2(1998):0-0)。
本发明的有益效果是：本发明降低了对硬件的高要求，解决了大量数据存储和传输的问题，提供一种低成本，效果优良的波前探测系统。
附录说明
图1为本发明的结构示意图；
图2为本发明流程示意图；
图3是波前斜率信号的稀疏字典Ψ图；
图4是表示原始的波前信号图；
图5表示复原的波前信号图；
图6表示二者的误差值图
其中，1、大气湍流，2、变形镜，3、波前传感器，4、计算机，5、采样矩阵Φ。
具体实施方式
结合图1的结构示意图和图2的流程示意图，一种基于压缩感知的波前探测器，包括变形镜2、波前传感器3和计算机4，变形镜1对大气湍流1导致的波前信息进行下采样然后通过波前传感器3得到压缩采样的波前斜率信息传输到计算机4，计算机4通过重构算法得到完整的波前信息，完整的波前信息计算步骤如下：
步骤一、建立直角坐标系，以平行于波前传感器中心轴线的方向为Z轴方向，以波前传感器3的微透镜阵列平面内的任意一条横向直线为X轴方向，以波前传感器3的微透镜阵列平面内的任意一条纵向直线为Y轴方向，射到波前传感器3的微透镜阵列平面内的波与Z轴的偏移量即与Z轴的夹角能通过该波与X方向夹角和该波与Y方向夹角确定，波前信号φ对应的X轴方向的斜率信号Gx，波前信号φ对应的Y轴方向的斜率信号Gy，波前传感器3的微透镜阵列平面上有256×256个子孔径，通过波前信息的特征建立起能够稀疏表示出波前斜率信号的稀疏字典Ψ，使得波前斜率信号等于极少量的系数乘以稀疏字典Ψ，即Gx＝Ψα,Gy＝Ψβ，α,β表示对应X轴方向和Y轴方向的斜率信号的系数，Gx、Gy、α、β、Ψ都为一个256×256的矩阵，
斜率信号Gx表示为
斜率信号Gy表示为
其中α、β中不为零的数据有3277个，
波前信号φ具有低频多高频少的特点，利用这一特点，通过黄金分割法(见李娟娟等."大气湍流畸变波前斜率的稀疏分解."光电工程2(2018))在频率域对频率分布进行分割，得到符合黄金分割法的频率序列F＝[2.55×10-5,6.08×10-5,...,1.02×103]，将F代入
得到大小为N×N的矩阵，将该矩阵分割为16×16的小块，得到256个块组成稀疏字典Ψ，其中m、n都为大于0小于等于256的正整数，u、v都为F中的某一个频率，j为虚数单位，图3是波前斜率信号的稀疏字典Ψ；
步骤二、变形镜2上有与波前传感器3的子孔径一一对应的变形镜3的子孔径，通过变形镜3的上加载的采样矩阵Φ5对波前信息进行下采样，采样矩阵Φ5为一个12×256的矩阵，矩阵中每个元素的取值为0或者1，其中1代表通，0代表断，波前传感器3得到波前斜率信号的测量信号，X轴方向的测量信号CSx＝ΦGx，Y轴方向的测量信号CSy＝ΦGy，CSx和CSy都为一个12×256矩阵；
步骤三、由公式CSx＝ΦΨα＝Aα和CSy＝ΦΨβ＝Aβ，通过计算机4运用迭代权重最小二乘法(见Chartrand,Rick,and W.Yin."Iteratively reweighted algorithms forcompressive sensing."(2008):3869-3872)，求解得到α、β，再通过Gx＝Ψα,Gy＝Ψβ获得波前信号φ对应的X轴方向的斜率信号Gx和波前信号φ对应的Y轴方向的斜率信号Gy；
步骤四、通过对波前信号φ对应的X轴方向的斜率信号Gx和波前信号φ对应的Y轴方向的斜率信号Gy使用区域法复原出波前信号φ("用于Hartmann-Shack波前探测器的区域法算法研究."强激光与粒子束10.2(1998):0-0)。图4表示原始的波前信号，图5表示复原的波前信号，图6表示二者的误差值。
Claims (1)
1.一种基于压缩感知的波前探测器，其特征在于，包括变形镜、波前传感器和计算机，变形镜对波前信息进行下采样然后通过波前传感器得到压缩采样的波前斜率信息传输到计算机，计算机通过重构算法得到完整的波前信息，计算机通过重构算法得到完整的波前信息的步骤如下：
步骤一、建立直角坐标系，以平行于波前传感器中心轴线的方向为Z轴方向，以波前传感器的微透镜阵列平面内的任意一条横向直线为X轴方向，以波前传感器的微透镜阵列平面内的任意一条纵向直线为Y轴方向，射到波前传感器的微透镜阵列平面内的波与Z轴的偏移量即与Z轴的夹角能通过该波与X方向夹角和该波与Y方向夹角确定，波前信号φ对应的X轴方向的斜率信号Gx，波前信号φ对应的Y轴方向的斜率信号Gy，波前传感器的微透镜阵列平面上有N×N个子孔径，通过波前信息的特征建立起能够稀疏表示出波前斜率信号的稀疏字典Ψ，使得波前斜率信号等于极少量的系数乘以稀疏字典Ψ，即Gx＝Ψα,Gy＝Ψβ，α,β表示对应X轴方向和Y轴方向的斜率信号的系数，Gx、Gy、α、β、Ψ都为一个N×N的矩阵，
斜率信号Gx表示为
斜率信号Gy表示为
其中α、β中不为零的数据有k个，k＜＜N×N，
波前信号φ具有低频多高频少的特点，利用这一特点，通过黄金分割法在频率域对频率分布进行分割，得到符合黄金分割法的频率序列F＝[f1,f2,f3,...,fn]，将F代入
得到大小为N×N的矩阵，将该矩阵分割为的小块，得到N个块组成稀疏字典Ψ，其中m、n都为大于0小于等于N的正整数，u、v都为F中的某一个频率，j为虚数单位；
步骤二、变形镜上有与波前传感器子孔径一一对应的变形镜子孔径，通过变形镜的上加载的采样矩阵Φ对波前信息进行下采样，采样矩阵Φ为一个M×N的矩阵，矩阵中每个元素的取值为0或者1，其中M＜N，其中1代表通，0代表断，波前传感器得到波前斜率信号的测量信号，X轴方向的测量信号CSx＝ΦGx，Y轴方向的测量信号CSy＝ΦGy，CSx和CSy都为一个M×N矩阵；
步骤三、由公式CSx＝ΦΨα＝Aα和CSy＝ΦΨβ＝Aβ，通过迭代权重最小二乘法，求解得到α、β，通过Gx＝Ψα,Gy＝Ψβ获得波前信号φ对应的X轴方向的斜率信号Gx和波前信号φ对应的Y轴方向的斜率信号Gy；
步骤四、通过对波前信号φ对应的X轴方向的斜率信号Gx和波前信号φ对应的Y轴方向的斜率信号Gy使用区域法复原出波前信号φ。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811029190.4A CN109115350B (zh) | 2018-09-04 | 2018-09-04 | 一种基于压缩感知的波前探测系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811029190.4A CN109115350B (zh) | 2018-09-04 | 2018-09-04 | 一种基于压缩感知的波前探测系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109115350A true CN109115350A (zh) | 2019-01-01 |
| CN109115350B CN109115350B (zh) | 2020-04-28 |
Family
ID=64858078
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811029190.4A Expired - Fee Related CN109115350B (zh) | 2018-09-04 | 2018-09-04 | 一种基于压缩感知的波前探测系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109115350B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113300769A (zh) * | 2021-04-01 | 2021-08-24 | 西安理工大学 | 一种基于压缩感知技术的自适应光学系统波前重构的方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104019899A (zh) * | 2014-05-28 | 2014-09-03 | 中国科学院空间科学与应用研究中心 | 一种超灵敏天文望远镜及其天文图像获取方法 |
| CN104915931A (zh) * | 2015-04-16 | 2015-09-16 | 南京信息工程大学 | 基于压缩感知理论的自适应光学波前信号去噪方法 |
| WO2017176901A1 (en) * | 2016-04-07 | 2017-10-12 | Arizona Board Of Regents On Behalf Of The University Of Arizona | Diffraction-based light beam scanner |
-
2018
- 2018-09-04 CN CN201811029190.4A patent/CN109115350B/zh not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104019899A (zh) * | 2014-05-28 | 2014-09-03 | 中国科学院空间科学与应用研究中心 | 一种超灵敏天文望远镜及其天文图像获取方法 |
| CN104915931A (zh) * | 2015-04-16 | 2015-09-16 | 南京信息工程大学 | 基于压缩感知理论的自适应光学波前信号去噪方法 |
| WO2017176901A1 (en) * | 2016-04-07 | 2017-10-12 | Arizona Board Of Regents On Behalf Of The University Of Arizona | Diffraction-based light beam scanner |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113300769A (zh) * | 2021-04-01 | 2021-08-24 | 西安理工大学 | 一种基于压缩感知技术的自适应光学系统波前重构的方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109115350B (zh) | 2020-04-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101893552B (zh) | 基于压缩感知的高光谱成像仪及成像方法 | |
| CN104121990B (zh) | 基于随机光栅的压缩感知宽波段高光谱成像系统 | |
| CN1844873B (zh) | 基于h-s波前传感器的超声速流场探测系统及其探测方法 | |
| CN100588920C (zh) | 变焦距全偏振光谱成像探测系统 | |
| CN103968945B (zh) | 基于二阶压缩感知的超灵敏光谱成像天文望远镜及方法 | |
| CN101598798B (zh) | 一种对高光谱干涉数据进行光谱重建的系统和方法 | |
| CN104019898B (zh) | 一种超灵敏光谱成像天文望远镜及天文光谱成像方法 | |
| CN202485671U (zh) | 基于压缩感知理论的航天相机成像装置 | |
| CN104992424B (zh) | 一种基于离散余弦变换的单像素快速主动成像系统 | |
| CN101546037A (zh) | 基于空间光调制器的无模型波前畸变校正系统 | |
| CN104019899B (zh) | 一种超灵敏天文望远镜及其天文图像获取方法 | |
| Mironov et al. | Statistical characterization of short wind waves from stereo images of the sea surface | |
| CN103234479B (zh) | 基于压缩感知脉冲编码调制的三维信息获取装置和方法 | |
| CN103472456A (zh) | 一种基于稀疏孔径压缩计算关联的主动成像系统及方法 | |
| CN102331336A (zh) | 长焦距大口径透镜的焦距测量方法及装置 | |
| CN103278143A (zh) | 基于机器视觉的零位走动量测试装置及方法 | |
| CN109115350A (zh) | 一种基于压缩感知的波前探测系统 | |
| CN110213467B (zh) | 一种多尺度调制压缩感知成像系统及其成像方法 | |
| CN101285712B (zh) | 基于分立光强测量器件的线性相位反演波前传感器 | |
| CN104391466A (zh) | 一种快速自适应光学控制器的设计方法 | |
| CN106872977A (zh) | 一种基于分段弱正交匹配追踪的层析sar三维成像方法 | |
| CN105973467B (zh) | 基于压缩感知的声光信号处理系统 | |
| CN102768070B (zh) | 一种互补测量的单光子计数成像系统及方法 | |
| CN104569976B (zh) | 基于稀疏测量的综合孔径辐射计遥感成像的方法和系统 | |
| CN106651928A (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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200428 |