[发明专利]一种漫射光相关谱血流测量的实时计算方法

说明书

本发明提供了一种漫射光相关谱血流测量的实时计算方法，包括S1、在归一化光场时间自相关函数g1(τ)曲线上选取n个点，产生新变量[x(i),y(i)],1≤i≤n；S2、执行关于[x(i),y(i)]的四则运算；S3、由步骤S2的结果获得一阶血流值Db1；S4、计算中间变量M1，更新变量y(i)的值并执行步骤S2的四则运算，获得二阶血流值Db2；S5、重复步骤S4，分别计算出三阶血流值Db3、四阶血流值Db4和五阶血流值Db5。本发明所述的漫射光相关谱血流测量的实时计算方法，克服了现有单纯性方法曲线拟合求解血流存在的不易计算、耗时长不足和稳定性差的问题。

技术领域

本发明属于生物医学工程中的组织血流测量技术领域，尤其是涉及从近红外漫射光相关谱数据中计算组织血流的一种方法。

背景技术

脑部血流的实时监测对于手术和重症病人的重要性不言而喻。例如，颈动脉翻瓣手术过程中需要实时监测由于血管钳夹引起的局部脑灌注不足。另外，其他脑部障碍(例如睡眠呼吸暂停或迷走神经性晕厥)也导致血流震荡或者不稳定。因此，脑部血流的测量对于多种疾病的诊断和治疗非常重要，尤其是脑部局部缺血、中风、睡眠中呼吸暂停、迷走神经性晕厥等重要疾病。除了脑部血流监测，骨骼肌血流也是评估肌肉功能和外周血管疾病、老年性肌纤维疼痛等病症的重要生理参数之一。

目前医院里使用的脑血流设备主要是超声多普勒(又称为经颅多普勒)和灌注磁共振(ASL-MRI)技术,但是，这两种技术都存在局限性。超声多普勒测量的是脑主干血管，但是很多脑疾病与局部微血管循环有关；ASL-MRI可以得到很高空间分辨率,但是测量成本昂贵，也不容易实现动态和连续测量。近红外漫射光谱(NIRS)是检测局部组织血氧的传统技术，近年来，一种新的NIRS技术得到了发展，称为漫射光相关谱技术(DCS)；DCS是利用时间自相关方法来快速估计红细胞的运动引起的光场扰动，是一种直接测量组织血流的新技术。该技术具有无创、低成本、便携等优点，近些年来发展很快。

在DCS数据采集过程中，将光源和探测器放置于被测组织(例如头部前额)的同侧，相距几毫米至几厘米，光源发射的近红外光子进入脑皮层后被运动的红细胞多次散射，造成光场扰动；少部分光子从表层逸出并被探测器收集；数字相关器将探测器采集的光子数(即光强)进行时间自相关计算，获得归一化光强自相关函数g2(τ)，并由西格特关系获得归一化光场时间自相关函数g1(τ)，这里τ为延迟时间。由于非归一化光场时间自相关函数G1(τ)满足扩散偏微分方程，因此目前的求解方法是将微分方程的解析表达式与g1(τ)曲线做拟合；曲线拟合的基本过程是：先将一个假定的血流值代入解析表达式，得到不同延迟时间τ的时间自相关函数g1(τ)的理论值，再将所有延迟时间的g1(τ)理论值与测量值求残差平方和；然后将假定的血流值迭代变化以寻求最小的残差平方和；这个过程主要用的是非线性的单纯形方法寻求最优解。但是，这种非线性拟合的方法不是直接的数据四则运算(即数据的加减乘除)，迭代变化的耗时长，不易实现实时的血流计算。另外，此方法容易受到噪声点的影响使得求解不稳定。

发明内容

有鉴于此，本发明旨在提出一种漫射光相关谱血流测量的实时计算方法，以克服现有单纯性方法曲线拟合求解血流存在的不易计算、耗时长不足和稳定性的问题。

为达到上述目的，本发明的技术方案是这样实现的：

一种漫射光相关谱血流测量的实时计算方法，包括以下步骤：

S1、在归一化光场时间自相关函数g1(τ)曲线上选取n个点，产生新变量[x(i),y(i)],1≤i≤n；

S2、执行关于[x(i),y(i)]的运算；

S3、由步骤S2的结果获得一阶血流值Db1；

S4、计算中间变量M1，更新变量y(i)的值并执行步骤S2的四则运算，获得二阶血流值Db2；

S5、计算中间变量M2，更新变量y(i)的值并执行步骤S2的四则运算，获得三阶血流值Db3；