CN102961129A - 一种远程医疗的异常心电张量分析方法 - Google Patents
一种远程医疗的异常心电张量分析方法 Download PDFInfo
- Publication number
- CN102961129A CN102961129A CN2012104169310A CN201210416931A CN102961129A CN 102961129 A CN102961129 A CN 102961129A CN 2012104169310 A CN2012104169310 A CN 2012104169310A CN 201210416931 A CN201210416931 A CN 201210416931A CN 102961129 A CN102961129 A CN 102961129A
- Authority
- CN
- China
- Prior art keywords
- tensor
- sigma
- data
- electrocardiogram
- ripple
- 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.)
- Granted
Links
Images
Landscapes
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明公布了一种远程医疗的异常心电张量分析方法,首先通过远程方式采集大量的标准12导联心电数据,然后通过短时傅立叶变换将心电转换为高维度的张量心电数据。然后以高维张量心电数据直接作为特征,使用直接以张量数据直接作为输入的特征抽取和特征降维的算法提取出直接用来分类的心电特征。由于这种方法是基于TTV变换法则的,所以最终可以得到基于向量存储的特征,然后使用SVM分类方法对这些向量特征进行分类。这种方法以张量心电数据直接作为输入,充分利用了心电的多导联心电的结构信息,消除了原先单导联心电单独分析带来的不精准缺陷,实现了心电分析的有效性。
Description
技术领域
本发明涉及的是一种针对远程心电诊断平台,多导联张量心电数据的特征抽取和辅助分类方法。
背景技术
本系统基于物联网技术,利用远程心电诊断平台,将上海市三级甲等医院的优质医疗服务延伸至远程的县市级甚至社区级医疗中心。设计了个性化服务系统,提高了远程医疗与健康监护平台的可用性及使用效率。利用基于物联网的构建远程医疗云服务平台,并建立基层医疗单位的远程心电诊断服务示范。构建远程医疗云服务主要包括医疗信息采集、传输、数据处理、反馈诊断等关键技术。当前心电信号采集设备和医疗数据远程传输技术已经成熟,本项目核心技术主要集中在医疗数据挖掘与诊断辅助决策支持平台方面,使得能为医疗中心医生实现快速准确的诊断支持。该医疗服务平台逐步实现基层医疗单位、家庭的两种远程医疗与健康监护服务模式,形成物联网下的新型远程医疗与健康监护医疗服务模式。该项目研发医疗技术云服务平台主要创新之处包括三个方面:远程医疗诊断服务模式创新、医疗数据挖掘技术平台创新和医疗数据服务创新。
基于远程心电云服务平台,通过实现心电分析的辅助决策支持,极大提高了诊断中心医生的诊断速度和精度,使得远程心电诊断新型服务模式成为可能。在远程心电智能诊断技术方面,主要在以下几个方面提供辅助决策支持:通过对心电数据分析,提供准确的心电分析基本参数;采用结合机器学习方法和智能模式识别方法,通过对心电信号的有效表征,利用最新模式识别方法(如核方法)对心电特征进行模式分类,实现对心电数据的预诊,云服务平台能即时推荐给熟悉相应病症的医生进行诊断,提高了诊断的效率。本系统通过网格化的数据存储很好的为非结构化的海量数据提供高效存储管理和查询服务功能。这样一套完备的分析数据库,为以后的医学和工程研究提供一个大型的完备数据库
在科学计算领域,基于心电数据进行特征抽取和降维并且进行分类分析的方法比较有名的就是C. Saritha, V. Sukanya, and Y. Narasimha Murthy, “ECG Signal Analysis Using WaveletTransforms,” Bulgarian Journal of Physics, vol. 35,pp. 68-77, 2008. (作者:C. Saritha, V. Sukanya, Y.Narasimha Murthy, 题目:基于小波变换的心电信号分析,杂志:保加利亚物理杂志,2008年35卷,68-77页) 和Kuo-Kuang Jen,and Yean-Ren Hwang, “ECGFeature Extraction and Classification Using spectrum and Neural Networks,” Journal of Medical and Biological Engineering, vol. 28, no. 1, 2008. (作者:Kuo-Kuang Jen, Yean-Ren Hwang, 题目:基于频谱和神经网络的心电特征抽取和分类,杂志:医学和生物工程杂志,2008年28卷,第一期)对于第一篇文章,首先对信号进行不同小波基底的小波分解,然后通过检测其系数来提取有效的针对不同疾病波形的特征进行分析。第二篇文章是将原来的心电转换成频谱,然后通过提取心电频谱上的特征来抽取有效的针对不同疾病的特征。最后使用神经网络来对原始的心电进行分类。需要指出的是,所有这些方法都是针对单导联心电的,而实际使用的心电诊断数据往往是多导联的。一般目前在医学上临床使用的心电是12标准导联的。这样的话针对单一导联的算法就不能完全满足需要,即使通过将多个导联心电拼接成单个导联的方法,多个导联的结构信息被破坏,有用信息有所减少并不是非常有效的方法。因此,需要研发新的将多导联心电张量心电数据作为直接输入的特征抽取和分析方法,来处理上述论文中方法无法直接处理的问题。
发明内容
本发明的目的在于克服现有分散式的心电诊断系统的非集中化带来的效率低下的不足,同时又针对原先基于单导联心电分析不够准确的缺陷,提出了一种针对多导联张量心电数据的心电分析方法。这种方法以张量心电数据直接作为输入,充分利用了心电的多导联心电的结构信息,消除了原先单导联心电单独分析带来的不精准缺陷,实现了心电分析的有效性。
本发明是通过以下技术方案实现的,其具体步骤如下:
一种远程医疗的异常心电张量分析方法,包括下述步骤:
(1)构造张量数据:
a.心电数据采集:
采集标准12导联心电数据;
b.数据预处理和去燥:
对采集的12导联心电数据首先对信号通过50hz的陷波滤波器进行滤波处理,然后对数据进行DB6小波的小波变换分解,然后去除其中最高频的和最低频的信号成份;
c.波形检测:
再次对小波进行DB6小波分解,然后找寻其中level2小波系数,然后采用过零点检测方法检测心电的QRS波尖峰,然后依次去除R波后检测P波和T波尖峰,然后采用detrend算法计算出波形的基线,通过计算各个波形与基线的交点,确定P波QRS波T波的起始和结束,也就是onset和offset;
d.逐跳切割:
一次采集心电大约包含20秒的心电数据,也就是25跳左右的心电,对于心电一次beat也就是一个P波一个QRS波一个T波进行切割;
e. R波对齐
对每一个beat的心电针对R波的峰值进行对齐,并且切割成统一的长度;
f. 短时fourier变换
为了有效抽取时频域的心电诊断特征,通过使用窗口为128点的短时fourier变换对心电进行转换,最后心电被转换为128×128×12的时频空的3阶张量;这里的空指的是导联位置就是指导联轴;对于12导联(lead×time)ECG信号,s[l,n]表示离散信号在时间点n对于导联l;在时间点 nΔt 和 频率 f 的短时fourier变换如下式:
这里 w[n] 是一个窗函数,变换完成后数据变成一个3阶张量;
(2)基于TTV映射的张量特征抽取:
a. 根据原始张量数据的判别性计算有效的投影张量:
b. 根据优化张量可分性特点计算更为优越的投影张量:
c. 对原始数据去除已抽取投影张量的维度,调整原始数据的结构:
d. 对原始数据计算投影张量进行加权处理:
通过调整每个张量的权重,来规避不合理张量的影响,等式如下:
最简单的方式就是取张量间距离的倒数distance (),
w(dxy)=dxy -n
或者如下式定义:
组合两种形式得到:
e. 算法整体迭代收敛过程如下:
整个算法过程,是一个严格的单调收敛过程,逐次迭代目标值关系如下式所示:
T->无穷时,算法收敛于最终目标极限收敛值;
f. 判断计算过程结束终止条件:
使用如下方法来判断算法是否收敛,并且判断算法什么时候应该停止;误差值小于一定阈值则判断算法停止:
||Fk-Fk-1||Fro≤ε
采用这种方法来判断收敛与否和终止整个算法;
(3)选择合理初值:
求解一个最为近似的张量:
a.无约束情况张量算法:
交替最小二乘方法求解目标等式如下:
其展开等式如下式:
张量梯度下降方法:
可以将上述目标函数进行展开,写成如下的形式:
第一项没有涉及变量,所以:
这里 0 表示一个 0 向量,长度为 In ,第二项犹如如下式子:
第二项求导之后得到如下式:
第三项如下:
因此
综合以上三项就可以得到最终结果;
b. 有约束张量情况:
带约束非线性最小二乘:
带约束优化问题,通过转换对求解的约束条件有所放松,再去求解就可以得到更加优越的计算结果,如下式:
将左边的优化问题转换为右边的优化问题,然后将结果代入,使用标准SQP算法求解,来求得最终的优化结果;
(4)分类比对:
最后使用SVM来对抽取出来的以vector方式存储的向量特征进行分类,
这直接通过求解如下优化主问题:
Subject to yi(wTφ(xi)+b)≥1-ξn,ξn≥0,i=1,2…,n
这里的参数 C>0 在松弛变量和惩罚因子间的平衡,而他的Lagrangian乘子变换成如下等式:
这里 {an≥0} 而且 { μn≥0} 是Lagrangian乘子, 而对偶Lagrangian 问题如下式:
本发明的有益效果是:
(1)充分利用远程心电集中式的优点,充分利用云技术和物联网技术。有效提升心电诊断的高度集中化统一调度诊断。使得处理速度非常快,计算效率高,可以做到实时处理。
(2)系统算法针对12导联医用标准数据库,比起原始但导联的算法具有更大的优势。
(3)提出直接针对张量心电的特征抽取和特征降维,使得高阶张量心电可以保留完整的结构信息不丢失,使诊断分类可以更加有效的进行。
(4)采用ECG多导联信号的时频域特性作为诊断系数,比起单纯时域的特征具有更好的性能。
附图说明
图 1 为本发明硬件架构图。
图 2 为本发明数据流程图。
图 3 为本发明方法实现流程图。
图 4 为本发明GTR1DA方法进行特征抽取后降维到3维和的分类效果图。
图 5 为ica方法进行特征抽取后降维到3维和的分类效果图。
图 6 为lda方法进行特征抽取后降维到3维和的分类效果图。
图 7 为pca方法进行特征抽取后降维到3维和的分类效果图。
图 8 为TR1DA方法进行特征抽取后降维到3维和的分类效果图。
图 9 为UMPCA方法进行特征抽取后降维到3维和的分类效果图。
具体实施方式
我们首先通过远程方式采集大量的标准12导联心电数据,然后通过短时傅立叶变换(Short-time Fourier Transform)将心电转换为高维度的张量心电数据(通常128×128×12)。然后以高维张量心电数据直接作为特征,使用直接以张量数据直接作为输入的特征抽取和特征降维的Generalized Tensor Rank One Discriminant Analysis算法提取出直接用来分类的心电特征。由于这种方法是基于TTV变换法则的,所以最终可以得到基于向量存储的特征,然后使用SVM分类方法对这些向量特征进行分类。提出一种以张量作为直接输入的特征降维和特征抽取算法来直接对于张量的心电数据进行处理是本发明的核心创新点。
为了更好地描述本发明的内容,首先描述一下远程心电系统的架构。
心电的远程医疗系统主要实现心电信号的远程采集、传输、诊断功能,框架上结合了物联网和分布式的优势特点,总共分为采集端、服务器端、诊断端三大模块,如图1所示。
患有心脏疾病的患者不需要前往大型医院排队,而是由社区医院的医生携带轻便的采集设备前往居民区,在现场对病人进行心电数据的采集;采集结束后,采集医生通过蓝牙、wifi等方式将心电数据从采集设备传输至数据接收前置机,该前置机可以是手机、平板电脑或PC等;前置机将心电数据和病人基本信息按照特定的网络传输协议进行打包,传输至服务器端的数据格式转换服务器。
数据格式转换服务器在接收到数据包之后,按照约定的传输协议进行解包,生成格式化的心电数据类,发送至数据库服务器,将心电数据保存至数据仓库,其中,由于各采集中心的配套设备存在差异,因此打包传输的网络传输协议可能有所不同,该数据格式转换服务器的作用及时形成统一接口,产生统一的数据类结构;数据库服务器对数据进行保存和判断,如果是新采集的心电数据,则立即推送至诊断中心服务器。
诊断中心服务器接收到新的推送数据(诊断请求),会根据各诊断医生的工作量动态地给各个诊断端推送诊断数据;诊断端是一款带有GUI界面的心电在线诊断软件,诊断医生根据波形显示和统计参数给出诊断结论并生成诊断报告,并将诊断报告发送至数据库服务器。数据库服务器一旦接收到新的诊断报告过来,会把该报告显示给该心电数据所属的病人患者及其采集医生。
整个数据流程图如图2所示。
为了更好地说明本系统的核心算法,给出如下定义:
如图3所示本发明包括以下步骤:
1.构造张量数据:
我们的系统使用的标准12导联心电数据,而不是单一导联的心电数据。而心电的有效特征又往往不是仅仅存在于时域信号中,在频域上往往也存在着有效的分类特征。所以,我们的方法首先通过预处理,去除心电波形中的各种干扰和噪声,然后通过短时Fourier变换,将心电转换为张量时频域的表征形式。
具体来说,包含如下步骤:
(1)数据预处理和去燥
对于心电来说,必须要去除包括工频噪声,肌肉电干扰,基线漂移等多类噪声。首先对信号做工频频率的滤波处理,使信号通过50hz的陷波滤波器。然后对数据进行DB6小波的小波变换分解,然后去除其中最高频的和最低频的信号成份。因为低频部分包含了基线漂移,而最高频部分包含了肌肉点干扰。
(2)波形检测
再次对小波进行DB6小波分解,然后找寻其中level2小波系数,然后采用过零点检测方法检测心电的QRS波尖峰。然后依次去除R波后检测P波和T波尖峰。然后采用detrend算法计算出波形的基线,通过计算各个波形与基线的交点,确定P波QRS波T波的起始和结束。也就是onset和offset。
(3)逐跳切割
一次采集心电大约包含20秒的心电数据,也就是25跳左右的心电。我们对于心电一次beat也就是一个P波一个QRS波一个T波进行切割。
(4)R波对齐
出去诊断有效性的考虑,我们对波形进行了对齐,也就是对每一个beat的心电针对R波的峰值进行对齐,并且切割成统一的长度。
(5)短时fourier变换
为了有效抽取时频域的心电诊断特征,我们通过使用窗口为128点的短时fourier变换对心电进行转换,最后心电被转换为128×128×12的时频空的3阶张量。这里的空指的是导联位置就是指导联轴。对于12导联(lead×time)ECG信号,s[l,n] 表示离散信号在时间点n 对于导联l。 在时间点 nΔt 和 频率 f 的短时fourier变换如下式:
这里 w[n] 是一个窗函数,变换完成后数据变成一个3阶张量。
2.基于TTV映射的张量特征抽取
(1)根据原始张量数据的判别性计算有效的投影张量
(2)根据优化张量可分性特点计算更为优越的投影张量
(3)对原始数据去除已抽取投影张量的维度,调整原始数据的结构
通过这个处理之后,后续计算过程就不用考虑,直接计算出来的投影张量的影响,完全可以独立进行运算。
(4)对原始数据计算投影张量进行加权处理。
通过调整每个张量的权重,来规避不合理张量的影响,等式如下:
w(dxy)=dxy -n
或者如下式定义:
组合两种形式得到:
(5)算法整体迭代收敛过程如下:
整个算法过程,是一个严格的单调收敛过程,逐次迭代目标值关系如下式所示:
T->无穷时,算法收敛于最终目标极限收敛值.
(6)判断计算过程结束终止条件
(7)我们使用了如下方法来算法是否收敛,并且判断算法什麽时候应该停止。误差值小于一定阈值:
||Fk-Fk-1||Fro≤ε
采用这种方法来判断收敛与否和终止整个算法。
3.选择合理初值
我们方法初值选择是否合理,严重影响结果的好坏,和算法的收敛情况。所以,这里使用了一种张量逼近算法,来提供的一个有效的初始张量来提高张量算法的计算结果。基本思路就是求解一个最为近似的张量:
(1)无约束情况张量算法
a.交替最小二乘方法
这种方法的求解目标等式如下:
我们可以展开这个等式如下式:
b.张量梯度下降方法
可以将上述目标函数进行展开,写成如下的形式:
第一项没有涉及变量,所以:
这里 0 表示一个 0 向量,长度为In,第二项犹如如下式子:
第二项求导之后得到如下式:
第三项如下:
因此,
综合以上三项就可以得到最终结果。
(2)有约束张量情况
a)带约束非线性最小二乘
带约束优化问题,通过转换对求解的约束条件有所放松,再去求解就可以得到更加优越的计算结果。如下式:
将左边的优化问题转换为右边的优化问题,然后将结果代入,使用标准SQP算法求解,来求得最终的优化结果。最终结果可以得到明显改善。
4.分类比对
我们的方法最后使用SVM来对抽取出来的以vector方式存储的向量特征进行分类。
这直接通过求解如下优化主问题:
Subject to yi(wTφ(xi)+b)≥1-ξn, ξn,≥0,i=1,2…,n
这里的参数 C>0 在松弛变量和惩罚因子间的平衡。而他的Lagrangian乘子变换成如下等式:
这里 {a n≥0} 而且 { μn≥0} 是Lagrangian乘子. 而对偶Lagrangian 问题如下式:
在实际应用中采用了医院方面提供的大约10万条心电作为样本,对本发明的张量心电的方法做了实验评测。该方法和基于向量的方法PCA,ICA和LDA做了比较,同时又和基于张量的方法UMPCA,TR1DA方法做了实验比较。图4-9为不同方法进行特征抽取后降维到3维和的分类效果。可以发现,我们的方法所展示的效果可分性明显优于其它方法。实验分类精度如下所示:图7:PCA(83.4841%),图5: ICA(81.0724%),图6:LDA(77.3902%),图9:UMPCA(86.6387%),图8:TR1DA(98.4388%),图4:GTR1DA(99.2140%). 实验结果解释我们的方法GTR1DA比其它基于向量的方法明显优越,而对于其它几种基于张量的算法也有不同程度的提升。
Claims (1)
1.一种远程医疗的异常心电张量分析方法,包括下述步骤:
(1)构造张量数据:
a.心电数据采集:
采集标准12导联心电数据;
b.数据预处理和去燥:
对采集的12导联心电数据首先对信号通过50hz的陷波滤波器进行滤波处理,然后对数据进行DB6小波的小波变换分解,然后去除其中最高频的和最低频的信号成份;
c.波形检测:
再次对小波进行DB6小波分解,然后找寻其中level2小波系数,然后采用过零点检测方法检测心电的QRS波尖峰,然后依次去除R波后检测P波和T波尖峰,然后采用detrend算法计算出波形的基线,通过计算各个波形与基线的交点,确定P波QRS波T波的起始和结束,也就是onset和offset;
d.逐跳切割:
一次采集心电大约包含20秒的心电数据,也就是25跳左右的心电,对于心电一次beat也就是一个P波一个QRS波一个T波进行切割;
e. R波对齐:
对每一个beat的心电针对R波的峰值进行对齐,并且切割成统一的长度;
f. 短时fourier变换
为了有效抽取时频域的心电诊断特征,通过使用窗口为128点的短时fourier变换对心电进行转换,最后心电被转换为128×128×12的时频空的3阶张量;这里的空指的是导联位置就是指导联轴;对于12导联(lead×time)ECG信号,s[l,n]表示在时间点n,对于第l个导联轴的离散信号值; 在时间点 nΔt 和 频率 f 的短时fourier变换如下式:
这里 w[n] 是一个窗函数,变换完成后数据变成一个3阶张量;
(2)基于TTV映射的张量特征抽取:
a. 根据原始张量数据的判别性计算有效的投影张量:
b. 根据优化张量可分性特点计算更为优越的投影张量:
对原始数据去除已抽取投影张量的维度,调整原始数据的结构:
c.对原始数据计算投影张量进行加权处理:
通过调整每个张量的权重,来规避不合理张量的影响,等式如下:
w(dxy)=dxy -n
或者如下式定义:
或者组合两种形式:
d. 算法整体迭代收敛过程如下:
整个算法过程,是一个严格的单调收敛过程,逐次迭代目标值关系如下式所示:
T->无穷时,算法收敛于最终目标极限收敛值;
e. 判断计算过程结束终止条件:
使用如下方法来判断算法是否收敛,并且判断算法什么时候应该停止;误差值小于一定阈值则判断算法停止:
||Fk-Fk-1||Fro≤ε
采用这种方法来判断收敛与否和终止整个算法;
(3)选择合理初值:
求解一个最为近似的张量:
a.无约束情况张量算法:
交替最小二乘方法求解目标等式如下:
其展开等式如下式:
这里 表示 kronecker product ,而 Z(n)表示按第n mode转换张量Z一个矩阵;这个问题的解犹如下式:
张量梯度下降方法:
可以将上述目标函数进行展开,写成如下的形式:
第一项没有涉及变量,所以:
这里 0 表示一个 0 向量,长度为 In,第二项犹如如下式子:
第二项求导之后得到如下式:
第三项如下:
因此
综合以上三项就可以得到最终结果;
b. 有约束张量情况:
带约束非线性最小二乘:
带约束优化问题,通过转换对求解的约束条件有所放松,再去求解就可以得到更加优越的计算结果,如下式:
将左边的优化问题转换为右边的优化问题,然后将结果代入,使用标准SQP算法求解,来求得最终的优化结果;
(4)分类比对:
最后使用SVM来对抽取出来的以vector方式存储的向量特征进行分类,
这直接通过求解如下优化主问题:
Subject to yi(wTφ(xi)+b)≥1-ξn, ξn,≥0,i=1,2…,n
这里的参数 C>0 在松弛变量和惩罚因子间的平衡,而他的Lagrangian乘子变换成如下等式:
这里 { an≥0} 而且 { μn≥0} 是Lagrangian乘子, 而对偶Lagrangian 问题如下式:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210416931.0A CN102961129B (zh) | 2012-10-26 | 2012-10-26 | 一种远程医疗的异常心电张量分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210416931.0A CN102961129B (zh) | 2012-10-26 | 2012-10-26 | 一种远程医疗的异常心电张量分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102961129A true CN102961129A (zh) | 2013-03-13 |
CN102961129B CN102961129B (zh) | 2015-11-25 |
Family
ID=47791767
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210416931.0A Expired - Fee Related CN102961129B (zh) | 2012-10-26 | 2012-10-26 | 一种远程医疗的异常心电张量分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102961129B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014169595A1 (zh) * | 2013-04-18 | 2014-10-23 | 深圳市科曼医疗设备有限公司 | 心律失常分析方法和系统 |
CN107669262A (zh) * | 2017-10-31 | 2018-02-09 | 王量弘 | 基于svm与wlt的多导联远程心电诊断与监护系统及方法 |
CN107887010A (zh) * | 2017-11-29 | 2018-04-06 | 天津中科爱乐芙医疗科技有限公司 | 一种心血管疾病数据采集与分诊平台 |
CN108090509A (zh) * | 2017-12-13 | 2018-05-29 | 四川大学 | 一种数据长度自适应的心电图分类方法 |
CN110085108A (zh) * | 2019-03-08 | 2019-08-02 | 台州学院 | 一种心电信号仿真系统 |
CN110897627A (zh) * | 2018-09-14 | 2020-03-24 | 杭州脉流科技有限公司 | 心电图信号特征提取方法、装置、设备、系统和存储介质 |
CN110974211A (zh) * | 2019-12-09 | 2020-04-10 | 上海数创医疗科技有限公司 | 高阶多项式激活函数的st段分类神经网络及其应用 |
CN112545528A (zh) * | 2020-12-28 | 2021-03-26 | 北京理工大学 | 基于分数阶傅里叶变换和张量分解的心电t波特征提取方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2004114193A2 (en) * | 2003-05-15 | 2004-12-29 | Widemed Ltd. | Adaptive prediction of changes of physiological/pathological states using processing of biomedical signals |
US20060100534A1 (en) * | 2004-07-09 | 2006-05-11 | Ansar, Inc. | Methods for real-time autonomic nervous system monitoring using total heart rate variability, and notched windowing |
US20080058659A1 (en) * | 2006-08-29 | 2008-03-06 | Board Of Regents, The University Of Texas System | Processing and analyzing physiological signals to detect a health condition |
CN101199416A (zh) * | 2006-09-06 | 2008-06-18 | 韦伯斯特生物官能公司 | 心电图与体表测量的相关性 |
CN101690659A (zh) * | 2009-09-29 | 2010-04-07 | 华东理工大学 | 脑电波分析方法 |
US20120179059A1 (en) * | 2010-02-03 | 2012-07-12 | Zhijun Gu | Analyzing Electrocardiogram Signals |
-
2012
- 2012-10-26 CN CN201210416931.0A patent/CN102961129B/zh not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2004114193A2 (en) * | 2003-05-15 | 2004-12-29 | Widemed Ltd. | Adaptive prediction of changes of physiological/pathological states using processing of biomedical signals |
US20060100534A1 (en) * | 2004-07-09 | 2006-05-11 | Ansar, Inc. | Methods for real-time autonomic nervous system monitoring using total heart rate variability, and notched windowing |
US7529579B2 (en) * | 2004-07-09 | 2009-05-05 | Ansar, Inc. | Methods for real-time autonomic nervous system monitoring using total heart rate variability, and notched windowing |
US20080058659A1 (en) * | 2006-08-29 | 2008-03-06 | Board Of Regents, The University Of Texas System | Processing and analyzing physiological signals to detect a health condition |
CN101199416A (zh) * | 2006-09-06 | 2008-06-18 | 韦伯斯特生物官能公司 | 心电图与体表测量的相关性 |
CN101690659A (zh) * | 2009-09-29 | 2010-04-07 | 华东理工大学 | 脑电波分析方法 |
US20120179059A1 (en) * | 2010-02-03 | 2012-07-12 | Zhijun Gu | Analyzing Electrocardiogram Signals |
Non-Patent Citations (2)
Title |
---|
ÜBEYLI, ELIF DERYA.: "ECG beats classification using multiclass support vector machines with error correcting output codes", 《DIGITAL SIGNAL PROCESSING》 * |
ZHAO Q 等: "ECG feature extraction and classification using wavelet transform and support vector machines[", 《NEURAL NETWORKS AND BRAIN》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014169595A1 (zh) * | 2013-04-18 | 2014-10-23 | 深圳市科曼医疗设备有限公司 | 心律失常分析方法和系统 |
CN107669262A (zh) * | 2017-10-31 | 2018-02-09 | 王量弘 | 基于svm与wlt的多导联远程心电诊断与监护系统及方法 |
CN107887010A (zh) * | 2017-11-29 | 2018-04-06 | 天津中科爱乐芙医疗科技有限公司 | 一种心血管疾病数据采集与分诊平台 |
CN108090509A (zh) * | 2017-12-13 | 2018-05-29 | 四川大学 | 一种数据长度自适应的心电图分类方法 |
CN108090509B (zh) * | 2017-12-13 | 2021-10-08 | 四川大学 | 一种数据长度自适应的心电图分类方法 |
CN110897627A (zh) * | 2018-09-14 | 2020-03-24 | 杭州脉流科技有限公司 | 心电图信号特征提取方法、装置、设备、系统和存储介质 |
CN110085108A (zh) * | 2019-03-08 | 2019-08-02 | 台州学院 | 一种心电信号仿真系统 |
CN110085108B (zh) * | 2019-03-08 | 2021-10-29 | 台州学院 | 一种心电信号仿真系统 |
CN110974211A (zh) * | 2019-12-09 | 2020-04-10 | 上海数创医疗科技有限公司 | 高阶多项式激活函数的st段分类神经网络及其应用 |
CN112545528A (zh) * | 2020-12-28 | 2021-03-26 | 北京理工大学 | 基于分数阶傅里叶变换和张量分解的心电t波特征提取方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102961129B (zh) | 2015-11-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102961129B (zh) | 一种远程医疗的异常心电张量分析方法 | |
CN105997064B (zh) | 一种用于人体下肢表面肌电信号的辨识方法 | |
CN109171707A (zh) | 一种智能心电图分类方法 | |
CN107945817A (zh) | 心肺音信号分类方法、检测方法、装置、介质和计算机设备 | |
CN102835954B (zh) | 一种心拍波形模板生成方法及模块 | |
CN109165556A (zh) | 一种基于grnn身份识别方法 | |
EP3167799A1 (en) | Movement recognition method and apparatus based on surface electromyographic signal | |
CN103750844B (zh) | 一种基于脑电相位同步的身份识别方法 | |
CN106803081A (zh) | 一种基于多分类器集成的脑电分类方法 | |
CN109077715A (zh) | 一种基于单导联的心电信号自动分类方法 | |
CN105411565A (zh) | 基于广义尺度小波熵的心率变异性特征分类方法 | |
CN105320969A (zh) | 基于多尺度Renyi熵的心率变异性特征分类方法 | |
CN101785670A (zh) | 模糊心电图智能在线分析仪系统 | |
CN110353649B (zh) | 一种心率检测方法 | |
CN103970975B (zh) | 心电数据处理方法和系统 | |
CN106264499A (zh) | 一种量化心肺系统交互作用的分析方法 | |
CN109497997A (zh) | 基于多数据采集的癫痫发作检测设备及预警系统 | |
CN109512390B (zh) | 基于eeg时域多维度特征及m-wsvm的睡眠分期方法及可穿戴装置 | |
CN104887263A (zh) | 一种基于心音多维特征提取的身份识别算法及其系统 | |
CN105286860A (zh) | 一种基于双树复小波能量差的运动想象脑电信号识别方法 | |
CN112957056B (zh) | 利用协同网络的肌肉疲劳等级特征的提取方法及系统 | |
CN108647584A (zh) | 基于稀疏表示和神经网络的心律不齐识别分类方法 | |
CN106951852A (zh) | 一种有效的人体行为识别系统 | |
CN115081486A (zh) | 一种癫痫发作前期颅内脑电网络定位癫痫灶系统及方法 | |
CN104757966A (zh) | 胎儿心电信号的提取方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20151125 Termination date: 20161026 |