CN109253997B - 一种基于频率调制和空间编码的拉曼断层成像系统 - Google Patents
一种基于频率调制和空间编码的拉曼断层成像系统 Download PDFInfo
- Publication number
- CN109253997B CN109253997B CN201811262326.6A CN201811262326A CN109253997B CN 109253997 B CN109253997 B CN 109253997B CN 201811262326 A CN201811262326 A CN 201811262326A CN 109253997 B CN109253997 B CN 109253997B
- Authority
- CN
- China
- Prior art keywords
- raman
- sample
- image
- frequency modulation
- spatial
- 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
- 238000001069 Raman spectroscopy Methods 0.000 title claims abstract description 164
- 238000003325 tomography Methods 0.000 title claims abstract description 33
- 230000005284 excitation Effects 0.000 claims abstract description 72
- 238000000034 method Methods 0.000 claims abstract description 39
- 238000003384 imaging method Methods 0.000 claims abstract description 32
- 238000012545 processing Methods 0.000 claims abstract description 16
- 238000004364 calculation method Methods 0.000 claims abstract description 15
- 230000003287 optical effect Effects 0.000 claims abstract description 12
- 238000001914 filtration Methods 0.000 claims abstract description 6
- 239000011159 matrix material Substances 0.000 claims description 24
- 108010076504 Protein Sorting Signals Proteins 0.000 claims description 18
- 230000008859 change Effects 0.000 claims description 12
- 238000013178 mathematical model Methods 0.000 claims description 4
- 238000005457 optimization Methods 0.000 claims description 4
- 230000008569 process Effects 0.000 claims description 4
- 238000010276 construction Methods 0.000 claims description 2
- 238000011084 recovery Methods 0.000 claims 7
- 238000009792 diffusion process Methods 0.000 claims 2
- 230000035945 sensitivity Effects 0.000 abstract description 12
- 230000008901 benefit Effects 0.000 abstract description 4
- 239000000126 substance Substances 0.000 abstract description 2
- 239000000835 fiber Substances 0.000 description 16
- 238000001514 detection method Methods 0.000 description 15
- 238000013461 design Methods 0.000 description 13
- 238000013480 data collection Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000003321 amplification Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000001678 irradiating effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/62—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
- G01N21/63—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
- G01N21/65—Raman scattering
Landscapes
- Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
Abstract
本发明属于光学三维成像技术领域，公开了一种基于频率调制和空间编码的拉曼断层成像系统；包括激发源模块、成像体控制模块、信号收集模块、控制与计算模块。所述信号收集模块包括空间编码单元，滤波单元、大数值孔径透镜、高灵敏度单点微弱信号探测器。所述控制与计算模块包括计算机控制单元和计算机处理单元。本发明通过对激发光的频率调制实现多束激发光的同时照射样本，激发拉曼散射信号；通过对出射拉曼散射信号的空间编码，实现用高灵敏度单点微弱信号探测器对微弱拉曼散射信号的高灵敏度收集。本发明对已有拉曼断层成像系统的数据获取方式进行了改进，可以快速、高灵敏度地收集拉曼散射信号，具有能够快速、高灵敏度获取样本三维体积内化学组分定量分布的优点。
Description
技术领域
本发明属于光学三维成像技术领域，尤其涉及一种基于频率调制和空间编码的拉曼断层成像系统。
背景技术
目前，业内常用的现有技术是这样的：拉曼断层成像技术可获取三维体积内的化学成分及其空间分布信息。作为一种漫射光学三维成像技术，拉曼断层成像的三维重建具有严重的病态性。通常可通过获取更多、更全面的数据来减少三维重建问题的病态性。目前，拉曼断层成像技术的数据收集方式可归结为三大类：第一类称为“单点激发-单点探测”的方式；采用一周光纤围绕成像体，成像时使激光器发出的激光入射到其中一束光纤，然后选择其他的某根光纤将拉曼散射信号投射到探测器上。在数据收集方式下，为了获得更多的数据量，需要通过旋转样本或者变换激发光纤和收集光纤的位置来改变激发点和收集点的位置，获得样本不同位置出射的拉曼散射信号，但是不断改变激发点和收集点的位置获取数据会消耗大量的时间。第二类称为“单点激发-多点探测”的方式。这种方式采用一根激发光纤传导激光，然后利用一组收集光纤以不同角度围绕在样本周围，光纤排列成线性阵列，将拉曼散射信号投射到探测器上。利用这种方式收集数据，可以收集样本在同一个激发点下多个角度的拉曼散射信号；然而，数据量仍不足以满足断层成像三维重建需求，需要变化激发光束所在光纤获取不同的激发点位置；这种情况下，需要改变激发光通道和收集光通道来增大数据量，同样地，收集数据会消耗大量时间。第三类称为“多点激发-多点探测”的方式。这种方式采用并行探测的方式收集数据，多根收集光纤排成阵列将拉曼散射信号投射到探测器上；前端同样采用多根光纤传导激发光束，使其入射到样本上。这种方式的优点在于能够提供所测体积内很大动态范围的信息，但多根光纤携带激光束同时照射样本，在探测光纤端并不能区分收集的拉曼散射信号来自于哪个激发点，而是所有激发点的叠加信号。这种情况下收集拉曼散射信号，需要从和信号中区分拉曼散射信号来自于哪个激发点并将其分解，因此不能增大采集的数据量。综上所述，为了增大数据采集量，改善拉曼断层成像技术三维重建的图像质量，已有三种数据收集方式均需要变换激发点位置、分时或并行获取样本表面出射的拉曼散射信号，这就造成需要花费大量时间，从而影响了拉曼断层成像技术的成像速度。另一方面，拉曼散射截面很小，因此样本发射的拉曼散射光强度很弱。已有拉曼断层成像系统中均采用CCD相机或电子增强CCD相机作为探测器收集样本出射的拉曼散射光信号。由于单点PMT比CCD相机具有更高的信号放大倍数和检测灵敏度，因此采用单点PMT代替CCD相机有可能提供更好的检测灵敏度。
综上所述，现有技术存在的问题是：
(1)现有的单点激发-单点探测的方式为了获得更多的数据量，需要通过旋转样本或者变换激发光纤和收集光纤的位置来改变激发点和收集点的位置，获得样本不同位置出射的拉曼散射信号。
(2)现有的单点激发-多点探测的方式数据量仍不足以满足断层成像三维重建需求，需要变化激发光束所在光纤获取不同的激发点位置；这种情况下，需要改变激发光通道和收集光通道来增大数据量。
(3)现有多点激发-多点探测的方式在探测光纤端并不能区分收集的拉曼散射信号来自于哪个激发点，而是所有激发点的叠加信号，不能增大采集的数据量。
解决上述技术问题的难度和意义：
难度：如何改善现有拉曼断层成像技术中的数据收集速度，以及提高拉曼断层成像系统的成像灵敏度。
意义：采用基于频率调制和空间编码的拉曼断层成像系统，可以增大拉曼散射信号数据量的同时，极大地加快成像速度；同时，高灵敏度单点探测器的使用，可提高系统的检测灵敏度，可以实现快速、高灵敏度地拉曼断层成像。
发明内容
针对现有技术存在的问题，本发明提供了一种基于频率调制和空间编码的拉曼断层成像系统。
本发明是这样实现的，一种基于频率调制和空间编码的拉曼断层成像系统，所述基于频率调制和空间编码的拉曼断层成像系统包括：
激发源模块，用于将激发光分成多束；
成像体控制模块，用于调整样本位置以及固定样本，处于成像视野的位置；
信号收集模块，样本被激发后发射拉曼散射光，经过滤波单元去除瑞利散射光，经过大数值孔径收集透镜收集聚焦后由高灵敏度单点微弱信号探测器接收，将光信号转化为电信号后，传输至控制与计算模块进行存储和后续处理；
控制与计算模块，用于控制多束激发光束产生单元、频率调制单元、成像体控制模块、空间编码单元和高灵敏度单点微弱信号探测器；计算机处理单元包括基于频率调制和空间编码的拉曼断层成像图像恢复方法，用于设计频率调制方式、空间编码模式以及对采集的序列拉曼散射光信号进行图像恢复。
进一步，所述激发源模块的频率调制单元对照射到样本不同位置的多束激发光束进行不同的频率调制；频率调制单元包括空间光调制器，用于对多束激光束进行光强度的频率调制，确保照射到样本不同位置的激光束强度具有不同的调制频率，则样本中激发的拉曼散射光携带有不同的调制频率。
进一步，所述信号收集模块包括空间编码单元，滤波单元、大数值孔径透镜、高灵敏度单点微弱信号探测器。
进一步，所述控制与计算模块包括计算机控制单元和计算机处理单元；
计算机控制单元包括图像采集卡，用于控制多束激发光束产生单元、频率调制单元、成像体控制模块、空间编码单元和高灵敏度单点微弱信号探测器；
计算机处理单元包括基于频率调制和空间编码的拉曼断层成像图像恢复方法，用于设计频率调制方式、空间编码模式以及对采集的序列拉曼散射光信号进行图像恢复。
进一步，所述激发源模块包括连续波激光器、分光镜；
所述分光镜一端连接于激光器一端。
本发明的另一目的在于提供一种实现所述基于频率调制和空间编码的拉曼断层成像系统的基于频率调制和空间编码的拉曼断层成像图像恢复方法，所述基于频率调制和空间编码的拉曼断层成像图像恢复方法包括以下步骤：
步骤一，输出并传输光源；连续波激光器发射的激光经分光镜分为多束激发光，并传输到频率调制单元；
步骤二，设计频率调制模式、激发拉曼散射光；设计空间光调制器的编码模式，使空间光调制器每个像元通过的光具有不同的强度调制频率，且随着时间的流逝，空间光调制器的编码模式不同；
步骤三，设计空间编码模式；拉曼散射光经漫射传输到样本表面，设计空间编码的模式，使样本表面不同位置的拉曼散射光带上不同的位置信息，通过变换空间编码模式，高灵敏度单点微弱信号探测器收集多个序列数据；
步骤四，收集并存储信号；利用信号收集模块采集单点序列数据，并传递到控制与计算模块进行存储和后续处理；
步骤五，建立空间编码和频率调制模型；利用空间编码的物理过程，建立空间编码的数学模型，并构造空间编码矩阵S，建立高灵敏度单点微弱信号探测器采集的拉曼散射信号序列R与待恢复样本表面拉曼图像D之间的数学关系式；
步骤六，恢复图像；基于稀疏正则化策略建立目标函数，采用优化方法对不同激发点对应的样本表面拉曼图像进行恢复。
进一步，所述步骤二的空间光调制器每个像元的调制频率由单点采集时间和待恢复图像像素个数按如下原则确定：
(1)由单点采集时间t确定最高调制频率fmax，
(2)由待恢复图像像素个数N确定时间序列采集点数，
(3)由时间序列采集点数n确定最低调制频率fmin，
(4)由最高调制频率、最低调制频率和待恢复图像像素个数N，确定调制频率分辨率以及每个像元的调制频率。
进一步，所述步骤五建立高灵敏度单点微弱信号探测器采集的拉曼散射信号序列R与待恢复样本表面拉曼图像D之间的数学关系式：
R＝SD；
其中，空间编码矩阵S、高灵敏度单点微弱信号探测器采集的拉曼散射信号序列与待恢复样本表面拉曼图像D之间的数学关系式；构造方法包括：
步骤一，经过频率调制的激光束与样本内拉曼信号发射源相互作用产生拉曼散射光，经漫射传输到样本表面，记为待恢复样本表面拉曼图像；将待恢复样本表面拉曼图像标记为二维矩阵形式，并按照规律转换为列向量D；
步骤二，设计空间编码的模式，使样本表面不同位置的拉曼散射光信号带上独特的位置信息；将空间编码模式以权重标记为二维矩阵形式，并按照规律转换为行向量S；
步骤三，建立高灵敏度单点微弱信号探测器信号强度与空间编码模式、样本表面拉曼图像之间的数学关系式：
式中，j表示第j个空间编码模式，M表示待恢复样本表面拉曼图像的像素点个数；
步骤四，按照一定规律改变空间编码模式，S向量发生变化；空间编码模式共变化m次，有：
步骤五，建立高灵敏度单点微弱信号探测器采集的拉曼散射信号序列与待恢复样本表面拉曼图像之间的数学关系式：
R＝SD；
式中，R为空间编码模式变化m次对应的高灵敏度单点微弱信号探测器采集的拉曼散射信号序列；S称为空间编码矩阵，其中每一行代表一种编码模式；D表示样本表面拉曼图像序列；
进一步，所述步骤二对多束激光频率调制方法包括：
步骤一，将第i个激发点对应的待恢复样本表面拉曼图像记录为二维矩阵，并按照规律转换为列向量Xi，激发点的数量记为N；
步骤二，将样本表面拉曼图像标记为二维矩阵形式，并按照规律转换为列向量D，列向量D中的每个像素点是所有激发点产生的拉曼散射信号的线性组合；
步骤三，设计空间光调制器每个像元的调制频率，随着时间流逝，空间光调制器形成的编码模式发生变化，从而导致频率调制编码矩阵发生变化，有：
步骤四，建立固定激发点对应的待恢复样本拉曼散射图像与样本表面出射总拉曼图像之间的数学关系式：
D＝WX；
式中，X表示固定激发点对应的待恢复样本拉曼散射图像；
步骤五，建立高灵敏度单点微弱信号探测器采集的拉曼散射信号序列与待恢复的不同激发点对应的样本表面拉曼图像之间的数学关系：
R＝SWX；
记为Y＝AX，其中Y＝R，A＝SW。
进一步，所述步骤六恢复图像目标函数为：
min arg||Y-AX||2+β|X|1；
其中，β是正则化因子。
本发明的另一目的在于提供一种应用所述基于频率调制和空间编码的拉曼断层成像图像恢复方法的光学三维成像系统。
综上所述，本发明的优点及积极效果为：成像速度快，采用频率调制方式，可以多束激发光同时照射样本，并能准确区分开收集的信号来自于哪个激发点。由于采用频率调制和空间编码的方式，成像的数据量的要求为已有方法的一半，即成像速度至少可以提高一倍；系统灵敏度高，由于采用了高灵敏度单点探测器，同时结合空间编码方式，可以通过单点探测方式获取样本表面的出射拉曼图像。EMCCD的增益可达10^4倍，而单点微弱信号探测器可以达到更高的增益，例如PMT的增益可达10^6以上，即对信号的放大倍数PMT约是EMCCD的100倍；灵敏度可以用输入信号与输出信号的比值表示，即使用单点PMT的灵敏度约可以提高100倍，所以本发明系统具有更高的灵敏度。
附图说明
图1是本发明实施例提供的基于频率调制和空间编码的拉曼断层成像系统结构示意图；
图2是本发明实施例提供的基于频率调制和空间编码的拉曼断层成像系统的原理示意图；
图中：1、激发源模块；1-1、连续光激光器；1-2、分光镜；1-3、空间光调制器；2、成像体控制模块；3、信号收集模块；3-1、空间编码单元；3-2、滤波单元；3-3、大数值孔径透镜；3-4、单点PMT；4、控制与计算模块；4-1、计算机控制单元；4-2、计算机处理单元。
图3是本发明实施例提供的基于频率调制和空间编码的拉曼断层成像图像恢复方法流程图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
本发明采用高灵敏度单点探测器，同时结合空间编码方式，可以通过单点探测方式获取样本表面的出射拉曼图像；单点PMT比EMCCD具有更高的放大倍数和检测灵敏度，具有更高的灵敏度，本发明实施例以单点PMT为例进行说明。
下面结合附图对本发明的应用原理作详细的描述。
如图1所示，本发明的基于频率调制和空间编码的拉曼断层成像系统包括：激发源模块1、成像体控制模块2、信号收集模块3、控制与计算模块4。
激发源模块1，用于将激发光分成多束；频率调制单元对照射到样本不同位置的多束激发光束进行不同的频率调制；频率调制单元包括空间光调制器，用于对多束激光束进行光强度的频率调制，确保照射到样本不同位置的激光束强度具有不同的调制频率，则样本中激发的拉曼散射光携带有不同的调制频率；
成像体控制模块2，用于调整样本位置以及固定样本，使其处于成像视野的最佳位置。
信号收集模块3包括空间编码单元3-1，滤波单元3-2、大数值孔径透镜3-3、单点PMT3-4；样本被激发后发射拉曼散射光，经过滤波单元去除瑞利散射光，之后经过大数值孔径收集透镜收集聚焦后由单点PMT接收，将光信号转化为电信号后，传输至控制与计算模块进行存储和后续处理；
控制与计算模块4，包括计算机控制单元4-1和计算机处理单元4-2；计算机控制单元包括图像采集卡，用于控制多束激发光束产生单元、频率调制单元、成像体控制模块、空间编码单元和单点PMT；计算机处理单元包括基于频率调制和空间编码的拉曼断层成像图像恢复方法，用于设计频率调制方式、空间编码模式以及对采集的序列拉曼散射光信号进行图像恢复。
如图2所示，激发源模块1包括连续波激光器1-1、分光镜1-2，所述分光镜一端连接于激光器一端。
如图3所示，本发明实施例提供的基于频率调制和空间编码的拉曼断层成像图像恢复方法包括以下步骤：
S301：输出并传输光源；连续波激光器发射的激光经分光镜分为多束激发光，并传输到频率调制单元；
S302：设计频率调制模式、激发拉曼散射光；设计空间光调制器的编码模式，使空间光调制器每个像元通过的光具有不同的强度调制频率，且随着时间的流逝，空间光调制器的编码模式不同；
S303：设计空间编码模式；拉曼散射光经漫射传输到样本表面，设计空间编码的模式，使样本表面不同位置的拉曼散射光带上不同的位置信息，通过变换空间编码模式，单点PMT可以收集多个序列数据；
S304：收集并存储信号；利用信号收集模块采集单点序列数据，并传递到控制与计算模块进行存储和后续处理；
S305：建立空间编码和频率调制模型；利用空间编码的物理过程，建立空间编码的数学模型，并构造空间编码矩阵S，建立高灵敏度单点微弱信号探测器采集的拉曼散射信号序列R与待恢复样本表面拉曼图像D之间的数学关系式；
S306：恢复图像；基于稀疏正则化策略建立目标函数，采用合适的优化方法对不同激发点对应的样本表面拉曼图像进行恢复。
在本发明的优选实施例中：步骤S302：空间光调制器每个像元的调制频率由单点采集时间和待恢复图像像素个数按如下原则确定：
(1)由单点采集时间t确定最高调制频率fmax，即
(2)由待恢复图像像素个数N确定时间序列采集点数，即
(3)由时间序列采集点数n确定最低调制频率fmin，即
(4)由最高调制频率、最低调制频率和待恢复图像像素个数N，确定调制频率分辨率以及每个像元的调制频率。
经过频率调制的激发光携带有不同的强度调制频率，同时照射样本后，与样本内的拉曼信号发射源相互作用激发拉曼散射光；
在本发明的优选实施例中：步骤S305：建立高灵敏度单点微弱信号探测器采集的拉曼散射信号序列R与待恢复样本表面拉曼图像D之间的数学关系式：
R＝SD；
其中，空间编码矩阵S、高灵敏度单点微弱信号探测器采集的拉曼散射信号序列与待恢复样本表面拉曼图像D之间的数学关系式的构造包括：
步骤一，经过频率调制的激光束与样本内拉曼信号发射源相互作用产生拉曼散射光，经漫射传输到样本表面，记为待恢复样本表面拉曼图像；将待恢复样本表面拉曼图像标记为二维矩阵形式，并按照规律转换为列向量D；
步骤二，设计空间编码的模式，使样本表面不同位置的拉曼散射光信号带上独特的位置信息；将空间编码模式以权重标记为二维矩阵形式，并按照规律转换为行向量S；
步骤三，建立高灵敏度单点微弱信号探测器信号强度与空间编码模式、样本表面拉曼图像之间的数学关系式：
式中，j表示第j个空间编码模式，M表示待恢复样本表面拉曼图像的像素点个数；
步骤四，按照一定规律改变空间编码模式，S向量发生变化；空间编码模式共变化m次，有：
步骤五，建立高灵敏度单点微弱信号探测器采集的拉曼散射信号序列与待恢复样本表面拉曼图像之间的数学关系式：
R＝SD；
式中，R为空间编码模式变化m次对应的高灵敏度单点微弱信号探测器采集的拉曼散射信号序列；S称为空间编码矩阵，其中每一行代表一种编码模式；D表示样本表面拉曼图像序列。
利用频率调制的物理过程，建立频率调制的数学模型，则对多束激光频率调制方法包括：
步骤一，将第i个激发点对应的待恢复样本表面拉曼图像记录为二维矩阵，并按照规律转换为列向量Xi，激发点的数量记为N；
步骤二，将样本表面拉曼图像标记为二维矩阵形式，并按照规律转换为列向量D，列向量D中的每个像素点是所有激发点产生的拉曼散射信号的线性组合；
步骤三，设计空间光调制器每个像元的调制频率，随着时间流逝，空间光调制器形成的编码模式发生变化，从而导致频率调制编码矩阵发生变化，有：
步骤四，建立固定激发点对应的待恢复样本拉曼散射图像与样本表面出射总拉曼图像之间的数学关系式：
D＝WX；
式中，X表示固定激发点对应的待恢复样本拉曼散射图像。
步骤五，建立高灵敏度单点微弱信号探测器采集的拉曼散射信号序列与待恢复的不同激发点对应的样本表面拉曼图像之间的数学关系：
R＝SWX；
记为Y＝AX，其中Y＝R，A＝SW。
在本发明的优选实施例中：步骤S306：恢复图像；基于稀疏正则化策略建立目标函数，采用合适的优化方法对不同激发点对应的样本表面拉曼图像进行恢复。其目标函数为：
min arg|Y-AX||2+β|X|1；
其中，β是正则化因子。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
Claims (6)
1.一种基于频率调制和空间编码的拉曼断层成像系统，其特征在于，所述基于频率调制和空间编码的拉曼断层成像系统包括：
激发源模块，用于将激发光分成多束；
成像体控制模块，用于调整样本位置以及固定样本，处于成像视野的位置；
信号收集模块，样本被激发后发射拉曼散射光，经过滤波单元去除瑞利散射光，经过大数值孔径收集透镜收集聚焦后由高灵敏度单点微弱信号探测器接收，将光信号转化为电信号后，传输至控制与计算模块进行存储和后续处理；
控制与计算模块，用于控制多束激发光束产生单元、频率调制单元、成像体控制模块、空间编码单元和高灵敏度单点微弱信号探测器；计算机处理单元包括基于频率调制和空间编码的拉曼断层成像图像恢复方法，用于设计频率调制方式、空间编码模式以及对采集的序列拉曼散射光信号进行图像恢复；
所述激发源模块的频率调制单元对照射到样本不同位置的多束激发光束进行不同的频率调制；频率调制单元包括空间光调制器，用于对多束激光束进行光强度的频率调制，确保照射到样本不同位置的激光束强度具有不同的调制频率，则样本中激发的拉曼散射光携带有不同的调制频率；
所述信号收集模块包括空间编码单元，滤波单元、大数值孔径透镜、高灵敏度单点微弱信号探测器；
所述控制与计算模块包括计算机控制单元和计算机处理单元；
计算机控制单元包括图像采集卡，用于控制多束激发光束产生单元、频率调制单元、成像体控制模块、空间编码单元和高灵敏度单点微弱信号探测器；
计算机处理单元包括基于频率调制和空间编码的拉曼断层成像图像恢复方法，用于设计频率调制方式、空间编码模式以及对采集的序列拉曼散射光信号进行图像恢复；
所述激发源模块包括连续波激光器、分光镜；
所述分光镜一端连接于激光器一端。
2.一种实现权利要求1所述基于频率调制和空间编码的拉曼断层成像系统的基于频率调制和空间编码的拉曼断层成像图像恢复方法，其特征在于，所述基于频率调制和空间编码的拉曼断层成像图像恢复方法包括以下步骤：
步骤一，输出并传输光源；连续波激光器发射的激光经分光镜分为多束激发光，并传输到频率调制单元；
步骤二，设计频率调制模式、激发拉曼散射光；设计空间光调制器的编码模式，使空间光调制器每个像元通过的光具有不同的强度调制频率，且随着时间的流逝，空间光调制器的编码模式不同；
步骤三，设计空间编码模式；拉曼散射光经漫射传输到样本表面，设计空间编码的模式，使样本表面不同位置的拉曼散射光带上不同的位置信息，通过变换空间编码模式，高灵敏度单点微弱信号探测器收集多个序列数据；
步骤四，收集并存储信号；利用信号收集模块采集单点序列数据，并传递到控制与计算模块进行存储和后续处理；
步骤五，建立空间编码和频率调制模型；利用空间编码的物理过程，建立空间编码的数学模型，并构造空间编码矩阵S，建立高灵敏度单点微弱信号探测器采集的拉曼散射信号序列R与待恢复样本表面拉曼图像D之间的数学关系式；
步骤六，恢复图像；基于稀疏正则化策略建立目标函数，采用优化方法对不同激发点对应的样本表面拉曼图像进行恢复。
3.如权利要求2所述的基于频率调制和空间编码的拉曼断层成像图像恢复方法，其特征在于，所述步骤二的空间光调制器每个像元的调制频率由单点采集时间和待恢复图像像素个数按如下原则确定：
(1)由单点采集时间t确定最高调制频率fmax，
(2)由待恢复图像像素个数N确定时间序列采集点数，
(3)由时间序列采集点数n确定最低调制频率fmin，
(4)由最高调制频率、最低调制频率和待恢复图像像素个数N，确定调制频率分辨率以及每个像元的调制频率。
4.如权利要求2所述的基于频率调制和空间编码的拉曼断层成像图像恢复方法，其特征在于，所述步骤五建立高灵敏度单点微弱信号探测器采集的拉曼散射信号序列R与待恢复样本表面拉曼图像D之间的数学关系式：
R＝SD；
其中，空间编码矩阵S、高灵敏度单点微弱信号探测器采集的拉曼散射信号序列与待恢复样本表面拉曼图像D之间的数学关系式；构造方法包括：
步骤一，经过频率调制的激光束与样本内拉曼信号发射源相互作用产生拉曼散射光，经漫射传输到样本表面，记为待恢复样本表面拉曼图像；将待恢复样本表面拉曼图像标记为二维矩阵形式，并按照规律转换为列向量D；
步骤二，设计空间编码的模式，使样本表面不同位置的拉曼散射光信号带上独特的位置信息；将空间编码模式以权重标记为二维矩阵形式，并按照规律转换为行向量S；
步骤三，建立高灵敏度单点微弱信号探测器信号强度与空间编码模式、样本表面拉曼图像之间的数学关系式：
式中，j表示第j个空间编码模式，M表示待恢复样本表面拉曼图像的像素点个数；
步骤四，按照一定规律改变空间编码模式，S向量发生变化；空间编码模式共变化m次，有：
步骤五，建立高灵敏度单点微弱信号探测器采集的拉曼散射信号序列与待恢复样本表面拉曼图像之间的数学关系式：
R＝SD；
式中，R为空间编码模式变化m次对应的高灵敏度单点微弱信号探测器采集的拉曼散射信号序列；S称为空间编码矩阵，其中每一行代表一种编码模式；D表示样本表面拉曼图像序列。
5.如权利要求2所述的基于频率调制和空间编码的拉曼断层成像图像恢复方法，其特征在于，所述步骤二对多束激光频率调制方法包括：
步骤一，将第i个激发点对应的待恢复样本表面拉曼图像记录为二维矩阵，并按照规律转换为列向量Xi，激发点的数量记为N；
步骤二，将样本表面拉曼图像标记为二维矩阵形式，并按照规律转换为列向量D，列向量D中的每个像素点是所有激发点产生的拉曼散射信号的线性组合；
步骤三，设计空间光调制器每个像元的调制频率，随着时间流逝，空间光调制器形成的编码模式发生变化，从而导致频率调制编码矩阵发生变化，有：
步骤四，建立固定激发点对应的待恢复样本拉曼散射图像与样本表面出射总拉曼图像之间的数学关系式：
D＝WX；
式中，X表示固定激发点对应的待恢复样本拉曼散射图像；
步骤五，建立高灵敏度单点微弱信号探测器采集的拉曼散射信号序列与待恢复的不同激发点对应的样本表面拉曼图像之间的数学关系：
R＝SWX；
记为Y＝AX，其中Y＝R，A＝SW；
所述步骤六恢复图像目标函数为：
min arg||Y-AX||2+β|X|1；
其中，β是正则化因子。
6.一种应用权利要求1所述基于频率调制和空间编码的拉曼断层成像系统。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811262326.6A CN109253997B (zh) | 2018-10-27 | 2018-10-27 | 一种基于频率调制和空间编码的拉曼断层成像系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811262326.6A CN109253997B (zh) | 2018-10-27 | 2018-10-27 | 一种基于频率调制和空间编码的拉曼断层成像系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109253997A CN109253997A (zh) | 2019-01-22 |
| CN109253997B true CN109253997B (zh) | 2020-12-01 |
Family
ID=65045023
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811262326.6A Active CN109253997B (zh) | 2018-10-27 | 2018-10-27 | 一种基于频率调制和空间编码的拉曼断层成像系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109253997B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111650180B (zh) * | 2020-06-08 | 2021-07-16 | 西安电子科技大学 | 一种基于信号编码和空间压缩的拉曼光谱成像系统 |
| CN112957011B (zh) * | 2021-02-01 | 2022-06-24 | 西安电子科技大学 | 高灵敏度微弱荧光信号探测系统、方法、存储介质及应用 |
| CN114608701B (zh) * | 2022-04-21 | 2024-08-16 | 西北核技术研究所 | 一种激光光强三维空间分布的非接触测量系统及方法 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101692971A (zh) * | 2009-10-13 | 2010-04-14 | 西安电子科技大学 | 非接触式光学断层成像方法 |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7394858B2 (en) * | 2003-08-08 | 2008-07-01 | Intel Corporation | Systems and methods for adaptive bit loading in a multiple antenna orthogonal frequency division multiplexed communication system |
| EP2452311A1 (en) * | 2009-07-08 | 2012-05-16 | Technion Research And Development Foundation Ltd. | Method and system for super-resolution signal reconstruction |
| US9029800B2 (en) * | 2011-08-09 | 2015-05-12 | Palo Alto Research Center Incorporated | Compact analyzer with spatial modulation and multiple intensity modulated excitation sources |
| CN102324101B (zh) * | 2011-10-09 | 2013-04-03 | 西安电子科技大学 | 基于光学投影断层成像系统的被测物图像拼接方法 |
| US9357928B2 (en) * | 2012-01-31 | 2016-06-07 | Purdue Research Foundation | Vibrational photoacoustic tomography using Raman laser |
| CN103234479B (zh) * | 2013-04-10 | 2015-07-08 | 浙江大学 | 基于压缩感知脉冲编码调制的三维信息获取装置和方法 |
| KR101646072B1 (ko) * | 2015-04-22 | 2016-08-05 | 유승진 | 동영상 데이터 암호화 장치 및 방법 |
| WO2017053825A1 (en) * | 2015-09-25 | 2017-03-30 | The Regents Of The University Of Michigan | Biopsy device for coherent raman imaging |
-
2018
- 2018-10-27 CN CN201811262326.6A patent/CN109253997B/zh active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101692971A (zh) * | 2009-10-13 | 2010-04-14 | 西安电子科技大学 | 非接触式光学断层成像方法 |
Non-Patent Citations (1)
| Title |
|---|
| "基于随机变量交替方向乘子法的荧光分子断层成像";侯榆青;《光学学报》;20170731;第37卷(第7期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109253997A (zh) | 2019-01-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11867610B2 (en) | Dynamic high-speed high-sensitivity imaging device and imaging method | |
| CN110793633B (zh) | 基于集束光纤的单像元多光谱计算成像系统及成像方法 | |
| CN109253997B (zh) | 一种基于频率调制和空间编码的拉曼断层成像系统 | |
| Coles et al. | Simultaneous measurements of angular scattering and intensity scintillation in the atmosphere | |
| US11727294B2 (en) | Method and system for quantum information processing and computation | |
| CN114813673A (zh) | 多通道超分辨基因检测仪及其检测方法 | |
| CN105044897A (zh) | 基于稀疏约束的快速随机光学重构成像系统及方法 | |
| CN103674264A (zh) | 一种基于周期衍射关联成像的图像融合装置及其方法 | |
| CN111579486B (zh) | 基于低功率受激发射损耗的超分辨成像方法及成像系统 | |
| CN114966739B (zh) | 基于光纤激光相控阵的目标探测成像系统及方法 | |
| CN102692394A (zh) | 一种基于热透镜效应的二维成像方法及装置 | |
| CN113785187B (zh) | 使用射频多路复用激发数据对荧光成像进行光谱解混 | |
| CN104568683A (zh) | 一种纳米颗粒粒度测量装置及方法 | |
| CN109900355B (zh) | 成像方法及装置 | |
| CN107621463A (zh) | 图像重建方法、装置及显微成像装置 | |
| CN109900356B (zh) | 关联成像方法及装置 | |
| CN109557070B (zh) | 一种基于空间编码光的拉曼成像系统 | |
| US20130250088A1 (en) | Multi-color confocal microscope and imaging methods | |
| US10663393B2 (en) | Spectrum inspecting apparatus | |
| Anas et al. | Towards a fast and safe LED-based photoacoustic imaging using deep convolutional neural network | |
| CN106841013B (zh) | 一种基于量子关联成像技术辨别微观生物体的方法 | |
| US8169695B2 (en) | Optical sub-wavelength super resolution imaging system producing nanometric spatial resolution of a scanned subject | |
| RU2579640C1 (ru) | Конфокальный спектроанализатор изображений | |
| CN116338724A (zh) | 一种基于感兴趣区的opa激光雷达三维图像传感方法 | |
| CN116125492A (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 |