一种实时的单导联心电R波检测方法与流程

一种实时的单导联心电r波检测方法
技术领域
1.本发明涉及健康检测技术领域，特别地是一种实时的单导联心电r波检测方法。


背景技术：

2.心血管疾病作为目前最难以治愈的疾病之一，有其不易察觉，不易治疗的特点。而心电图的判读诊断，在无创的心血管疾病诊断中，比起b超等有着更安全，更易用的优势。而i导联的心电设备，比起多导联的设备，更有全天候的使用的潜力。而在单导联设备中心电图更容易受到噪声干扰，所以在单导联设备上的r波检测方法需要有更高的鲁棒性。
3.现有的qrs波检测大多基于离线分析的方式，实时性不好，且针对于单导联设备的心电信号容易出现错检和漏检等情况。


技术实现要素：

4.本发明的目的在于提供一种具有良好的实时性且鲁棒性好的实时单导联心电r波检测方法。
5.本发明通过以下技术方案实现的：
6.一种实时的单导联心电r波检测方法，包括以下步骤：
7.步骤s1，原始心电数据预处理；
8.步骤s2，将预处理信号以数据段形式输入r波检测模块；
9.步骤s3，r波检测模块对数据段信号分析，计算得到r波位置信息。
10.进一步地，所述步骤s1包括以下步骤：
11.步骤s11，去除原始心电信号的基线漂移；
12.步骤s12，去除原始心电信号的工频干扰；
13.步骤s12，去除原始心电信号的高频毛刺噪声。
14.进一步地，所述步骤s3包括以下步骤：
15.步骤s31，采用缓存机制存储上一个数据段内容；
16.步骤s32，分离qrs波高频分量；
17.步骤s34，基于qrs波高频分量分析信号的r波信息。
18.进一步地，所述步骤31包括以下步骤：
19.步骤s311，算法初始化时创建一个两倍数据段长度的队列，队列中数据置零；步骤s312，有数据段输入r波检测模块时将其添加到队列，并同时从队列中清除一个数据段。
20.进一步地，所述步骤s32包括以下步骤：
21.步骤s321，使用队列中第一个数据段的后端数据对第二个数据段的前端进行扩展以形成新数据段；
22.步骤s322，基于特定的fir低通滤波器滤除新数据段的高频分量；
23.步骤s323，将了滤波前的数据段与滤波后的低频数据段做差，得到qrs波的高频分量；
24.步骤s324，重复步骤s311和步骤s312，再次创建并维护一个队列，用来存储当前时段与上一时段的qrs波高频分量。
25.进一步地，所述步骤s3还包括步骤s33，所述步骤s33包括以下步骤：
26.步骤s331，遍历qrs波的高频分量，当为信号段起始点时，判断标志位为1，分情况处理；
27.步骤s332，若标志位为0，进行一般的r波检测操作，判断该点幅值和与上一个r波的距离是否符合标准，若该点符合标准，暂停遍历，将该点设为候选点，进行下一步操作；
28.步骤s333，计算该点当前位置之后的0.6倍平均rr间距是否超出该高频分量的范围，若超出计算当前位置到末尾的最大值，并记录坐标a，并设定标志位为1，分情况处理；
29.步骤s334，若该点当前位置之后的0.6倍平均rr间距未超出该高频分量数据段的范围，计算该范围内的最大值坐标，判断是否与当前位置一致；分情况处理；
30.步骤s335，若当前位置与最大值坐标一致，该点为r波位置，更新自适应系数；
31.步骤s336，在当前检测出r波位置的基础上后移上一次0.6倍平均rr间距的距离继续遍历。
32.进一步地，所述步骤s331包括以下步骤：
33.步骤s3311，判断标志位为1；
34.步骤s3312，将qrs波高频分量队列中的新旧两个数据段拼接，检测从上段数据记录坐标a到之后0.6倍平均rr间距的范围内的最大值；
35.步骤s3313，跳转到步骤s335，标志位设0。
36.进一步地，所述步骤s333包括以下步骤：
37.步骤s3311，当前位置之后的0.6倍平均rr间距是超出该高频分量的范围，记录坐标a；
38.步骤s3312，标志位设为1；
39.步骤s3313，将当前位置设为数据段末尾，结束当前数据段遍历。
40.进一步地，所述步骤s334包括以下步骤：
41.步骤s3341，判断当前位置与最大值位置不一致；
42.步骤s3342，基于两个位置生成qrs模板；
43.步骤s3343，分别计算这两个qrs模板与自适应qrs模板的误差；
44.步骤s3344，若当前位置的qrs模板误差小于最大值位置的误差，且rr间隔更接近平均的rr间隔，则该点为r峰的位置，更新自适应系数，跳转到步骤s336；
45.步骤s3345，若当前位置的qrs模板误差大于最大值位置的误差，设定当前位置为最大值位置，跳转到步骤s33。
46.进一步地，所述步骤s335包括以下步骤：
47.步骤s3351，判断当前位置与最大值位置一致；
48.步骤s3352，将符合一定标准的rr间隔加入自适应rr间隔列表；
49.步骤s3353，将符合一定标准的r波峰值加入自适应r波峰值列表；
50.步骤s3354，将该点所对应的qrs模板以0.1倍加到自适应模板上，并归一化自适应模板。
51.本发明的有益效果：
52.本发明通过分时段的心电信号处理方式以实现实时功能，相比与离线r波检测方法更具实时性。为了保证检测的连续性，算法中使用缓存机制临时存储上一个时段的输入信号段。该方法首先基于数字信号处理技术分别对信号进行去除基线漂移、工频干扰记忆高频毛刺噪声的处理，之后由数据段形式交由r检测主模块中的一个特定的fir低通滤波器提取信号中的低频分量并与原信号做差。该差值中含有清晰的qrs波分量且噪声干扰少，相比于直接基于原信号的r波检测算法，基于该差值检测方法更强的抗噪声能力，之后基于次差值进行r波检测。该模块中同时三个自适应系数，包括自适应qrs模板、自适应rr间隔列表，自适应r波峰值列表，通过对这三个系数的使用以有效避免漏检和错减的情况。除滤除“基线漂移”、“工频干扰”以及“高频毛刺噪声”的处理操作之外，该r波检测方法的输出延迟时间约等于输入的信号段的时间长度，具有良好的实时性且鲁棒性好。
附图说明
53.图1为本发明实施例一种实时的单导联心电r波检测方法的流程图；
54.图2为本发明实施例r波检测模块的流程图。
具体实施方式
55.下面将结合附图以及具体实施例来详细说明本发明，在此以本发明的示意下面将结合附图以及具体实施例来详细说明本发明，在此以本发明的示意性实施例及说明用来解释本发明，但并不作为对本发明的限定。
56.需要说明，在本发明中如涉及“第一”、“第二”的描述仅用于描述目的，而不能理解为指示或暗示其相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。另外，各个实施例之间的技术方案可以相互结合，但是必须是以本领域普通技术人员能够实现为基础，当技术方案的结合出现相互矛盾或无法实现时应当认为这种技术方案的结合不存在，也不在本发明要求的保护范围之内。
57.在本发明中，除非另有明确的规定和限定，术语“连接”应做广义理解，例如，“连接”可以是固定连接，也可以是可拆卸连接，或成一体；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通或两个元件的相互作用关系，除非另有明确的限定。对于本领域的普通技术人员而言，可以根据具体情况理解上述术语在本发明中的具体含义。
58.如图1和图2所示，一种实时的单导联心电r波检测方法，包括以下步骤：
59.步骤s1，原始心电数据预处理；
60.步骤s2，将预处理信号以数据段形式输入r波检测模块；
61.步骤s3，r波检测模块对数据段信号分析，计算得到r波位置信息。
62.具体的，本实施例方案中，所述步骤s1包括以下步骤：
63.步骤s11，去除原始心电信号的基线漂移；
64.步骤s12，去除原始心电信号的工频干扰；
65.步骤s12，去除原始心电信号的高频毛刺噪声。
66.具体的，本实施例方案中，所述步骤s3包括以下步骤：
67.步骤s31，采用缓存机制存储上一个数据段内容；
68.步骤s32，分离qrs波高频分量；
69.步骤s34，基于qrs波高频分量分析信号的r波信息。
70.具体的，本实施例方案中，所述步骤31包括以下步骤：
71.步骤s311，算法初始化时创建一个两倍数据段长度的队列，队列中数据置零；
72.步骤s312，有数据段输入r波检测模块时将其添加到队列，并同时从队列中清除一个数据段。
73.具体的，本实施例方案中，所述步骤s32包括以下步骤：
74.步骤s321，使用队列中第一个数据段的后端数据对第二个数据段的前端进行扩展以形成新数据段；
75.步骤s322，基于特定的fir低通滤波器滤除新数据段的高频分量；
76.步骤s323，将了滤波前的数据段与滤波后的低频数据段做差，得到qrs波的高频分量；
77.步骤s324，重复步骤s311和步骤s312，再次创建并维护一个队列，用来存储当前时段与上一时段的qrs波高频分量。
78.具体的，本实施例方案中，所述步骤s3还包括步骤s33，所述步骤s33包括以下步骤：
79.步骤s331，遍历qrs波的高频分量，当为信号段起始点时，判断标志位为1，分情况处理；
80.步骤s332，若标志位为0，进行一般的r波检测操作，判断该点幅值和与上一个r波的距离是否符合标准，若该点符合标准，暂停遍历，将该点设为候选点，进行下一步操作；
81.步骤s333，计算该点当前位置之后的0.6倍平均rr间距是否超出该高频分量的范围，若超出计算当前位置到末尾的最大值，并记录坐标a，并设定标志位为1，分情况处理；
82.步骤s334，若该点当前位置之后的0.6倍平均rr间距未超出该高频分量数据段的范围，计算该范围内的最大值坐标，判断是否与当前位置一致；分情况处理；
83.步骤s335，若当前位置与最大值坐标一致，该点为r波位置，更新自适应系数；
84.步骤s336，在当前检测出r波位置的基础上后移上一次0.6倍平均rr间距的距离继续遍历。
85.具体的，本实施例方案中，所述步骤s331包括以下步骤：
86.步骤s3311，判断标志位为1；
87.步骤s3312，将qrs波高频分量队列中的新旧两个数据段拼接，检测从上段数据记录坐标a到之后0.6倍平均rr间距的范围内的最大值；
88.步骤s3313，跳转到步骤s335，标志位设0。
89.具体的，本实施例方案中，所述步骤s333包括以下步骤：
90.步骤s3311，当前位置之后的0.6倍平均rr间距是超出该高频分量的范围，记录坐标a；
91.步骤s3312，标志位设为1；
92.步骤s3313，将当前位置设为数据段末尾，结束当前数据段遍历。
93.具体的，本实施例方案中，所述步骤s334包括以下步骤：
94.步骤s3341，判断当前位置与最大值位置不一致；
95.步骤s3342，基于两个位置生成qrs模板；
96.步骤s3343，分别计算这两个qrs模板与自适应qrs模板的误差；
97.步骤s3344，若当前位置的qrs模板误差小于最大值位置的误差，且rr间隔更接近平均的rr间隔，则该点为r峰的位置，更新自适应系数，跳转到步骤s336；
98.步骤s3345，若当前位置的qrs模板误差大于最大值位置的误差，设定当前位置为最大值位置，跳转到步骤s33。
99.具体的，本实施例方案中，所述步骤s335包括以下步骤：
100.步骤s3351，判断当前位置与最大值位置一致；
101.步骤s3352，将符合一定标准的rr间隔加入自适应rr间隔列表；
102.步骤s3353，将符合一定标准的r波峰值加入自适应r波峰值列表；
103.步骤s3354，将该点所对应的qrs模板以0.1倍加到自适应模板上，更新后的qrs模板再除其数据内部的最大值进行标准化。或将该点所对应的qrs模板以0.1倍加到自适应模板上，并归一化自适应模板。
104.本发明通过分时段的心电信号处理方式以实现实时功能，相比与离线r波检测方法更具实时性。为了保证检测的连续性，算法中使用缓存机制临时存储上一个时段的输入信号段。该方法首先基于数字信号处理技术分别对信号进行去除基线漂移、工频干扰记忆高频毛刺噪声的处理，之后由数据段形式交由r检测主模块中的一个特定的fir低通滤波器提取信号中的低频分量并与原信号做差。该差值中含有清晰的qrs波分量且噪声干扰少，相比于直接基于原信号的r波检测算法，基于该差值检测方法更强的抗噪声能力，后基于次差值进行r波检测。该模块中同时三个自适应系数，包括自适应qrs模板、自适应rr间隔列表，自适应r波峰值列表，通过对这三个系数的使用以有效避免漏检和错减的情况。除滤除“基线漂移”、“工频干扰”以及“高频毛刺噪声”的处理操作之外，该r波检测方法的输出延迟时间约等于输入的信号段的时间长度，具有良好的实时性且鲁棒性好。
105.以上对本发明实施例所提供的技术方案进行了详细介绍，本文中应用了具体个例对本发明实施例的原理以及实施方式进行了阐述，以上实施例的说明只适用于帮助理解本发明实施例的原理；同时，对于本领域的一般技术人员，依据本发明实施例，在具体实施方式以及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。