CN113349754A - 一种基于便携式动态心电监护仪的心率计算方法 - Google Patents
一种基于便携式动态心电监护仪的心率计算方法 Download PDFInfo
- Publication number
- CN113349754A CN113349754A CN202110816029.7A CN202110816029A CN113349754A CN 113349754 A CN113349754 A CN 113349754A CN 202110816029 A CN202110816029 A CN 202110816029A CN 113349754 A CN113349754 A CN 113349754A
- Authority
- CN
- China
- Prior art keywords
- data
- heart rate
- peak
- steps
- portable dynamic
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 27
- 238000000034 method Methods 0.000 claims abstract description 45
- 238000012545 processing Methods 0.000 claims abstract description 33
- 238000001914 filtration Methods 0.000 claims abstract description 27
- 238000007781 pre-processing Methods 0.000 claims abstract description 12
- 230000002159 abnormal effect Effects 0.000 claims abstract description 10
- 238000000718 qrs complex Methods 0.000 claims description 12
- 238000001514 detection method Methods 0.000 claims description 8
- 230000002457 bidirectional effect Effects 0.000 claims description 7
- 238000005259 measurement Methods 0.000 claims description 6
- 238000012935 Averaging Methods 0.000 claims description 4
- 230000008030 elimination Effects 0.000 claims description 4
- 238000003379 elimination reaction Methods 0.000 claims description 4
- 238000012544 monitoring process Methods 0.000 description 5
- 230000005189 cardiac health Effects 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 230000036541 health Effects 0.000 description 3
- 238000013186 photoplethysmography Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 208000024172 Cardiovascular disease Diseases 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000007774 longterm Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005485 electric heating Methods 0.000 description 1
- 230000005802 health problem Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 230000033764 rhythmic process Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/024—Detecting, measuring or recording pulse rate or heart rate
- A61B5/02438—Detecting, measuring or recording pulse rate or heart rate with portable devices, e.g. worn by the patient
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
- A61B5/332—Portable devices specially adapted therefor
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
- A61B5/346—Analysis of electrocardiograms
- A61B5/349—Detecting specific parameters of the electrocardiograph cycle
- A61B5/352—Detecting R peaks, e.g. for synchronising diagnostic apparatus; Estimating R-R interval
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
- A61B5/346—Analysis of electrocardiograms
- A61B5/349—Detecting specific parameters of the electrocardiograph cycle
- A61B5/366—Detecting abnormal QRS complex, e.g. widening
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7225—Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/725—Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
Abstract
本发明公开了一种基于便携式动态心电监护仪的心率计算方法,包括以下步骤:步骤S1:步骤S2:数据接收;将接收到的数据记为S;步骤S3:数据预处理;对步骤S2中的数据S进行预处理,去除其中的异常数据点,得到数据S1;步骤S4:数据滤波处理;通过信号滤波算法对步骤S3中的数据进行滤波处理,去除其中的噪声干扰,得到数据S2;步骤S5:更新缓存,利用缓存数据计算心率;步骤S6:数据正则化;步骤S7:心率及心电图波形展示。本发明使用用户的心电数据计算心率,可信度更高。在计算心率时提前进行了数据预处理与信号滤波,同时利用缓存池的方法来计算心率,使心率计算更加稳定。在进行心电图数据展示时进行了正则化处理,数据的可读性更好。
Description
技术领域
本发明属于心率计算技术领域,具体为一种基于便携式动态心电监护仪的心率计算方法。
背景技术
心脏健康属于人体最重要的几项健康指标之一,根据世界卫生组织统计,目前心血管疾病的死亡率居全球首位。对于心血管疾病高发人群,长时间的持续心脏监测可以帮助提前发现潜在危险,提醒患者及时就医治疗。在心脏的各项指标中,心率可以最直接地反映人体心脏的健康状态,它代表心脏每分钟跳动的次数,心率过快、过慢或心率不齐都可能反映出心脏存在潜在的健康问题。
在医院中,最普遍的检查方法是使用动态心电图机进行心脏检查,但是其测量时间较短,无法实现长期监测;且携带起来很不方便,容易影响人们正常的生活状态。便携式动态心电监护仪则可以有效解决此问题,其体积小巧,佩戴方便,能够在不影响人们正常生活状态的情况下持续记录佩戴者的心电数据。
数字信号处理可以利用数字方法对信号进行分析、变换、滤波、检测、调制、解调等操作,从而获取到人们所需要的特定属性或知识。借助数字信号处理技术,可以对采集到的心电数据进行分析处理,让用户知晓自己的实时心率,掌握自身心脏健康状态。
目前市面上存在多种多样的可穿戴智能设备,可以对用户心率进行监测,如手环、手表等设备。但是这类设备通常是利用光电容积脉搏波描记法(PPG)来计算心率的,通过光源照射皮肤,监测反射光或者折射光信号的周期性来对心率进行计算。
在同类型的动态心电监护仪产品中,一般是直接在动态心电监护仪采集信号时由硬件计算心率,并展示原始心电图波形。但这种方法计算出来的心率不稳定,波动幅度较大;其展示的心电图波形由于未经过处理,可读性也较差。
发明内容
本发明的目的在于提供一种基于便携式动态心电监护仪的心率计算方法,以解决背景技术中提出的现有技术计算出来的心率不稳定,波动幅度较大,展示的心电图波形可读性较差的问题。
为解决上述技术问题,本发明所采用的技术方案是:
一种基于便携式动态心电监护仪的心率计算方法,包括以下步骤:
步骤S1:将便携式动态心电监护仪贴在人体表面;
步骤S2:数据接收;通过移动终端与便携式动态心电监护仪无线连接,接收便携式动态心电监护仪所检测的心电信号数据,将接收到的数据记为S;
步骤S3:数据预处理;对步骤S2中的数据S进行预处理,去除其中的异常数据点,得到数据S1;
步骤S4:数据滤波处理;通过信号滤波算法对步骤S3中的数据进行滤波处理,去除其中的噪声干扰,得到滤波后的数据S2;
步骤S5:更新缓存,利用缓存数据计算心率;
步骤S6:数据正则化;
步骤S7:心率及心电图波形展示;
步骤S8:重复步骤S2-S7直至测量结束。
根据上述技术方案,所述步骤S3中的具体处理为:
步骤S31:确定S的取值范围为-1mV到1mV:
步骤S32:将S中大于1mV的值置为1mV,小于-1mV的值置为-1mV,得到处理后的数据S1,处理后的S1即为去除掉异常数据点的数据。
根据上述技术方案,所述步骤S4中的具体处理为:
步骤S41:使用0.5Hz-80Hz的带通滤波器进行滤波,去除小于0.5Hz的低频噪声及大于80Hz的高频噪声干扰,得到数据S11;
步骤S42:首先对S11使用1Hz低通滤波器得到S11wander,再使用S11减去S11wander,得到去除了基线漂移影响的S12;
步骤S43:对步骤S42得到的S12使用48Hz-52Hz的带阻滤波器,去除50Hz左右的工频干扰,经过以上滤波处理,得到处理后的数据S2。
根据上述技术方案,所述步骤S5中的具体处理为:
步骤S51:在接收数据的移动终端建立一个数据缓存区;
步骤S52:在数据缓存区建立一个双向队列;双向队列用于记录当前时刻前4秒的所有ECG数据为S-4s;
步骤S53:当有新的数据S2到来时,进行入队列操作,更新缓存数据S-4s,并使用数据S-4s计算心率。
根据上述技术方案,通过数据S-4s计算心率的方法为:
步骤S54:首先使用XQRS检测算法检测缓存数据S-4s中的所有R峰位置[r1,r2,...,rn];
步骤S55:根据步骤S54中检测到的所有R峰位置计算所有相邻R峰之间的RR间期,得到RR间期序列RRs=[r2-r1,r3-r2,...,rn-rn-1];其中,每个RR间期的单位为秒;
步骤S56:最后根据RRs计算平均心率HR,计算方法为HR=60/mean(RRs),其中mean为取均值函数。
根据上述技术方案,所述心率的计算步骤具体为:
步骤S57:将缓存数据S-4s通过5Hz-20Hz的带通滤波器进行过滤,得到滤波后的数据S-4s-f;
步骤S58:将数据S-4s-f通过Ricker小波变换进行处理,提取QRS波群信息,得到S-4s-i;
步骤S59:从S-4s-i中寻找所有局部最大值点组成R峰候选项;
步骤S510:设置QRS波群检测阈值相关参数,参数包括:R峰处的幅值应大于0.13mV、该QRS波群的R峰位置距离上一个QRS波群的R峰位置应大于0.2秒;
步骤S511:对每一个R峰候选项按时间顺序进行遍历,根据阈值条件判断其是否属于真正的QRS波群,同时根据当前波形中的R峰幅值与位置,更新步骤S510中的R峰幅值与上一个R峰位置阈值参数;
步骤S512:重复步骤S511直到遍历完所有R峰候选项,得到所有真正的R峰位置;
步骤S513:根据RR间期计算该段数据的平均心率;
步骤S514:输出心率结果。
根据上述技术方案,所述步骤S6的具体处理为:
步骤S61:首先采用S3=(S2-μ)/σ的方法对数据S2进行标准化,其中μ为S2的均值μ=mean(S2),σ为S2的标准差σ=std(S2),mean为取均值函数,std为取标准差的函数;
步骤S62:采用S4=S3/max(abs(S3))的方法将取值范围固定到[-1,1]之间,其中max为取最大值函数,abs为取绝对值函数。
根据上述技术方案,所述步骤S7的具体处理为:
通过将心率HR和心电图数据S4在APP上进行可视化展示,让用户可以观察到自己的心电图波形和实时心率。
与现有技术相比,本发明具有以下有益效果:
本发明基于便携式动态心电监护仪采集的心电数据设计处理方法,用于用户进行长期自我监测的场景,用户测量时可以实时查看自己的心率及心电图波形。其计算的心率和展示的心电图波形相对稳定,抗干扰能力强,能够从一定程度上正确反映用户的心脏健康状态,帮助用户进行自我监测。
本发明针对便携式动态心电监护仪的应用场景,设计了APP端实时心率计算方法和心电图展示方法。不同于常见的利用PPG计算心率,本发明使用用户的心电数据计算心率,可信度更高。在同类型的使用心电数据计算心率的方法中,不同于其他方法直接根据设备采集数据进行计算,本发明在计算心率时提前进行了数据预处理与信号滤波,同时利用缓存池的方法来计算心率,使心率计算更加稳定。在进行心电图数据展示时不同于其他的直接展示方法,本发明进行了正则化处理,数据的可读性更好。
附图说明
图1为基于便携式动态心电监护仪的心率计算方法整体流程图;
图2为基于便携式动态心电监护仪的心率计算方法计算流程流程图;
图3为基于便携式动态心电监护仪的心率计算方法计算心率时不同信号的示例图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例一
请参阅图1-图2,一种基于便携式动态心电监护仪的心率计算方法,包括以下步骤:
步骤S1:将便携式动态心电监护仪贴在人体表面;
步骤S2:数据接收;通过移动终端与便携式动态心电监护仪无线连接,接收便携式动态心电监护仪所检测的心电信号数据,将接收到的数据记为S;
步骤S3:数据预处理;对步骤S2中的数据S进行预处理,去除其中的异常数据点,得到数据S1;
步骤S4:数据滤波处理;通过信号滤波算法对步骤S3中的数据进行滤波处理,去除其中的噪声干扰,得到滤波后的数据S2;
步骤S5:更新缓存,利用缓存数据计算心率;
步骤S6:数据正则化;
步骤S7:心率及心电图波形展示;
步骤S8:重复步骤S2-S7直至测量结束。
所述步骤S3中的具体处理为:
步骤S31:确定S的取值范围为-1mV到1mV:
步骤S32:将S中大于1mV的值置为1mV,小于-1mV的值置为-1mV,得到处理后的数据S1,处理后的S1即为去除掉异常数据点的数据。
所述步骤S4中的具体处理为:
步骤S41:使用0.5Hz-80Hz的带通滤波器进行滤波,去除小于0.5Hz的低频噪声及大于80Hz的高频噪声干扰,得到数据S11;
步骤S42:首先对S11使用1Hz低通滤波器得到S11wander,再使用S11减去S11wander,得到去除了基线漂移影响的S12;
步骤S43:对步骤S42得到的S12使用48Hz-52Hz的带阻滤波器,去除50Hz左右的工频干扰,经过以上滤波处理,得到处理后的数据S2。
所述步骤S5中的具体处理为:
步骤S51:在接收数据的移动终端建立一个数据缓存区;
步骤S52:在数据缓存区建立一个双向队列;双向队列用于记录当前时刻前4秒的所有ECG数据为S-4s;
步骤S53:当有新的数据S2到来时,进行入队列操作,更新缓存数据S-4s,并使用数据S-4s计算心率。
通过数据S-4s计算心率的方法为:
步骤S54:首先使用中的XQRS检测算法检测缓存数据S-4s中的所有R峰位置[r1,r2,...,rn];
步骤S55:根据步骤S54中检测到的所有R峰位置计算所有相邻R峰之间的RR间期,得到RR间期序列RRs=[r2-r1,r3-r2,...,rn-rn-1];其中,每个RR间期的单位为秒;
步骤S56:最后根据RRs计算平均心率HR,计算方法为HR=60/mean(RRs),其中mean为取均值函数。
所述心率的计算步骤具体为:
步骤S57:将缓存数据S-4s通过5Hz-20Hz的带通滤波器进行过滤,得到滤波后的数据S-4s-f;
步骤S58:将数据S-4s-f通过Ricker小波变换进行处理,提取QRS波群信息,得到S-4s-i;
步骤S59:从S-4s-i中寻找所有局部最大值点组成R峰候选项;
步骤S510:设置QRS波群检测阈值相关参数,参数包括:R峰处的幅值应大于0.13mV、该QRS波群的R峰位置距离上一个QRS波群的R峰位置应大于0.2秒;
步骤S511:对每一个R峰候选项按时间顺序进行遍历,根据阈值条件判断其是否属于真正的QRS波群,同时根据当前波形中的R峰幅值与位置,更新步骤S510中的R峰幅值与上一个R峰位置阈值参数;
步骤S512:重复步骤S511直到遍历完所有R峰候选项,得到所有真正的R峰位置;
步骤S513:根据RR间期计算该段数据的平均心率;
步骤S514:输出心率结果。
所述步骤S6的具体处理为:
步骤S61:首先采用S3=(S2-μ)/σ的方法对数据S2进行标准化,其中μ为S2的均值μ=mean(S2),σ为S2的标准差σ=std(S2),mean为取均值函数,std为取标准差的函数;
步骤S62:采用S4=S3/max(abs(S3))的方法将取值范围固定到[-1,1]之间,其中max为取最大值函数,abs为取绝对值函数。
所述步骤S7的具体处理为:
通过将心率HR和心电图数据S4在APP上进行可视化展示,让用户可以观察到自己的心电图波形和实时心率。
实施例二
本发明主要用于对便携式动态心电监护仪采集的心电信号数据进行处理,在手机APP端计算用户的实时心率,并展示实时心电图波形,对用户的心脏健康状态进行反馈。其整体结构如图1所示,其中便携式动态心电监护仪属于数据采集设备,在测量时贴在人体上持续采集用户的心电信号;APP端与便携式动态心电监护仪通过蓝牙连接,持续接收用户心电信号,并利用数字信号处理技术实现实时心率计算,并将用户实时心率及实时心电图波形进行展示。
原理说明:
APP端实时心率计算方法如图1所示,其具体流程如下所示:
1.数据接收。手机APP与便携式动态心电监护仪通过蓝牙连接传输数据,在接收到足够的数据点用于一屏展示时,存储该段数据为S,进入下一步操作。
2.数据预处理。对S进行预处理,去除其中的异常数据点。本方法主要针对过大过小异常数据进行处理。根据正常心电信号的取值范围,及相应便携式动态心电监护仪采集信号的先验知识,确定正常取值范围为-1mV到1mV之间。因此将大于1mV的值置为1mV,小于-1mV的值置为-1mV,减小异常点对后续处理的影响,得到处理后的数据S1。
3.信号滤波算法。该算法主要采用频率域滤波器对ECG信号进行滤波处理,去除其中的噪声干扰。针对ECG信号中可能存在的低频/高频干扰、基线漂移、国内50Hz工频干扰等影响,本步骤使用巴特沃斯滤波器和切比雪夫滤波器两种滤波器来对信号进行处理。主要包括:a)使用0.5Hz-80Hz的巴特沃斯带通滤波器进行滤波,去除低频噪声及高频噪声干扰;b)使用全通滤波器结果减去巴特沃斯1Hz低通滤波器结果,去除基线漂移影响;c)使用切比雪夫48Hz-52Hz的带阻滤波器,去除50Hz左右的工频干扰。经过以上滤波处理,得到处理后的数据S2。
4.更新缓存,利用缓存数据计算心率。由于数据段S2的长度通常为1秒左右(S1、S同理),根据此数据计算的心率很不稳定,因此在缓存区建立一个双向队列,用来记录当前时刻前4秒的所有ECG数据为S-4s,每当由新数据S2到来时,进行入队列操作,更新缓存数据S-4s,之后使用缓存数据S-4s计算心率。
心率计算的具体计算流程如附图2所示。首先使用XQRS检测算法检测缓存数据S-4s中的所有R峰位置;之后根据检测到的所有R峰位置计算RR间期序列RRs;最后根据RRs计算平均心率HR,计算方法为HR=60/mean(RRs),其中mean为取均值函数。
5.数据正则化。主要是对数据S2进行处理用于APP展示,包含两个操作:首先采用S3=(S2-μ)/σ的方法对数据S进行标准化,其中μ为均值,σ为标准差;之后采用S4=S3/max(abs(S3))的方法将取值范围固定到[-1,1]之间,其中max为取最大值函数,abs为取绝对值函数。
6.心率及心电图波形展示。将心率HR和心电图数据S4在APP上进行可视化展示,让用户可以观察到自己的心电图波形和实时心率。
实施例三
本实施例为实施例一的进一步细化。心率计算方面,本发明设置了心电信号缓冲池,用于记录当前时刻前四秒的心电数据,之后利用此段数据计算平均心率,提高心率的稳定性。值得注意的是在该方法中,四秒只是一个可选值,该值的选择是可变的,更大的值意味着心率更稳定,但是实时性更差;更小的值意味着实时性更好,但是波动浮动比较大。与此同时,在进行心电图波形展示时,事先采取了正则化方法,使得图形能够完整地呈现在APP界面中,可读性更强。
需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。
最后应说明的是:以上所述仅为本发明的优选实施例而已,并不用于限制本发明,尽管参照前述实施例对本发明进行了详细的说明,对于本领域的技术人员来说,其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (8)
1.一种基于便携式动态心电监护仪的心率计算方法,其特征在于:包括以下步骤:
步骤S1:将便携式动态心电监护仪贴在人体表面;
步骤S2:数据接收;通过移动终端与便携式动态心电监护仪无线连接,接收便携式动态心电监护仪所检测的心电信号数据,将接收到的数据记为S;
步骤S3:数据预处理;对步骤S2中的数据S进行预处理,去除其中的异常数据点,得到数据S1;
步骤S4:数据滤波处理;通过信号滤波算法对步骤S3中的数据进行滤波处理,去除其中的噪声干扰,得到滤波后的数据S2;
步骤S5:更新缓存,利用缓存数据计算心率;
步骤S6:数据正则化;
步骤S7:心率及心电图波形展示;
步骤S8:重复步骤S2-S7直至测量结束。
2.根据权利要求1所述的一种基于便携式动态心电监护仪的心率计算方法,其特征在于:所述步骤S3中的具体处理为:
步骤S31:确定S的取值范围为-1mV到1mV:
步骤S32:将S中大于1mV的值置为1mV,小于-1mV的值置为-1mV,得到处理后的数据S1,处理后的S1即为去除掉异常数据点的数据。
3.根据权利要求1所述的一种基于便携式动态心电监护仪的心率计算方法,其特征在于:所述步骤S4中的具体处理为:
步骤S41:使用0.5Hz-80Hz的带通滤波器进行滤波,去除小于0.5Hz的低频噪声及大于80Hz的高频噪声干扰,得到数据S11;
步骤S42:首先对S11使用1Hz低通滤波器得到S11wander,再使用S11减去S11wander,得到去除了基线漂移影响的S12;
步骤S43:对步骤S42得到的S12使用48Hz-52Hz的带阻滤波器,去除50Hz左右的工频干扰,经过以上滤波处理,得到处理后的数据S2。
4.根据权利要求1所述的一种基于便携式动态心电监护仪的心率计算方法,其特征在于:所述步骤S5中的具体处理为:
步骤S51:在接收数据的移动终端建立一个数据缓存区;
步骤S52:在数据缓存区建立一个双向队列;双向队列用于记录当前时刻前4秒的所有ECG数据为S-4s;
步骤S53:当有新的数据S2到来时,进行入队列操作,更新缓存数据S-4s,并使用数据S-4s计算心率。
5.根据权利要求4所述的一种基于便携式动态心电监护仪的心率计算方法,其特征在于:通过数据S-4s计算心率的方法为:
步骤S54:首先使用中的XQRS检测算法检测缓存数据S-4s中的所有R峰位置[r1,r2,...,rn];
步骤S55:根据步骤S54中检测到的所有R峰位置计算所有相邻R峰之间的RR间期,得到RR间期序列RRs=[r2-r1,r3-r2,...,rn-rn-1];其中,每个RR间期的单位为秒;
步骤S56:最后根据RRs计算平均心率HR,计算方法为HR=60/mean(RRs),其中mean为取均值函数。
6.根据权利要求4所述的一种基于便携式动态心电监护仪的心率计算方法,其特征在于:所述心率的计算步骤具体为:
步骤S57:将缓存数据S-4s通过5Hz-20Hz的带通滤波器进行过滤,得到滤波后的数据S-4s-f;
步骤S58:将数据S-4s-f通过Ricker小波变换进行处理,提取QRS波群信息,得到S-4s-i;
步骤S59:从S-4s-i中寻找所有局部最大值点组成R峰候选项;
步骤S510:设置QRS波群检测阈值相关参数,参数包括:R峰处的幅值应大于0.13mV、该QRS波群的R峰位置距离上一个QRS波群的R峰位置应大于0.2秒;
步骤S511:对每一个R峰候选项按时间顺序进行遍历,根据阈值条件判断其是否属于真正的QRS波群,同时根据当前波形中的R峰幅值与位置,更新步骤S510中的R峰幅值与上一个R峰位置阈值参数;
步骤S512:重复步骤S511直到遍历完所有R峰候选项,得到所有真正的R峰位置;
步骤S513:根据RR间期计算该段数据的平均心率;
步骤S514:输出心率结果。
7.根据权利要求1所述的一种基于便携式动态心电监护仪的心率计算方法,其特征在于:所述步骤S6的具体处理为:
步骤S61:首先采用S3=(S2-μ)/σ的方法对数据S2进行标准化,其中μ为S2的均值μ=mean(S2),σ为S2的标准差σ=std(S2),mean为取均值函数,std为取标准差的函数;
步骤S62:采用S4=S3/max(abs(S3))的方法将取值范围固定到[-1,1]之间,其中max为取最大值函数,abs为取绝对值函数。
8.根据权利要求1所述的一种基于便携式动态心电监护仪的心率计算方法,其特征在于:所述步骤S7的具体处理为:
通过将心率HR和心电图数据S4在APP上进行可视化展示,让用户可以观察到自己的心电图波形和实时心率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110816029.7A CN113349754A (zh) | 2021-07-20 | 2021-07-20 | 一种基于便携式动态心电监护仪的心率计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110816029.7A CN113349754A (zh) | 2021-07-20 | 2021-07-20 | 一种基于便携式动态心电监护仪的心率计算方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113349754A true CN113349754A (zh) | 2021-09-07 |
Family
ID=77539916
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110816029.7A Pending CN113349754A (zh) | 2021-07-20 | 2021-07-20 | 一种基于便携式动态心电监护仪的心率计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113349754A (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105342601A (zh) * | 2015-12-21 | 2016-02-24 | 重庆大学 | 一种家用心电检测仪 |
CN105982661A (zh) * | 2015-02-02 | 2016-10-05 | 四川理工学院 | 一种便携式心电监护系统及其数据处理方法 |
CN106991291A (zh) * | 2017-04-14 | 2017-07-28 | 北京工业大学 | 胎儿心电信号实时转换为胎儿胎音的计算机仿真转换方法 |
CN107361764A (zh) * | 2017-06-16 | 2017-11-21 | 华南理工大学 | 一种心电信号特征波形r波的快速提取方法 |
CN107616782A (zh) * | 2017-09-06 | 2018-01-23 | 郑州云海信息技术有限公司 | 一种心电信号质量检测方法与装置 |
US20210106240A1 (en) * | 2019-10-07 | 2021-04-15 | Owlet Baby Care, Inc. | Fetal Heart Rate Prediction from Electrocardiogram |
-
2021
- 2021-07-20 CN CN202110816029.7A patent/CN113349754A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105982661A (zh) * | 2015-02-02 | 2016-10-05 | 四川理工学院 | 一种便携式心电监护系统及其数据处理方法 |
CN105342601A (zh) * | 2015-12-21 | 2016-02-24 | 重庆大学 | 一种家用心电检测仪 |
CN106991291A (zh) * | 2017-04-14 | 2017-07-28 | 北京工业大学 | 胎儿心电信号实时转换为胎儿胎音的计算机仿真转换方法 |
CN107361764A (zh) * | 2017-06-16 | 2017-11-21 | 华南理工大学 | 一种心电信号特征波形r波的快速提取方法 |
CN107616782A (zh) * | 2017-09-06 | 2018-01-23 | 郑州云海信息技术有限公司 | 一种心电信号质量检测方法与装置 |
US20210106240A1 (en) * | 2019-10-07 | 2021-04-15 | Owlet Baby Care, Inc. | Fetal Heart Rate Prediction from Electrocardiogram |
Non-Patent Citations (2)
Title |
---|
王余涛: "基于嵌入式系统的便携式心电监护系统的研制", 《中国优秀硕士学位论文全文数据库 (工程科技Ⅱ辑)》 * |
龚著琳等: "《生物医学数据挖掘》", 31 August 2011, 上海科学技术出版社 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
AU2016293110B2 (en) | Processing biological data | |
US11234658B2 (en) | Photoplethysmogram data analysis and presentation | |
US4938228A (en) | Wrist worn heart rate monitor | |
CN100413461C (zh) | 一种获取呼吸暂停事件和睡眠结构图信息的方法 | |
US8172764B2 (en) | Method and system for high-resolution extraction of quasi-periodic signals | |
EP3773185A1 (en) | Photoplethysmogram data analysis and presentation | |
CN109222961A (zh) | 一种便携式睡眠监测系统及相关睡眠监测方法 | |
US20140148714A1 (en) | Automatic online delineation of a multi-lead electrocardiogram bio signal | |
CN106539580B (zh) | 一种自主神经系统动态变化的连续监测方法 | |
CN110507317B (zh) | 一种心电信号r波的自适应ca-cfar定位方法 | |
Paliakaitė et al. | Modeling of artifacts in the wrist photoplethysmogram: Application to the detection of life-threatening arrhythmias | |
CN111297340A (zh) | 基于gps与心率结合的运动状态监测方法 | |
Bashar et al. | Smartwatch based atrial fibrillation detection from photoplethysmography signals | |
Zhang et al. | Hear the heart: Daily cardiac health monitoring using Ear-ECG and machine learning | |
De Giovanni et al. | Ultra-low power estimation of heart rate under physical activity using a wearable photoplethysmographic system | |
KR20160008368A (ko) | 재귀 추정을 이용한 생체 신호 처리 방법 및 그 장치 | |
Modak et al. | A novel adaptive multilevel thresholding based algorithm for QRS detection | |
CN113349754A (zh) | 一种基于便携式动态心电监护仪的心率计算方法 | |
Tun et al. | Analysis of computer aided identification system for ECG characteristic points | |
CN113349753A (zh) | 一种基于便携式动态心电监护仪的心律失常检测方法 | |
Wu et al. | Development of full-featured ECG system for visual stress induced heart rate variability (HRV) assessment | |
Kouekeu et al. | Embedded QRS complex detection based on ECG signal strength and trend | |
Fallet et al. | An adaptive organization index to characterize atrial fibrillation using wrist-type photoplethysmographic signals | |
Karnewar et al. | Preprocessing ECG signal by eliminating various noises using denoising methods | |
Jahan et al. | A low cost stethoscopic system for real time auscultation of heart sound |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210907 |
|
RJ01 | Rejection of invention patent application after publication |