CN112763190B - 一种oct系统灵敏度测量方法以及测量装置 - Google Patents
一种oct系统灵敏度测量方法以及测量装置 Download PDFInfo
- Publication number
- CN112763190B CN112763190B CN202011564335.8A CN202011564335A CN112763190B CN 112763190 B CN112763190 B CN 112763190B CN 202011564335 A CN202011564335 A CN 202011564335A CN 112763190 B CN112763190 B CN 112763190B
- Authority
- CN
- China
- Prior art keywords
- signal
- sampling
- optical signal
- oct system
- attenuation
- 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
- 230000035945 sensitivity Effects 0.000 title claims abstract description 84
- 238000000034 method Methods 0.000 title claims abstract description 31
- 230000003287 optical effect Effects 0.000 claims abstract description 190
- 238000005070 sampling Methods 0.000 claims abstract description 140
- 238000001228 spectrum Methods 0.000 claims abstract description 37
- 239000000835 fiber Substances 0.000 claims description 26
- 238000003384 imaging method Methods 0.000 claims description 26
- 230000001427 coherent effect Effects 0.000 claims description 15
- 238000006243 chemical reaction Methods 0.000 claims description 14
- 238000000691 measurement method Methods 0.000 claims description 14
- 238000002310 reflectometry Methods 0.000 claims description 10
- 230000003595 spectral effect Effects 0.000 claims description 2
- 230000009466 transformation Effects 0.000 claims 1
- 230000002238 attenuated effect Effects 0.000 abstract description 9
- 239000013307 optical fiber Substances 0.000 description 12
- 230000000052 comparative effect Effects 0.000 description 11
- 238000005259 measurement Methods 0.000 description 11
- 238000003325 tomography Methods 0.000 description 9
- 230000005540 biological transmission Effects 0.000 description 8
- 230000007246 mechanism Effects 0.000 description 6
- 238000001514 detection method Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 3
- 238000004891 communication Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000033001 locomotion Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000035559 beat frequency Effects 0.000 description 1
- 238000003759 clinical diagnosis Methods 0.000 description 1
- 238000012733 comparative method Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000002526 effect on cardiovascular system Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011503 in vivo imaging Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 230000009545 invasion Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M11/00—Testing of optical apparatus; Testing structures by optical methods not otherwise provided for
- G01M11/02—Testing optical properties
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M11/00—Testing of optical apparatus; Testing structures by optical methods not otherwise provided for
- G01M11/02—Testing optical properties
- G01M11/0207—Details of measuring devices
Landscapes
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本申请公开了一种OCT系统灵敏度测量方法，包括如下步骤：扫描光源发出光信号；将光信号分成初始采样光信号和初始参考光信号；初始参考光信号进入参考臂；初始采样光信号进入采样臂，采样光信号在采样臂内发生2N+1次衰减；从参考臂反射出来的输出参考光信号与从采样臂反射出来的输出采样光信号发生干涉形成干涉光信号；将干涉光信号转换为电信号；将电信号转换为FFT频谱，通过FFT频谱得到信噪比，从而计算出OCT系统灵敏度；OCT系统灵敏度＝采样光信号的衰减量+信噪比，OCT系统灵敏度为S，采样光信号的衰减量为FA，信噪比为SNR，即S＝FA+SNR；采样光信号的衰减量＝2N×第一衰减量+第二衰减量；N≥1。本申请还提供一种OCT系统灵敏度测量装置。
Description
技术领域
本申请涉及光学干涉断层成像技术领域，具体涉及一种OCT系统灵敏度测量方法以及测量装置。
背景技术
光学干涉断层成像(OCT)是一种类似超声的生物医学成像技术。它通过测量样品背向散射光或反射光来获得生物组织的横切面图像和三维图像，具有非侵入、分辨率高、可在体成像等优点，目前已广泛应用于眼科、皮肤科和心血管内科等的临床诊断和研究上。
早期的OCT系统主要以时域OCT为主。为了提高OCT系统的成像速度，研究人员提出了频域OCT系统。与时域OCT相比，频域OCT具有更高的灵敏度。灵敏度反映了OCT系统对弱信号的探测能力。灵敏度越高，系统探测弱信号的能力越强，样品的成像深度越大，从而能够获得更多的样品结构信息。因此，对于OCT系统灵敏度的测试显得尤为重要。
光学干涉断层成像系统的灵敏度是指光学干涉断层成像系统的样品臂能允许的最大的对光信号衰减。影响光学干涉断层成像系统灵敏度的就是接收光信号的强度和系统噪声，接收光信号强度取决于光源输出光功率和光路衰减，系统噪声来自光源、无源光路和光接收。一般成像系统主要依赖于成像透镜的光学结构，而内窥式OCT是通过光纤透镜采集的光信号返回干涉仪与参考光进行干涉拍频实现高灵敏度光相干接收，采用平衡探测去除直流信号提高光电转换效率，完成光电转换再重建图像，光传输性能对OCT成像质量有直接影响。和光通信传输链路一样，OCT光传输链路的噪声来自3方面：光源噪声(量子噪声和模式噪声RIN、光路噪声MPI、探测器热噪声，扫频OCT光源采用的是极低噪声的单模激光器，模式噪声极低RIN<120dBc/Hz可以忽略，只有散粒噪声(量子噪声)，通过控制好光路的连接点降低MPI噪声、优质的光相干接收和平衡探测，能把MPI噪声和接收机热噪声也降到低于光源量子噪声，实现量子极限传输。在1310nm波长，SS-OCT系统能实现接近理论量子极限传输120dB左右灵敏度。
由于光接收机的有效动态范围只有30～40dB，测量120dB的灵敏度范围，必须引入固定衰减量，把光信号预先衰竭到接收机的最佳接收动态范围内，然后测量接收到的信号的信噪比，固定衰减量+信噪比就是该系统的光灵敏度。但是，引入固定衰减量之后，如何准确测量信噪比、如何得到整个系统的灵敏度，学术界和工业界还没有共识。目前最常用的方法是采用光谱分析仪OSA(optical spectrum analyzer)测量进入接收机的光信号的SNR、或者采用光束分析仪测量接收光的点扩展函数PSF(point spread function)，PSF的峰值与噪声的差值就是信噪比，灵敏度就是固定衰减量+信噪比。这个方法的缺点是OSA和光束分析仪里面的光电转换其实不同于OCT系统里的光电接收机，因此得到的信噪比并不代表整个OCT系统的噪声性能。
发明内容
有鉴于此，本申请提供一种OCT系统灵敏度测量方法，该方法以及装置既可以方便地调整采样臂内的衰减量，又保证了采样光信号的高质量，还去可以直接利用OCT系统本身来测量系统噪声，能够简便、准确地得到整个OCT系统的灵敏度，既不需要昂贵的光束分析仪和光谱分析仪OSA，又包括了OCT系统本身的光接收机噪声，非常适合在实际生产中使用。
为实现上述目的，根据本申请的一个方面，提供了一种OCT系统灵敏度测量方法。
本申请提供如下技术方案。
1、一种OCT系统灵敏度测量方法，其特征在于，包括如下步骤：
扫描光源发出光信号；
将所述光信号分成初始采样光信号和初始参考光信号；
初始参考光信号进入参考臂；
初始采样光信号进入采样臂后发生2N+1次衰减；
从参考臂反射出来的输出参考光信号与从采样臂反射出来的输出采样光信号发生干涉形成干涉光信号；
将所述干涉光信号转换为电信号；
将所述电信号转换为FFT频谱，通过所述FFT频谱得到信噪比，从而计算出OCT系统灵敏度；
所述OCT系统灵敏度＝采样光信号的衰减量+信噪比，所述OCT系统灵敏度为S，所述采样光信号的衰减量为FA，所述信噪比为SNR，即S＝FA+SNR；
所述采样光信号的衰减量＝2N×第一衰减量+第二衰减量；N≥1。
2、根据项1所述的一种OCT系统灵敏度测量方法，其特征在于，所述OCT系统包括扫描光源、分光镜、参考臂、成像系统接口、干涉模块、平衡探测器、FFT变换模块以及显示模块，所述成像系统接口连接采样臂，所述采样臂的一侧设置有采样物。
3、根据项2所述的一种OCT系统灵敏度测量方法，其特征在于，所述扫描光源发射弱相干光信号，所述分光镜将所述弱相干光信号分成初始采样光信号和初始参考光信号，通过所述参考臂以及采样臂反射出来的输出参考光信号与输出采样光信号在干涉模块内发生干涉，形成干涉光信号，所述干涉光信号通过平衡探测器转换为电信号，所述电信号通过FFT变换模块以及显示模块转变为FFT频谱显示在显示模块的屏幕上。
4、根据项3所述的一种OCT系统灵敏度测量方法，其特征在于，所述采样臂包括可变光衰减器和光纤透镜，初始采样光信号进入所述采样臂内，首先依次经过可变光衰减器、光纤透镜照射到采样物上，采样光信号在采样物上发射反射后又依次经过光纤透镜和可变光衰减器输出；
通过所述可变光衰减器来调整初始采样光信号的衰减量，从而改变所述采样臂输出的采样光信号。
5、根据项1-4任一项所述的一种OCT系统灵敏度测量方法，其特征在于，所述采样光信号2N次经过可变光衰减器发生的衰减量为2N×第一衰减量；采样光信号在采样物上发射后的衰减量为第二衰减量；
所述第一衰减量为可变光衰减器的固定衰减量；
所述采样物的反射率为R；
所述第二衰减量为Δx＝10lg(R)。
6、根据项4所述的一种OCT系统灵敏度测量方法，其特征在于，所述可变光衰减器的固定衰减量为0～50dB。
7、根据项4所述的一种OCT系统灵敏度测量方法，其特征在于，所述采样物为平面镜，所述平面镜的反射率为30％～100％，优选为80％～100％；所述采样物的平面和所述光纤透镜的出射光束垂直设置。
8、根据项1所述的一种OCT系统灵敏度测量方法，其特征在于，所述信噪比＝FFT谱峰值高度-噪声平均值-标准偏差；
所述FFT谱峰值高度为FFT频谱峰值达到最大时的峰值高度。
9、根据项2所述的一种OCT系统灵敏度测量方法，其特征在于，所述扫频光源发出的弱相干光信号的波长为820nm～1500nm，带宽为50nm～260nm，扫描频率为5kHz～10MHz。
10、根据项2所述的一种OCT系统灵敏度测量方法，其特征在于，所述参考臂内设置有平面反射镜，所述初始参考光信号进入所述参考臂内，经过所述平面反射镜反射后射出，形成所述输出参考光信号。
11、一种OCT系统灵敏度测量装置，其特征在于，包括OCT系统和测量机构，所述测量机构包括采样臂和采样物，所述采样臂包括可变光衰减器和光纤透镜，且所述可变光衰减器、光纤透镜以及采样物同轴，所述OCT系统与所述可变光衰减器的第一端连接，所述可变光衰减器的第二端与所述光纤透镜连接。
12、根据项11所述的一种OCT系统灵敏度测量装置，其特征在于，所述OCT系统设置有成像系统接口，所述可变光衰减器的第一端设置有光输入接口，所述OCT系统的成像系统接口与所述可变光衰减器的第一端上的光输入接口连接。
13、根据项11所述的一种OCT系统灵敏度测量装置，其特征在于，所述可变光衰减器的第二端与所述光线透镜通过光导纤维连接。
14、根据项11所述的一种OCT系统灵敏度测量装置，其特征在于，所述采样物为平面镜，所述平面镜的镜面和所述光纤透镜的出射光束垂直设置，以使从所述光纤透镜射出来的光线垂直射在所述平面镜上；
所述平面镜的反射率为30％～100％，优选为80％～100％。
15、根据项11所述的一种OCT系统灵敏度测量装置，其特征在于，所述OCT系统包括扫描光源、分光镜、参考臂、干涉模块、平衡探测器、FFT变换模块以及显示模块，所述扫描光源发射弱相干光信号，所述分光镜将所述弱相干光信号分成初始采样光信号和初始参考光信号，通过所述参考臂以及采样臂反射出来的输出参考光信号与输出采样光信号在干涉模块内发生干涉，形成干涉光信号，所述干涉光信号通过平衡探测器转换为电信号，所述电信号通过FFT变换模块以及显示模块转变为FFT频谱显示在显示模块的屏幕上。
16、根据项11所述的一种OCT系统灵敏度测量装置，其特征在于，所述采样臂的长度与所述OCT系统的参考臂的长度一致。
17、根据项11所述的一种OCT系统灵敏度测量装置，其特征在于，所述光纤透镜与所述采样物之间的距离不超过1.5mm。
本申请提供的OCT系统灵敏度测量方法，通过所述扫描光源发射光信号，所述光信号分为初始采样光信号和初始参考光信号，通过输出参考光信号与输出采样光信号两路光的干涉、光电转换、FFT变化、图像采样，获得FFT频谱图。通过移动采样物远离光纤透镜，在FFT频谱图上可获得采样物对应位置的干涉光强。本测量方法可用于不同波长、频率的OCT系统，具有普遍适用性。
本申请提供的OCT系统灵敏度测量方法或装置，可以通过逐渐增加可变光衰减器的固定衰减量，获得对应衰减量的信噪比，最终获得OCT系统极限灵敏度，因此本申请的测量方法更便捷、有效，低成本，同时便于操作。本申请的装置结构简单，能够简便、准确地得到整个OCT系统的灵敏度，既不需要昂贵的光束分析仪和光谱分析仪OSA，又包括了OCT系统本身的光接收机噪声，非常适合在实际生产中使用。
附图说明
附图用于更好地理解本申请，不构成对本申请的不当限定。其中：
图1是根据本申请的OCT系统灵敏度测量装置的结构示意图。
图2是根据本申请的OCT系统灵敏度测量装置的结构示意图。
附图标记列表
1-OCT系统，2-光导纤维，3-可变光衰减器，4-光纤透镜，5-采样物。
具体实施方式
以下对本申请的示范性实施例做出说明，其中包括本申请实施例的各种细节以助于理解，应当将它们认为仅仅是示范性的。因此，本领域普通技术人员应当认识到，可以对这里描述的实施例做出各种改变和修改，而不会背离本申请的范围和精神。同样，为了清楚和简明，以下的描述中省略了对公知功能和结构的描述。
由于扫频光学干涉断层成像系统的灵敏度只取决于扫频光源的光源输出光功率和光路衰减，系统噪声来自光源和光接收，与成像导管本身的机械移动和导管设计没有任何关系，而扫频光学干涉断层成像系统的灵敏度是90～120分贝，为了准确测量，必须排除机械移动和导管设计等所有干扰，因此在测试中对于光路做了调整，用光纤透镜+平面镜来替代成像导管，这个调整完全不会改变系统本身的灵敏度，仅用来优化信号、排除干扰、提高测量精度。
本申请所述的测量方法是从OCT系统里直接得到信噪比SNR，既不需要昂贵的光束分析仪和光谱分析仪OSA，又包括了OCT系统本身的光接收机噪声。OCT系统是一个点对点的光通信传输链路：用光源发出的激光从样品臂照射被检测物体，物体返回的反射/散射光与参考臂干涉，干涉光信号进入平衡探测器进行光电转化，然后对得到的电信号采样，FFT图像处理，在OCT系统显示屏幕上可以直接显示出FFT光强的频谱。这个FFT光强频谱就是系统接收到的光信号的信噪比。不同于来自光束分析仪和光谱分析仪OSA的SNR，这个FFT谱的SNR来自系统本身的光电探测器，因此包含了系统本身的光电探测噪声，并且OCT重建图像就是基于这个FFT信号，因此该信号的信噪比直接决定重建图像的质量。
因为扫频光学干涉断层成像系统灵敏度非常高，如果靠直接在样品臂做光衰减，观察系统FFT信号峰值降至噪声相当来检测，很难得到稳定可靠的测量结果。因此，合理的测量方式应该是光衰减和信噪比相结合：首先，把光路衰减到一定程度，让光接收机得到最好信噪比，然后准确测量系统信噪比(信噪比＝峰值高度-噪声平均值-标准偏差)，衰减量+信噪比＝灵敏度。
本申请提供一种OCT系统灵敏度测量方法，包括如下步骤：
步骤一：扫描光源发出光信号；
步骤二：将所述光信号分成初始采样光信号和初始参考光信号；
步骤三：初始参考光信号进入参考臂；
步骤四：初始采样光信号进入采样臂后发生2N+1次衰减；
步骤五：从参考臂反射出来的输出参考光信号与从采样臂反射出来的输出采样光信号发生干涉形成干涉光信号；
步骤六：将所述干涉光信号转换为电信号；
步骤七：将所述电信号转换为FFT频谱，通过所述FFT频谱得到信噪比，从而计算出OCT系统灵敏度；
所述OCT系统灵敏度＝采样光信号的衰减量+信噪比，所述OCT系统灵敏度为S，所述采样光信号的衰减量为FA，所述信噪比为SNR，即S＝FA+SNR。
所述采样光信号的衰减量＝2N×第一衰减量+第二衰减量，N≥1。
本申请提供的OCT系统灵敏度测量方法，通过所述扫描光源发射光信号，所述光信号分为初始采样光信号和初始参考光信号，所述初始采样光信号在经过采样臂发生衰减后输出的是输出采样光信号，所述初始参考光信号经过参考臂输出的是输出参考光信息，所述输出采样光信号与输出参考光信号发生干涉、光电转换、FFT变化、图像采样，获得FFT频谱。从所述FFT频谱上可以读出信噪比。
在本申请中，所述OCT系统包括扫描光源、分光镜、参考臂、成像系统接口、干涉模块、平衡探测器、FFT变换模块以及显示模块，所述成像系统接口连接有测量机构，所述测量机构包括采样臂和采样物，所述成像系统接口与所述采样臂的第一端连接，所述采样物靠近所述采样臂的第二端设置，且所述采样物与所述采用壁同轴且间隔设置。
所述扫描光源用来提供光信号，所述分光镜用来将所述光信号分为初始参考光信号和初始采样光信号，所述参考臂用来将初始参考光信号通过反射输出，所述采样臂用来将初始采样光信号光量衰减，并反射输出，所述干涉模块用来为所述输出采样光信号和输出参考光信号发生干涉提供条件，所述平衡探测器用来将干涉光信号转化为电信号，FFT变换模块用来将所述电信号转换为FFT频谱，所述显示模块用来显示所述FFT频谱。
所述参考臂与所述采样臂的长度相同，所述输出参考光信号与所述输出采样光信号同时输出发生干涉。
在本申请中，所述扫描光源发射弱相干光信号，所述分光镜将所述弱相干光信号分成初始采样光信号和初始参考光信号，通过所述参考臂以及采样臂反射出来的输出参考光信号与输出采样光信号在干涉模块内发生干涉，形成干涉光信号，所述干涉光信号通过平衡探测器转换为电信号，所述电信号通过FFT变换模块以及显示模块转变为FFT频谱显示在显示模块的屏幕上。
在本申请中，所述采样臂包括可变光衰减器、光纤透镜以及采样物，初始采样光信号进入所述采样臂内，首先依次经过可变光衰减器、光纤透镜照射到采样物上，采样光信号在采样物上发射反射后又依次经过光纤透镜和可变光衰减器输出；
通过所述可变光衰减器来调整初始采样光信号的衰减量，从而改变所述采样臂输出的采样光信号。
所述可变光衰减器是光纤通信中一种重要的光无源器件，通过衰减传输光功率来实现对信号的实时控制。
所述初始采样光信号在所述采样臂内发生2N+1次衰减，其中2N次经过可变光衰减器发生衰减，且2N次的衰减量相同，衰减量由所述可变光衰减器的参数决定。即所述采样光信号两次经过可变光衰减器发生的衰减量为2N×第一衰减量；采样光信号在所述采样物上发射后的衰减量为第二衰减量；
所述第一衰减量为可变光衰减器的固定衰减量；
所述采样物的反射率为R；
所述第二衰减量为Δx＝10lg(R)。
所述采样臂内可以设置一个或多个可变光衰减器，具体根据实际情况来确定。当采样臂内的可变光衰减器为1个时，所述初始采样光信号在所述采样臂内发生3次衰减，所述可变光衰减器的数量为N。
在本申请中，所述可变光衰减器的固定衰减量为0-50dB。
所述可变光衰减器的固定衰减量可以为0、5dB、10dB、15dB、20dB、25dB、30dB、35dB、40dB、45dB、50dB。
在本申请中，所述采样物为平面镜，所述平面镜的反射率为30％～100％，优选80％～100％；所述采样物的平面和所述光纤透镜的出射光束垂直设置，以减少光的损失。
所述平面镜的反射率可以为30％、40％、50％、60％、70％、80％、90％、100％。
在本申请中，所述信噪比＝FFT谱峰值高度-噪声平均值-标准偏差；
所述FFT谱峰值高度为FFT频谱峰值达到最大时的峰值高度。其中FFT谱峰值高度、噪声平均值、标准偏差可以通过OCT系统得到。
在本申请中，所述扫频光源发出的弱相干光信号的波长为820nm～1500nm，带宽为50nm～260nm，扫描频率为5kHz～10MHz。
在本申请中，所述参考臂内设置有平面反射镜，所述初始参考光信号进入所述参考臂内，经过所述平面反射镜反射后射出，形成所述输出参考光信号。
如图1和图2所示，本申请还提供一种OCT系统灵敏度测量装置，包括OCT系统1和测量机构，所述测量机构包括采样臂和采样物5，所述采样臂包括可变光衰减器3和光纤透镜4，且所述可变光衰减器3、光纤透镜4以及采样物5同轴，所述OCT系统1与所述可变光衰减器的第一端连接，所述可变光衰减器的第二端与所述光纤透镜4连接。
初始采样光信号进入所述可变光衰减器3后发生衰减，衰减后的采样光信号穿过光纤透镜4，且在光纤透镜4中不发生衰减，采样光信号射在所述采样物5上，然后采样光信号在采样物5上发生反射，反射后的采样光信号依次射入光纤透镜4和可变光衰减器3，进入OCT系统1中。
所述OCT系统1可以为光学干涉断层成像系统。
前述OCT系统灵敏度测量方法是通过OCT系统灵敏度测量装置来进行测量。
在本申请中，所述OCT系统1设置有成像系统接口，所述可变光衰减器3的第一端设置有光输入接口，所述OCT系统1的成像系统接口与所述可变光衰减器3的第一端上的光输入接口连接。所述成像系统接口与所述光输入接口通过光导纤维2连通。
在本申请中，所述可变光衰减器3的第二端与所述光线透镜通过光导纤维2连接。
在本申请中，所述采样物5为平面镜，所述平面镜的镜面正对所述光纤透镜4设置，以使从所述光纤透镜4射出来的光线垂直射在所述平面镜上。
在本申请中，所述OCT系统1包括扫描光源、分光镜、参考臂、干涉模块、平衡探测器、FFT变换模块以及显示模块，所述扫描光源发射弱相干光信号，所述分光镜将所述弱相干光信号分成初始采样光信号和初始参考光信号，通过所述参考臂以及采样臂反射出来的输出参考光信号与输出采样光信号在干涉模块内发生干涉，形成干涉光信号，所述干涉光信号通过平衡探测器转换为电信号，所述电信号通过FFT变换模块以及显示模块转变为FFT频谱显示在显示模块的屏幕上。
在本申请中，所述采样臂的长度与所述OCT系统1的参考臂的长度一致，以使从参考臂反射出来的光信号和从采样臂输出的光信号能够发生干涉。
在本申请中，所述光纤透镜4与所述采样物5之间的距离不超过1.5mm，以减小采样光信号的衰减量的误差，从而提高检测装置的测量准确的，减小误差。
实施例1
一种OCT系统灵敏度测量方法，包括如下步骤：
步骤一：扫描光源发出光信号；
步骤二：将所述光信号分成初始采样光信号和初始参考光信号；
步骤三：初始参考光信号进入参考臂；
步骤四：初始采样光信号进入采样臂，采样光信号在所述采样臂经过两次固定衰减(固定衰减量为第一衰减量)，两次固定衰减的衰减量为2×35dB，即第一衰减量为35dB，同时经过光纤透镜到达平面镜后反射回光纤透镜，其平面镜发射率为80％；
步骤五：从参考臂反射出来的输出参考光信号与从采样臂反射出来的输出采样光信号发生干涉形成干涉光信号；
步骤六：将所述干涉光信号转换为电信号；
步骤七：将所述电信号转换为FFT频谱，通过所述FFT频谱得到信噪比，从而计算出OCT系统灵敏度。
实施例2和3与实施例1不同之处在于，实施例2以及实施例3中的第一衰减量和平面镜发射率不同。实施例2的第一衰减量为37.5dB，平面镜反射率为90％。实施例3的第一衰减量为40dB，平面镜的反射率为100％。
对比例1
OCT系统灵敏度测量方法，包括如下步骤：
步骤一：扫描光源发出光信号；
步骤二：将所述光信号分成初始采样光信号和初始参考光信号；
步骤三：初始参考光信号进入参考臂；
步骤四：初始采样光信号进入采样臂，采样光信号在所述采样臂经过2*30dB固定衰减，同时经过连接的成像导管到达平面镜后发射回成像导管，其平面镜发射率为0.1％(由于成像导管出光角度和平面镜的相对位置的限制，经反射镜返回的光信号很弱)；
步骤五：从参考臂反射出来的输出参考光信号与从采样臂反射出来的输出采样光信号发生干涉形成干涉光信号；
步骤六：将所述干涉光信号转换为电信号；
步骤七：将所述电信号转换为FFT频谱，观察系统FFT信号峰值和噪声相当，从而计算出OCT系统灵敏度。
对比例2和3与对比例1不同之处在于，对比例2以及对比例3中的第一衰减量和成像导管经平面镜发射率不同。对比例2的第一衰减量为35dB，平面镜反射率为0.2％。对比例3的第一衰减量为40dB，平面镜的反射率为0.4％。
表1为各实施例以及对比例的方法检测的灵敏度以及其他各项参数
| 信噪比/dB | 平面镜反射率 | 第一衰减量/dB | 第二衰减量/dB | 灵敏度/dB | |
| 实施例1 | 48.16 | 80％ | 35 | 0.97 | 119.13 |
| 实施例2 | 43.81 | 90％ | 37.5 | 0.46 | 119.27 |
| 实施例3 | 39.35 | 100％ | 40 | 0 | 119.35 |
| 对比例1 | 0 | 0.1％ | 30 | 30 | 90 |
| 对比例2 | 0 | 0.2％ | 35 | 26.99 | 96.99 |
| 对比例3 | 0 | 0.4％ | 40 | 23.98 | 103.98 |
小结：由上表可见，采用本申请所述的方法，检测到的灵敏度接近OCT系统的极限灵敏度(120dB)，而对比例1-3中由于由于成像导管出光角度和平面镜的相对位置的限制，经反射镜返回的光信号很弱，因此信噪比为0，检测到的灵敏度与极限灵敏度相差较大，而且测量结果波动较大。因此采用本申请所述的方法通过调整采样臂内的衰减量，同时经过平面镜的衰减很小，保证了采样光信号的高质量，而且还避免了成像导管引入的不确定衰减，还可以直接利用OCT系统本身测量系统的噪声和干涉信号获得信噪比，从而能够简便、准确地得到整个OCT系统的灵敏度。
尽管以上结合对本申请的实施方案进行了描述，但本申请并不局限于上述的具体实施方案和应用领域，上述的具体实施方案仅仅是示意性的、指导性的，而不是限制性的。本领域的普通技术人员在本说明书的启示下和在不脱离本申请权利要求所保护的范围的情况下，还可以做出很多种的形式，这些均属于本申请保护之列。
Claims (11)
1.一种OCT系统灵敏度测量方法，其特征在于，包括如下步骤：
扫描光源发出光信号；
将所述光信号分成初始采样光信号和初始参考光信号；
初始参考光信号进入参考臂；
初始采样光信号进入采样臂后发生2N+1次衰减；
从参考臂反射出来的输出参考光信号与从采样臂反射出来的输出采样光信号发生干涉形成干涉光信号；
将所述干涉光信号转换为电信号；
将所述电信号转换为FFT频谱，通过所述FFT频谱得到信噪比，从而计算出OCT系统灵敏度；
所述OCT系统灵敏度=采样光信号的衰减量+信噪比，所述OCT系统灵敏度为S，所述采样光信号的衰减量为FA，所述信噪比为SNR，即S=FA+SNR；
所述采样光信号的衰减量=2N×第一衰减量+第二衰减量；N≥1。
2.根据权利要求1所述的一种OCT系统灵敏度测量方法，其特征在于，所述OCT系统包括扫描光源、分光镜、参考臂、成像系统接口、干涉模块、平衡探测器、FFT变换模块以及显示模块，所述成像系统接口连接采样臂，所述采样臂的一侧设置有采样物。
3.根据权利要求2所述的一种OCT系统灵敏度测量方法，其特征在于，所述扫描光源发射弱相干光信号，所述分光镜将所述弱相干光信号分成初始采样光信号和初始参考光信号，通过所述参考臂以及采样臂反射出来的输出参考光信号与输出采样光信号在干涉模块内发生干涉，形成干涉光信号，所述干涉光信号通过平衡探测器转换为电信号，所述电信号通过FFT变换模块以及显示模块转变为FFT频谱显示在显示模块的屏幕上。
4.根据权利要求3所述的一种OCT系统灵敏度测量方法，其特征在于，所述采样臂包括可变光衰减器和光纤透镜，初始采样光信号进入所述采样臂内，首先依次经过可变光衰减器、光纤透镜照射到采样物上，采样光信号在采样物上发射反射后又依次经过光纤透镜和可变光衰减器输出；
通过所述可变光衰减器来调整初始采样光信号的衰减量，从而改变所述采样臂输出的采样光信号。
5.根据权利要求1-4任一项所述的一种OCT系统灵敏度测量方法，其特征在于，所述采样光信号2N次经过可变光衰减器发生的衰减量为2N×第一衰减量；采样光信号在采样物上发射后的衰减量为第二衰减量；
所述第一衰减量为可变光衰减器的固定衰减量；
所述采样物的反射率为R;
所述第二衰减量为∆x=10lg(R)。
6.根据权利要求4所述的一种OCT系统灵敏度测量方法，其特征在于，所述可变光衰减器的固定衰减量为0~50dB。
7.根据权利要求4所述的一种OCT系统灵敏度测量方法，其特征在于，所述采样物为平面镜，所述平面镜的反射率为30%~100%；所述采样物的平面和所述光纤透镜的出射光束垂直设置。
8.根据权利要求7所述的一种OCT系统灵敏度测量方法，其特征在于，所述平面镜的反射率为80%~100%。
9.根据权利要求1所述的一种OCT系统灵敏度测量方法，其特征在于，所述信噪比=FFT谱峰值高度-噪声平均值-标准偏差；
所述FFT谱峰值高度为FFT频谱峰值达到最大时的峰值高度。
10.根据权利要求2所述的一种OCT系统灵敏度测量方法，其特征在于，所述扫描光源发出的弱相干光信号的波长为820nm~1500nm，带宽为50nm~260nm，扫描频率为5kHz ~10MHz。
11.根据权利要求2所述的一种OCT系统灵敏度测量方法，其特征在于，所述参考臂内设置有平面反射镜，所述初始参考光信号进入所述参考臂内，经过所述平面反射镜反射后射出，形成所述输出参考光信号。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211040958.4A CN115307877A (zh) | 2019-12-25 | 2020-12-25 | 一种oct系统灵敏度测量装置 |
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911356322 | 2019-12-25 | ||
| CN2019113563229 | 2019-12-25 |
Related Child Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211040958.4A Division CN115307877A (zh) | 2019-12-25 | 2020-12-25 | 一种oct系统灵敏度测量装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112763190A CN112763190A (zh) | 2021-05-07 |
| CN112763190B true CN112763190B (zh) | 2022-08-05 |
Family
ID=75694613
Family Applications (3)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211040958.4A Pending CN115307877A (zh) | 2019-12-25 | 2020-12-25 | 一种oct系统灵敏度测量装置 |
| CN202023196483.7U Active CN213957153U (zh) | 2019-12-25 | 2020-12-25 | 一种oct系统灵敏度测量装置 |
| CN202011564335.8A Active CN112763190B (zh) | 2019-12-25 | 2020-12-25 | 一种oct系统灵敏度测量方法以及测量装置 |
Family Applications Before (2)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211040958.4A Pending CN115307877A (zh) | 2019-12-25 | 2020-12-25 | 一种oct系统灵敏度测量装置 |
| CN202023196483.7U Active CN213957153U (zh) | 2019-12-25 | 2020-12-25 | 一种oct系统灵敏度测量装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (3) | CN115307877A (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115307877A (zh) * | 2019-12-25 | 2022-11-08 | 武汉阿格斯科技有限公司 | 一种oct系统灵敏度测量装置 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101617214A (zh) * | 2007-02-21 | 2009-12-30 | 爱克发医疗保健公司 | 用于光学相干断层摄影的系统 |
| CN107389367A (zh) * | 2017-07-24 | 2017-11-24 | 安徽容知日新科技股份有限公司 | 一种基于最优信噪比的信号自适应解包络方法及计算设备 |
| CN109171659A (zh) * | 2018-09-28 | 2019-01-11 | 南京航空航天大学 | 基于琼斯矩阵的光纤型扫频偏振敏感oct成像方法及系统 |
| CN213957153U (zh) * | 2019-12-25 | 2021-08-13 | 武汉阿格斯科技有限公司 | 一种oct系统灵敏度测量装置 |
Family Cites Families (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101793823A (zh) * | 2002-01-24 | 2010-08-04 | 通用医疗公司 | Lci和oct信号的测距和降噪的装置和方法 |
| US7145661B2 (en) * | 2003-12-31 | 2006-12-05 | Carl Zeiss Meditec, Inc. | Efficient optical coherence tomography (OCT) system and method for rapid imaging in three dimensions |
| US8553219B2 (en) * | 2011-01-24 | 2013-10-08 | Vanderbilt University | Common detector for combined raman spectroscopy-optical coherence tomography |
| CN102499648B (zh) * | 2011-11-16 | 2013-11-06 | 清华大学 | 一种基于菲涅尔光谱仪的频域光学相干层析成像系统 |
| CN102835948B (zh) * | 2012-09-12 | 2016-03-30 | 无锡微奥科技有限公司 | 一种扫频光源oct实时图像显示方法及其系统 |
| JP6262762B2 (ja) * | 2012-12-06 | 2018-01-17 | リーハイ・ユニバーシティー | 空間分割多重光コヒーレンストモグラフィー装置 |
| CN104027073B (zh) * | 2014-06-11 | 2016-06-22 | 无锡微奥科技有限公司 | 基于扫频光源的共路光学相干层析成像系统及方法 |
| CN105476605B (zh) * | 2015-12-31 | 2020-04-14 | 东莞理工学院 | 高速光学相干层析成像系统和方法 |
| CN105748041B (zh) * | 2016-02-15 | 2018-07-03 | 苏州大学 | 光学相干断层扫描成像中散斑噪声的抑制系统和方法 |
-
2020
- 2020-12-25 CN CN202211040958.4A patent/CN115307877A/zh active Pending
- 2020-12-25 CN CN202023196483.7U patent/CN213957153U/zh active Active
- 2020-12-25 CN CN202011564335.8A patent/CN112763190B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101617214A (zh) * | 2007-02-21 | 2009-12-30 | 爱克发医疗保健公司 | 用于光学相干断层摄影的系统 |
| CN107389367A (zh) * | 2017-07-24 | 2017-11-24 | 安徽容知日新科技股份有限公司 | 一种基于最优信噪比的信号自适应解包络方法及计算设备 |
| CN109171659A (zh) * | 2018-09-28 | 2019-01-11 | 南京航空航天大学 | 基于琼斯矩阵的光纤型扫频偏振敏感oct成像方法及系统 |
| CN213957153U (zh) * | 2019-12-25 | 2021-08-13 | 武汉阿格斯科技有限公司 | 一种oct系统灵敏度测量装置 |
Non-Patent Citations (3)
| Title |
|---|
| 一种提高扫频光学相干层析成像系统灵敏度的方法;王瑄等;《中国激光》;20170831;第44卷(第8期);第261-267页 * |
| 对超声多普勒胎儿心率仪综合灵敏度检测的探讨;官辉煜;《中国医疗器械信息》;20091231;第15卷(第11期);第63-65页 * |
| 扫频OCT系统关键技术研究;曹蛟;《中国优秀博硕士学位论文全文数据库(硕士)信息科技辑》;20180215(第02期);第I138-1474页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112763190A (zh) | 2021-05-07 |
| CN213957153U (zh) | 2021-08-13 |
| CN115307877A (zh) | 2022-11-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4423450B2 (ja) | 新規の干渉計を用いた光学コヒーレンス断層撮影 | |
| EP1782020B1 (en) | Process, system and software arrangement for determining at least one location in a sample using an optical coherence tomography | |
| CA2421113C (en) | Optical amplification in coherence reflectometry | |
| CA2594311C (en) | Interferometric system for complex image extraction | |
| EP2278265A2 (en) | Common-Path Interferometer for Endoscopic OCT | |
| CN108572161B (zh) | 基于分波阵面干涉仪的光学相干层析成像装置 | |
| JP5120509B1 (ja) | 光プローブおよび光学的測定方法 | |
| CN101214145A (zh) | 大探测深度的频域光学相干层析成像方法及系统 | |
| CN106361279A (zh) | 光学相干层析成像系统全探测深度色散补偿方法 | |
| CN112763190B (zh) | 一种oct系统灵敏度测量方法以及测量装置 | |
| CN106092968B (zh) | 光学检测装置及方法 | |
| CN118356145B (zh) | 一种基于光子集成芯片的时域相干断层扫描系统及方法 | |
| CN106580258B (zh) | 一种宽谱高分辨率光学相干层析成像系统 | |
| CN116849625B (zh) | 一种基于光计算的强度调制式光学相干层析成像系统 | |
| CN111756436B (zh) | 光纤通信中分布式插损回损测量装置及方法 | |
| CN211213098U (zh) | 一种单光谱仪多光束光学相干层析成像仪 | |
| KR20040039139A (ko) | 미첼슨 간섭계와 이를 이용한 광학적 간섭 단층 촬영 시스템 | |
| CN116849626B (zh) | 一种基于光计算的相位调制式光学相干层析成像系统 | |
| CN111637906B (zh) | 基于自差相干的光纤光栅解调装置和方法 | |
| JP5700165B2 (ja) | 光学的測定装置 | |
| Kang et al. | Fourier domain common-path fiber OCT with tunable reference: analysis and optimization | |
| KR101134506B1 (ko) | 고분자 ｗｉｎｃ를 이용한 고해상도ｏｃｔ시스템 | |
| JP2014025953A (ja) | 光学的測定方法 | |
| JP2013165954A (ja) | 光プローブおよび光学的測定方法 |
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 | ||
| REG | Reference to a national code |
Ref country code: HK Ref legal event code: DE Ref document number: 40050041 Country of ref document: HK |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |