CN105342568B

CN105342568B - 联合相位和幅值的光学相干造影方法及系统

Info

Publication number: CN105342568B
Application number: CN201510843502.5A
Authority: CN
Inventors: 周传清; 秦林
Original assignee: Shanghai Jiaotong University
Current assignee: Shanghai Jiaotong University
Priority date: 2015-11-26
Filing date: 2015-11-26
Publication date: 2018-02-02
Anticipated expiration: 2035-11-26
Also published as: CN105342568A

Abstract

一种OCT技术领域的联合相位和幅值的光学相干造影方法和系统，包括以下步骤：扫描样品，采集干涉信号得到样品的体数据；对样品体数据进行快速傅里叶变换重建信号得到样品组织信息；提取样品组织的相位信息，结合信号幅值对组织运动引起的伪信号进行相位补偿，得到多普勒血流速度断面图像和多普勒血流图像；对多普勒血流断面图像进行直方图分析，估计样品中静止组织的相位分布；根据相位分布，在快速扫描方向上进行高通滤波，分离出血流与静止组织信号；将分离出的血流信号与多普勒血流速度断面图像结合并进行处理得到血流成像。本发明可对生物组织血流进行高分辨率成像，便于血流分布、血管形态以及血液动力学变化的检测。

Description

联合相位和幅值的光学相干造影方法及系统

技术领域

本发明涉及一种OCT(Optical Coherence Tomography，光学相干层析成像)领域的技术，具体是一种联合相位和幅值的光学相干造影方法和系统。

背景技术

生物组织光学不均匀性导致背向散射光光强变化，OCT技术通过探测光强变化来获得生物组织的结构图像，实现非侵入、高分辨率的三维生物组织成像。在临床上，尤其是在眼科的疾病诊断和治疗中具有广泛的应用。大多数情况下，尤其是在疾病的早期阶段，传统的血管造影技术不仅需要注射造影剂，可能会引起不良反应，而且无法得到更多的血流功能信息，而利用OCT技术则可以非侵入、无标记地发掘更多的生理信息。

目前，OCT技术大致可以分为两大类，一类基于散斑效应，另一类基于多普勒效应。基于散斑效应的OCT技术利用血流流动引起反射信号散斑图样的变化进行成像，对于血流而言，反射信号强度随时间发生变化，而组织信号由于保持静止，故信号变化较小，从而能够分离得到血管的形态结构图像。基于散斑效应的造影方法具有较高的灵敏度和对比度，且不受到相位不稳定的影响，但需要对同一样品位置进行多次扫描采集数据，对于临床应用的OCT系统而言，扫描速度较慢，不适于进行多次重复扫描，且不宜进行过采样采集数据。而基于多普勒效应的OCT技术，由于组织中的运动物体，如红细胞等，使得反射回来的信号存在多普勒频移，而且多普勒频移与血流速度有关，其中组织与血流分别对应信号频谱的低频和高频部分，因此利用信号处理的方法提取血流的多普勒频移进行血流与组织信号的分离，能够提供血流的速度方向等信息，但是由于眼底血管分布复杂，赖以获取血流速度的多普勒角难以准确确定，因此多普勒OCT获取的血流速度不准确，临床意义有限。

经过对现有技术的检索发现，美国专利文献号US8718743B2，公开(公告)日2014.05.06，公开了一种基于改进的希尔伯特变换算法的单通道容积双向血流的谱域OCT造影技术，进行过采样采集数据，并设定阈值分离血流和组织。但是由于血流不断变化，且相邻的血流之间也存在差异，利用固定的阈值分离血流和组织的准确率不高，误差较大。

中国专利文献号CN102835947A，公开(公告)日2012.12.26，公开了一种基于散斑相关度的OCT图像分析方法，首先，获取同一成像对象在不同波长段的相对应的若干幅OCT图像；或者采用同一光源获得同一成像对象在不同角度的相对应的若干幅OCT图像；对上述不同的OCT图像，在每幅OCT图像中的同一位置分别取相同尺寸的处理窗口；然后，计算处理窗口内像素强度并通过相关函数的计算进行相关性分析；同步移动处理窗口，直至遍历整个OCT的图像；根据若干幅OCT图像中各个部分相关度的分析，重建散斑分析后的图像，获取成像对象的微结构信息。但是该发明依赖于信号的散斑效应，需要对眼动及其他非血流运动进行精确的消除；同时无法获取血流速度，且并不适用于低速扫描的OCT。

发明内容

本发明针对现有技术存在的上述不足，提出了一种联合相位和幅值的光学相干造影方法及系统，能够自适应的根据当前帧的血流分布，分离出血流和组织，并适用于扫描速度较低的OCT系统，实现了无标记、高分辨的血流成像。

本发明是通过以下技术方案实现的，

本发明涉及一种联合相位和幅值的光学相干造影方法，包括以下步骤：

步骤S1，光谱OCT系统沿着x-y轴扫描样品，采集干涉信号得到样品干涉信号的体数据I(x,y,λ)；

步骤S2，对光谱OCT系统采集的干涉信号在λ方向上重建信号，得到样品组织信息A(x,y,z)；

步骤S3，对样品组织信息A(x,y,z)进行相位提取得到相位结合信号幅值|A(x,y,z)|对组织运动引起的伪信号进行相位补偿，得到多普勒血流速度断面图像

步骤S4，对多普勒血流速度断面图像进行直方图分析，估计样品中静止组织的相位分布F_cut-off；

步骤S5，根据估计的静止组织相位分布F_cut-off，在x轴方向上进行高通滤波A_blood(x,z)＝A(x,z)*h_high-pass(F_cut-off)，分离得到血流信号A_blood(x,z)；

步骤S6，将分离出的血流A_blood(x,z)和多普勒血流速度断面图像结合，得到每个断面的血流成像结合y方向上每个断面的血流成像I_blood(x,z)得到血流三维体数据I_blood(x,y,z)，并对I_blood(x,y,z)在深度z方向进行投影，得到对样品无标记的血流成像。

所述的相位补偿是指：完成，其中：表示第i条A-line的相位，I_i(z)表示第i条A-line的信号幅值，表示补偿后的相位。

本发明涉及一种实现上述方法的光学造影系统，包括：低相干光源、分束器、参考臂、探测臂和光信号分析支路，其中：探测臂包括：依次设置的透镜、X-Y扫描振镜、透镜和物镜；低相干光源产生的低相干光由分束器分束，分别入射到参考臂和探测臂，参考臂接受低相干光并反射至分束器，进入探测臂的低相干光经过透镜、X-Y扫描振镜、透镜和物镜入射到样品表面，样品接收低相干光并发出背向散射光，背向散射光依次经过物镜、透镜、X-Y扫描振镜和透镜至分束器；参考臂反射光与背向散射光在分束器中相互干涉，分束器将干涉信号传输至光信号分析支路进行分析，得到样品的血流成像。

技术效果

与现有技术相比，本发明基于生物组织血流成像的频谱分析，对生物组织血流进行高分辨率成像。实现生物组织血流分布、血管形态以及血液动力学变化的检测，可广泛用于生理和疾病状态下的眼底和脑皮层局部血流成像。

附图说明

图1为本发明的方法流程图；

图2为本发明的系统结构示意图；

图3为现有OCT技术得到的人眼眼底结构投影图；

图4为本发明得到的人眼视网膜眼底血管投影图；

图中：低相干光源1、分束器2、参考臂3、准直透镜3a、透镜3b、平面镜3c、探测臂4、透镜4a、X-Y扫描振镜4b、物镜4c、光信号分析支路5、准直透镜5a、光栅5b、透镜5c，电耦合装置5d、处理系统5e、样品6。

具体实施方式

下面对本发明的实施例作详细说明，本实施例在以本发明技术方案为前提下进行实施，给出了详细的实施方式和具体的操作过程，但本发明的保护范围不限于下述的实施例。

实施例1

如图1所示，本实施例涉及利用上述系统进行造影的方法，包括以下步骤：

步骤S1，光谱OCT系统沿着x-y轴扫描样品6，采集干涉信号得到干涉信号的体数据I(x,y,λ)；所述的体数据I(x,y,λ)是指在y方向上的每一个断层图像I(x,λ)排列起来得到的三维数据；

步骤S2，对光谱OCT系统采集的干涉信号在λ方向上进行快速傅里叶变换重建信号：FT_λ→z[I(x,y,λ)]＝A(x,y,z)，得到样品6的组织信息A(x,y,z)；

步骤S3，重建后样品6的组织信息可表示为对样品6的组织信息进行相位提取得到相位结合信号幅值|A(x,y,z)|对组织运动引起的伪信号进行相位补偿，得到多普勒血流速度断面图像

步骤S4，通过对多普勒血流速度断面图像进行直方图分析：估计样品6中静止组织的相位分布F_cut-off；

步骤S6，将分离出的血流A_blood(x,z)和多普勒血流速度断面图像结合，得到每个断面的血流成像结合y方向上每个断面的血流成像I_blood(x,z)得到血流三维体数据I_blood(x,y,z)，对I_blood(x,y,z)在深度z方向进行投影，得到对样品6无标记的血流成像。

所述的相位补偿通过公式：完成，其中：表示第i条A-line的相位，I_i(z)表示第i条A-line的信号幅值，表示补偿后的相位。

如图2所示，本实施例涉及实现上述方法的光学造影系统，包括：低相干光源1、分束器2、参考臂3、探测臂4和光信号分析支路5，其中：探测臂4包括：依次设置的透镜4a、X-Y扫描振镜4b、透镜4a和物镜4c；低相干光源1产生的低相干光由分束器2分束，分别入射到参考臂3和探测臂4，参考臂3接受低相干光并反射至分束器2，进入探测臂4的低相干光经过透镜4a、X-Y扫描振镜4b、透镜4a和物镜4c入射到样品6的表面，样品6接收低相干光并发出背向散射光，背向散射光依次经过物镜4c、透镜4a、X-Y扫描振镜4b和透镜4a至分束器2；参考臂反射光与背向散射光在分束器2中相互干涉，分束器2将干涉信号传输至光信号分析支路5进行分析，得到样品6的血流成像。

所述的光信号分析支路5包括：依次设置的准直透镜5a、光栅5b、透镜5c、电耦合装置5d和处理系统5e。

如图3和图4所示，即为现有的图像与经本发明处理后所得图像的比较。根据图4可见，采用本方法可以明显提高图像清晰度。

Claims

1.一种联合相位和幅值的光学相干造影方法，包括以下步骤：

步骤S1，光谱OCT系统沿着x-y轴扫描样品，采集干涉信号得到干涉信号的体数据I(x,y,λ)；

其特征在于，所述方法还包括：

步骤S3，对重建后的样品组织信息A(x,y,z)进行相位提取得到相位结合信号幅值|A(x,y,z)|对组织运动引起的伪信号进行相位补偿，得到多普勒血流速度断面图像

步骤S6，将分离出的血流信号A_blood(x,z)与多普勒血流速度断面图像结合，得到每个断面的血流成像结合y方向上每个断面的血流成像I_blood(x,z)得到血流三维体数据I_blood(x,y,z)，并对I_blood(x,y,z)在深度z方向进行投影，得到对样品无标记的血流成像。

2.根据权利要求1所述的联合相位和幅值的光学相干造影方法，其特征是，所述的相位补偿通过公式：完成，其中：表示第i条A-line的相位，I_i(z)表示第i条A-line的信号幅值，表示补偿后的相位。

3.根据权利要求1所述的联合相位和幅值的光学相干造影方法，其特征是，所述的重建信号通过快速傅里叶变换FT_λ→z[I(x,y,λ)]＝A(x,y,z)得到。

4.根据权利要求1所述的联合相位和幅值的光学相干造影方法，其特征是，所述的重建后样品的组织信息表示为

5.根据权利要求1所述的联合相位和幅值的光学相干造影方法，其特征是，所述的直方图分析表示为

6.根据权利要求1所述的联合相位和幅值的光学相干造影方法，其特征是，所述的体数据I(x,y,λ)为处理轴向扫描信号得到的横向扫描断层图像的集合。

7.根据权利要求1所述的联合相位和幅值的光学相干造影方法，其特征是，所述的血流成像包括：血管结构图像与血流分布图像。

8.一种实现上述任一权利要求所述方法的光学造影系统，包括：低相干光源、分束器、参考臂、探测臂和光信号分析支路，其中：探测臂包括：依次设置的透镜、X-Y扫描振镜、透镜和物镜；低相干光源产生的低相干光由分束器分束，分别入射到参考臂和探测臂，参考臂接受低相干光并反射至分束器，进入探测臂的低相干光经过透镜、X-Y扫描振镜、透镜和物镜入射到样品表面，样品接收低相干光并发出背向散射光，背向散射光依次经过物镜、透镜、X-Y扫描振镜和透镜至分束器；参考臂反射光与背向散射光在分束器中相互干涉，分束器将干涉信号传输至光信号分析支路进行分析，得到样品的血流成像。