CN104570000B - 一种基于混沌压缩编码的光学综合孔径成像系统及方法 - Google Patents
一种基于混沌压缩编码的光学综合孔径成像系统及方法 Download PDFInfo
- Publication number
- CN104570000B CN104570000B CN201510006433.2A CN201510006433A CN104570000B CN 104570000 B CN104570000 B CN 104570000B CN 201510006433 A CN201510006433 A CN 201510006433A CN 104570000 B CN104570000 B CN 104570000B
- Authority
- CN
- China
- Prior art keywords
- chaos
- matrix
- optical
- dimension
- light
- 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.)
- Expired - Fee Related
Links
- 230000003287 optical effect Effects 0.000 title claims abstract description 96
- 230000000739 chaotic effect Effects 0.000 title claims abstract description 52
- 238000003384 imaging method Methods 0.000 title claims abstract description 37
- 238000000034 method Methods 0.000 title claims abstract description 19
- 238000005070 sampling Methods 0.000 claims abstract description 73
- 230000005855 radiation Effects 0.000 claims abstract description 6
- 239000011159 matrix material Substances 0.000 claims description 100
- 230000015572 biosynthetic process Effects 0.000 claims description 53
- 238000003786 synthesis reaction Methods 0.000 claims description 53
- 238000010276 construction Methods 0.000 claims description 21
- 230000009466 transformation Effects 0.000 claims description 18
- 230000005540 biological transmission Effects 0.000 claims description 6
- 238000001514 detection method Methods 0.000 claims description 6
- 238000007639 printing Methods 0.000 claims description 5
- 230000008569 process Effects 0.000 claims description 5
- 230000008859 change Effects 0.000 claims description 4
- 238000009826 distribution Methods 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 239000011800 void material Substances 0.000 claims 2
- 241000208340 Araliaceae Species 0.000 claims 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims 1
- 235000003140 Panax quinquefolius Nutrition 0.000 claims 1
- 235000008434 ginseng Nutrition 0.000 claims 1
- 238000012544 monitoring process Methods 0.000 abstract description 2
- 238000012545 processing Methods 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 6
- 239000004973 liquid crystal related substance Substances 0.000 description 4
- 239000000523 sample Substances 0.000 description 4
- 230000006835 compression Effects 0.000 description 3
- 238000007906 compression Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000013500 data storage Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000000903 blocking effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/88—Lidar systems specially adapted for specific applications
- G01S17/89—Lidar systems specially adapted for specific applications for mapping or imaging
- G01S17/90—Lidar systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Electromagnetism (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Studio Devices (AREA)
Abstract
本发明涉及光学综合孔径成像领域,特别涉及一种基于混沌压缩编码的光学综合孔径成像系统及方法。包括光学综合孔径子系统、混沌压缩编码子系统和探测成像子系统,光学综合孔径子系统负责接收目标辐射光波,混沌压缩编码子系统对目标辐射光波进行混沌相位调制和随机采样,调制采样后的光波最终被探测成像子系统探测并输入到计算机中进行高分辨率图像重构和显示。本发明可用于机载和星载或其他载体平台,实现对地侦察、监视、遥感和监测等各种应用。
Description
技术领域
本发明涉及光学综合孔径成像领域,特别涉及一种基于混沌压缩编码的光学综合孔径成像系统及方法,属于被动辐射探测、遥感和侦察领域。
背景技术
随着航天遥感技术的不断发展,高分辨率光学遥感器已成为各国争先研制的目标,以满足军事侦查详查等需要。而根据角分辨率公式:(为中心波长,为孔径直径),当波长一定时,分辨率越高,要求空间光学系统的口径越大,但在实际应用中口径的增大受到材料、制造工艺、机械结构、有效载荷舱体积以及成本等诸多因素的限制。光学综合孔径成像系统采用易制造的小口径光学元件或光学系统来合成大口径系统,是实现超大口径、 甚高分辨率空间光学遥感器的有效途径和发展方向之一。
压缩感知(Compressed Sensing,CS)理论是一种利用信号稀疏性或可压缩性的全新的信号采样理论。该理论指出,对于稀疏或可压缩信号,可以以低于甚至远低于奈奎斯特采样速率对其采样而不损失信号的细节。从应用角度考虑来看,CS 理论将图像获取端的成本转移到了后端图像重构端,即以软件效率换取硬件成本。将 CS理论引入光学综合孔径遥感成像系统中,能够用比香农/奈奎斯特采样定理少得多的采样数据完全重构出信号,因此,可以减少采样时间,减少信号处理时间,降低计算成本,降低数据存储和传输代价,提高成像速度;同时,降低对后端光电探测器的分辨率要求,从而降低成本。
CS理论的光学实现方法有两种, 一种是顺序处理,如美国 Rice 大学 Baraniuk等人研制出的单像素相机(参见《Single pixel imaging via compressive sampling.IEEE Signal Processing Magazine, 2008, 25(2): 83~91》),该系统采用数字微镜器件DMD对成像目标进行编码,但是需要连续采样投影M次才能获得重构图像所需的全部测量值,成像时间长,不适于遥感应用。另一种是并行处理,如基于编码孔径的压缩成像(参见《Compressive coded aperture superresolution image reconstruction. IEEEInternational Conference on Acoustics, Speech and Signal Processing, 2008,833~836》),基于随机卷积编码的压缩成像(参见《CMOS compressed sensing by randomconvolution. IEEE International Conference on Acoustics, Speech and SignalProcessing, 2009, 19~24》),并行处理的方法采用编码的方式实现图像的压缩采样,只需一次即可得到全部测量值,更适于对动目标进行快速成像。现有专利(一种基于压缩感知的稀疏孔径成像系统及方法,201310329810.7)将压缩感知理论应用于稀疏孔径成像系统中,采用空间光调制器对光场做随机调制后会聚到点探测器,多次映射获得全部测量值,属于顺序处理。目前基于压缩感知并行处理的光学综合孔径成像系统尚未有文献报道。
发明内容
本发明所要解决的技术问题是:如何实现快速、低成本的光学综合孔径成像系统中的压缩成像
本发明所采用的技术方案是:一种基于混沌压缩编码的光学综合孔径成像系统,包括光学综合孔径子系统、混沌压缩编码子系统和探测成像子系统,光学综合孔径子系统负责接收目标辐射光波,混沌压缩编码子系统对目标辐射光波进行混沌相位调制和随机采样,调制采样后的光波最终被探测成像子系统探测并输入到计算机中进行高分辨率图像重构和显示。
多个小口径望远镜组成的望远镜阵列接收目标辐射的光波,光束变换单元调整光路使光波会聚在光束合成透镜上,用于实现光学综合孔径的直接成像,光束合成透镜焦平面上放置空间光调制器,空间调制器根据混沌相位编码调制矩阵对光场做混沌相位调制,调制后的光线经过一个傅里叶变换透镜后,再通过混沌稀疏采样掩膜板进行光线的随机采样,最终在光电探测器上记录出射光场。将光电探测器采集到的光场强度信号与构造混沌相位编码调制矩阵和混沌稀疏采样矩阵的参数送入计算机进行高分辨率图像重构和显示。
本发明提供的基于混沌压缩编码的光学综合孔径成像系统包括光学综合孔径子系统、混沌压缩编码子系统和探测成像子系统;各所述子系统工作如下:光学综合孔径子系统负责接收目标辐射光波,混沌压缩编码子系统对目标辐射光波进行混沌相位调制和随机采样,调制采样后的光波最终被探测成像子系统探测并输入到计算机中进行高分辨率图像重构和显示。
所述光学综合孔径子系统由望远镜阵列、光束变换单元、光束合成透镜组成;望远镜阵列接收携带目标信息的辐射光波,每个子望远镜的主镜和次镜共焦,使出射光波为平行光,各平行光波经过光束变换单元后分别与原光束平行,构成与原光束等比列缩小的平行光束,最后会聚到光束合成透镜上。
所述望远镜阵列包括至少三个望远镜;所述光束变换单元包括至少三个反射镜组,每个反射镜组包含四个全反镜;一子望远镜,一反射镜组构成一条光路,各条光路上入射的光信号分别入射到所述光束合成透镜上。
所述混沌压缩编码子系统由空间光调制器、傅里叶变换透镜和混沌稀疏采样掩膜板组成;空间光调制器根据混沌相位编码调制矩阵对光束合成透镜焦平面的入射光场做混沌相位调制,调制后的光线经过傅里叶变换透镜后,照射到混沌稀疏采样掩膜板,获得经过随机采样的透射光场。
所述探测成像子系统由光电探测器和计算机组成;光电探测器记录混沌压缩编码子系统的出射光场,将光电探测器采集到的光场强度信号与构造混沌相位编码调制矩阵和混沌稀疏采样矩阵的参数送入计算机进行高分辨率图像重构和显示。
所述光学综合孔径子系统中所述望远镜阵列结构包括:Golay结构、三臂结构、环形结构、环面结构;
所述子望远镜结构包括:折射式望远镜、反射式望远镜;
所述空间光调制器采用透射型数字微透镜阵列或透射型液晶光阀中的任意一种;其分辨率由期望的高分辨率图像的分辨率决定,若高分辨率图像的分辨率为n×n,则所述空间光调制器的分辨率为n×n。
所述空间光调制器加载的混沌相位编码调制矩阵按如下方式产生:
1)设所述空间光调制器(4)的分辨率为n×n,构造1×N维(N=n×n)混沌序列,为混沌序列h的元素。混沌序列h可以采用如下任意一种:Logistic混沌序列、Tent混沌序列;
2)利用混沌序列h按如下方式构造n×n维矩阵T:
其中,矩阵中的各元素取自混沌序列h中的元素;
3)构造n×n维混沌相位编码调制矩阵G,该矩阵的元素为,其中j为虚数单位,表示矩阵T第k行第l列的元素,,方括号[ ]表示取整。矩阵G的元素为1或0,分别对应于空间光调制器独立单元的开和闭。
所述混沌稀疏采样掩膜板对应的矩阵如下设计:
由期望的高分辨率图像的分辨率确定掩膜板的大小,若高分辨率图像的分辨率为n×n,则混沌稀疏采样掩膜板对应的矩阵H为n×n维;
构造1×N维(N=n×n)混沌序列,为混沌序列h的元素。混沌序列h可以采用如下任意一种:Logistic混沌序列、Tent混沌序列;
将1×N维混沌序列h变换成n×n维混沌矩阵H0,则按如下方式构造混沌稀疏采样掩膜板对应的n×n维混沌稀疏采样矩阵H,该矩阵的元素为1或0,对应于掩膜板相应位置的透光和不透光。
其中,sgn()表述符号函数,表述混沌矩阵H0第k行第l列的元素。
所述光电探测器采用CCD或CMOS,且分辨率小于期望的高分辨率图像的分辨率。
本发明提供的一种基于混沌压缩编码的光学综合孔径成像方法,包括以下步骤:
预先构造混沌相位编码调制矩阵和混沌稀疏采样掩膜板对应的矩阵
设所述空间光调制器的分辨率为n×n,则混沌相位编码调制矩阵为n×n维;构造1×N维(N=n×n)混沌序列,为混沌序列h的元素。混沌序列h可以采用如下任意一种:Logistic混沌序列、Tent混沌序列;
利用混沌序列h按如下方式构造n×n维矩阵T:
其中,矩阵中的各元素取自混沌序列h中的元素;
构造n×n维混沌相位编码调制矩阵G,该矩阵的元素为,其中j为虚数单位,表示矩阵T第k行第l列的元素,,方括号[ ]表示取整。矩阵G的元素为1或0,分别对应于空间光调制器独立单元的开和闭;
将1×N维混沌序列h变换成n×n维混沌矩阵H0,则按如下方式构造混沌稀疏采样掩膜板对应的n×n维混沌稀疏采样矩阵H,该矩阵的元素为1或0,对应于掩膜板相应位置的透光和不透光。
其中,sgn(·)表示符号函数,表示混沌矩阵H0第k行第l列的元素。
2)混沌相位调制步骤
光学综合孔径入射的成像光信号X传输到空间光调制器上,所述空间光调制器根据加载在其上的混沌相位编码调制矩阵G对入射光波进行混沌相位调制。
3)混沌稀疏采样步骤
调制后的光波经过一个傅里叶变换透镜后,得到光信号z=F-1(GF(X)),其中X为光学综合孔径入射的成像光信号,F(·) 表示傅里叶变换运算,F-1(·)表示逆傅里叶变换运算,表示任意两个相同维数矩阵相应元素的乘积运算,即光信号z为入射的成像光信号X作傅里叶变换后与混沌相位编码调制矩阵G的各元素相乘,然后作逆傅里叶变换后得到的结果。该信号经过混沌稀疏采样掩膜板,随机选择部分光线通过,得到稀疏采样光信号,其中,H表示混沌稀疏采样矩阵。最后,利用光电探测器记录稀疏采样光信号,由于所述光电探测器的分辨率比所述空间光调制器分辨率低,因此光电探测器会对稀疏采样光信号分块求和,实现降采样,从而获得最终的低维测量值Y=D,其中D为降采样因子,用来表示光电探测器对稀疏采样光信号的降采样过程。
成像重构步骤
将光电探测器的输出Y以及构造所述混沌相位编码调制矩阵和混沌稀疏采样矩阵的参数输入计算机,选取合适的稀疏基来稀疏表示光学综合孔径入射的成像光信号X,通过压缩感知重构算法进行重构,最终实现光学综合孔径的压缩成像。
本发明中n、N、i、k、l为自然数。
本发明的有益效果是:本发明提供的一种基于混沌压缩编码的光学综合孔径成像系统和方法,实现了压缩感知理论在光学综合孔径领域的应用,同时由于本发明光信号的采样和压缩以低速率进行,从而降低数据量,降低计算成本,降低数据存储和传输代价,提高成像速度;降低对后端光电探测器的分辨率要求,或者在光电探测器分辨率相同条件下,可以获得比传统光学综合孔径成像系统更高的成像性能。此外,本发明采用混沌序列构造混沌相位编码调制矩阵和混沌稀疏采样掩膜板,相比传统的随机信号构造方法,无需存储整个矩阵,而只需存储混沌序列的参数,从而节省存储空间。本发明可用于机载和星载或其他载体平台,实现对地侦察、监视、遥感和监测等各种应用。
附图说明
图1为本发明基于混沌压缩编码的光学综合孔径成像系统图;
其中:1、望远镜阵列,2、光束变换单元,3、光束合成透镜,4、空间光调制器,5、傅里叶变换透镜;6、混沌稀疏采样掩膜板,7、光电探测器,8、计算机,201-204、全反镜。
具体实施方式
现结合附图对本发明作进一步的描述。
图1为本发明的基于混沌压缩编码的光学综合孔径成像系统在一个实施例中的结构示意图,该系统包括望远镜阵列1、由4个全反镜201-204构成的光束变换单元2、光束合成透镜3、空间光调制器4、傅里叶变换透镜5、混沌稀疏采样掩膜板6、光电探测器7和计算机8。空间光调制器4位于光束合成透镜3的后焦面,傅里叶变换透镜5的前焦面;混沌稀疏采样掩膜板6位于傅里叶变换透镜5的后焦面。
所述望远镜阵列接收携带目标信息的辐射光波,每个子望远镜的主镜和次镜共焦,使出射光波为平行光,各平行光波经过光束变换单元(2)后分别与原光束平行,构成与原光束等比列缩小的平行光束,最后会聚到光束合成透镜(3)上。
所述望远镜阵列(1)包括至少三个子望远镜;每个子望远镜的主镜和次镜共焦,使出射光波为平行光;所述光束变换单元(2)包括至少三个反射镜组,每个反射镜组包含四个全反镜;一子望远镜,一反射镜组构成一条光路,各条光路的光波为平行光波,与子望远镜出射的原光束平行,且构成与原光束等比列缩小的平行光束,最后会聚到光束合成透镜(3)上。
所述望远镜阵列(1)结构包括:Golay结构、三臂结构、环形结构、环面结构;所述子望远镜结构包括:折射式望远镜、反射式望远镜。
本实施例中,望远镜阵列(1)采用三个卡塞格林望远镜构成Golay3结构(望远镜中心构成等边三角形)。
所述空间光调制器(4)采用透射型数字微透镜阵列或透射型液晶光阀中的任意一种;其分辨率由期望的高分辨率图像的分辨率决定,即系统最终的分辨率由空间光调制器(4)决定。
本实施例中,设期望的高分辨率图像的分辨率为n×n,则空间光调制器(4)的分辨率为n×n。空间光调制器(4)采用透射型电寻址液晶光阀,该液晶光阀的每个单元在随时间变化的电驱动信号的控制下独立地开和闭,从而控制光线的透过与遮挡,该液晶光阀根据加载在其上的混沌相位编码调制矩阵来控制光线的透过与遮挡。
所述空间光调制器(4)加载的混沌相位编码调制矩阵利用混沌序列产生,混沌序列可以采用如下任意一种:Logistic混沌序列、Tent混沌序列。
所述混沌稀疏采样掩膜板(6)对应的矩阵利用混沌序列产生,混沌序列可以采用如下任意一种:Logistic混沌序列、Tent混沌序列。
所述空间光调制器(4)加载的混沌相位编码调制矩阵和混沌稀疏采样掩膜板(6)可以采用不同的混沌序列产生。
本实施例中,利用Logistic混沌序列构造混沌相位编码调制矩阵,并加载到空间光调制器(4)上。混沌稀疏采样掩膜板对应的矩阵也由Logistic混沌序列产生。
所述光电探测器(7)采用CCD或CMOS,且分辨率小于期望的高分辨率图像的分辨率。
本实施例中,光电探测器(7)采用电荷耦合器件CCD,且分辨率是空间光调制器4的四分之一,即空间光调制器4的分辨率为n×n,则光电探测器(7)的分辨率为n/2×n/2。
本实施例中,系统按如下方式工作:望远镜阵列(1)接收目标辐射的光波,光束变换单元(2)调整光路使光波会聚在光束合成透镜(3)上,用于实现光学综合孔径的直接成像,空间调制器(4)根据加载在其上的混沌相位编码调制矩阵对光场做混沌相位调制,调制后的光线经过一个傅里叶变换透镜(5)后,再通过混沌稀疏采样掩膜板(6)进行光线的随机采样,最终在光电探测器(7)上记录出射光场。将光电探测器(7)采集到的光场强度信号与构造混沌相位编码调制矩阵和混沌稀疏采样矩阵的参数送入计算机(8),然后根据压缩感知重构算法进行高分辨率图像重构和显示。
本发明提供的基于混沌压缩编码的光学综合孔径成像方法包括以下步骤:
预先构造混沌相位编码调制矩阵和混沌稀疏采样掩膜板
设所述空间光调制器(4)的分辨率为n×n,则混沌相位编码调制矩阵为n×n维;按如下方式构造1×N维(N=n×n)Logistic混沌序列,为混沌序列h的元素:
其中,为实值控制参数。本实施例中,取,初始值。
利用混沌序列h按如下方式构造n×n维矩阵T:
其中,矩阵中的各元素取自混沌序列h中的元素;
构造n×n维混沌相位编码调制矩阵G,该矩阵的元素为,其中j为虚数单位,表示矩阵T第k行第l列的元素,,方括号[ ]表示取整。矩阵G的元素为1或0,分别对应于空间光调制器4独立单元的开和闭;
将1×N维混沌序列h变换成n×n维混沌矩阵H0,则按如下方式构造混沌稀疏采样掩膜板(6)对应的n×n维混沌稀疏采样矩阵H,该矩阵的元素为1或0,对应于掩膜板相应位置的透光和不透光。
其中,sgn(·)表示符号函数,表示混沌矩阵H0第k行第l列的元素。
混沌相位调制步骤
光学综合孔径入射的成像光信号X传输到空间光调制器(4)上,所述空间光调制器(4)根据加载在其上的混沌相位编码调制矩阵G对入射光波进行混沌相位调制。
混沌稀疏采样步骤
调制后的光波经过一个傅里叶变换透镜(5)后,得到光信号z=F-1(GF(X)),其中X为光学综合孔径入射的成像光信号,F(·) 表示傅里叶变换运算,F-1(·)表示逆傅里叶变换运算,表示任意两个相同维数矩阵相应元素的乘积运算,即光信号z为入射的成像光信号X作傅里叶变换后与混沌相位编码调制矩阵G的各元素相乘,然后作逆傅里叶变换后得到的结果。该信号经过混沌稀疏采样掩膜板(6),随机选择部分光线通过,得到稀疏采样光信号,其中,H表示混沌稀疏采样矩阵。最后,利用光电探测器(7)记录稀疏采样光信号,由于所述光电探测器(7)的分辨率比所述空间光调制器4分辨率低,因此光电探测器7会对稀疏采样光信号分块求和,实现降采样,从而获得最终的低维测量值Y=D,其中D为降采样因子,用来表示光电探测器7对稀疏采样光信号的降采样过程。
4)成像重构步骤
将光电探测器(7)的输出Y以及构造所述混沌相位编码调制矩阵和混沌稀疏采样矩阵的参数输入计算机,选取合适的稀疏基来稀疏表示光学综合孔径入射的成像光信号X,通过压缩感知重构算法进行重构,最终实现光学综合孔径的压缩成像。
以上所述的具体描述,对发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (5)
1.一种基于混沌压缩编码的光学综合孔径成像系统,其特征在于,包括光学综合孔径子系统、混沌压缩编码子系统和探测成像子系统,光学综合孔径子系统负责接收目标辐射光波,混沌压缩编码子系统对目标辐射光波进行混沌相位调制和随机采样,调制采样后的光波最终被探测成像子系统探测并输入到计算机中进行高分辨率图像重构和显示,光学综合孔径子系统包括望远镜阵列、光束变换单元、光束合成透镜,望远镜阵列接收携带目标信息的辐射光波,望远镜阵列包括至少三个子望远镜,光束变换单元包括至少三个反射镜组,每个反射镜组包含四个全反镜,望远镜阵列的每个子望远镜的主镜和次镜共焦,使出射光波为平行光,各平行光波经过光束变换单元后分别与原光束平行,构成与原光束等比列缩小的平行光束,入射到光束合成透镜上,混沌压缩编码子系统包括空间光调制器、傅里叶变换透镜和混沌稀疏采样掩膜板,空间光调制器根据混沌相位编码调制矩阵对光束合成透镜焦平面的入射光场做混沌相位调制,调制后的光线经过傅里叶变换透镜后,照射到混沌稀疏采样掩膜板,获得经过随机采样的透射光场。
2.根据权利要求1所述一种基于混沌压缩编码的光学综合孔径成像系统,其特征在于:探测成像子系统包括光电探测器和计算机,光电探测器记录混沌压缩编码子系统的透射光场,将光电探测器采集到的光场强度信号与构造混沌相位编码调制矩阵和混沌稀疏采样矩阵的参数送入计算机进行高分辨率图像重构和显示。
3.根据权利要求2所述一种基于混沌压缩编码的光学综合孔径成像系统,其特征在于:空间光调制器加载的混沌相位编码调制矩阵按如下方式产生:
空间光调制器的分辨率为n×n,构造1×N维(N=n×n)混沌序列,为混沌序列h的元素,混沌序列h可以采用Logistic混沌序列、Tent混沌序列任意一种;
利用混沌序列h按如下方式构造n×n维矩阵T:
其中,矩阵中的各元素取自混沌序列h中的元素;
构造n×n维混沌相位编码调制矩阵G,该矩阵的元素为,其中j为虚数单位,表示矩阵T第k行第l列的元素,,方括号[ ]表示取整,矩阵G的元素为1或0,分别对应于空间光调制器独立单元的开和闭;
其中n、N、i、k、l为自然数。
4.根据权利要求3所述一种基于混沌压缩编码的光学综合孔径成像系统,其特征在于:所述混沌稀疏采样掩膜板对应的矩阵如下设计:
由期望的高分辨率图像的分辨率确定掩膜板的大小,若高分辨率图像的分辨率为n×n,则混沌稀疏采样掩膜板对应的矩阵H为n×n维;
构造1×N维(N=n×n)混沌序列,为混沌序列h的元素,混沌序列h可以采用如下任意一种:Logistic混沌序列、Tent混沌序列;
将1×N 维混沌序列h变换成n×n维混沌矩阵H0,则按如下方式构造混沌稀疏采样掩膜板对应的n×n维混沌稀疏采样矩阵H,该矩阵的元素为1或0,对应于掩膜板相应位置的透光和不透光,
其中,sgn()表述符号函数,表述混沌矩阵H0第k行第l列的元素,n、N、k、l为自然数。
5.一种利用权利要求4成像系统成像方法,其特征在于,该方法包括以下步骤:
步骤一、预先构造混沌相位编码调制矩阵和混沌稀疏采样掩膜板对应的矩阵
设所述空间光调制器的分辨率为n×n,则混沌相位编码调制矩阵为n×n维;构造1×N维(N=n×n)混沌序列,为混沌序列h的元素,混沌序列h可以采用如下任意一种:Logistic混沌序列、Tent混沌序列;
利用混沌序列h按如下方式构造n×n维矩阵T:
其中,矩阵中的各元素取自混沌序列h中的元素;
构造n×n维混沌相位编码调制矩阵G,该矩阵的元素为,其中j为虚数单位,表示矩阵T第k行第l列的元素,,方括号[ ]表示取整;矩阵G的元素为1或0,分别对应于空间光调制器独立单元的开和闭;
将1×N维混沌序列h变换成n×n维混沌矩阵H0,则按如下方式构造混沌稀疏采样掩膜板对应的n×n维混沌稀疏采样矩阵H,该矩阵的元素为1或0,对应于掩膜板相应位置的透光和不透光:
其中,sgn(·)表示符号函数,表示混沌矩阵H0第k行第l列的元素;
步骤二、混沌相位调制步骤
光学综合孔径入射的成像光信号X传输到空间光调制器上,所述空间光调制器根据加载在其上的混沌相位编码调制矩阵G对入射光波进行混沌相位调制;
步骤三、混沌稀疏采样步骤
调制后的光波经过一个傅里叶变换透镜后,得到光信号z=F-1(GF(X)),其中X为光学综合孔径入射的成像光信号,F(·) 表示傅里叶变换运算,F-1(·)表示逆傅里叶变换运算,表示任意两个相同维数矩阵相应元素的乘积运算,即光信号z为入射的成像光信号X作傅里叶变换后与混沌相位编码调制矩阵G的各元素相乘,然后作逆傅里叶变换后得到的结果;
该信号经过混沌稀疏采样掩膜板,随机选择部分光线通过,得到稀疏采样光信号,其中,H表示混沌稀疏采样矩阵;最后,利用光电探测器记录稀疏采样光信号,由于所述光电探测器的分辨率比所述空间光调制器分辨率低,因此光电探测器会对稀疏采样光信号分块求和,实现降采样,从而获得最终的低维测量值Y=D,其中D为降采样因子,用来表示光电探测器对稀疏采样光信号的降采样过程;
步骤四、成像重构步骤
将光电探测器的输出Y以及构造所述混沌相位编码调制矩阵和混沌稀疏采样矩阵的参数输入计算机,选取合适的稀疏基来稀疏表示光学综合孔径入射的成像光信号X,通过压缩感知重构算法进行重构,最终实现光学综合孔径的压缩成像;
其中:n、N、i、k、l为自然数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510006433.2A CN104570000B (zh) | 2015-01-07 | 2015-01-07 | 一种基于混沌压缩编码的光学综合孔径成像系统及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510006433.2A CN104570000B (zh) | 2015-01-07 | 2015-01-07 | 一种基于混沌压缩编码的光学综合孔径成像系统及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104570000A CN104570000A (zh) | 2015-04-29 |
CN104570000B true CN104570000B (zh) | 2017-02-22 |
Family
ID=53086516
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510006433.2A Expired - Fee Related CN104570000B (zh) | 2015-01-07 | 2015-01-07 | 一种基于混沌压缩编码的光学综合孔径成像系统及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104570000B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104898126A (zh) * | 2015-05-19 | 2015-09-09 | 哈尔滨工业大学 | 一种基于微小卫星平台的合成孔径光学成像系统及方法 |
EP3350998A1 (en) * | 2015-09-17 | 2018-07-25 | Thomson Licensing | An apparatus and a method for encoding an image captured by an optical acquisition system |
WO2017062834A1 (en) * | 2015-10-09 | 2017-04-13 | Visby Camera Corporation | Holographic light field imaging device and method of using the same |
CN105741336B (zh) * | 2016-03-07 | 2019-04-30 | 中国科学院上海高等研究院 | 一种基于随机采样的分块压缩感知成像方法及系统 |
CN106301383B (zh) * | 2016-08-25 | 2019-05-21 | 东北大学 | 一种基于压缩感知的信号处理方法 |
CN107580164A (zh) * | 2017-09-11 | 2018-01-12 | 陕西美星恒祺新能源科技有限公司 | 一种单次测量的压缩感知超分辨率成像系统及其成像方法 |
CN108833919B (zh) * | 2018-06-29 | 2020-02-14 | 东北大学 | 基于随机循环矩阵的彩色单像素成像方法及系统 |
CN110213467B (zh) * | 2019-05-27 | 2020-10-23 | 中国科学院国家空间科学中心 | 一种多尺度调制压缩感知成像系统及其成像方法 |
JP7472136B2 (ja) | 2019-07-12 | 2024-04-22 | 株式会社小糸製作所 | イメージング装置およびその照明装置、車両、車両用灯具 |
CN110794576A (zh) * | 2019-11-01 | 2020-02-14 | 中国科学院光电技术研究所 | 一种基于相位调制的光学合成孔径成像望远镜阵列偏心误差探测方法 |
CN111177742B (zh) * | 2019-11-29 | 2023-11-24 | 上海理工大学 | 基于光场原理和混沌系统的多图像加密方法和解密方法 |
CN111999616B (zh) * | 2020-08-31 | 2023-03-10 | 南京工程学院 | 一种基于空间编码掩模板的单像素日盲紫外成像系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2003010917A1 (en) * | 2001-07-24 | 2003-02-06 | Atlinks Usa, Inc. | Narrow band chaotic frequency modulation using skewed tent maps |
CN102087411A (zh) * | 2010-12-02 | 2011-06-08 | 上海电机学院 | 量子成像方法及量子成像系统 |
CN102608617A (zh) * | 2012-03-09 | 2012-07-25 | 太原理工大学 | 一种基于混沌激光的超宽带穿墙雷达探测装置 |
CN104144343A (zh) * | 2014-07-11 | 2014-11-12 | 东北大学 | 一种数字图像压缩加密联合编码方法 |
-
2015
- 2015-01-07 CN CN201510006433.2A patent/CN104570000B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2003010917A1 (en) * | 2001-07-24 | 2003-02-06 | Atlinks Usa, Inc. | Narrow band chaotic frequency modulation using skewed tent maps |
CN102087411A (zh) * | 2010-12-02 | 2011-06-08 | 上海电机学院 | 量子成像方法及量子成像系统 |
CN102608617A (zh) * | 2012-03-09 | 2012-07-25 | 太原理工大学 | 一种基于混沌激光的超宽带穿墙雷达探测装置 |
CN104144343A (zh) * | 2014-07-11 | 2014-11-12 | 东北大学 | 一种数字图像压缩加密联合编码方法 |
Non-Patent Citations (3)
Title |
---|
一种基于编码孔径的压缩传感图像重构算法的改进;张智诠,丁晟;《光学技术》;20140331;第40卷(第2期);第134页 * |
光学综合孔径成像技术实验研究;王长伟 等;《应用光学》;20100131;第31卷(第1期);第9-10页 * |
基于混沌信号的SAR抗距离向相干干扰方法;谢春健 等;《计算机仿真》;20101130;第27卷(第11期);第133页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104570000A (zh) | 2015-04-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104570000B (zh) | 一种基于混沌压缩编码的光学综合孔径成像系统及方法 | |
Mait et al. | Computational imaging | |
CN103968945B (zh) | 基于二阶压缩感知的超灵敏光谱成像天文望远镜及方法 | |
Bian et al. | Efficient single pixel imaging in Fourier space | |
CN101893552B (zh) | 基于压缩感知的高光谱成像仪及成像方法 | |
Ragazzoni et al. | Adaptive optics for 100-m-class telescopes: new challenges require new solutions | |
CN103471715B (zh) | 一种共光路组合式光场光谱成像方法及装置 | |
CN104019899B (zh) | 一种超灵敏天文望远镜及其天文图像获取方法 | |
CN105589210B (zh) | 一种基于光瞳调制的数字化合成孔径成像方法 | |
CN104019898A (zh) | 一种超灵敏光谱成像天文望远镜及天文光谱成像方法 | |
CN103398729B (zh) | 一种基于压缩感知的稀疏孔径成像系统及方法 | |
CN103913228B (zh) | 编码模板多目标超分辨率时间飞行成像系统及方法 | |
CN101228460A (zh) | 编码孔径成像系统 | |
Ferreira et al. | Real-time end-to-end AO simulations at ELT scale on multiple GPUs with the COMPASS platform | |
CN109828285B (zh) | 一种双波段时域压缩感知高速成像方法及装置 | |
WO2020124992A1 (zh) | 一种基于透射式双缝孔径编码成像系统及其超分辨方法 | |
CN105116542A (zh) | 一种双视场计算关联成像系统及方法 | |
CN116245726A (zh) | 基于深度学习框架的压缩感知偏振超分辨成像方法 | |
Zhang et al. | MEMS-based super-resolution remote sensing system using compressive sensing | |
CN116609942B (zh) | 一种分孔径压缩感知偏振超分辨率成像方法 | |
CN105120141B (zh) | 一种压缩感知光电成像方法及装置 | |
CN102768070B (zh) | 一种互补测量的单光子计数成像系统及方法 | |
CN109781260B (zh) | 超紧凑型快照式偏振光谱成像探测装置及探测方法 | |
CN106991637A (zh) | 一种利用gpu并行计算实现多分辨率光场分解的方法 | |
CN107580164A (zh) | 一种单次测量的压缩感知超分辨率成像系统及其成像方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170222 |
|
CF01 | Termination of patent right due to non-payment of annual fee |