[发明专利]基于横向流速模型的血流流速成像方法

说明书

技术领域

本发明涉及血流成像技术领域，特别涉及一种基于横向流速模型的血流流速成像方法。

背景技术

体表浅层微循环系统血流流速分布的快速成像，在临床上有广泛的应用，目前，非侵入性的血流成像技术包括：散斑技术、超声多普勒流速仪、激光多普勒流速仪、多普勒光学相干层析、基于高速CMOS相机的激光多普勒血流成像技术。散斑技术不能进行定量测量，也不适用于较快的流速，超声多普勒流速仪只能用于大血管，无法用于微循环系统。多普勒光学相干层析和激光多普勒流速仪都是单点测量，需要机械扫描实现较大面积成像，成像速度慢。基于高速CMOS相机的激光多普勒血流成像技术是对激光多普勒流速仪的改进，和激光多普勒流速仪一样都是基于多普勒效应，使用帧频率为20KHZ左右的高速CMOS相机，代替机械扫描，实现快速成像。但是这种高速CMOS相机及所需要的巨大的数据传输处理速度，使这种成像技术的成本很高，同时血流流速结果受多普勒角不确定性影响较大。

发明内容

本发明目的在于提供一种基于横向流速模型的血流流速成像方法，它能通过红细胞断续进出探测光焦斑导致的背向散射光强信号的低频涨落去计算血流流速，以代替目前使用的由多普勒效应导致的高频涨落计算血流流速的方法。

本发明所提出的技术解决方案是这样的：

一种基于横向流速模型的血流流速成像方法，该成像方法包括如下步骤：用经扩束的激光5照射体表6，用CMOS相机7对被照区域连续成像，CMOS相机7每个像素记录的光强信号形成独立的信号序列，计算每个信号序列的归一化自相关函数，归一化自相关函数的第一个零点之前部分线性区域的斜率为红细胞1横向穿过探测光焦斑2的渡越时间                                               ，红细胞的横向速度为，其中，为探测光焦斑的横向宽度。

本发明的原理如下：经扩束的激光5以和体表6接近垂直的方向照射体表6，用CMOS相机7对被照区域进行连续成像，CMOS相机7的整个像素矩阵等价于并行的探测器阵列，每一个像素即是一个独立的探测器，记录对应焦点处随时间变化的背向散射光强信号。假定某个像素对应的焦斑处的血流情况如图1（a）所示，图中黑点表示流动的红细胞，点线矩形框表示探测光的焦斑2，当红细胞没有进入探测光焦斑2时，反射回去的信号很弱，当红细胞1进入探测光焦斑2时，会有较强的反射信号，因此受流动的红细胞1断续进出探测光焦斑2的调制，背向散射光强信号随时间形成一个一个的脉冲，当红细胞1移动速度较慢时，穿过探测光焦斑2的时间较长，产生的脉冲宽度较大，当红细胞1移动速度较快时，穿过探测光焦斑2的时间较短，产生的脉冲宽度较窄，脉冲的时间宽度对应红细胞1穿过探测光焦斑2的时间，即渡越时间。图1（a）中，在非聚焦区域穿过探测光束的红细胞3、4，由于这些红细胞3、4处于焦斑2之外，对背向散射光强信号的影响较小，因此，焦斑2外的红细胞的影响可以忽略。图1（b）中矩形框表示放大的焦斑2，红细胞以速度流过焦斑2，探测光沿竖直方向，多普勒角为，的纵向分量导致多普勒频移，假定为；红细胞以横向速度穿过探测光焦斑2，则背向散射光可表示为，假定血管周围静止组织的背向散射光为，这两部分光都进入探测器进行干涉（即外差法），则探测到的光为，因为远小于，则探测到的光强信号可以简化为，如图2（a）所示，其中的高频振荡来源于多普勒效应，即，所以，基于多普勒效应的血流成像方法都是分析高频振动信号计算流速，对相机的帧频要求比较高，一般在20KHZ。图2（a）中的低频包络（如图2（b）所示）对应于，图2（b）中的脉冲就是红细胞穿过探测光焦斑所形成，因此脉冲的平均时间宽度就是红细胞穿过探测光焦斑的时间，即渡越时间，假定探测光焦斑的横向宽度为，则红细胞的横向速度。由于图2（b）中的低频脉冲的频率较低，因此对相机的帧频要求不低于500HZ即可。当相机的曝光时间大于多普勒效应导致的高频振荡的几个周期时，则多普勒效应导致的高频振荡就被平均，信号中只有红细胞断续进出探测光焦斑导致的低频脉冲，如图2（b）所示。

虽然横向流速模型只能得到流速的横向分量，即，为实际流速，为多普勒角，由于体表浅层血管大部分都是和体表接近平行分布，当成像系统和照射激光和体表接近垂直时，多普勒角的范围在90°～80°之间，则流速横向分量的范围在之间，非常接近实际流速，所以用横向流速可以代表实际流速，不用考虑未知的多普勒角的影响，而目前的基于多普勒效应的血流成像方法得到的是流速的纵向分量，即，当多普勒角的范围在度之间时，结果为，和实际流速相差较大，不能反映真实流速。