一种应用于光学相干断层扫描血管成像方法及OCT系统

摘要

本发明涉及一种应用于光学相干断层扫描血管成像方法及OCT系统，本发明基于分频的思想，将OCT干涉条纹分解为几个波数带，减少被拍摄组织运动产生的噪声。在获取分频的强度图像后，并进一步采用改进CM方法，计算相邻连续扫描的断层图像之间的采用皮尔森相关系数计算强度的变化程度，增强血管的信号。并结合探测点相邻点的信息，提高对血管探测的敏感性，同时减少对眼动的敏感性，适用于对散射较强的生物组织的血管成像。在眼科中可以用于眼前节巩膜或虹膜等组织的血管成像，也可以应用人体其他部分的微血管成像。

说明书

技术领域

本发明涉及一种应用于光学相干断层扫描血管成像方法及OCT系统。

背景技术

微血管形态成像技术在医学成像中具有重要的应用。微血管形态评价对疾病，尤其是血管性疾病的诊断、监测与治疗评价的有重要的价值。离体、侵入性的成像方法可以获得高分辨率三维形态成像，但目前这类成像技术速度慢及对人体具有一定损伤的，限制了其在临床诊断的应用，因此基于光学的成像方法是非侵入性的观察在体血管心态的重要手段，在医学临床中具有重要的应用价值。

目前已有多种基于光学成像技术用于生物组织的血管形态的成像。荧光血管造影技术利用造影剂增强血管与组织的对比度，目前广泛应用于眼科临床，是部分的眼底疾病诊断的金标准。但作为一种具有侵入性的检查方法，部分人可能会对造影剂产生严重的过敏反应，从而限制了造影检查的临床应用。激光多普勒血流仪利用光的多普勒现象，分析运动红细胞与静态组织之间的拍频获取速度，可以获取眼底血流动力学的信息。视网膜功能成像仪利用高速眼底相机在短时间内连续拍摄眼底图像，通过分析图像之间的由于运动红细胞产生的微小差异，增强获取眼底血管的状态。近年来，发生的功能自适应光学成像技术，也可以获取眼底的微血管心态成像。同时这些仪器多受限于二维的平面成像，无法提供三维的深度上的信息。

光学相干断层成像(OCT)技术一种高分辨率、非侵入性、深度分辨的成像技术，其最大的优势在于可以获得轴向高分辨结构图像。已经在许多领域获得广泛应用，尤其是在眼科，已经成为一种不可替代的检查仪器。目前已经发展到第二代的傅立叶OCT(FD-OCT)系统，主要分为谱域OCT(SD-OCT)系统和扫频光源OCT(SS-OCT)系统。由于第二代的FD-OCT成像分辨率与成像速度都得到了极大提高，因此依赖于高速及高分辨率的功能OCT成像技术也得到快速发展。多普勒OCT成像技术作为一种功能OCT成像技术，其利用光的多普勒现象，用于探测血管的流速，但绝对流速的获取依赖于多普勒角度的探测，从而限制了其临床的应用。

与多普勒OCT相比，OCT血管造影技术舍弃了血流流速的信息，增强对微血管形态的成像的灵敏度，目前已经发展有几种OCT血管造影技术。总体上分为与相位相关和基于强度两类血管造影方法。光学微血管造影方法是相位相关的算法，通过使用改良的希尔伯特变化，区分运动与静态的物体，可以实现高灵敏度的微血管成像效果。但基于相位相关的算法依赖于系统的相位的稳定性。与SS-OCT系统相比，SD-OCT系统具有更好的相位稳定性，但是扫描速度相对较慢，随着探测深度的增加系统的信噪比也出现明显的下降，另外对于流速相对较快的血管会出现无法探测到干涉信号的现象。因此单纯基于强度信息的血管造影算法显示出更好的稳定性。近年来美国俄勒冈健康科技大学的大卫黄小组提出频谱分离幅度去相关血管造影(SSADA)在基于强度的血管造影算法的基础上加入了分频方法，减少了深度方向上运动导致的血管造影噪声。在眼科中临床中得到广泛应用。另一种相关成像(CM)的方法，则在相关的成像的基础上采用计算相邻A扫描之间的OCT信号的变化，增强对血管成像的敏感性，可以对微小的血管的进行成像，由于高的敏感性，在拍摄过程组织微小的移动也会带入明显的噪声。在高色散的生物组织中，如眼前的巩膜组织，SSADA和CM方法仍存在一定的问题，SSADA的方法在高色散的组织中对微小的血管成像不敏感，无法很好的成像，很高的敏感性，但是其容易对眼动敏感，而产生明显的噪声。因此OCT血管造影方法在高色散的生物组织中，如眼前节的血管成像中，仍存在一定的困难。

发明内容

本发明针对上述问题提供一种应用于光学相干断层扫描血管成像方法及OCT系统。

本发明所采取的技术方案如下：一种应用于光学相干断层扫描血管成像方法，包括以下步骤：

1)在OCT系统中，采用MB扫描模式采集血管的OCT干涉信号，即在同一个位置B进行N次扫描后，再移动至下一个位置；

2)通过分频将获取得到的干涉信号分解为M个干涉信号；

3)将步骤2)分解后的干涉信号转化为强度图像；

4)获取完整的干涉信号的结构图像，利用得到的在同一位置B扫描的结构图像之间的差异，通过相位相关的方法在获取相邻两幅图片之间的位移(Δx,Δy)，矫正位移；

5)将步骤4)矫正位移后的图像平均，获得血管成像信号；

6)将血管成像信号投影成像。

步骤5)中，图像平均的计算公式如下：

$\left(\begin{array}{c}\bar{D}(x,y)\\ =1\\ -(\frac{1}{N-1}\underset{n=1}{\overset{N-1}{\Sigma }}\underset{m=1}{\overset{M}{\Sigma }}\frac{{\Sigma }_{p=0}^P{\Sigma }_{q=0}^Q[I_{n,m}(x+p,y+q)-\overline{I_{n,m}(x,y)}][I_{n+1,m}(x+p,y+q)-\overline{I_{n+1,m}(x,y)}]}{\sqrt{{\Sigma }_{p=0}^P{\Sigma }_{q=0}^Q{\left(I_{n,m}\right(x+p,y+q)-\overline{I_{n,m}(x,y)})}^2}\sqrt{{\Sigma }_{p=0}^P{\Sigma }_{q=0}^Q{\left(I_{n+1,m}\right(x+p,y+q)-\overline{I_{n+1,m}(x,y)})}^2}}\\ +1)/2\end{array}\right)$

其中N为在同一个位置B扫描的数目，M为分频的数目，则为获取的血管成像信号。

步骤5)中，对于得到的血管成像信号，统计信号的平均值和标准差s，进行阈值处理，对于低于血管成像信号值将其置零，同时调整对比度。

步骤6)中，投影成像前，探测表面边界，并对边界的以下的部位进行部分投影增强血管成像的对比度。

一种应用于上述的应用于光学相干断层扫描血管成像方法的OCT系统，所述OCT系统为扫频光源OCT系统，所述扫频光源OCT系统包括扫频光源、平衡探测器、第一光纤耦合器、第二光纤耦合器、计算机、色散补偿、反射镜、振镜、参考臂、样品臂，所述扫频光源连接第一光纤耦合器的输入端，所述第一光纤耦合器的输出端连接参考臂和样品臂，参考臂和样品臂连接第二光纤耦合器，所述第二光纤耦合器连接平衡探测器，所述扫频光源和平衡探测器连接计算机，所述参考臂设有色散补偿、反射镜，所述样品臂设有振镜。

所述参考臂、样品臂均设有偏振控制器调整偏振。

所述第一光纤耦合器为80:20，其中20％的光进入样品臂，80％的光进入参考臂。

本发明的有益效果如下：本发明基于分频的思想，将OCT干涉条纹分解为几个波数带，减少被拍摄组织运动产生的噪声。在获取分频的强度图像后，并进一步采用改进CM方法，计算相邻连续扫描的断层图像之间的采用皮尔森相关系数计算强度的变化程度，增强血管的信号。并结合探测点相邻点的信息，提高对血管探测的敏感性，同时减少对眼动的敏感性，适用于对散射较强的生物组织的血管成像。在眼科中可以用于眼前节巩膜或虹膜等组织的血管成像，也可以应用人体其他部分的微血管成像。

附图说明

图1用于血管成像的扫频光源OCT系统。

图2血管成像算法流程图。

图3基于强度的边界探测，其中(a)为一个断面的图像，(b)为对(a)中蓝色竖直线的强度分析。

图4巩膜及角膜缘处血管成像结果，其中(a)为对强度图像直接做投影的结果，(b)则为利用本文提出的方法获取的血管成像图像，(c)和(d)分别为(b)中红线和蓝线对应的横断面B扫描图像。

图中，1，扫频光源；2，平衡探测器；3，第一光纤耦合器；4.第二光纤耦合器；5，计算机；6，偏振控制器；7，色散补偿；8，反射镜；9，振镜；10，检测部位；11，参考臂；12，样品臂。

具体实施方式

下面结合附图以及具体实施方式，可更好地说明本发明。

如图1所示，扫频光源1发出的光经过80:20的第一光纤耦合器3，其中20％的光进入样品臂12，80％的光进入参考臂11。在样品臂12和参考臂11均有偏振控制器6调整偏振。进入样品臂12的光，通过准直镜和聚焦透镜照射于人眼前节部位，其中通过振镜9的摆动实现三维的数据采集获取。80％的光进入参考臂，并经过色散补偿7和反射镜8返回。返回的参考臂和样品臂的光进入50:50的第二光纤耦合器，并在平衡探测器2处获取干涉的信号。系统的时钟信号(Clk)和触发信号(Tri in)都由光源产生，产生的干涉信号由通过通道A(Ch A)由采集卡获取。扫描振镜触发信号则通过采集卡产生的触发信号(Tri out)同步，最终将光源，采集卡和扫描振镜同步。本文的采用MB扫描模式，即在同一个B扫描的位置进行多次(N次)扫描后，再移动至下一个位置。

如图2所示，采集到OCT干涉信号后，对获取得到的干涉信号已经是按照波数采集的信号，通过分频将获取得到的干涉信号分解为M个干涉信号，即对获取得到的完整的干涉信号做高斯滤波，不同的高斯函数的宽度相同，中心位置在采集得到的干涉信号上等间隔分布。从干涉信号到强度图像的主要的处理步骤如下：去除固定噪声、傅里叶变换、逆傅里叶变换及数值色散补偿及再次的傅立叶变换。最终通过舍弃获得的复数信号的相位部分，得到强度图像(I_n,m)。

在获取不同干涉信号上强度图像的同时，获取完整的干涉信号的结构图像，对同一个位置的B扫描通过相位相关的方法在获取相邻两幅图片之间的位移(Δx,Δy)，通过矫正位移，可以消除眼动的影响。对于相邻位移较大的图像，及位移量大于一定的阈值，则获取得到强度图像后，对于已经做位移矫正的相邻的两幅幅图像(I_n,m和I_n+1,m)，进行皮尔森相关运算，具体如下式所示：

$\left(\begin{array}{c}{C}_n(x,y)\\ =\frac{{\Sigma }_{p=0}^P{\Sigma }_{q=0}^Q[I_{n,m}(x+p,y+q)-\overline{I_{n,m}(x,y)}][I_{n+1,m}(x+p,y+q)-\overline{I_{n+1,m}(x,y)}]}{\sqrt{{\Sigma }_{p=0}^P{\Sigma }_{q=0}^Q{\left(I_{n,m}\right(x+p,y+q)-\overline{I_{n,m}(x,y)})}^2}\sqrt{{\Sigma }_{p=0}^P{\Sigma }_{q=0}^Q{\left(I_{n+1,m}\right(x+p,y+q)-\overline{I_{n+1,m}(x,y)})}^2}}\end{array}\right)$

其中P和Q为选取相邻位置区块的大小，则对原图像做模板大小为P×Q均值滤波后的结果，通过与原位置强度相减可以对血管进行增强，在本文中采用的P和Q均为3，最后获取的相邻两幅图像的相关值C_n(x,y)，其范围在-1到1之间。在相关计算的基础上，通过分频及多次扫描的方法进一步增强血管，公式可变为下式所示：

$\left(\begin{array}{c}\bar{D}(x,y)\\ =1\\ -(\frac{1}{N-1}\underset{n=1}{\overset{N-1}{\Sigma }}\underset{m=1}{\overset{M}{\Sigma }}\frac{{\Sigma }_{p=0}^P{\Sigma }_{q=0}^Q[I_{n,m}(x+p,y+q)-\overline{I_{n,m}(x,y)}][I_{n+1,m}(x+p,y+q)-\overline{I_{n+1,m}(x,y)}]}{\sqrt{{\Sigma }_{p=0}^P{\Sigma }_{q=0}^Q{\left(I_{n,m}\right(x+p,y+q)-\overline{I_{n,m}(x,y)})}^2}\sqrt{{\Sigma }_{p=0}^P{\Sigma }_{q=0}^Q{\left(I_{n+1,m}\right(x+p,y+q)-\overline{I_{n+1,m}(x,y)})}^2}}\\ +1)/2\end{array}\right)$

其中N为在同一个位置B扫描的数目，M为分频的数目，则为获取的血管成像信号。对每一幅B扫描的图像进行处理，最终可以得到三维的血管成像信号。

为增强血管的成像效果，对于三维血管成像信号，统计信号的平均值和标准差s，进行阈值处理，认为低于血管成像信号值不是血管信号，将其置零，同时调整对比度。

在对血管的图像进行投影前，根据眼前节角巩膜部位的组织结构特点，在结膜表面没有血管，但表面的泪液产生明显的强反射会引起较大的噪声。本文提出通过探测表面边界，并对边界的以下的部位进行部分投影增强血管成像的对比度。如上图3所示，通过对每一条A扫描(图3a中蓝色竖直线)强度分析，分别设置强度曲线(图3b中蓝色线)和一阶导数曲线(图3b中红色线)阈值，判断角巩膜部位第一条边界位置，并向下取150个像素点做投影的图像。图3a中的红线和绿线显示了探测上界和下界，两条线之间的区域为投影区域。

图4显示角结膜缘处血管成像图像，其中(a)为对强度图像直接做投影的结果，无法观察到清晰的血管心态，(b)则为利用本文提出的方法获取的血管成像图像，可以看到本文血管成像的方法可以清晰显示角结膜缘处的微血管网络。其中(c)和(d)分别为(b)中红线和蓝线对应的横断面B扫描图像，可以看大在箭头所指的部位，本文提出的血管成像的方法可以对巩膜内深部的血管进行清晰的成像。

其中使用扫频光源OCT系统用于人眼前节的微血管成像，但本发明所述的血管成像方法可以应用于其他类型的OCT系统，如谱域OCT系统，拍摄部位可以是生物组织的其他部位，如眼底视网膜，表面皮肤等。

以上所述仅为本发明的一种实施例，并非用来限制本发明的保护范围；本发明的保护范围由权利要求书中的权利要求限定，并且凡是依发明所作的等效变化与修改，都在本发明专利的保护范围之内。