一种心电特征提取方法、装置及存储介质

1.本发明涉及心电监测技术领域，具体而言，涉及一种心电特征提取方法、装置及存储介质。


背景技术：

2.心电图是临床上用于诊断心血管疾病的技术手段，可辅助医护人员辨识心率失常和心脏病变等的性质和严重程度。其中，心电信号属于微弱信号，容易受到50hz工频和外界其它噪声的干扰。
3.为了降低干扰，从心电信号中提取出准确的r波等心电特征，进而提高心率的计算准确度。目前，常通过改善模拟前端电路提高前级输入阻抗的方式来处理噪声，配合主动电极可使输入阻抗达到400gω。但是这种方式可能会造成对心电信号的过度处理，带来明显的基线偏移等问题，导致心电信号失真，无法有效提取到准确的心电特征，进而无法准确体现出心电信号中的特征和异常信息等。


技术实现要素：

4.本发明解决的问题是现有技术中提取心电特征时，需要对心电信号进行深度滤波以滤除噪声，会导致心电信号失真，进而无法提取到准确的心电特征的问题。
5.为解决上述问题，本发明提供一种心电特征提取方法、装置及存储介质。
6.第一方面，本发明提供了一种心电特征提取方法，包括：
7.获取心电信号，将所述心电信号中的心电数据依次存入缓冲区，并确定所述缓冲区内所述心电数据的第一平均值；
8.采用预设长度的窗口滑动截取所述心电信号，确定所述窗口内所有所述心电数据的第二平均值；
9.判断所述第一平均值和所述第二平均值是否满足预设规则，若是，则在所述窗口中确定所述心电数据的局部极大值；
10.当确定的所述局部极大值达到预设数量时，在所有所述局部极大值中确定最大值所在的波形为r波。
11.可选地，所述确定所述缓冲区内所述心电数据的第一平均值包括：
12.当所述心电信号中的第一个所述心电数据存入所述缓冲区时，确定所述缓冲区内所有所述心电数据的所述第一平均值为初始平均值；
13.当所述缓冲区内所述心电数据的数据长度达到预设阈值时，采用第一公式对所述第一平均值进行更新，所述第一公式包括：
[0014][0015]
其中，mb为更新后的第一平均值，lb为所述缓冲区内所述心电数据的数据长度，m0为所述初始平均值。
[0016]
可选地，所述预设规则采用第二公式表示，所述第二公式包括：
[0017][0018]
其中，mb为所述更新后的第一平均值，m0为所述初始平均值，λ为权重因子，i表示所述心电信号中的第i个所述心电数据，表示与第i个所述心电数据对应的窗口内所有所述心电数据的第二平均值，为与第i-1个所述心电数据对应的窗口内所有所述心电数据的第二平均值。
[0019]
可选地，所述在所有所述局部极大值中确定最大值所在的波形为r波之后，还包括：
[0020]
以n个所述局部极大值为周期采用第三公式更新所述权重因子，所述第三公式包括：
[0021][0022]
其中，λ为所述权重因子，为第k个所述局部极大值，m0为所述初始平均值。
[0023]
可选地，所述在所有所述局部极大值中确定最大值所在的波形为r波之后，还包括：
[0024]
根据确定的所述r波采用第四公式计算心率，所述第四公式包括：
[0025][0026]
其中，hr表示所述心率，nr表示检测到的所述心电信号中的r波数，t
scale
表示所述心电信号的时间长度。
[0027]
可选地，所述将所述心电信号中的心电数据依次存入缓冲区之前，还包括：对所述心电信号进行数字滤波。
[0028]
可选地，所述对所述心电信号进行数字滤波包括：
[0029]
将所述心电信号输入无限冲激响应数字滤波器，通过所述无限冲激响应数字滤波器消除所述心电信号中的基线漂移，获得第一滤波信号；
[0030]
将所述第一滤波信号输入陷波器，通过所述陷波器滤除工频噪声，获得第二滤波信号；
[0031]
将所述第二滤波信号输入二阶高通滤波器，通过所述二阶高通滤波器滤除低频噪声，获得滤波后的心电信号。
[0032]
第二方面，本发明提供了一种心电特征提取装置，包括：
[0033]
获取模块，用于获取心电信号，将所述心电信号中的心电数据依次存入缓冲区，并确定所述缓冲区内所述心电数据的第一平均值；
[0034]
处理模块，用于采用预设长度的窗口滑动截取所述心电信号，确定所述窗口内所有所述心电数据的第二平均值；
[0035]
判断模块，用于判断所述第一平均值和所述第二平均值是否满足预设规则，若是，则在所述窗口中确定所述心电数据的局部极大值；
[0036]
提取模块，用于当确定的所述局部极大值达到预设数量时，在所有所述局部极大值中确定最大值所在的波形为r波。
[0037]
第三方面，本发明提供了一种电子设备，包括存储器和处理器；
[0038]
所述存储器，用于存储计算机程序；
[0039]
所述处理器，用于当执行所述计算机程序时，实现如第一方面任一项所述的心电特征提取方法。
[0040]
第四方面，本发明提供了一种计算机可读存储介质，所述存储介质上存储有计算机程序，当所述计算机程序被处理器执行时，实现如第一方面任一项所述的心电特征提取方法。
[0041]
本发明的心电特征提取方法、装置及存储介质的有益效果是：将获取的心电信号中的心电数据依次存入缓冲区，确定缓冲区内心电数据的第一平均值。采用预设长度的窗口对心电信号进行滑动处理，确定窗口内所有心电数据的第二平均值，根据预设规则比对第一平均值和第二平均值，根据比对结果确定心电信号中心电数据的局部极大值。r波为心电信号中的正向波，通过确定局部极大值缩小r波所在的范围，然后在预设数量的局部极大值中确定最大值所在的波形就是r波。本发明的技术方案对心电信号中的心电数据进行分析，准确定位心电信号中的r波，降低了对心电信号滤波过程的要求，避免滤波过程中对心电信号进行过度处理导致心电信号失真，提高了提取的心电特征的准确性。
附图说明
[0042]
图1为本发明实施例的一种心电特征提取方法的流程示意图；
[0043]
图2为本发明另一实施例的一种心电特征提取装置的结构示意图。
具体实施方式
[0044]
为使本发明的上述目的、特征和优点能够更为明显易懂，下面结合附图对本发明的具体实施例做详细的说明。
[0045]
需要说明的是，本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换，以便这里描述的本发明的实施例能够以除了在这里图示或描述的那些以外的顺序实施。
[0046]
如图1所示，本发明实施例提供的一种心电特征提取方法，包括：
[0047]
步骤s110，获取心电信号，将所述心电信号中的心电数据依次存入缓冲区，并确定所述缓冲区内所述心电数据的第一平均值。
[0048]
步骤s120，采用预设长度的窗口滑动截取所述心电信号，确定所述窗口内所有所述心电数据的第二平均值。
[0049]
具体地，可采用滑动窗口滤波方法对心电信号进行平滑处理，能够使得心电信号的波动更加平滑，对于一个心电数据，采用该心电数据对应的窗口内所有心电数据的平均值进行替代。窗口的预设长度可为其中，fs为心电信号的频率。
[0050]
步骤s130，判断所述第一平均值和所述第二平均值是否满足预设规则，若是，则在
所述窗口中确定所述心电数据的局部极大值。
[0051]
步骤s140，当确定的所述局部极大值达到预设数量时，在所有所述局部极大值中确定最大值所在的波形为r波。
[0052]
具体地，根据滑动窗口和缓冲区存储的方式依次对心电信号中的各个心电数据进行处理，获得多个局部极大值。当局部极大值的数量达到预设数量时，预设数量可优选为8，确定局部极大值中最大值所在的波形为r波。
[0053]
本实施例中，将获取的心电信号中的心电数据依次存入缓冲区，确定缓冲区内心电数据的第一平均值。采用预设长度的窗口对心电信号进行滑动处理，确定窗口内所有心电数据的第二平均值，根据预设规则比对第一平均值和第二平均值，根据比对结果确定心电信号中心电数据的局部极大值。r波为心电信号中的正向波，通过确定局部极大值缩小r波所在的范围，然后在预设数量的局部极大值中确定最大值所在的波形就是r波。本发明的技术方案对心电信号中的心电数据进行分析，准确定位心电信号中的r波，降低了对心电信号滤波过程的要求，避免滤波过程中对心电信号进行过度处理导致心电信号失真，提高了提取的心电特征的准确性。
[0054]
可选地，所述确定所述缓冲区内所述心电数据的第一平均值包括：
[0055]
当所述心电信号中的第一个所述心电数据存入所述缓冲区时，确定所述缓冲区内所有所述心电数据的所述第一平均值为初始平均值；
[0056]
当所述缓冲区内所述心电数据的数据长度达到预设阈值时，采用第一公式对所述第一平均值进行更新，所述第一公式包括：
[0057][0058]
其中，mb为更新后的第一平均值，lb为所述缓冲区内所述心电数据的数据长度，m0为所述初始平均值。
[0059]
具体地，当缓冲区内心电数据的数据长度lb≥2fs时，fs表示心电信号的频率，即经过上述步骤处理后，已处理的心电数据的长度超过预设阈值的两倍，假定到达r波出现阶段，缩小r波的出现范围，便于进一步判断以确定r波的精确位置。
[0060]
可选地，所述预设规则采用第二公式表示，所述第二公式包括：
[0061][0062]
其中，mb为所述更新后的第一平均值，m0为所述初始平均值，λ为权重因子，其随着心电信号动态调整，i表示所述心电信号中的第i个所述心电数据，表示与第i个所述心电数据对应的窗口内所有所述心电数据的第二平均值，为与第i-1个所述心电数据对应的窗口内所有所述心电数据的第二平均值。
[0063]
具体地，当第一平均值和第二平均值满足第二公式时，确定第i-1个所述心电数据为局部极大值，可记为该局部极大值所在位置为可能的r波位置。通过确定心电信号中的局部极大值，进一步缩小r波位置的范围，便于后续精确定位r波的位置。
[0064]
可选地，所述在所有所述局部极大值中确定最大值所在的波形为r波之后，还包
括：
[0065]
以n个所述局部极大值为周期采用第三公式更新所述权重因子，所述第三公式包括：
[0066][0067]
其中，λ为所述权重因子，为第k个所述局部极大值，m0为所述初始平均值。
[0068]
具体地，n优选为8，即第三公式优选为：
[0069][0070]
本可选的实施例中，根据心电信号中心电数据的变化动态调整权重因子的数值，使权重因子适应于心电信号的变化，提高局部极大值的判断精度。
[0071]
可选地，所述在所有所述局部极大值中确定最大值所在的波形为r波之后，还包括：
[0072]
根据确定的所述r波采用第四公式计算心率，所述第四公式包括：
[0073][0074]
其中，hr表示所述心率，nr表示检测到的所述心电信号中的r波数，t
scale
表示所述心电信号的时间长度，需将其换算为分钟进行表示，因此需要除以60s。
[0075]
本可选的实施例中，通过本发明实施例的心电特征提取方法能够准确定位心电信号中的r波，不需要预先采用复杂滤波电路对心电信号进行深入滤波，能够防止滤波过程中造成的信号失真，提高了r波的提取准确度，进而提高了通过r波计算心率的准确度。
[0076]
可选地，所述将所述心电信号中的心电数据依次存入缓冲区之前，还包括：对所述心电信号进行数字滤波。
[0077]
具体地，虽然采用复杂滤波电路对心电信号进行滤波，可能会造成心电信号失真，但是若心电信号中的噪声较多，也会影响r波的提取精度，因此需要在保证不造成信号失真，同时尽可能降低噪声对r波提取的影响的前提下，对心电信号进行合理滤波。
[0078]
可选地，所述对所述心电信号进行数字滤波包括：
[0079]
将所述心电信号输入无限冲激响应数字滤波器，通过所述无限冲激响应数字滤波器消除所述心电信号中的基线漂移，获得第一滤波信号；
[0080]
将所述第一滤波信号输入陷波器，通过所述陷波器滤除工频噪声，获得第二滤波信号；
[0081]
将所述第二滤波信号输入二阶高通滤波器，通过所述二阶高通滤波器滤除低频噪声，获得滤波后的心电信号。
[0082]
具体地，通过对心电信号波形进行分析，可以发现，除了r波清晰以外，其它特征波形几乎都被噪声淹没。对噪声进行功率谱分析，可以发现，心电信号中的噪声重要来自50hz的工频干扰及其高次谐波干扰，其次还有电路引入的噪声和人体抖动(例如呼吸和轻微颤抖等)引起的低频噪声。因此，首先可通过无限冲激响应数字滤波器消除呼吸等人体抖动引起的基线漂移，可通过设定两个截止频率来截取原始心电数据中有用的信号，设两个截止
频率分别为f1和f2，则基于巴特沃斯的带通滤波器可采用下列公式表示：
[0083][0084]
其中，h(ω)表示滤波后的信号，ω表示输入的心电信号的频率，ωc＝[f1,f2]为截止频率，通过上述公式可以确定巴特沃斯数字滤波器的滤波系数。由于心电信号的动态性，当获得巴特沃斯滤波器系数后，通常使用零相位数字滤波器对心电信号进行降噪。
[0085]
其次，通过一个50hz的陷波器，对工频噪声进行滤除；最后，通过一个截止频率为0.5hz的二阶高通滤波器，滤除心电信号中的低频噪声。
[0086]
本可选的实施例中，对心电信号进行数字滤波，能够针对性地滤除心电信号中的各种噪声，滤波效果好，同时避免对心电信号的过度处理，不会造成心电信号失真，影响提取的心电特征的准确性。
[0087]
可通过第五公式对滤波后的心电信号的信噪比进行量化评估，第五公式包括：
[0088][0089]
其中，snr
ecg
为信噪比，代表滤波后的心电信号中qrs波群在该时间段内的电压有效值，表示在一个心电周期末的t波和下一个心电周期的p波之间的，且与qrs波群等长度的电压有效值。
[0090]
如图2所示，本发明另一实施例提供的一种心电特征提取装置，包括：
[0091]
获取模块，用于获取心电信号，将所述心电信号中的心电数据依次存入缓冲区，并确定所述缓冲区内所述心电数据的第一平均值；
[0092]
处理模块，用于采用预设长度的窗口滑动截取所述心电信号，确定所述窗口内所有所述心电数据的第二平均值；
[0093]
判断模块，用于判断所述第一平均值和所述第二平均值是否满足预设规则，若是，则在所述窗口中确定所述心电数据的局部极大值；
[0094]
提取模块，用于当确定的所述局部极大值达到预设数量时，在所有所述局部极大值中确定最大值所在的波形为r波。
[0095]
本发明又一实施例提供的一种电子设备包括存储器和处理器；所述存储器，用于存储计算机程序；所述处理器，用于当执行所述计算机程序时，实现如上所述的心电特征提取方法。
[0096]
本发明又一实施例提供的一种计算机可读存储介质上存储有计算机程序，当所述计算机程序被处理器执行时，实现如上任一项所述的心电特征提取方法。
[0097]
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的程序可存储于一计算机可读取存储介质中，该程序在执行时，可包括如上述各方法的实施例的流程。其中，所述的存储介质可为磁碟、光盘、只读存储记忆体(read-only memory，rom)或随机存储记忆体(random access memory，ram)等。在本技术中，所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现
本发明实施例方案的目的。另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以是两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。
[0098]
虽然本发明公开披露如上，但本发明公开的保护范围并非仅限于此。本领域技术人员在不脱离本发明公开的精神和范围的前提下，可进行各种变更与修改，这些变更与修改均将落入本发明的保护范围。