一种基于电阻抗断层成像识别无效吸气努力的方法

1.本发明涉及一种基于电阻抗断层成像的机械通气无效吸气努力检测方法，属于医学信号处理领域。


背景技术：

2.在使用呼吸机对重症病人实施有创机械通气期间，常会发生人机不同步，当呼吸神经的吸气时间与机械通气时间不匹配时容易发生人机不同步。人机不同步可能造成机械通气时间延长，患者预后变差等结果，因此准确地识别出人机不同步，进而修正人机不同步具有重要的医学研究价值。
3.无效吸气努力是人机不同步中最常见的类型之一。无效吸气努力原理是由于在呼气阶段，患者本能的呼吸肌肉收缩想要吸气，但吸气的幅度比较弱，气道压力未达到呼吸机的送气阈值。所以检测无效吸气努力只需要判断呼气阶段患者是否存在微弱的吸气努力。
4.目前，尽管医学领域对无效吸气努力的机制和识别进行了很多研究，但很少注意不同步发生时肺部局部潮气量分布的变化。医生通常通过呼吸机波形来间接判断是否发生无效吸气努力，但有的时候，患者有吸气努力却在呼吸机波形上没有明显变化，导致部分无效吸气努力被遗漏。部分医生通过食道压或膈肌电监测等方式来识别吸气努力，尽管这些方法被认为是金标准，但方法需要操作过于麻烦，不适用于临床日常工作。


技术实现要素：

5.为了克服现有技术的不足，本发明提出一种基于电阻抗断层成像(electrical impedance tomography,eit)的方法来判断患者是否存在吸气努力的动作，eit是监测局部肺通气的有效工具，基于eit可以根据局部潮气量变化来分析无效吸气努力的发生。本发明通过eit计算得到局部密度分数(regional intensity fraction,rif)曲线，进而根据rif来鉴别无效吸气努力，为识别无效吸气努力提供一种更精准合理的途径。
6.本发明解决其技术问题所采用的技术方案是：
7.一种基于电阻抗断层成像识别无效吸气努力的方法，包括以下步骤：
8.a.获取单次呼吸的肺部eit图像采样数据，采样频率ωhz；
9.b.将每张肺部eit图像水平划分为大小相等两个区域，分别为腹侧区域和背侧区域。将每个区域的像素值分别累加，记为i
v
(t)和i
d
(t)；
10.c.对i
v
(t)和i
d
(t)求一阶差分，得到i
v
′
(t)和i
d
′
(t)，进而将差分结果进行归一化处理，归一化到γ1到γ2之间；归一化的公式如下：
[0011][0012]
其中，max是i
′
(t)的最大值，min是i
′
(t)的最小值，norm.i
′
(t)表示对i
′
(t)进行归一化处理，归一化到γ1到γ2之间；
[0013]
d.把c步骤得到的背部和腹部区域结果相除，我们称得到的结果为区域密度分数
曲线rif，其中每次呼吸的rif计算公式如下：
[0014][0015]
e.定义呼气阶段rif的峰值为rif
max
，呼气开始时刻一直到幅值达到rif
max
时刻期间，rif的最小值为rif
min
，定义δrif＝rif
max
‑
rif
min
，若δrif>α则该次呼吸判断为无效吸气努力，否则判断为非无效吸气努力；
[0016]
f.将上述步骤应用于呼吸机导出的连续的肺部eit图像采样数据，即能自动检测出符合要求的无效吸气努力。
[0017]
进一步，所述步骤a中，用呼吸机导出的连续的肺部eit图像采样数据，采样频率应高于20hz。
[0018]
再进一步，所述步骤b中，eit图像划分两个区域，分割线是图像正中间的水平线，划分为上半部分为背侧区域，下半部分为腹侧区域。图像每帧像素会随时间变化，所以累加值i
v
(t)和i
d
(t)都是随时间动态变化。
[0019]
所述步骤c中，做一阶差分可以得到变化趋势，归一化是为了防止d步骤的除法运算出现零的情况。
[0020]
更进一步，所述步骤d中，每次呼吸都存在对应的rif(t)，具体的计算过程已经在公式中展示。
[0021]
所述步骤e中，在呼气阶段，可能由于干扰因素存在若干波峰，取最大的波峰为rif
max
，该波峰之前的最小值为rif
min
，两值之差δrif>α则认为是无效吸气努力，如果不存在波峰或者δrif很小等其它情况都是因为肺部潮气量变化幅度不明显，也可以说患者没有产生自发的吸气动作，属于非无效吸气努力。
[0022]
本发明的有益效果主要表现在：通过肺部eit图像准确识别出无效吸气努力的方法，通过分析rif形态特征准确识别患者是否存在自发的吸气努力，甚至检测出从流速压力波形上不明显的无效吸气努力。本发明为医生提供一种检测无效吸气努力的新方法。由于方法是从患者肺部气体分布变化来直接分析患者的吸气动作，而不是从压力流速波形上间接分析，更不需要食道压或膈肌电等复杂昂贵的手段来检测，本发明可以为医生提供一种检测无效吸气努力的新思路。
附图说明
[0023]
图1为本发明的流程示意图。
[0024]
图2为本发明在肺部eit图划分感兴趣区域的示意图。
[0025]
图3为本发明根据rif判断无效吸气努力的示意图。
具体实施方式
[0026]
为了更好的解释本发明，以便于理解，下面结合附图，通过具体实施方式，对本发明作详细描述。
[0027]
参照图1～图3，一种基于电阻抗断层成像识别无效吸气努力的方法，通过连续eit图像数据采样，划分为两个感兴趣区域，再累加区域像素值后差分归一化，最后这两个区域的结果相除得到rif，用rif大小来判断是否为无效吸气努力，该方法是一种更精准更合理
更直接的途径，包括以下步骤：
[0028]
a.获取单次呼吸的肺部eit图像采样数据，采样频率20hz；
[0029]
b.图2给每张肺部eit图像划分为大小相等两个感兴趣区域roi
v
和roi
d
，分别表示背部区域和腹部区域，每个区域的像素值累加。也就是说，每帧都有一张eit图像，对应两个累加值依次为i
v
(t)和i
d
(t)；
[0030]
c.对i
v
(t)和i
d
(t)都求一阶差分依次得到i
′
v
(t)和i
′
d
(t)，进而把差分结果做归一化处理，归一化到0.1到1之间；归一化的公式如下：
[0031][0032]
其中，max是f(t)中的最大值，min是f(t)中的最小值，norm.f(t)表示对f(t)进行归一化处理，归一化到0.1到1之间。把i
′
v
(t)和i
′
d
(t)代入f(t)就可以得到计算结果。
[0033]
d.把背部区域norm.i
′
v
(t)和腹部区域norm.i
′
d
(t)相除，我们称得到的商为区域密度分数曲线(rif)，其中每次呼吸的rif计算公式如下：
[0034][0035]
其中，f
′
(t)表示f(t)的一阶差分；
[0036]
e.如图3所示，绿色虚线是呼气起始位置，在呼气阶段，rif的峰值最大的波峰幅值为rif
max
，呼气开始时刻一直到幅值达到rif
max
时刻期间，rif的最小值为rif
min
。我们定义δrif＝rif
max
‑
rif
min
，若δrif>0.25则该次呼吸判断为无效吸气努力，否则判断为非无效吸气努力；
[0037]
f.将上述步骤应用于eit设备导出的连续的肺部eit图像采样数据，即能自动检测出符合要求的无效吸气努力。
[0038]
本实施例的步骤a中，用eit设备导出的连续的肺部eit图像采样数据，采样频率应在20hz。
[0039]
步骤b中，eit图像划分两个区域，分割线是图像正中间的水平线，划分为上半部分为背部区域，下半部分为腹部区域。图像每帧像素会随时间变化，所以累加值i
v
(t)和i
d
(t)都是随时间动态变化。
[0040]
步骤c中，做一阶差分可以得到变化趋势，归一化是为了防止d步骤的除法运算出现零的情况。
[0041]
步骤d中，每次呼吸都存在对应的rif(t)，具体的计算过程已经在公式中展示。
[0042]
步骤e中，在呼气阶段，可能由于干扰因素存在若干波峰，取最大的波峰为rif
max
，该波峰之前的最小值为rif
min
，两数之差δrif>0.25，则认为是无效吸气努力。如果不存在波峰或者δrif很小等等其他情况都是因为肺部潮气量变化幅度不明显，也可以说患者没有产生自发的吸气动作，属于非无效吸气努力。如图3所示，(a)属于非无效吸气努力，(b)和(c)都属于无效吸气努力，我们可以看出，更具本方法可以很好地区分无效吸气努力，而且(c)这种从流速波形看吸气表现不明显的情况，用rif方法也能精准识别出来。
[0043]
本发明中的采样频率ω根据呼吸机具体型号而定；γ1和γ2在实例中基于经验设定为0.1和1，其中γ1也可以设定为其他大于零的很小的数；α阈值设定范围0.2<α<0.3；在应用中，可根据实际情况分析后做适当调整。
[0044]
本发明是一种通过肺部eit图像准确识别出无效吸气努力的方法，通过分析rif形态特征准确识别患者是否存在自发的吸气努力，甚至检测出从流速压力波形上不明显的无效吸气努力。将本发明应用于可以导出肺部eit图像的呼吸机中，能自动精准的检测到无效吸气努力的结果，更能准确挖掘出真正的无效吸气努力。
[0045]
本发明为医生提供一种检测无效吸气努力的新方法，本方法不是从压力流速波形上间接分析，更不需要食道压或膈肌电等复杂昂贵的方式，所以本发明可以为医生提供一种检测无效吸气努力的新思路。自动快速精准的检测减少了医护人员的工作量。将来应用到呼吸机上可以实时的提醒医生，可以更好的指导进一步肺保护通气的实施。
[0046]
对6022次呼吸周期进行算法评估。算法检测出3865个无效吸气努力以及2157个非无效吸气努力。医生专家最后评估核实，算法检测出的无效吸气努力中有2945个正确，920个错误；而算法检测出的非无效吸气努力全部正确，表明算法具有很高的灵敏度和特异性，可以说，算法检测不会漏掉任何一个无效吸气努力，在医学上这种高精度预警指标是非常重要的。
[0047]
本说明书的实施例所述的内容仅仅是对发明构思的实现形式的列举，仅作说明用途。本发明的保护范围不应当被视为仅限于本实施例所陈述的具体形式，本发明的保护范围也及于本领域的普通技术人员根据本发明构思所能想到的等同技术手段。