CN113890624A - 一种频域鬼成像光谱检测方法及装置 - Google Patents
一种频域鬼成像光谱检测方法及装置 Download PDFInfo
- Publication number
- CN113890624A CN113890624A CN202111339117.9A CN202111339117A CN113890624A CN 113890624 A CN113890624 A CN 113890624A CN 202111339117 A CN202111339117 A CN 202111339117A CN 113890624 A CN113890624 A CN 113890624A
- Authority
- CN
- China
- Prior art keywords
- frequency domain
- coding
- frequency
- spectrum
- ghost imaging
- 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
- 238000001228 spectrum Methods 0.000 title claims abstract description 58
- 238000001514 detection method Methods 0.000 title claims abstract description 39
- 238000003384 imaging method Methods 0.000 title claims abstract description 38
- 230000003287 optical effect Effects 0.000 claims abstract description 50
- 239000006185 dispersion Substances 0.000 claims abstract description 31
- 238000000034 method Methods 0.000 claims abstract description 10
- 230000003595 spectral effect Effects 0.000 claims description 24
- 239000000835 fiber Substances 0.000 claims description 7
- 239000011159 matrix material Substances 0.000 description 13
- 230000000875 corresponding effect Effects 0.000 description 10
- 238000013507 mapping Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 238000011084 recovery Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B10/00—Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
- H04B10/50—Transmitters
- H04B10/572—Wavelength control
Landscapes
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
- Optical Modulation, Optical Deflection, Nonlinear Optics, Optical Demodulation, Optical Logic Elements (AREA)
Abstract
本发明公开了一种频域鬼成像光谱检测方法，采集待测目标在频域预编码结构光场中的光强信息，然后通过将采集到的光强信息与所述频域预编码结构光场作互相关运算以重建出待测目标频谱；所述频域预编码结构光场的生成方法如下：利用色散傅里叶变换对光脉冲进行时域拉伸，然后用预设的伪随机编码对时域拉伸后的光脉冲进行时域上的强度调制。本发明还公开了一种频域鬼成像光谱检测装置。相比现有技术，本发明通过色散拉伸结合强度调制的方式，以低成本手段实现对频域光场的伪随机编码，无需使用可编程光滤波器，从而在频域鬼成像无需高成本探测器的基础上，进一步降低系统的调制成本，实现对频谱目标的高分辨低成本探测。
Description
技术领域
本发明涉及一种光谱检测方法，尤其涉及一种频域鬼成像光谱检测方法。
背景技术
鬼成像技术是一种通过光场强度的关联测量非局域地获得物体目标的新型成像手段。通过对区别于有无经过目标的两束光路得到的光场分布和积分结果进行互相关计算，可以以超出系统带宽限制的频率传输信号，具有可成像波段范围广，能抵抗恶劣环境干扰，实现超高分辨成像等优点。由于时空频之间的对偶关系，近年来鬼成像逐渐扩展到时频域，能够在单臂架构中，通过预设对应域的结构光场编码，经过探测目标之后，利用单个慢速探测器的积分测量结果与编码光场之间的互相关计算实现对应域目标信号恢复，为实现时频域目标的低成本高分辨探测成像提供了全新、可行的解决方案。
频域鬼成像通过对光场频谱强度进行关联测量实现目标频谱信息的获取。其中，对光场的频谱强度分布进行高分辨的编码是实现频域鬼成像光谱检测的关键。目前已提出的频域鬼成像编码方案包括：(1)利用超连续谱光源自身的光谱随机强度波动作为编码信息；(2)通过可编程光滤波器对宽谱光源进行频谱强度的伪随机调制编码。然而，前者的编码信息未知，需要设立单独的参考臂，通过高速探测同步获取对应频谱的强度波动，再与积分得到的强度作互相关运算实现目标频谱的检测；而后者受限于可编程光滤波器的分辨率和调制速度，不仅难以实现快速、高精度的频谱强度编码，也难以降低系统检测成本，限制了频域鬼成像光谱检测方案的发展与应用。
发明内容
本发明所要解决的技术问题在于克服现有技术不足，提供一种频域鬼成像光谱检测方法，采用色散拉伸傅里叶变换技术以低成本的方式实现频域鬼成像中光场频谱强度的伪随机编码，从而实现对频域目标的低成本高分辨率探测。
本发明具体采用以下技术方案解决上述技术问题：
一种频域鬼成像光谱检测方法，采集待测目标在频域预编码结构光场中的光强信息，然后通过将采集到的光强信息与所述频域预编码结构光场作互相关运算以重建出待测目标频谱；所述频域预编码结构光场的生成方法如下：利用色散傅里叶变换对光脉冲进行时域拉伸，然后用预设的伪随机编码对时域拉伸后的光脉冲进行强度调制。
优选地，通过大色散啁啾光纤光栅实现所述时域拉伸。
优选地，使用马赫曾德尔调制器实现所述强度调制。
优选地，使用锁模激光器生成所述光脉冲。
优选地，时域拉伸后的相邻光脉冲之间互不重叠。
基于同一发明构思还可以得到以下技术方案：
一种频域鬼成像光谱检测装置，包括：用于生成频域预编码结构光场的结构光生成模块，用于采集待测目标在所述频域预编码结构光场中的光强信息的光强采集模块，用于通过将采集到的光强信息与所述频域预编码结构光场作互相关运算以重建出待测目标频谱的互相关运算模块；所述结构光生成模块包括：
色散拉伸模块，用于利用色散傅里叶变换对光脉冲进行时域拉伸；
强度调制器，用于用预设的伪随机编码对时域拉伸后的光脉冲进行强度调制。
优选地，所述色散拉伸模块为大色散啁啾光纤光栅。
优选地，所述强度调制器为马赫曾德尔调制器。
优选地，使用锁模激光器生成所述光脉冲。
优选地，时域拉伸后的相邻光脉冲之间互不重叠。
相比现有技术，本发明技术方案具有以下有益效果：
本发明通过色散拉伸结合强度调制的方式，以低成本手段实现对光场频谱强度的伪随机编码，无需使用可编程光滤波器，从而在频域鬼成像无需高成本探测器的基础上，进一步降低系统的调制成本，实现对频谱目标的高分辨低成本探测。
附图说明
图1是本发明频域鬼成像光谱检测装置的一种具体结构示意图；
图2为光脉冲经过色散拉伸模块前后的时间波形以及脉冲频谱波形对比；
图3为强度调制后得到的对应编码变化的频谱形状示意图；
图4为基于Hadamard矩阵的伪随机编码波形，其中的(a)表示64*64的Hadamard矩阵形式，(b)为Hadamard矩阵做自相关的结果，表示为HTH，其中上标T表示转置。
具体实施方式
现有的频域鬼成像光谱检测架构通过可编程光滤波器实现对宽谱光源的强度调制，由于调制信息可控，作用到频谱目标上的光场分布I(λ)已知，因此通过积分探测器采集总光强的响应为：
Si＝∫Ii(λ)T(λ)dλ (1)
式中Ii(λ)为第i次采样的光场分布，T(λ)为目标对光场的调制系数。通过编码的不同模式所对应的N次采样后，重建的频谱目标可由二阶关联算法给出：
将该方程表示为矩阵的形式，H代表编码矩阵，T代表目标的矢量形式，G代表重建目标。则可以得到二阶关联算法的矩阵表示形式
然而现有的这种频域鬼成像光谱检测架构需要使用可编程光滤波器对光场频谱强度进行伪随机滤波调制，一方面导致系统成本过高，另一方面系统性能会受到可编程光滤波器的分辨率和模式切换速度限制。
为解决这一问题，本发明的思路是通过色散拉伸的方式利用色散傅里叶技术实现脉冲信号频谱到时域的映射，通过对时域波形强度的伪随机调制，实现对应频谱强度分布的伪随机编码，从而解决频域鬼成像中高精度频谱编码分辨率受限、调制速度受限的难题，以低成本方式实现对频谱目标的高精度检测。
本发明的解决方案具体如下：
一种频域鬼成像光谱检测方法，采集待测目标在频域预编码结构光场中的光强信息，然后通过将采集到的光强信息与所述频域预编码结构光场作互相关运算以重建出待测目标频谱；所述频域预编码结构光场的生成方法如下：利用色散傅里叶变换对光脉冲进行时域拉伸，然后用预设的伪随机编码对时域拉伸后的光脉冲进行强度调制。
一种频域鬼成像光谱检测装置，包括：用于生成频域预编码结构光场的结构光生成模块，用于采集待测目标在所述频域预编码结构光场中的光强信息的光强采集模块，用于通过将采集到的光强信息与所述频域预编码结构光场作互相关运算以重建出待测目标频谱的互相关运算模块；所述结构光生成模块包括：
色散拉伸模块，用于利用色散傅里叶变换对光脉冲进行时域拉伸；
强度调制器，用于用预设的伪随机编码对时域拉伸后的光脉冲进行强度调制。
所述色散拉伸模块优选为色散光纤或大色散啁啾光纤光栅；此外，为实现脉冲光谱与时域映射波形一致，最好保证时域拉伸后的相邻光脉冲之间互不重叠。
优选地，所述强度调制器为马赫曾德尔调制器。
优选地，使用锁模激光器生成所述光脉冲。
为了便于公众理解，下面通过一个优选实施例，并结合附图来对本发明的技术方案进行详细说明：
如图1所示，本实施例采用锁模激光器作为激光源产生重频稳定的光脉冲，然后通过大色散啁啾光纤光栅实现光脉冲在时域上的色散拉伸，同时将光脉冲频谱形状映射到时域上；随后利用马赫曾德尔调制器将微波源产生的一系列伪随机编码波形强度调制到拉伸后的时域光脉冲上，得到对应强度编码的脉冲频谱；将一系列强度编码的脉冲频谱通过频谱目标，通过光功率计实现对应光信号强度的探测；最终通过对探测得到的光强度与编码光谱之间做互相关运算，实现对待测目标频谱的高分辨还原。
图2为光脉冲经过色散拉伸模块前后时间波形以及脉冲频谱波形的对比。如图2所示：当光脉冲经过色散拉伸模块时，大色散啁啾光纤光栅为光脉冲的频谱引入二阶相位调制，输出光谱与输入光谱之间的关系满足：
而利用傅里叶逆变换可以得到输出光脉冲的时域表达式：
其中，z为光传播距离，β2为群速度色散系数。当公式满足时域远场衍射条件时(β2z→∞)，积分仅在ω＝T/(β2z)+ω0处取非零值，因此输出光脉冲在T时刻的振幅只与输入光谱中ω＝T/(β2z)+ω0的频谱分量有关，输出光功率表示为：
由此可知，色散傅里叶变换会使得输出光脉冲时域与输入光脉冲光谱之间呈现一一对应关系，且色散量越大，相邻波长分量在时间上分得越开。因此能够实现光脉冲“频域-时间”映射。
图3为利用马赫曾德尔调制器对色散拉伸的光脉冲进行时域调制后，得到的对应编码变化的频谱形状示意图。由于经过色散拉伸后的光脉冲实现了频谱到时域的映射，因此，对时间波形强度进行时域上的伪随机调制，将会一一映射到对应的频谱位置，从而实现对脉冲频谱的伪随机编码。使马赫-曾德尔调制器工作在推挽模式，偏置在正交点，则电信号会改变内部折射率从而影响光信号的相位分布，进而影响强度分布，实现电光强度调制。需要注意的是，由于时域光脉冲呈周期排列，因此色散拉伸后的频谱也会时间映射到时域波形的多个周期。在进行强度调制时，伪随机信号的时钟周期需要和光脉冲保持一致，可以通过时钟同步产生伪随机信号的微波源与锁模激光器实现，调制速度能够达到μs级别。
所述预设伪随机编码可使用现有技术所使用的各种编码方案，优选采用基于Hadamard矩阵、傅里叶基、小波基等正交编码方案，相比非正交的伪随机编码方式，能够以更少的调制次数和数据量级实现对信号的高精度还原。图4即显示了基于Hadamrd矩阵的伪随机编码波形。如其中的(a)所示，该64*64的Hadamard矩阵由+1，-1构成，其中浅色的像素点代表“+1”，深色的像素点代表“-1”。为了用更适合于实际物理模型的1，0来表示对时间强度波形的调制，优化后的Hadamard矩阵可以表示为
其中的(b)是为了验证Hadamard矩阵调制的正交性给出的Hadamard矩阵自相关结果，从图中可以看出，HTH是一个标量矩阵，因此代入回式(3)中可以得到
从而能够以最小调制次数和数据量级恢复时域信号。
而考虑到光脉冲频谱自身强度起伏对目标探测造成影响，设脉冲频谱在频域上的分布可用矢量A表示，则只需要通过光谱仪预先获取光脉冲功率谱的包络形状，在目标的恢复重建过程中乘上A-1，表示为：
即可消除脉冲频谱自身强度起伏对目标探测的影响。
综上所述，通过本发明技术方案，可在现有频域鬼成像技术的基础上突破可编程光滤波器对频谱编码分辨率和调制速度的限制，通过频时映射的方式实现对脉冲光谱强度的高精度编码。具有成本低、结构简单，高分辨率、高速采集的优势，解决了频域高精度调制探测的难题，可实现光谱的高精度检测。
Claims (10)
1.一种频域鬼成像光谱检测方法，采集待测目标在频域预编码结构光场中的光强信息，然后通过将采集到的光强信息与所述频域预编码结构光场作互相关运算以重建出待测目标频谱；其特征在于，所述频域预编码结构光场的生成方法如下：利用色散傅里叶变换对光脉冲进行时域拉伸，然后用预设的伪随机编码对时域拉伸后的光脉冲进行强度调制。
2.如权利要求1所述频域鬼成像光谱检测方法，其特征在于，通过大色散啁啾光纤光栅实现所述时域拉伸。
3.如权利要求1所述频域鬼成像光谱检测方法，其特征在于，使用马赫曾德尔调制器实现所述强度调制。
4.如权利要求1所述频域鬼成像光谱检测方法，其特征在于，使用锁模激光器生成所述光脉冲。
5.如权利要求1所述频域鬼成像光谱检测方法，其特征在于，时域拉伸后的相邻光脉冲之间互不重叠。
6.一种频域鬼成像光谱检测装置，包括：用于生成频域预编码结构光场的结构光生成模块，用于采集待测目标在所述频域预编码结构光场中的光强信息的光强采集模块，用于通过将采集到的光强信息与所述频域预编码结构光场作互相关运算以重建出待测目标频谱的互相关运算模块；其特征在于，所述结构光生成模块包括：
色散拉伸模块，用于利用色散傅里叶变换对光脉冲进行时域拉伸；
强度调制器，用于用预设的伪随机编码对时域拉伸后的光脉冲进行强度调制。
7.如权利要求6所述频域鬼成像光谱检测装置，其特征在于，所述色散拉伸模块为大色散啁啾光纤光栅。
8.如权利要求6所述频域鬼成像光谱检测装置，其特征在于，所述强度调制器为马赫曾德尔调制器。
9.如权利要求6所述频域鬼成像光谱检测装置，其特征在于，使用锁模激光器生成所述光脉冲。
10.如权利要求6所述频域鬼成像光谱检测装置，其特征在于，时域拉伸后的相邻光脉冲之间互不重叠。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111339117.9A CN113890624B (zh) | 2021-11-12 | 2021-11-12 | 一种频域鬼成像光谱检测方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111339117.9A CN113890624B (zh) | 2021-11-12 | 2021-11-12 | 一种频域鬼成像光谱检测方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113890624A true CN113890624A (zh) | 2022-01-04 |
| CN113890624B CN113890624B (zh) | 2024-02-23 |
Family
ID=79017952
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111339117.9A Active CN113890624B (zh) | 2021-11-12 | 2021-11-12 | 一种频域鬼成像光谱检测方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113890624B (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114609645A (zh) * | 2022-02-21 | 2022-06-10 | 复旦大学 | 一种兼具异步通信功能的关联成像系统 |
| CN114690236A (zh) * | 2022-06-01 | 2022-07-01 | 成都理工大学 | 针对束流强度分布的快速鬼成像方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| AU2009273954A1 (en) * | 2008-07-24 | 2010-01-28 | The Regents Of The University Of California | Apparatus and method for dispersive Fourier-transform imaging |
| CN106371201A (zh) * | 2016-11-03 | 2017-02-01 | 清华大学 | 基于计算鬼成像的傅里叶重叠关联成像系统及方法 |
| CN108683458A (zh) * | 2018-05-22 | 2018-10-19 | 华中科技大学 | 一种基于正交频分复用的时域鬼成像装置及方法 |
-
2021
- 2021-11-12 CN CN202111339117.9A patent/CN113890624B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| AU2009273954A1 (en) * | 2008-07-24 | 2010-01-28 | The Regents Of The University Of California | Apparatus and method for dispersive Fourier-transform imaging |
| CN106371201A (zh) * | 2016-11-03 | 2017-02-01 | 清华大学 | 基于计算鬼成像的傅里叶重叠关联成像系统及方法 |
| CN108683458A (zh) * | 2018-05-22 | 2018-10-19 | 华中科技大学 | 一种基于正交频分复用的时域鬼成像装置及方法 |
Non-Patent Citations (1)
| Title |
|---|
| 赵生妹等: "鬼成像原理与进展研究", Retrieved from the Internet <URL:nyzr.njupt.edu.cn/ch/reader/create_pdf.aspx?file_no=202101009&flag=1&journal_id=nyzr&year_id=2021> * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114609645A (zh) * | 2022-02-21 | 2022-06-10 | 复旦大学 | 一种兼具异步通信功能的关联成像系统 |
| CN114609645B (zh) * | 2022-02-21 | 2023-11-24 | 复旦大学 | 一种兼具异步通信功能的关联成像系统 |
| CN114690236A (zh) * | 2022-06-01 | 2022-07-01 | 成都理工大学 | 针对束流强度分布的快速鬼成像方法 |
| CN114690236B (zh) * | 2022-06-01 | 2022-08-02 | 成都理工大学 | 针对束流强度分布的快速鬼成像方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113890624B (zh) | 2024-02-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113890624B (zh) | 一种频域鬼成像光谱检测方法及装置 | |
| Wang et al. | Large time-bandwidth product microwave arbitrary waveform generation using a spatially discrete chirped fiber Bragg grating | |
| Nazarathy et al. | Spectral analysis of optical mixing measurements | |
| US6724783B2 (en) | Method and apparatus for arbitrary waveform generation using photonics | |
| Xie et al. | STFT based on bandwidth-scaled microwave photonics | |
| CN107219002B (zh) | 一种超高分辨率光谱测量方法及系统 | |
| CN111609875A (zh) | 一种基于啁啾连续光的数字域可调分布式光纤传感系统及其传感方法 | |
| CN107064919B (zh) | 基于光电复合的超宽带信号参数估计方法 | |
| CN107135005B (zh) | 基于光电复合的超宽带信号多路并行压缩采样方法 | |
| CN113359370A (zh) | 一种光学数模转换方法及装置 | |
| CN209265419U (zh) | 一种基于光子自发辐射的高速量子随机数发生器 | |
| CN108566244A (zh) | 多通道并行的光器件频谱响应测量方法及装置 | |
| CN104699986A (zh) | 基于有限新息率采样理论的太赫兹脉冲信号采样方法 | |
| CN109981198B (zh) | 一种基于多模光纤的宽带射频谱感知方法及系统 | |
| US20230011819A1 (en) | Intelligent real-time full-field measurement method and system for high-repetition-rate femtosecond pulse | |
| CN113091795B (zh) | 光电器件与信道的测量方法及系统、装置、介质 | |
| Mididoddi et al. | All-optical random sequence generation for compressive sensing detection of RF signals | |
| CN115941053A (zh) | 一种基于色散拉伸的频域鬼成像光谱测量方法及装置 | |
| Scott et al. | Continuous, real-time, full-field waveform measurements via spectral slicing and parallel digital coherent detection | |
| Luo et al. | Broadband optical chaos generation by constructing a simple hybrid feedback loop | |
| Ding et al. | Photonic compressive sampling of sparse broadband rf signals using a multimode fiber | |
| Scott et al. | Frequency-to-time-assisted interferometry for full-field optical waveform measurements with picosecond resolution and microsecond record lengths | |
| CN113660036B (zh) | 微波信号时频特性测量方法及装置 | |
| Zhou et al. | Optical Signal Sampling Based on Compressive Sensing with Adjustable Compression Ratio | |
| Azaña et al. | Study of optical pulses-fiber gratings interaction by means of joint time-frequency signal representations |
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 |