CN106405546B - 基于压缩散斑的快速关联成像***及方法 - Google Patents
基于压缩散斑的快速关联成像***及方法 Download PDFInfo
- Publication number
- CN106405546B CN106405546B CN201610979352.5A CN201610979352A CN106405546B CN 106405546 B CN106405546 B CN 106405546B CN 201610979352 A CN201610979352 A CN 201610979352A CN 106405546 B CN106405546 B CN 106405546B
- Authority
- CN
- China
- Prior art keywords
- compressed
- scene
- information
- matrix
- speckles
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 43
- 238000000034 method Methods 0.000 title claims abstract description 25
- 238000007906 compression Methods 0.000 claims abstract description 45
- 230000006835 compression Effects 0.000 claims abstract description 45
- 238000005070 sampling Methods 0.000 claims abstract description 34
- 238000005516 engineering process Methods 0.000 claims abstract description 10
- 239000011159 matrix material Substances 0.000 claims description 52
- 230000002596 correlated effect Effects 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 6
- 230000008707 rearrangement Effects 0.000 claims description 6
- 238000005286 illumination Methods 0.000 claims description 4
- 238000013144 data compression Methods 0.000 claims description 3
- 238000009499 grossing Methods 0.000 claims description 3
- 230000001678 irradiating effect Effects 0.000 claims description 3
- 238000002310 reflectometry Methods 0.000 claims description 3
- 230000009471 action Effects 0.000 claims description 2
- 230000003287 optical effect Effects 0.000 claims 3
- 238000011084 recovery Methods 0.000 claims 1
- 238000001514 detection method Methods 0.000 abstract description 3
- 238000002474 experimental method Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 3
- 101100277917 Caenorhabditis elegans dmd-3 gene Proteins 0.000 description 2
- 230000003760 hair shine Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Image Processing (AREA)
- Studio Devices (AREA)
Abstract
本发明公开了一种基于压缩散斑的快速关联成像***及方法,结合压缩采样技术、关联成像技术和稀疏信息复原技术,利用压缩散斑对场景信息进行压缩采样,利用压缩散斑和探测强度复原场景压缩信息。通过压缩散斑降低了数据量,提高了关联成像***的成像速度。最后利用场景压缩信息通过压缩感知实现了对完整准确场景信息的复原,该专利在关联成像领域具有广泛的应用价值。
Description
技术领域
本发明涉及关联成像技术领域,具体是一种基于压缩散斑的快速关联成像***及方法。
背景技术
近年来关联成像技术一直被人们研究,传统关联成像***利用散斑光对场景进行照射,采用单像素探测采集***获取场景回波信号,利用散斑光分布信息和探测强度信息进行相关运算获取场景图像信息。该类***正朝着工程实际应用发展,但***需要不同散斑光的很多次照射采集才能有效获取场景图像信息,因此***实时成像能力受到了很大约束。目前,提高关联成像***的成像速度和质量是该领域研究方向之一。利用绝大部分自然场景图像信息具有稀疏性的特点,本项发明利用压缩散斑对场景信息进行压缩采样,降低了关联成像***数据量,可以提高关联成像***运算速度,在关联成像领域具有广阔的应用前景。
发明内容
本发明的目的是提供一种基于压缩散斑的快速关联成像***及方法,以提高关联成像***的成像速度。
为了达到上述目的,本发明所采用的技术方案为:
基于压缩散斑的快速关联成像***及方法,其特征在于:包括光源、扩束器、光调制器DMD、投影透镜、汇聚透镜、单像素探测器和数据采集***,所述光源发出的光经扩束器后照射到光调制器DMD上,光调制器DMD对光进行调制,使得调制光场产生的散斑具有压缩稀疏形式,能够实现对场景的压缩采样,调制光场产生的压缩散斑经投影透镜照射场景,场景反射光经汇聚透镜汇聚到单像素探测器上,单像素探测器与数据采集***连接,由数据采集***通过单像素探测器获取场景总的反射光强信息。
所述的基于压缩散斑的快速关联成像***及方法,其特征在于:光调制器DMD可以实现对光场的调制,能够产生压缩散斑,且该散斑分布形式能够通过调制信息准确获知,单像素探测器可以对总反射光强信息进行探测。
所述的基于压缩散斑的快速关联成像***及方法,其特征在于:采用压缩散斑对物体信息进行压缩采样有效地降低了关联成像中的数据量,实现快速关联成像。
所述***的快速关联成像方法,其特征在于:包括以下步骤:
(1)、首先固定压缩矩阵,然后依次在压缩矩阵的基础上进行不同的稀疏采样并加载到光调制器DMD上来实现压缩散斑,利用压缩散斑对场景进行照射,单像素探测器对场景反射信号进行探测;采集***对很多个不同压缩散斑照射场景的反射信号进行采集并保存;设置压缩矩阵B,其为M×M大小的含0,1的矩阵,其中有效采样1的点数为N×N,满足条件M>N,数据压缩率γ=N2/M2;假设第j次利用稀疏采样矩阵Sj与压缩矩阵B进行作用后的压缩散斑对场景进行照射,采集***获取的光强信息为i,则该过程可以表示为如公式(1)所示:
i=∑M×MBSjR, (1),
这里,R表示场景反射率信息,∑M×M表示对M×M矩阵的元素进行求和;
(2)、利用压缩矩阵中的有效采样位置索引对照射场景的散斑进行矩阵调整,剔除数据中的未采样点即压缩矩阵中0值点,数据量大大降低,能够有效提高关联成像的运算速度,进行算法运算后可以获得压缩场景信息;采用压缩矩阵B对散斑进行了压缩,只选出采样点为1位置,按照矩阵N×N对上式(1)进行改写为如公式(2)所示:
i=∑N×NSj′R′, (2),
其中,Sj′为BSj经过重新排列形式,R′为R的采样点重新排列形式,大小为N×N,即为原信息的压缩形式;按照公式(2)形式对k1次散斑照明采集的过程用矩阵表示为如公式(3)所示:
其中,为k1次采集的光强信息形式,大小为k1×1,为R′的列向量表示形式,为k1次不同压缩散斑的矩阵表示形式,大小为k1×N2;通过对公式的分析可以看出,需要处理的数据矩阵大小为N2×k1,数据的压缩比为γ;与此同时由于需要求解R′未知数的个数相对于R的要少很多,因此可以减小散斑照射次数,实际***中数据量的压缩程度比γ还要小;
(3)、利用上述步骤获取的压缩场景信息,按照有效采样位置索引进行排列获取与原场景大小一致稀疏信息,采用稀疏信息复原技术可以获得完整准确的场景信息;根据压缩矩阵的有效索引信息把调整为大小为M×M的稀疏场景信息RB;利用压缩感知稀疏信息复原技术,由RB还原出完整准确的场景信息R,求解最优公式为:
argmin[||RB-BR||2+λΦTV(R)] (4),
公式(4)中,ΦTV代表总变分模型函数,λ表示平滑系数。
本发明中光调制器DMD按照压缩散斑形式对光进行调制,产生压缩散斑,能够实现对场景的压缩采样。与传统关联成像利用散斑实现场景信息的全采样不同,本发明实现了场景信息的压缩采样,数据量大幅减小,可以提高关联成像的成像速度。
本发明的优点是:本发明结合压缩采样技术、关联成像技术和稀疏信息复原技术,利用压缩散斑降低数据量,利用信息相关性实现场景压缩信息的快速获取,提高关联成像***的成像速度,最后通过压缩感知实现完整准确场景信息的复原,该方法在关联成像领域具有广泛的应用价值。
附图说明
图1为本发明的结构示意图。
图2***实验结果图,其中:
图2a为具体实施实验中要进行成像的物体图像,图2b为利用压缩矩阵信息获取的物体稀疏信息图,图2c为采用压缩感知对稀疏信息进行复原得到的信息图。
具体实施方式
如图1所示,基于压缩散斑的快速关联成像***,其装置包括光源1、扩束器2、光调制器DMD3、投影透镜4、汇聚透镜5、单像素探测器6和数据采集处理***7;
光源1发光,经扩束器2后具有较大的截面分布,然后光束照射到光调制器DMD3上产生压缩散斑经投影透镜4对场景进行照射,场景反射光经汇聚透镜5到达单像素探测器6上,然后经数据采集处理***7进行数据采集、保存和处理。
首先固定压缩矩阵,然后依次在压缩矩阵的基础上进行不同的稀疏采样来实现压缩稀疏散斑,利用压缩稀疏散斑对场景进行照射,单像素探测器对场景反射信号进行探测。采集***对很多次不同压缩稀疏散斑照射场景的反射信号进行采集并保存。设置压缩矩阵B,其为M×M大小的含0,1的矩阵,其中有效采样1的点数为N×N,满足条件(M>N),数据压缩率γ=N2/M2。假设第j次利用稀疏采样矩阵Sj与压缩矩阵B进行作用后的压缩散斑对场景进行照射,采集***获取的光强信息为i,则该过程可以表示为:
i=∑M×MBSjR, (1),
这里,R表示场景反射率信息,∑M×M表示对M×M矩阵的元素进行求和。
利用压缩矩阵中的有效采样位置索引对照射场景的相关压缩散斑进行矩阵调整,剔除数据中的未采样点即压缩矩阵中0值点,数据量大大降低,能够有效提高关联成像的运算速度,进行算法运算后可以获取压缩场景信息。采用压缩矩阵B对散斑进行了压缩,只选出采样点为1位置,按照矩阵N×N对上式(1)进行改写为:
i=∑N×NSj′R′, (2),
其中,Sj′为BSj经过重新排列形式,R′为R的采样点重新排列形式,大小为N×N,即为原信息的压缩形式。按照公式(2)形式对k1次散斑照明采集的过程用矩阵表示为:
其中,为k1次采集的光强信息形式,大小为k1×1,为R′的列向量表示形式,为k1次不同压缩稀疏散斑的矩阵表示形式,大小为k1×N2。通过对公式的分析可以看出,采用传统的方法进行算法求解时需要处理的数据矩阵大小为M2×k1,而利用本专利需要处理的数据矩阵大小为N2×k1,数据的压缩比为γ。与此同时由于需要求解R′未知数的个数相对于R的要少很多,因此可以减小散斑照射次数,实际***中数据量的压缩程度比γ还要小。
利用上述步骤获取的压缩场景信息,按照有效采样位置索引进行排列获取与原场景大小一致稀疏信息,采用稀疏信息复原技术可以获得完整准确的场景信息。根据压缩矩阵的有效索引信息把调整为大小为M×M的稀疏场景信息RB。利用压缩感知稀疏信息复原技术,由RB还原出完整准确的场景信息R,求解最优公式为:
argmin[||RB-BR||2+λΦTV(R)] (4),
这里,ΦTV代表总变分(TV)模型函数,λ表示平滑系数。
显然,本领域的技术人员可以通过对本发明所涉及的压缩散斑关联成像***进行改动和变型而不脱离本发明的精神和范围。这样,如果这些修改和变动属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些修改和变型在内。
实验结果
为了验证本发明的可行性,设计***并进行了一次实验,实验结果如图2所示。实验中要进行成像的物体如图2a所示,成像分辨率大小为80×80。使用压缩矩阵有效采样数为4096,即压缩比为0.64。利用压缩稀疏散斑对物体进行3400次照射,利用压缩矩阵信息获取的物体稀疏信息如图2b所示。最后采用压缩感知对稀疏信息进行复原得到的信息如图2c所示,通过结果可以看出该方法可以有效获取物体信息。
Claims (1)
1.基于压缩散斑的快速关联成像方法,其特征在于:该方法由基于压缩散斑的快速关联成像***实现,所述的***包括光源、扩束器、光调制器DMD、投影透镜、汇聚透镜、单像素探测器和数据采集***,所述光源发出的光经扩束器后照射到光调制器DMD上,光调制器DMD对光进行调制,使得调制光场产生的散斑具有压缩稀疏形式,能够实现对场景的压缩采样,调制光场产生的压缩散斑经投影透镜照射场景,场景反射光经汇聚透镜汇聚到单像素探测器上,单像素探测器与数据采集***连接,由数据采集***通过单像素探测器获取场景总的反射光强信息;
光调制器DMD可以实现对光场的调制,能够产生压缩散斑,且该散斑分布形式能够通过调制信息准确获知,单像素探测器可以对总反射光强信息进行探测;
采用压缩散斑对物体信息进行压缩采样有效地降低了关联成像中的数据量,实现快速关联成像;
所述基于压缩散斑的快速关联成像方法,包括以下步骤:
(1)、首先固定压缩矩阵,然后依次在压缩矩阵的基础上进行不同的稀疏采样并加载到光调制器DMD上来实现压缩散斑,利用压缩散斑对场景进行照射,单像素探测器对场景反射信号进行探测;采集***对很多个不同压缩散斑照射场景的反射信号进行采集并保存;设置压缩矩阵B,其为M×M大小的含0,1的矩阵,其中有效采样1的点数为N×N,满足条件M>N,数据压缩率γ=N2/M2;假设第j次利用稀疏采样矩阵Sj与压缩矩阵B进行作用后的压缩散斑对场景进行照射,采集***获取的光强信息为i,则该过程可以表示为如公式(1)所示:
i=∑M×MBSjR, (1),
这里,R表示场景反射率信息,∑M×M表示对M×M矩阵的元素进行求和;
(2)、利用压缩矩阵中的有效采样位置索引对照射场景的压缩散斑进行矩阵调整,剔除数据中的未采样点即压缩矩阵中0值点,数据量大大降低,能够有效提高关联成像的运算速度,进行算法运算后可以获得压缩场景信息;采用压缩矩阵B对散斑进行了压缩,只选出采样点为1位置,按照矩阵N×N对上式(1)进行改写为如公式(2)所示:
i=∑N×NS′jR′, (2),
其中,S′j为BSj经过重新排列形式,R′为R的采样点重新排列形式,大小为N×N,即为原信息的压缩形式;按照公式(2)形式对k1次散斑照明采集的过程用矩阵表示为如公式(3)所示:
其中,为k1次采集的光强按行排列的形式,大小为k1×1,为R′的列向量表示形式,为k1次不同压缩散斑的矩阵表示形式,大小为k1×N2;通过对公式的分析可以看出,需要处理的数据矩阵大小为N2×k1,数据的压缩比为γ;与此同时由于需要求解R′未知数的个数相对于R的要少很多,因此可以减小散斑照射次数,实际***中数据量的压缩程度比γ还要小;
(3)、利用上述步骤获取的压缩场景信息,按照有效采样位置索引进行排列获取与原场景大小一致稀疏信息,采用稀疏信息复原技术可以获得完整准确的场景信息;根据压缩矩阵的有效索引信息把调整为大小为M×M的稀疏场景信息RB;利用稀疏信息复原技术,由RB还原出完整准确的场景信息R,求解最优公式为:
argmin[||RB-BR||2+λΦTV(R)] (4),
公式(4)中,ΦTV代表总变分模型函数,λ表示平滑系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610979352.5A CN106405546B (zh) | 2016-11-08 | 2016-11-08 | 基于压缩散斑的快速关联成像***及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610979352.5A CN106405546B (zh) | 2016-11-08 | 2016-11-08 | 基于压缩散斑的快速关联成像***及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106405546A CN106405546A (zh) | 2017-02-15 |
CN106405546B true CN106405546B (zh) | 2021-01-15 |
Family
ID=58015383
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610979352.5A Active CN106405546B (zh) | 2016-11-08 | 2016-11-08 | 基于压缩散斑的快速关联成像***及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106405546B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108957448B (zh) * | 2018-06-06 | 2022-10-28 | 西安电子科技大学 | 一种基于广义全变差正则化的雷达关联成像方法 |
CN108895985B (zh) * | 2018-06-19 | 2020-05-15 | 中国科学院合肥物质科学研究院 | 一种基于单像素探测器的物体定位方法 |
CN109814128B (zh) * | 2019-01-23 | 2020-08-11 | 北京理工大学 | 时间飞行与关联成像相结合的高分辨快速成像***及方法 |
CN110111271B (zh) * | 2019-04-24 | 2021-06-22 | 北京理工大学 | 一种基于侧抑制网络的单像素成像方法 |
CN110389440A (zh) * | 2019-07-22 | 2019-10-29 | 上海理工大学 | 基于关联成像和改进光纤束的内窥镜成像***及方法 |
CN110646810B (zh) * | 2019-09-27 | 2021-08-06 | 北京理工大学 | 一种散斑优化压缩感知鬼成像方法及*** |
CN111707413B (zh) * | 2020-05-15 | 2021-11-26 | 中国科学院合肥物质科学研究院 | 一种基于单像素探测器的质心探测方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103363924A (zh) * | 2013-07-15 | 2013-10-23 | 中国科学院空间科学与应用研究中心 | 一种压缩的三维计算鬼成像***及方法 |
CN105892061A (zh) * | 2016-06-24 | 2016-08-24 | 北京国承万通信息科技有限公司 | 显示装置及方法 |
-
2016
- 2016-11-08 CN CN201610979352.5A patent/CN106405546B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103363924A (zh) * | 2013-07-15 | 2013-10-23 | 中国科学院空间科学与应用研究中心 | 一种压缩的三维计算鬼成像***及方法 |
CN105892061A (zh) * | 2016-06-24 | 2016-08-24 | 北京国承万通信息科技有限公司 | 显示装置及方法 |
Non-Patent Citations (2)
Title |
---|
《L1范数的总变分正则化超分辨率图像重建》;刘志文 等;《微处理机》;20120630(第3期);第37-39页 * |
《基于稀疏测量的超分辨率鬼成像研究》;李龙珍;《中国优秀硕士学位论文全文数据库 信息科技辑》;20160215;第I138-1336页 * |
Also Published As
Publication number | Publication date |
---|---|
CN106405546A (zh) | 2017-02-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106405546B (zh) | 基于压缩散斑的快速关联成像***及方法 | |
CN103323396B (zh) | 一种基于符合测量的二维压缩鬼成像***及方法 | |
Jiang et al. | Advancing opportunistic sensing in hydrology: A novel approach to measuring rainfall with ordinary surveillance cameras | |
US20140327944A1 (en) | Hologram processing method and system | |
CN108895985B (zh) | 一种基于单像素探测器的物体定位方法 | |
CN105931196A (zh) | 基于傅里叶光学建模的编码光圈相机图像恢复方法 | |
JP2009536748A (ja) | 位相検索および位相ホログラムの合成 | |
US20120330162A1 (en) | Modulated aperture imaging for automatic moving target detection | |
CN114721006B (zh) | 高精度单像素成像方法及*** | |
JP2023522442A (ja) | モデル生成方法、画像透視図の決定方法、装置、ディバイス及び媒体 | |
CN116893429A (zh) | 基于圆谐傅里叶光场的单像素成像方法和目标识别方法 | |
CN116157867A (zh) | Lfa测试条的神经网络分析 | |
CN112634380A (zh) | 一种单帧超光学记忆效应的多目标彩色散射成像方法 | |
CN111123284B (zh) | 探测方法及探测设备 | |
CN205038150U (zh) | 一种太赫兹成像装置 | |
CN113163082B (zh) | 基于散射介质和空间光调制器共同调制的单像素成像*** | |
CN110398477B (zh) | 材质分类及识别方法 | |
CN105241816B (zh) | 一种太赫兹成像装置和方法 | |
López-Cepero | AstroDART: Astronomical data analysis and recovery from tracklets | |
Hu et al. | Fourier ptychographic layer-based imaging of hazy environments | |
CN116148197B (zh) | 一种基于时空调制的非重复光谱高速测量***和方法 | |
CN117538900A (zh) | 一种林下遮蔽目标试验模拟成像装置及方法 | |
CN105654441B (zh) | 基于空间二阶相关的抑制压缩成像中动态噪声的方法 | |
CN113971642A (zh) | 微光环境下单像素成像***一维信号的降噪方法 | |
Song et al. | The application of compressed sensing algorithm based on total variation method into ghost image reconstruction |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |