CN109856112B - 一种基于cars的断层扫描成像装置、检测系统及方法 - Google Patents

一种基于cars的断层扫描成像装置、检测系统及方法 Download PDF

Info

Publication number
CN109856112B
CN109856112B CN201711242905.XA CN201711242905A CN109856112B CN 109856112 B CN109856112 B CN 109856112B CN 201711242905 A CN201711242905 A CN 201711242905A CN 109856112 B CN109856112 B CN 109856112B
Authority
CN
China
Prior art keywords
sample
detected
light
laser beam
stokes
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
Application number
CN201711242905.XA
Other languages
English (en)
Other versions
CN109856112A (zh
Inventor
夏栩笙
刘金波
郭敬为
袁红
桑凤亭
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Dalian Institute of Chemical Physics of CAS
Original Assignee
Dalian Institute of Chemical Physics of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Dalian Institute of Chemical Physics of CAS filed Critical Dalian Institute of Chemical Physics of CAS
Priority to CN201711242905.XA priority Critical patent/CN109856112B/zh
Publication of CN109856112A publication Critical patent/CN109856112A/zh
Application granted granted Critical
Publication of CN109856112B publication Critical patent/CN109856112B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

本发明公开了一种基于CARS的断层扫描成像装置、检测系统及方法，其包括：探测光源，其输出由泵浦光及斯托克斯光所合成的探测激光光束；待测样品承载台；光学机械装置，其使得光束穿过待测样品时改变光束到座标原点距离l和入射角θ；滤光片，以滤除泵浦光以及斯托克斯光；以及探测器以对待测样品受激所产生的相干反斯托克斯拉曼散射光进行检测；其中，光束探测激光光束穿过待测样品，基于待测样品的三阶非线性效应所产生相干反斯托克斯拉曼散射光，经由滤光片滤除混杂于相干反斯托克斯拉曼散射光中的泵浦光以及斯托克斯光后，由探测器进行检测。本发明解决了已有技术不能满足对稳定分布样品浓度、温度场等信息进行大范围断层扫描等使用需求的缺陷。

Description

一种基于CARS的断层扫描成像装置、检测系统及方法
技术领域
本发明属于激光光谱技术领域，具体涉及一种基于CARS的断层扫描成像装置、检测系统及方法。
背景技术
CARS(相干反斯托克斯拉曼散射)光谱技术是一种三阶非线性激光光谱技术，利用分子的拉曼活性共振获取气态、液态和固态中分子的光谱信号。与普通的拉曼光谱技术相比，具有SN比高、信号不易受到荧光信号干扰的优点。但是目前的CARS技术的应用多注重一维、二维的成像。例如，探测均匀气体、液体中的样品浓度；或者利用CARS显微镜技术，通过扫描振镜的二维扫描快速取得样品的二维样品浓度分布；或者通过扫描获得火焰的二维温度场。但是上述技术尚不能满足对稳定分布样品浓度、温度场等信息进行断层扫描等使用需求。
发明内容
鉴于已有技术不能满足对稳定分布样品浓度、温度场等信息进行断层扫描等使用需求的缺陷，本发明的目的是要提供一种基于CARS的断层扫描成像装置，其结合了相干反斯托克斯拉曼光谱技术和Radon变换技术来对样品进行断层扫描，以得到待测样品浓度和其温度的三维分布。
为了实现上述目的，本发明的技术方案：
一种基于CARS的断层扫描成像装置，其特征在于，包括：
探测光源，其用以输出由泵浦光以及斯托克斯光所合成的探测激光光束；
待测样品承载台；
光学机械装置，其用以使得所述探测激光光束穿过承载于其上的待测样品时能够改变光束到座标原点距离l和入射角θ；
滤光片，其用以滤除泵浦光以及斯托克斯光；
以及探测器，其用以对待测样品受激所产生的相干反斯托克斯拉曼散射光进行检测；
其中，所述探测激光光束穿过待测样品，基于待测样品的三阶非线性效应所产生相干反斯托克斯拉曼散射光，经由滤光片滤除混杂于相干反斯托克斯拉曼散射光中的泵浦光以及斯托克斯光后，由探测器对待测样品受激所产生的相干反斯托克斯拉曼散射光进行检测。
基于上述方案，进一步优选的，
所述探测激光光束依据待测样品所对应的拉曼频移设置泵浦光以及斯托克斯光之间的频率差。
基于上述方案，进一步优选的，
所述装置依据待测样品检测需求设置探测激光光束以及探测器类型；若检测温度断层扫描参数则所述探测激光光束选择宽带激光光束，同时所述探测器选择光谱仪光纤；若检测样品浓度分布扫描参数则所述探测激光光束选择窄带激光光束，同时所述探测器选择Pin管。
基于上述方案，进一步优选的，
所述装置还包括数据分析器，其用以依据探测器所检测的相干反斯托克斯拉曼散射光信号分析获得待测样品各断层所对应的样品浓度分布扫描参数或者温度断层扫描参数。
基于上述方案，进一步优选的，
所述样品浓度分布扫描参数的分析过程为：若假定所述Pin管所检测到相干反斯托克斯拉曼散射光信号强度为I，泵浦光信号强度为Ip，斯托克斯光信号强度为Is；待测样品的浓度为f，则依据非线性光学的理论，上述物理量对应的关系为
Figure GDA0002431808650000021
其中，比例系数k通过测已知浓度的被测样品进行标定，对应的样品浓度分布扫描参数即样品浓度场为非均匀场f(x,y)，对应的样品中各处产生相干反斯托克斯拉曼散射光信号强度也是非均匀场，表示为
Figure GDA0002431808650000022
自上述公式(1)可知，鉴于探测器所检测的相干反斯托克斯拉曼散射光信号的电流强度反映的是I(x,y)沿着光线的积分值，则引入Radon变换进行求解，即若探测激光光束的入射光线到光轴座标原点距离和入射角分别为l和θ，入射光线所在的直线为L，则I(x,y)对应的Radon变换式如下式(2)所示函数
Figure GDA0002431808650000023
Figure GDA0002431808650000024
则公式(2)对应的傅里叶变换与和I(x,y)的傅里叶变换存在如下关系：
Figure GDA0002431808650000031
式中
Figure GDA0002431808650000032
为对入射光线到座标原点距离l作傅立叶变换，
Figure GDA0002431808650000033
表示I(x,y)的傅里叶变换
Figure GDA0002431808650000034
最后通过光学机械装置不断改变光束到座标原点距离l和入射角θ以获得全部
Figure GDA0002431808650000035
对所述
Figure GDA0002431808650000036
作逆傅里叶变换求出对应的I(x,y)后，基于
Figure GDA0002431808650000037
求出待测样品各断层所对应的样品浓度分布扫描参数f(x,y)。
基于上述方案，进一步优选的，
所述温度断层扫描参数的分析过程为：
首先通过光谱仪光纤探测待测样品中多个转动能级的粒子分布构成的连续转动光谱P支光谱以及Q支光谱所构成的函数Iλ(l,θ)＝∫LI(x,y)dL，并通过光学机械装置遍历l和θ；其中下标λ表示波长为λ的光谱线在光谱仪上的积分；
其次依据所述P支光谱和Q支光谱计算对应的温度，所述计算过程为
对Iλ(l,θ)中的l作傅里叶变换，得到
Figure GDA0002431808650000038
再对kx,ky分别作傅里叶逆变换，得到
Figure GDA0002431808650000039
最后根据光谱分布函数Iλ(x,y)可以求解出温度分布函数T(x,y)。
发明的另一目的是要提供一种基于CARS的断层扫描成像的检测系统，其包括上述任意一项方案所述的CARS的断层扫描成像装置。
发明的另一目的是要提供一种基于CARS的断层扫描成像的检测方法，其基于上述任意一项方案所述的CARS的断层扫描成像装置进行检测，其具体包括如下步骤：
步骤1、将待测样品装入待测样品承载台中；
步骤2、基于待测样品所对应的拉曼频移数据设置探测光源；
步骤3、通过探测器、光学机械装置使得激光光束进入样品时在某一座标原点距离l和某一入射角θ条件下待测样品受激所产生的相干反斯托克斯拉曼散射光进行检测；
步骤4、通过光学机械装置依次改变光光束到座标原点距离和入射角获得全部检测数据。
所述方法还包括步骤5、依据探测器所检测的相干反斯托克斯拉曼散射光信号分析获得待测样品断层所对应的样品浓度分布扫描参数或者温度断层扫描参数。
基于上述方案，进一步优选的，
所述探测激光光束依据待测样品所对应的拉曼频移设置泵浦光以及斯托克斯光之间的频率差且依据待测样品检测需求设置探测激光光束以及探测器类型；若检测温度断层扫描参数则所述探测激光光束选择宽带激光光束，同时所述探测器选择光谱仪光纤；若检测样品浓度分布扫描参数则所述探测激光光束选择窄带激光光束，同时所述探测器选择Pin管。
基于上述方案，进一步优选的，
所述样品浓度分布扫描参数的分析过程为：若假定所述Pin管所检测到相干反斯托克斯拉曼散射光信号强度为I，泵浦光信号强度为Ip，斯托克斯光信号强度为Is；待测样品的浓度为f，则依据非线性光学的理论，上述物理量对应的关系为
Figure GDA0002431808650000041
其中，比例系数k通过测已知浓度的被测样品进行标定，对应的样品浓度分布扫描参数即样品浓度场为非均匀场f(x,y)，对应的样品中各处产生相干反斯托克斯拉曼散射光信号强度也是非均匀场，表示为
Figure GDA0002431808650000042
自上述公式(1)可知，鉴于探测器所检测的相干反斯托克斯拉曼散射光信号的电流强度反映的是I(x,y)沿着光线的积分值，则引入Radon变换进行求解，即若探测激光光束的入射光线到光轴座标原点距离和入射角分别为l和θ，入射光线所在的直线为L，则I(x,y)对应的Radon变换式如下式(2)所示函数
Figure GDA0002431808650000043
Figure GDA0002431808650000044
则公式(2)对应的傅里叶变换与和I(x,y)的傅里叶变换存在如下关系：
Figure GDA0002431808650000051
式中
Figure GDA0002431808650000052
为对入射光线到座标原点距离l作傅立叶变换，
Figure GDA0002431808650000053
表示I(x,y)的傅里叶变换
Figure GDA0002431808650000054
最后通过光学机械装置不断改变光束到座标原点距离l和入射角θ以获得全部
Figure GDA0002431808650000055
对所述
Figure GDA0002431808650000056
作逆傅里叶变换求出对应的I(x,y)后，基于
Figure GDA0002431808650000057
求出待测样品各断层所对应的样品浓度分布扫描参数f(x,y)；
所述温度断层扫描参数的分析过程为：
首先通过光谱仪光纤探测待测样品中多个转动能级的粒子分布构成的连续转动光谱即P支光谱以及Q支光谱；
其次依据所述P支光谱和Q支光谱计算对应的温度，所述计算过程为
基于Radon变换原理获得待测样品中断层中CARS光谱分布函数Iλ(x,y)，且由于Iλ表示波长为λ的光谱线在光谱仪上的积分，即
Figure GDA0002431808650000058
最后计算出温度断层扫描参数即温度分布T(x,y)。首先通过光谱仪光纤探测待测样品中多个转动能级的粒子分布构成的连续转动光谱P支光谱以及Q支光谱所构成的函数Iλ(l,θ)＝∫LI(x,y)dL，并通过光学机械装置遍历l和θ；其中下标λ表示波长为λ的光谱线在光谱仪上的积分；
其次依据所述P支光谱和Q支光谱计算对应的温度，所述计算过程为
对Iλ(l,θ)中的l作傅里叶变换，得到
Figure GDA0002431808650000059
再对kx,ky分别作傅里叶逆变换，得到
Figure GDA00024318086500000510
最后根据光谱分布函数Iλ(x,y)可以求解出温度分布函数T(x,y)。
与现有技术相比，本发明的有益效果：
为了解决前述问题，本发明结合了相干反斯托克斯拉曼光谱技术和Radon变换技术，实现了对待测样品进行断层扫描并得到待测物种浓度和其温度的三维分布；与普通的窄带CARS光谱技术相比，其实现了测量样品内部的物种浓度分布并作断层扫描；和共焦CARS显微镜技术相比，其实现了对体积比较大的样品进行扫描，并克服了共焦CARS显微镜技术探测深度所受到镜头焦深的限制；因而可以说本发明扩展了CARS光谱技术的应用领域。
附图说明
图1为本发明所述装置的原理结构示意图；
图2是本发明所述装置-用于测量静态样品的优化结构示意图；
图3是本发明所述装置-化工现场测量的优化结构示意图。
图中：1-探测器，2-待测样品承载台，3-滤波片，4-激光光束，5-激光光源，6.光电机械装置，7-反射镜，8-产生稳定的部分预混层流火焰的化工反应器，9、光轴。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
针对以上应用需求和目前技术的不足之处，本发明借鉴计算机断层扫描(CT)的概念，提出了一种基于CARS的断层扫描装置，以对具有稳定分布样品浓度、温度场等信息进行断层扫描，进而进行三维重建，所述断层扫描是利用探测激光光束与灵敏度极高的探测器以逐一改变方位即其与座标原点距离和入射角的方式对各个断面进行扫描；而这种三维重建在化工流体力学、航空发动机流场分析中十分具有应用潜力。
基于上述设计目的，如图1所示原理：一种基于CARS的断层扫描成像装置，其特征在于，包括：探测光源5，其用以输出由泵浦光以及斯托克斯光所合成的探测激光光束4；待测样品承载台2；光学机械装置6，其用以使得所述探测激光光束穿过承载于其上的待测样品时能够改变激光光束到光轴座标原点距离l和入射角θ；滤光片3，其用以滤除泵浦光以及斯托克斯光使得相干反斯托克斯拉曼散射光通过；以及探测器1，其用以对待测样品受激所产生的相干反斯托克斯拉曼散射光进行检测；其中，所述探测激光光束穿过待测样品，基于待测样品的三阶非线性效应所产生相干反斯托克斯拉曼散射光，经由滤光片滤除混杂于相干反斯托克斯拉曼散射光中的泵浦光以及斯托克斯光后，由探测器对待测样品受激所产生的相干反斯托克斯拉曼散射光进行检测。
基于上述方案，进一步的优选例1，所述探测器的选取根据具体要求选择光电探头或者光谱仪光纤头中任意一种。
基于上述优选例1，进一步的优选例2，所述探测激光光束依据待测样品所对应的拉曼频移设置泵浦光以及斯托克斯光即Stokes光之间的频率差，且所述泵浦光和Stokes光首先合为一束细而平行度高的激光束4后次穿过待测样品、滤光片照射到光电探测器(或光谱仪光纤头)上；在穿过待测样品时由于待测样品的三阶非线性效应，产生了相干反斯托克斯拉曼散射光即CARS光。
基于上述优选例2，进一步的优选例3，所述光学机械装置具有承载台且能够实现使得承载台进行平移或者旋转，进而改变激光光束到座标原点距离l和入射角θ。
基于上述优选例3，进一步的优选例4，所述所述装置还包括数据分析器，其用以依据探测器所检测的相干反斯托克斯拉曼散射光信号分析获得待测样品断层所对应的样品浓度分布扫描参数或者温度断层扫描参数。
基于上述优选例4，进一步的优选例5，本发明通过上述装置对样品作样品浓度分布扫描的原理是利用窄带CARS光谱和Radon变换；其中Radon变换成像原理和医学影像学CT类似，不同的是，医疗CT采用X光吸收光谱成像，而本发明采用CARS光谱成像；具体的所述样品浓度分布扫描参数的分析过程为：若假定所述探测器所检测到相干反斯托克斯拉曼散射光信号强度为I，泵浦光信号强度为Ip，斯托克斯光信号强度为Is；待测样品的浓度为f，则依据非线性光学的理论，上述物理量对应的关系为
Figure GDA0002431808650000071
其中，比例系数k通过测已知浓度的被测样品进行标定，对应的样品浓度分布扫描参数即样品浓度场为非均匀场f(x,y)，对应的样品中各处产生相干反斯托克斯拉曼散射光信号强度也是非均匀场，表示为
Figure GDA0002431808650000072
自上述公式(1)可知，鉴于探测器所检测的相干反斯托克斯拉曼散射光信号的电流强度反映的是I(x,y)沿着光线的积分值，则引入Radon变换进行求解，即若探测激光光束的入射光线到座标原点距离和入射角分别为l和θ，入射光线所在的直线为L，则I(x,y)对应的Radon变换式如下式(2)所示光强对于l和θ的函数
Figure GDA0002431808650000081
(l,θ)
Figure GDA0002431808650000082
自上式可得出探测器直接探测到的物理量就是
Figure GDA0002431808650000083
(l,θ)；
则公式(2)对应的傅里叶变换与和I(x,y)的傅里叶变换存在如下关系：
Figure GDA0002431808650000084
式中
Figure GDA0002431808650000085
为对入射光线到座标原点距离l作傅立叶变换，
Figure GDA0002431808650000086
表示I(x,y)的傅里叶变换
Figure GDA0002431808650000087
最后通过光学机械装置不断改变光束到座标原点距离l和入射角θ以获得全部
Figure GDA0002431808650000088
对所述
Figure GDA0002431808650000089
作逆傅里叶变换求出对应的I(x,y)后，基于
Figure GDA00024318086500000810
求出待测样品断层所对应的样品浓度分布扫描参数f(x,y)。
基于上述优选例5，进一步的优选例6，所述数据分析器还用以依据光电探测器所检测的相干反斯托克斯拉曼散射光信号分析获得温度断层扫描参数。
其对应的测量原理利用了宽带CARS光谱和Radon变换实现的；由于可以通过宽带CARS光谱测量相邻两个转动能级的粒子数之比，并利用Boltzman分布公式，即下式
Figure GDA00024318086500000811
来计算温度T，其中k为Boltzman常数，ΔE为两个转动能级的能量之差；且由于基于样品中单点的CARS光谱Iλ计算单点的温度T已经有成熟方法，则可以考虑利用Radon变换取得Iλ(x,y)，
基于上述原理，在本发明所述温度断层扫描参数的分析过程为：
首先通过光谱仪光纤探测待测样品中多个转动能级的粒子分布构成的连续转动光谱P支光谱以及Q支光谱所构成的函数Iλ(l,θ)＝∫LI(x,y)dL，并通过光学机械装置遍历l和θ；其中下标λ表示波长为λ的光谱线在光谱仪上的积分；
其次依据所述P支光谱和Q支光谱计算对应的温度，所述计算过程为
对Iλ(l,θ)中的l作傅里叶变换，得到
Figure GDA0002431808650000091
再对kx,ky分别作傅里叶逆变换，得到
Figure GDA0002431808650000092
最后根据光谱分布函数Iλ(x,y)可以求解出温度分布函数T(x,y)。
具体求解方法为，利用Boltzman公式计算不同温度下转动能级上的粒子数，通过粒子数计算理论上的光谱加宽线型，再根据实验测到的光谱的加宽线型进行拟合温度T。
本发明的另一目的是要提供一种基于CARS的断层扫描成像的检测系统，其包括上述任意一项方案所述的CARS的断层扫描成像装置。
发明的另一目的是要提供一种基于CARS的断层扫描成像的检测方法，其基于上述任意一项方案所述的CARS的断层扫描成像装置进行检测，其具体包括如下步骤：
步骤1、将待测样品装入待测样品承载台中；
步骤2、基于待测样品所对应的拉曼频移数据设置探测光源；
步骤3、通过探测器、光学机械装置使得激光光束进入样品时在某一座标原点距离l和某一入射角θ条件下待测样品受激所产生的相干反斯托克斯拉曼散射光进行检测；
步骤4、通过光学机械装置依次改变光光束到座标原点距离和入射角获得全部检测数据。
步骤5、依据探测器所检测的相干反斯托克斯拉曼散射光信号分析获得待测样品断层所对应的样品浓度分布扫描参数或者温度断层扫描参数。
步骤2中所述探测激光光束依据待测样品所对应的拉曼频移设置泵浦光以及斯托克斯光之间的频率差且依据待测样品检测需求设置探测激光光束以及探测器类型；若检测温度断层扫描参数则所述探测激光光束选择宽带激光光束，同时所述探测器选择光谱仪光纤；若检测样品浓度分布扫描参数则所述探测激光光束选择窄带激光光束，同时所述探测器选择Pin管。
步骤5中的所述样品浓度分布扫描参数的分析过程为：若假定所述Pin管所检测到相干反斯托克斯拉曼散射光信号强度为I，泵浦光信号强度为Ip，斯托克斯光信号强度为Is；待测样品的浓度为f，则依据非线性光学的理论，上述物理量对应的关系为
Figure GDA0002431808650000101
其中，比例系数k通过测已知浓度的被测样品进行标定，对应的样品浓度分布扫描参数即样品浓度场为非均匀场f(x,y)，对应的样品中各处产生相干反斯托克斯拉曼散射光信号强度也是非均匀场，表示为
Figure GDA0002431808650000102
自上述公式(1)可知，鉴于探测器所检测的相干反斯托克斯拉曼散射光信号的电流强度反映的是I(x,y)沿着光线的积分值，则引入Radon变换进行求解，即若探测激光光束的入射光线到光轴座标原点距离和入射角分别为l和θ，入射光线所在的直线为L，则I(x,y)对应的Radon变换式如下式(2)所示函数
Figure GDA0002431808650000103
Figure GDA0002431808650000104
则公式(2)对应的傅里叶变换与和I(x,y)的傅里叶变换存在如下关系：
Figure GDA0002431808650000105
式中
Figure GDA0002431808650000106
为对入射光线到座标原点距离l作傅立叶变换，
Figure GDA0002431808650000107
表示I(x,y)的傅里叶变换
Figure GDA0002431808650000108
最后通过光学机械装置不断改变光束到座标原点距离l和入射角θ以获得全部
Figure GDA0002431808650000109
对所述
Figure GDA00024318086500001010
作逆傅里叶变换求出对应的I(x,y)后，基于
Figure GDA00024318086500001011
求出待测样品各断层所对应的样品浓度分布扫描参数f(x,y)；
所述温度断层扫描参数的分析过程为：
首先通过光谱仪光纤探测待测样品中多个转动能级的粒子分布构成的连续转动光谱P支光谱以及Q支光谱所构成的函数Iλ(l,θ)＝∫LI(x,y)dL，并通过光学机械装置遍历l和θ；其中下标λ表示波长为λ的光谱线在光谱仪上的积分；
其次依据所述P支光谱和Q支光谱计算对应的温度，所述计算过程为
对Iλ(l,θ)中的l作傅里叶变换，得到
Figure GDA0002431808650000111
再对kx,ky分别作傅里叶逆变换，得到
Figure GDA0002431808650000112
最后根据光谱分布函数Iλ(x,y)可以求解出温度分布函数T(x,y)。
下述以具体实施例对本发明所述方案作进一步说明及论证：
实施例1：采用本发明所述装置对某种冷冻生物组织样品中NO的含量进行检测；其中NO是动物体内一种的小分子，在很多生理活动中起重要的信号传递作用，检测不同部位NO的浓度分布具有临床意义。
其实际实验装置如附图2所示；由于NO的拉曼频移为1874.35cm-1，则探测器(采用PIN管)，其采用波长为532nm的固体激光作为泵浦光(绿色)，则相应Stokes光的波长为590.92nm(黄色)；又由于本实验只需要测浓度，因此用可调谐窄带染料激光器即可，染料用RhB；将532nm的固体激光和590.92nm的染料激光合束之后，通过反射镜7射入样品中，进而产生波长为483.76nm的CARS光(应为蓝色)；同时选择选取截止波长为500nm的短通滤波片安装在探测器上；优选的，如图所述光电机械装置的载物台可以经由轴承、丝杠、电机等进行二维调整，即能够平移l和转动θ；随后遍历所有的l和θ，并记录PIN管响应I(θ,l)；最后通过对I(θ,l)中l作逆Radon变换，求出I(x,y)，进一步基于
Figure GDA0002431808650000113
求出断层NO的浓度分布f(x,y)。
实施例2：采用本发明所述装置对火焰温度进行断层扫描；由于稳定火焰温度是燃烧诊断研究中最常用到也是最重要的物理量，是对燃烧过程最直观的描述；则测量火焰温度就能够了解燃烧过程中发生的各种物理化学反应的机理，了解燃烧过程，评估燃烧能量利用，进而对火箭发动机、化工反应器总体评价。
所述装置如图3所示，其中8是能够产生稳定预混火焰的化工反应器；由于氮气的拉曼频移是2330cm-1左右则对应激光光源中Stokes波长中心为607.27nm；同时为了获得氮气在高温高压下加宽的光谱，需要采用宽带CARS技术，因此激光采用共线的532nm、线宽0.3cm-1窄带泵浦光和中心波长607nm、线宽5nm的宽带染料激光器，染料采用RhB；双波长激光射出，穿过预混火焰，达到探测器，其采用光谱仪探头采集光谱。其运动控制装置即光电机械装置，通过轴承、丝杠、电机控制光轴的旋转角θ和光轴9的偏移量l；遍历θ和l后得到光谱Iλ(θ,l)，对Iλ(θ,l)做逆Radon变换，就得到了断层中每一个点所对应的CARS光谱Iλ(x,y)。进而根据光谱的形状可以求出每一个点的氮气的温度。
综上所述，本发明实现了对待测样品进行断层扫描并得到待测样品浓度和其温度的三维分布情况，其扩展了CARS光谱技术的应用领域。
以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，根据本发明的技术方案及其发明构思加以等同替换或改变，都应涵盖在本发明的保护范围之内。

Claims (7)

1.一种基于CARS的断层扫描成像装置，其特征在于，包括：
探测光源，其用以输出由泵浦光以及斯托克斯光所合成的探测激光光束；
待测样品承载台；
光学机械装置，其用以使得所述探测激光光束穿过承载于其上的待测样品时能够改变光束到座标原点距离l和入射角θ；
滤光片，其用以滤除泵浦光以及斯托克斯光；
以及探测器，其用以对待测样品受激所产生的相干反斯托克斯拉曼散射光进行检测；
其中，所述探测激光光束穿过待测样品，基于待测样品的三阶非线性效应所产生相干反斯托克斯拉曼散射光，经由滤光片滤除混杂于相干反斯托克斯拉曼散射光中的泵浦光以及斯托克斯光后，由探测器对待测样品受激所产生的相干反斯托克斯拉曼散射光进行检测；
所述装置依据检测需求设置探测激光光束以及探测器类型；若检测温度断层扫描参数则所述探测激光光束选择宽带激光光束，同时所述探测器选择光谱仪光纤；若检测样品浓度分布扫描参数则所述探测激光光束选择窄带激光光束，同时所述探测器选择Pin管；
所述装置还包括数据分析器，其用以依据探测器所检测的相干反斯托克斯拉曼散射光信号分析获得待测样品各断层所对应的样品浓度分布扫描参数或者温度断层扫描参数；
所述样品浓度分布扫描参数的分析过程为：若假定所述Pin管所检测到相干反斯托克斯拉曼散射光信号强度为I，泵浦光信号强度为Ip，斯托克斯光信号强度为Is；待测样品的浓度为f，则依据非线性光学的理论，上述物理量对应的关系为
Figure FDA0002444712630000011
其中，比例系数k通过测已知浓度的被测样品进行标定，对应的样品浓度分布扫描参数即样品浓度场为非均匀场f(x,y)，对应的样品中各处产生相干反斯托克斯拉曼散射光信号强度也是非均匀场，表示为
Figure FDA0002444712630000012
自上述公式(1)可知，鉴于探测器所检测的相干反斯托克斯拉曼散射光信号的电流强度反映的是I(x，y)沿着光线的积分值，则引入Radon变换进行求解，即若探测激光光束的入射光线到光轴座标原点距离和入射角分别为l和θ，入射光线所在的直线为L，则I(x，y)对应的Radon变换式如下式(2)所示函数
Figure FDA0002444712630000021
Figure FDA0002444712630000022
则公式(2)对应的傅里叶变换与和I(x，y)的傅里叶变换存在如下关系：
Figure FDA0002444712630000023
式中
Figure FDA0002444712630000024
为对入射光线到座标原点距离l作傅立叶变换，
Figure FDA00024447126300000211
表示I(x，y)的傅里叶变换
Figure FDA0002444712630000025
最后通过光学机械装置不断改变光束到座标原点距离l和入射角θ，遍历所有可能的值，以获得全部
Figure FDA0002444712630000026
对所述
Figure FDA0002444712630000027
作逆傅里叶变换求出对应的I(x，y)后，基于
Figure FDA0002444712630000028
求出待测样品各断层所对应的样品浓度分布扫描参数f(x，y)。
2.根据权利要求1所述的基于CARS的断层扫描成像装置，其特征在于：
所述探测激光光束依据待测样品所对应的拉曼频移设置泵浦光以及斯托克斯光之间的频率差。
3.根据权利要求1所述的基于CARS的断层扫描成像装置，其特征在于：
所述温度断层扫描参数的分析过程为：
首先通过光谱仪光纤探测待测样品中多个转动能级的粒子分布构成的连续转动光谱P支光谱以及Q支光谱所构成的函数Iλ(l，θ)＝∫L I(x，y)dL，并通过光学机械装置遍历l和θ；其中下标λ表示波长为λ的光谱线在光谱仪上的积分；
其次依据所述P支光谱和Q支光谱计算对应的温度，所述计算过程为
对Iλ(l，θ)中的l作傅里叶变换，得到
Figure FDA0002444712630000029
再对kx，ky分别作傅里叶逆变换，得到
Figure FDA00024447126300000210
最后根据光谱分布函数Iλ(x，y)求解出温度分布函数T(x，y)。
4.一种基于CARS的断层扫描成像的检测系统，其包括权利要求1-3任意一项所述的断层扫描成像装置。
5.一种基于CARS的断层扫描成像的检测方法，其基于上述权利要求1-3任意一项所述的CARS的断层扫描成像装置进行检测，具体包括如下步骤：
步骤1、将待测样品装入待测样品承载台中；
步骤2、基于待测样品所对应的拉曼频移数据设置探测光源；
步骤3、通过探测器、光学机械装置使得激光光束进入样品时在某一座标原点距离l和某一入射角θ条件下待测样品受激所产生的相干反斯托克斯拉曼散射光进行检测；
步骤4、通过光学机械装置依次改变光光束到座标原点距离和入射角获得全部检测数据；
步骤5、依据探测器所检测的相干反斯托克斯拉曼散射光信号分析获得待测样品断层所对应的样品浓度分布扫描参数或者温度断层扫描参数。
6.根据权利要求5所述的基于CARS的断层扫描成像的检测方法，其特征在于：
所述探测激光光束依据待测样品所对应的拉曼频移设置泵浦光以及斯托克斯光之间的频率差且依据待测样品检测需求设置探测激光光束以及探测器类型；若检测温度断层扫描参数则所述探测激光光束选择宽带激光光束，同时所述探测器选择光谱仪光纤；若检测样品浓度分布扫描参数则所述探测激光光束选择窄带激光光束，同时所述探测器选择Pin管。
7.根据权利要求5所述的基于CARS的断层扫描成像的检测方法，其特征在于：
所述样品浓度分布扫描参数的分析过程为：若假定所述Pin管所检测到相干反斯托克斯拉曼散射光信号强度为I，泵浦光信号强度为Ip，斯托克斯光信号强度为Is；待测样品的浓度为f，则依据非线性光学的理论，上述物理量对应的关系为
Figure FDA0002444712630000031
其中，比例系数k通过测已知浓度的被测样品进行标定，对应的样品浓度分布扫描参数即样品浓度场为非均匀场f(x，y)，对应的样品中各处产生相干反斯托克斯拉曼散射光信号强度也是非均匀场，表示为
Figure FDA0002444712630000041
自上述公式(1)可知，鉴于探测器所检测的相干反斯托克斯拉曼散射光信号的电流强度反映的是I(x，y)沿着光线的积分值，则引入Radon变换进行求解，即若探测激光光束的入射光线到光轴座标原点距离和入射角分别为l和θ，入射光线所在的直线为L，则I(x，y)对应的Radon变换式如下式(2)所示函数
Figure FDA0002444712630000042
Figure FDA0002444712630000043
则公式(2)对应的傅里叶变换与和I(x，y)的傅里叶变换存在如下关系：
Figure FDA0002444712630000044
式中
Figure FDA0002444712630000045
为对入射光线到座标原点距离l作傅立叶变换，
Figure FDA00024447126300000412
表示I(x，y)的傅里叶变换
Figure FDA0002444712630000046
最后通过光学机械装置不断改变光束到座标原点距离l和入射角θ以获得全部
Figure FDA0002444712630000047
对所述
Figure FDA0002444712630000048
作逆傅里叶变换求出对应的I(x，y)后，基于
Figure FDA0002444712630000049
求出待测样品各断层所对应的样品浓度分布扫描参数f(x，y)；
所述温度断层扫描参数的分析过程为：
首先通过光谱仪光纤探测待测样品中多个转动能级的粒子分布构成的连续转动光谱P支光谱以及Q支光谱所构成的函数Iλ(l，θ)＝∫L I(x，y)dL，并通过光学机械装置遍历l和θ；其中下标λ表示波长为λ的光谱线在光谱仪上的积分；
其次依据所述P支光谱和Q支光谱计算对应的温度，所述计算过程为
对Iλ(l，θ)中的l作傅里叶变换，得到
Figure FDA00024447126300000410
再对kx，ky分别作傅里叶逆变换，得到
Figure FDA00024447126300000411
最后根据光谱分布函数Iλ(x，y)求解出温度分布函数T(x，y)。
CN201711242905.XA 2017-11-30 2017-11-30 一种基于cars的断层扫描成像装置、检测系统及方法 Active CN109856112B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711242905.XA CN109856112B (zh) 2017-11-30 2017-11-30 一种基于cars的断层扫描成像装置、检测系统及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711242905.XA CN109856112B (zh) 2017-11-30 2017-11-30 一种基于cars的断层扫描成像装置、检测系统及方法

Publications (2)

Publication Number Publication Date
CN109856112A CN109856112A (zh) 2019-06-07
CN109856112B true CN109856112B (zh) 2020-05-22

Family

ID=66888754

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711242905.XA Active CN109856112B (zh) 2017-11-30 2017-11-30 一种基于cars的断层扫描成像装置、检测系统及方法

Country Status (1)

Country Link
CN (1) CN109856112B (zh)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102116929A (zh) * 2011-01-30 2011-07-06 中国科学院西安光学精密机械研究所 高速宽视场相干反斯托克斯拉曼散射显微系统及方法
CN104713866A (zh) * 2013-12-15 2015-06-17 中国科学院大连化学物理研究所 一种宽带cars 探测1δ氧装置及其使用方法
CN105043559A (zh) * 2015-09-15 2015-11-11 哈尔滨工业大学 一种基于双焦透镜的cars光谱测温装置
CN105136329A (zh) * 2015-09-15 2015-12-09 哈尔滨工业大学 一种基于双焦透镜的cars光谱测温实验装置
CN105203222A (zh) * 2015-09-15 2015-12-30 哈尔滨工业大学 一种基于菲涅耳透镜和cars光谱对火焰一维扫描测温的装置
CN106226284A (zh) * 2016-08-23 2016-12-14 深圳大学 相干反斯托克斯拉曼散射光谱探测和显微成像系统及方法
CN106679810A (zh) * 2015-11-05 2017-05-17 中国科学院大连化学物理研究所 一种前向相干增强反斯托克斯拉曼光谱仪
CN106680259A (zh) * 2015-11-05 2017-05-17 中国科学院大连化学物理研究所 一种后向相干增强反斯托克斯拉曼光谱仪

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4618341B2 (ja) * 2008-06-26 2011-01-26 ソニー株式会社 コヒーレントアンチストークスラマン散乱光を利用した生体内物質量測定方法
WO2011099269A1 (ja) * 2010-02-10 2011-08-18 国立大学法人大阪大学 顕微鏡及び観察方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102116929A (zh) * 2011-01-30 2011-07-06 中国科学院西安光学精密机械研究所 高速宽视场相干反斯托克斯拉曼散射显微系统及方法
CN104713866A (zh) * 2013-12-15 2015-06-17 中国科学院大连化学物理研究所 一种宽带cars 探测1δ氧装置及其使用方法
CN105043559A (zh) * 2015-09-15 2015-11-11 哈尔滨工业大学 一种基于双焦透镜的cars光谱测温装置
CN105136329A (zh) * 2015-09-15 2015-12-09 哈尔滨工业大学 一种基于双焦透镜的cars光谱测温实验装置
CN105203222A (zh) * 2015-09-15 2015-12-30 哈尔滨工业大学 一种基于菲涅耳透镜和cars光谱对火焰一维扫描测温的装置
CN106679810A (zh) * 2015-11-05 2017-05-17 中国科学院大连化学物理研究所 一种前向相干增强反斯托克斯拉曼光谱仪
CN106680259A (zh) * 2015-11-05 2017-05-17 中国科学院大连化学物理研究所 一种后向相干增强反斯托克斯拉曼光谱仪
CN106226284A (zh) * 2016-08-23 2016-12-14 深圳大学 相干反斯托克斯拉曼散射光谱探测和显微成像系统及方法

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
Fourier-transform coherent anti-Stokes Raman scattering microscopy;Jennifer P. Ogilvie et al.;《OPTICS LETTERS》;20060215;第31卷(第4期);第480-482页 *
Phase-interfacial stimulated Raman scattering generated in strongly pumped water;HONG YUAN et al.;《Optics Letters》;20160715;第41卷(第14期);第3335-3338页 *
分时窄带ＣＡＲＳ光谱技术在测量烯烃浓度方面的应用;杨杰等;《光谱学与光谱分析》;20170430;第37卷(第4期);第1037-1041页 *
基于TDLAS二维温度场分布反演软件设计;陈福多;《中国优秀硕士学位论文全文数据库基础科学辑》;20170115;正文全文 *
基于相干反斯托克斯拉曼散射的二维温度场扫描测量;李国华等;《光学精密工程》;20160131;第24卷(第1期);第14-17页 *
激光吸收光谱技术测量非均匀燃烧流场研究进展;洪延姬等;《航空学报》;20150325;第36卷(第3期);第724-736页 *

Also Published As

Publication number Publication date
CN109856112A (zh) 2019-06-07

Similar Documents

Publication Publication Date Title
Marcelli et al. Biological applications of synchrotron radiation infrared spectromicroscopy
Reif et al. Analytical model of light reflectance for extraction of the optical properties in small volumes of turbid media
Swartling et al. Comparison of spatially and temporally resolved diffuse-reflectance measurement systems for determination of biomedical optical properties
Nothdurft et al. In vivo fluorescence lifetime tomography
Gulsen et al. Design and implementation of a multifrequency<? xpp qa?> near-infrared diffuse optical tomography system
Oldenburg et al. Spectroscopic optical coherence tomography and microscopy
Ozaki Recent advances in molecular spectroscopy of electronic and vibrational transitions in condensed phase and its application to chemistry
Patwardhan et al. Quantitative diffuse optical tomography for small animals using an ultrafast gated image intensifier
Ignatenko et al. Applications of Mueller matrix polarimetry to biological and agricultural diagnostics: a review
Zheng et al. Noninvasive blood glucose detection using a miniature wearable Raman spectroscopy system
Kirchner et al. Multiple illumination learned spectral decoloring for quantitative optoacoustic oximetry imaging
Leavesley et al. An excitation wavelength–scanning spectral imaging system for preclinical imaging
CN109856112B (zh) 一种基于cars的断层扫描成像装置、检测系统及方法
Matthews et al. Fourier domain multispectral multiple scattering low coherence interferometry
Wood et al. Molecular imaging of red blood cells by raman spectroscopy
US20070171433A1 (en) Systems and processes for providing endogenous molecular imaging with mid-infrared light
Chu et al. 2D full-field measurement of oxygen concentration based on the phase fluorometry technique that uses the four-frame integrating-bucket method
Sacco et al. Graphene edge method for three‐dimensional probing of Raman microscopes focal volumes
Morsin et al. Development of plasmonic sensor for detection of toxic materials
Pham et al. Quantifying the optical properties and chromophore concentrations of turbid media by chemometric analysis of hyperspectral diffuse reflectance data collected using a Fourier interferometric imaging system
Liger et al. Temperature Measurements by Wavelength Modulation Diode Laser Absorption Spectroscopy with Logarithmic Conversion and 1 f Signal Detection
Vishwanath et al. Fluorescence quenching by polystyrene microspheres in UV-visible and NIR tissue-simulating phantoms
Artemyev et al. Study of spurious optical signals in a fiber-optic Raman spectroscopy system
Wang et al. Advanced label-free laser scanning microscopy and its biological imaging application
Al-Qadi et al. Optical characterization and rotational dynamics observation of colloidal gold nanorods based on polarized light scattering microscopy

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