CN106108885B - 评估身体运动状态及心电干扰大小的方法及其仪器 - Google Patents
评估身体运动状态及心电干扰大小的方法及其仪器 Download PDFInfo
- Publication number
- CN106108885B CN106108885B CN201610420036.4A CN201610420036A CN106108885B CN 106108885 B CN106108885 B CN 106108885B CN 201610420036 A CN201610420036 A CN 201610420036A CN 106108885 B CN106108885 B CN 106108885B
- Authority
- CN
- China
- Prior art keywords
- electrocardio
- characteristic point
- data
- electrocardiogram
- section
- 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.)
- Active
Links
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/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]
-
- 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
- A61B5/7207—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal of noise induced by motion artifacts
- A61B5/721—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal of noise induced by motion artifacts using a separate sensor to detect motion or using motion information derived from signals other than the physiological signal to be measured
-
- 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
- A61B5/7217—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal of noise originating from a therapeutic or surgical apparatus, e.g. from a pacemaker
Abstract
本发明公开了一种评估身体运动状态及心电干扰大小的方法及其仪器。本发明采用了一种特殊的数据处理技术和算法,从中提取出活动因子,用来评估身体运动状态及心电干扰大小,使得仪器在记录动态心电图的同时,可以监测身体活动状态,拓展动态心电图的应用新领域。
Description
技术领域
本发明涉及心电数据处理算法技术,具体涉及一种采用心电数据来评估身体运动状态及心电干扰大小的方法及其仪器。
背景技术
心脏周围的组织和体液都能导电,因此可将人体看成为一个具有长、宽、厚三度空间的容积导体。心脏好比电源,无数心肌细胞动作电位变化的总和可以传导并反映到体表。在体表很多点之间存在着电位差,也有很多点彼此之间无电位差是等电的。
心电图是反映心脏兴奋的电活动过程,心电图上所记录的电位变化是一系列瞬间心电综合向量在不同导联轴上的反映,也就是平面向量环在有关导联轴上的再投影。心电图代表整个心脏电激动的综合过程,以一个个心肌细胞的电激动为基础,心肌激动时细胞内发生电传变化。
动态心电图是一种可以长时间连续记录并分析人体心脏在活动和安静状态下心电图变化的方法。
动态心电图技术于1957年由Holter首先应用于监测心脏电活动的研究,所以又称Holter监测心电图仪,目前已成为临床心血管领域中非创伤性检查的重要诊断方法之一。与普通心电图相比,动态心电图于24小时内可连续记录多达10万次左右的心电信号,这样可以提高对非持续性心律失常,尤其是对一过性心律失常及短暂的心肌缺血发作的检出率,因此扩大了心电图临床运用的范围。
心电图是在体表记录心脏的电活动,心电数据受到很多因素的干扰。干扰因素包括心电记录系统因素和非心电记录系统因素。其中,非心电记录系统因素主要有电极因素(电极与皮肤接触不良、放置部位不当、电极质量差等)、肌肉电因素(活动产生的肌肉电干扰)、其他因素(静电、电磁干扰)。
动态心电图由于受试者为正常的生活状态,身体定有不同程度的运动,使得所记录的运动心电信号中可能混入不同程度的干扰信号,包括呼吸引起的心电信号基线漂移,电极与皮肤间的摩擦噪声,体位移动所引发的扰动,电源噪声以及其它电磁噪声等,这些干扰信号的消除是不可能完全依靠数字信号处理技术来完成的。
到目前为止,由于身体活动带来心电图的变化都作为一种信号干扰来处理。
发明内容
针对心电图测量中身体活动对信号的干扰存在必然性,以及动态心电图机仅作为心脏检测的单一功能的问题,本发明的目的在于提供一种采用心电数据来评估身体运动状态及心电干扰大小的方法,并提供一种可以监测运动状态及其可记录动态心电图的仪器,拓展动态心电图仪的应用功能。
为了达到上述目的,本发明采用如下的技术方案:
一种评估身体运动状态及心电干扰大小的方法,其特征在于,
该方法具体包括如下步骤:
A、选取一段心电数据;
B、拾取心电波形的峰和谷点;
C、计算每个心电周期的特征点参数,包括相临R波的时间间隔RR、P波和R波的时间间隔PR、S波和T波的时间间隔ST、Q波和T波的时间间隔QT、R波峰面积ER、P波峰面积EP、T波峰面积ET;
D、统计特征点参数的相对误差;包括该段心电数据内所有RR特征点参数的相对误差RRrsd、该段心电数据内所有PR特征点参数的相对误差PRrsd、该段心电数据内所有ST特征点参数的相对误差STrsd、该段心电数据内所有QT特征点参数的相对误差QTrsd以及该段心电数据内所有ER特征点参数的相对误差ERrsd、该段心电数据内所有EP特征点参数的相对误差EPrsd、该段心电数据内所有ET特征点参数的相对误差ETrsd;
E、计算出活动因子:对该段心电数据内所有心电周期的特征点参数的相对误差值,使用加权平均的方法,计算出活动因子;
F、显示心电综合干扰程度和判断身体活动状态:通过活动因子来判断身体是处于安静状态还是和运动状态,活动阈值为通过试验确定的一个经验参数,当活动因子小于活动阈值,表示身体处于安静状态,反之表示身体处于运动状态;同时使用活动因子表示心电综合干扰程度并显示。
优选的,在步骤A中,心电数据的长度最好选在3到10之间。
优选的,在步骤F中,活动阈值为8-12。
通过上述的算法可知,本发明通过在心电数据数据分析中采用了一种特殊的数据处理技术和算法,引入活动因子和活动阈值概念,实现实时评估身体运动状态及心电干扰大小,使得仪器在记录动态心电图的同时,可以监测身体活动状态和反映心电信号的干扰程度,提高心电信号分析的可信度。
一种可以监测运动状态及其记录动态心电图的仪器,其特征在于,包括心电电极、心电采集盒、数据记录和处理设备;其中,心电图电极包括两个可导电的接触电极;所述心电电极通过导线与心电采集盒相连;所述心电采集盒和接收端相连;心电采集盒采集接触电极的心电电信号,并把数据发送到数据记录和处理设备;数据记录和处理设备可以是装有本发明算法的心电数据记录和处理软件的手机、计算机、平板电脑。
本发明的方法采用了一种特殊的数据处理技术和算法,从中提取出活动因子,用来评估身体运动状态及心电干扰大小,使得仪器在记录动态心电图的同时,可以监测身体活动状态,拓展动态心电图的应用新领域。
本发明提供的数据处理技术具有如下优点:
(1)动态心电仪能够实现自动记录身体的活动状态,开辟了动态心电图新的应用领域,可以实现用心电仪来进行被测试者的活动状态分析和睡眠质量的分析;
(2)心电仪能够实时得到心电综合干扰程度,用于判断心电数据的可靠性,以便正确地分析心电图;
(3)“活动因子”的提出,为心电数据的误差分析提供了一个新思路,为心电周期特征增加了一个新参数。
(4)本发明实现了仅用心电数据来监测心率和运动状态。
附图说明
以下结合附图和具体实施方式来进一步说明本发明。
图1为可以监测运动状态及其记录动态心电图的仪器的原理图;
图2为活动因子计算方法及身体活动状态判断流程图;
图3为心电周期的特征点示意图;
图4为样本数据心电图;
图5为静息状态实际测量心电图;
图6为活动状态实际测量心电图。
具体实施方式
为了使本发明实现的技术手段、创作特征、达成目的与功效易于明白了解,下面结合具体图示,进一步阐述本发明。
参见图1,其所示为可以监测运动状态及其记录动态心电图的仪器的原理图。由图可知,该仪器100主要包括心电电极101、心电采集盒102、数据记录和处理设备103三部分。心电电极101为心电传感器,包括两个可导电的接触电极;心电采集盒102由数据采集模块104和蓝牙通讯模块105组成,实现心电电信号采集,并通过蓝牙把数据发送到数据记录和处理设备;心电电极101和心电采集盒102之间由两根导电线连接;数据记录和处理设备103为心电数据记录和处理软件的手机、计算机或平板电脑;
参见图2,其所示为本发明的活动因子计算方法及身体活动状态判断流程图。由图可知,本发明算法步骤如下:
A、选取一段心电数据,如图4所示;
B、拾取心电波形的峰和谷点:寻峰参数可使用经验参数,保证峰和谷为有效的心电周期的特征点,特征点如图3所示;心电周期的特征点包括P波点、Q波点、R波点、S波点、T波点、U波点、tg点、ug点;
C、计算每个心电周期的特征点参数,包括相临R波的时间间隔RR、P波和R波的时间间隔PR、S波和T波的时间间隔ST、Q波和T波的时间间隔QT、R波峰面积ER、P波峰面积EP、T波峰面积ET;
D、统计特征点参数的相对误差;包括该段心电数据内所有RR特征点参数的相对误差RRrsd、该段心电数据内所有PR特征点参数的相对误差PRrsd、该段心电数据内所有ST特征点参数的相对误差STrsd、该段心电数据内所有QT特征点参数的相对误差QTrsd以及该段心电数据内所有ER特征点参数的相对误差ERrsd、该段心电数据内所有EP特征点参数的相对误差EPrsd、该段心电数据内所有ET特征点参数的相对误差ETrsd;
E、计算出活动因子:对该段心电数据内所有心电周期的特征点参数的相对误差值,使用加权平均的方法,计算出活动因子;公式为:
f = RRrsd *Krr + PRrsd *Kpr+ STrsd *Kst+ QTrsd *Kqt+ ERrsd *Ker +EPrsd *Kep+ ETrsd *Ket。
其中RRrsd、PRrsd、STrsd、QTrsd、ERrsd、EPrsd、ETrsd为特征点参数的相对误差;Krr、Kpr、Kst、Kqt、Ker、Kep、Ket为对应特征点参数相对误差权值,并且Krr+Kpr+Kst+Kqt+Ker+Kep+Ket = 1。
例如:
f = RRrsd *0.2 + PRrsd *0.1+ STrsd *0.2+ QTrsd *0.1+ ERrsd *0.15 +EPrsd *0.1+ ETrsd *0.15 。
参见图5,其所示为静息状态实际测量心电图。由图可知,“活动因子”f为3 。
参见图6,其所示为活动状态实际测量心电图。由图可知,“活动因子”f为15 。
F、活动阈值的选定:活动阈值用来划分身体是处于安静状态还是和运动状态。活动阈值为通过试验确定的一个经验参数,一般为10。
G、显示心电综合干扰程度和判断身体活动状态:通过活动因子来判断身体是处于安静状态还是和运动状态,当活动因子小于活动阈值,表示身体处于安静状态,反之表示身体处于运动状态;同时使用活动因子表示心电综合干扰程度并显示。取活动阈值为9,图5为活动因子f为3,小于活动阈值,运动状态可定为安静;图6 为活动因子f为15,大于活动阈值,运动状态可定为活动中。
其中,算法第一步所述的:选取一段心电数据,如图4所示,为样本数据心电图,这里样本数据长度为3秒,本方法选取一段心电数据的长度最好选在3到10之间,一方面数据量不至于过大而影响计算速度,一方面数据也能满足统计要求。
以图4的样本数据心电图为例,该段心电数据长度3秒,有三个心电周期,心电特征点参数有3组,各特征点参数的相对误差就是三个心电周期中各特征点参数和三个参数平均值之间的误差,再除以参数平均值就得到该参数的相对误差。
以上显示和描述了本发明的基本原理、主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。
Claims (4)
1.一种评估身体运动状态及心电干扰大小的方法,其特征在于,
该方法具体包括如下步骤:
A、选取一段心电数据;
B、拾取心电波形的峰和谷点;
C、计算每个心电周期的特征点参数,包括相临R波的时间间隔RR、P波和R波的时间间隔PR、S波和T波的时间间隔ST、Q波和T波的时间间隔QT、R波峰面积ER、P波峰面积EP、T波峰面积ET;
D、统计特征点参数的相对误差;包括该段心电数据内所有RR特征点参数的相对误差RRrsd、该段心电数据内所有PR特征点参数的相对误差PRrsd、该段心电数据内所有ST特征点参数的相对误差STrsd、该段心电数据内所有QT特征点参数的相对误差QTrsd以及该段心电数据内所有ER特征点参数的相对误差ERrsd、该段心电数据内所有EP特征点参数的相对误差EPrsd、该段心电数据内所有ET特征点参数的相对误差ETrsd;
特征点参数的相对误差的计算方法是心电周期中各特征点参数和参数平均值之间的误差,再除以参数平均值就得到该参数的相对误差;
E、计算出活动因子:对该段心电数据内所有心电周期的特征点参数的相对误差值,使用加权平均的方法,计算出活动因子;
F、显示心电综合干扰程度和判断身体活动状态:通过活动因子来判断身体是处于安静状态还是和运动状态,活动阈值为通过试验确定的一个经验参数,当活动因子小于活动阈值,表示身体处于安静状态,反之表示身体处于运动状态;同时使用活动因子表示心电综合干扰程度并显示。
2.根据权利要求1所述的评估身体运动状态及心电干扰大小的方法,其特征在于,在步骤A中,心电数据的长度在3到10之间。
3.根据权利要求1所述的评估身体运动状态及心电干扰大小的方法,其特征在于,在步骤F中,活动阈值为8-12。
4.一种可以监测运动状态及其记录动态心电图的仪器,其特征在于,包括心电电极、心电采集盒、数据记录和处理设备;其中,心电图电极包括两个可导电的接触电极;所述心电电极通过导线与心电采集盒相连;所述心电采集盒和接收端相连;心电采集盒采集接触电极的心电电信号,并把数据发送到数据记录和处理设备;数据记录和处理设备是采用如权利要求1-3任一项所述方法的心电数据记录和处理软件的手机、计算机、平板电脑。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610420036.4A CN106108885B (zh) | 2016-06-14 | 2016-06-14 | 评估身体运动状态及心电干扰大小的方法及其仪器 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610420036.4A CN106108885B (zh) | 2016-06-14 | 2016-06-14 | 评估身体运动状态及心电干扰大小的方法及其仪器 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106108885A CN106108885A (zh) | 2016-11-16 |
CN106108885B true CN106108885B (zh) | 2021-07-23 |
Family
ID=57270043
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610420036.4A Active CN106108885B (zh) | 2016-06-14 | 2016-06-14 | 评估身体运动状态及心电干扰大小的方法及其仪器 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106108885B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108198623A (zh) * | 2017-12-15 | 2018-06-22 | 东软集团股份有限公司 | 人体身体状态检测方法、装置、存储介质及电子设备 |
CN109512422A (zh) * | 2018-09-25 | 2019-03-26 | 维灵(杭州)信息技术有限公司 | 一种ecg波形混乱程度计算方法 |
CN111166354B (zh) * | 2020-01-23 | 2022-11-18 | 北京津发科技股份有限公司 | 影响情绪变化的因素的分析方法及电子设备 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100493714B1 (ko) * | 2004-11-26 | 2005-06-02 | 주식회사 메디코아 | 자율신경 검사장치 |
CN101061949A (zh) * | 2006-04-28 | 2007-10-31 | 万威科研有限公司 | 运动数据设备 |
CN101822539A (zh) * | 2009-03-02 | 2010-09-08 | 王建成 | 一种消除心电图运动干扰的方法 |
CN102397067A (zh) * | 2011-12-07 | 2012-04-04 | 秦皇岛市康泰医学系统有限公司 | 一种运动负荷心电波形叠加分析方法及其系统 |
CN103099615A (zh) * | 2013-01-23 | 2013-05-15 | 深圳市理邦精密仪器股份有限公司 | 一种消除运动心电信号干扰的方法和装置 |
US8513945B2 (en) * | 2009-04-03 | 2013-08-20 | New York University | System, method and computer-accessible medium for providing breath-hold multi-echo fast spin-echo pulse sequence for accurate R2 measurement |
CN103330550A (zh) * | 2013-03-04 | 2013-10-02 | 北京中医药大学 | Mems液压传触的脉诊三部九侯信息自动获取和识别装置及方法 |
CN104027109A (zh) * | 2013-03-07 | 2014-09-10 | 精工爱普生株式会社 | 心房颤动解析装置以及程序 |
CN104188652A (zh) * | 2014-09-09 | 2014-12-10 | 广东工业大学 | 心电图的数据质量实时控制方法及其系统 |
CN111956210A (zh) * | 2020-08-27 | 2020-11-20 | 重庆邮电大学 | 一种基于eemd和信号结构分析的心电信号r波识别方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2532450B (en) * | 2014-11-19 | 2019-05-15 | Suunto Oy | Wearable sports monitoring equipment with context determination capabilities and relating method |
-
2016
- 2016-06-14 CN CN201610420036.4A patent/CN106108885B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100493714B1 (ko) * | 2004-11-26 | 2005-06-02 | 주식회사 메디코아 | 자율신경 검사장치 |
CN101061949A (zh) * | 2006-04-28 | 2007-10-31 | 万威科研有限公司 | 运动数据设备 |
CN101822539A (zh) * | 2009-03-02 | 2010-09-08 | 王建成 | 一种消除心电图运动干扰的方法 |
US8513945B2 (en) * | 2009-04-03 | 2013-08-20 | New York University | System, method and computer-accessible medium for providing breath-hold multi-echo fast spin-echo pulse sequence for accurate R2 measurement |
CN102397067A (zh) * | 2011-12-07 | 2012-04-04 | 秦皇岛市康泰医学系统有限公司 | 一种运动负荷心电波形叠加分析方法及其系统 |
CN103099615A (zh) * | 2013-01-23 | 2013-05-15 | 深圳市理邦精密仪器股份有限公司 | 一种消除运动心电信号干扰的方法和装置 |
CN103330550A (zh) * | 2013-03-04 | 2013-10-02 | 北京中医药大学 | Mems液压传触的脉诊三部九侯信息自动获取和识别装置及方法 |
CN104027109A (zh) * | 2013-03-07 | 2014-09-10 | 精工爱普生株式会社 | 心房颤动解析装置以及程序 |
CN104188652A (zh) * | 2014-09-09 | 2014-12-10 | 广东工业大学 | 心电图的数据质量实时控制方法及其系统 |
CN111956210A (zh) * | 2020-08-27 | 2020-11-20 | 重庆邮电大学 | 一种基于eemd和信号结构分析的心电信号r波识别方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106108885A (zh) | 2016-11-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Satija et al. | Real-time signal quality-aware ECG telemetry system for IoT-based health care monitoring | |
Taji et al. | Impact of skin–electrode interface on electrocardiogram measurements using conductive textile electrodes | |
Phinyomark et al. | Fractal analysis features for weak and single-channel upper-limb EMG signals | |
CN102462494B (zh) | 一种新型智能心电测试保健器 | |
CN110831494A (zh) | 对来自便携式传感器装置的心音图数据和心电图数据进行分析 | |
CN106108885B (zh) | 评估身体运动状态及心电干扰大小的方法及其仪器 | |
JP2005512720A (ja) | 医療の電極システムと方法 | |
KR20120003662A (ko) | 생체 신호를 측정하는 장치 및 방법 | |
Khan et al. | Detection of myocardial infarction using pulse plethysmograph signals | |
CN108209912A (zh) | 一种肌电信号采集方法及装置 | |
Issac et al. | CUEDETA: A real time heart monitoring system using android smartphone | |
Kuzmin et al. | Device and software for mobile heart monitoring | |
CN113160921A (zh) | 一种基于血流动力学的数字人体心血管系统的构建方法和应用 | |
Burden | Surface electromyography | |
Popović-Maneski et al. | Properties of different types of dry electrodes for wearable smart monitoring devices | |
Grönlund et al. | On-line signal quality estimation of multichannel surface electromyograms | |
Fraser et al. | Biosignal quality analysis of surface EMG using a correlation coefficient test for normality | |
Zhang et al. | An improved real-time R-wave detection efficient algorithm in exercise ECG signal analysis | |
Adi et al. | ECG-LPWAN based for Real-time monitoring Patient’s Heart Beat Status | |
CN105232041A (zh) | 人体体表电特性分布的测量装置、数据终端及方法 | |
Gan et al. | A new method to reduce motion artifact in electrocardiogram based on an innovative skin-electrode impedance model | |
Plesnik et al. | ECG signal acquisition and analysis for telemonitoring | |
Alcaraz et al. | Application of non-linear methods in the study of atrial fibrillation organization | |
Anas et al. | On-line monitoring and analysis of bioelectrical signals | |
Kontunen et al. | Low-latency EMG onset and termination detection for facial pacing |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |