CN102793539A - 心电向量图检测分析方法及系统 - Google Patents
心电向量图检测分析方法及系统 Download PDFInfo
- Publication number
- CN102793539A CN102793539A CN201210317717XA CN201210317717A CN102793539A CN 102793539 A CN102793539 A CN 102793539A CN 201210317717X A CN201210317717X A CN 201210317717XA CN 201210317717 A CN201210317717 A CN 201210317717A CN 102793539 A CN102793539 A CN 102793539A
- Authority
- CN
- China
- Prior art keywords
- analysis
- lead
- vectorcardiogram
- check
- leads
- 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
本发明提供了一种心电向量图检测分析方法及系统,该心电向量图检测分析方法包括如下步骤:A.接收心电信号;B.判断导联体系类型;在所述步骤B中包括如下步骤:判断是否有预设信息,如有、那么根据预设信息判断导联体系类型,否则根据导联数目判断导联体系类型。本发明的有益效果是本发明公开的心电向量图检测分析方法及系统能够对不同应用场景下采集到的心电信号自动选择合适的导联体系判断方法判断当前所采用的导联体系,判断准确快速,为后续的工作做好准备。
Description
技术领域
本发明涉及心电信息处理领域,尤其涉及心电向量图检测分析方法及系统。
背景技术
心电向量图和常规心电图一样,都是记录心脏生物电的一种表现形式,两者都是从心脏电活动这一角度来反映心脏的生理和病理状态。不同的是,心电向量图记录心脏各瞬间所产生的电活动在空间的方向及大小,而心电图是心电向量图在不同导联上的二次投影,因此,心电向量图在判断心房心室肥大、束支传导阻滞、肺心病、心肌梗死、冠状动脉供血不足及预激综合征等方面有着较好的参考作用,具有更多的优越性,在临床中应用也越来越广泛。
类似于常规心电图的12导联系统,临床上多采用Frank导联系统来记录X、Y、Z导联的心电向量。但Frank导联系统与标准12导联系统在电极安放不同,因此,临床上想要同时获得一个病人的常规心电图和心电向量图,需要安装更多的电极和导联,而目前主流的心电图机装置(或系统),并非都能满足该要求。比如某医院临床使用了12道心电图机,其标配只有10个电极及导联,无法满足同时接两个系统,医生需要先后按不同的导联体系安放电极,分两次操作实现数据采集。而有些心电图机装置(或系统)又不支持Frank导联体系,此时无法通过采集获取心电向量图信息。为了解决此类问题,通常采用两种方法:1.采用12导联推导Frank X、Y、Z导联;2.更换现有心电图机装置(或系统),采用15导甚至更多导的心电图机装置(或系统)采集心电信号。此外,临床应用场景和需求不同,比如有些医生希望能在做常规心电图的获取心电向量图信息供诊断参考,有些医生对心电向量图有更高的要求,需要直接利用Frank导联体系直接采集,但目前仍没有一个操作更便捷、功能更完善的心电向量图系统适用于不同应用场景,满足不同医生的差异性需求。并且目前的设备无法自动判断和辨别所采用的是哪种导联体系,导联体系的类型、目前设备无法判断和辨别。
发明内容
为了解决现有技术中无法判断和辨别导联体系类型的问题,本发明提供了一种心电向量图检测分析方法。
本发明提供了一种心电向量图检测分析方法,包括如下步骤:
A. 接收心电信号;
B. 判断导联体系类型;
在所述步骤B中包括如下步骤:
判断是否有预设信息,如有、那么根据预设信息判断导联体系类型,否则根据导联数目判断导联体系类型。
作为本发明的进一步改进,该心电向量图检测分析方法还包括步骤:
C. 根据导联体系类型设定相应的分析导联信号。
作为本发明的进一步改进,在所述步骤C中包括如下步骤:
C1. 判断是否需要推导X、Y、Z导联,如需要推导、那么执行步骤C2,否则执行步骤C3;
C2. 标准12导联推导X、Y、Z导联;
C3. 设定分析导联信号。
作为本发明的进一步改进,该心电向量图检测分析方法还包括步骤:
D. 对各分析导联的信号进行自动检测分析,获取平均模板及特征点。
作为本发明的进一步改进,所述步骤D包括如下步骤:
首先判断是否含标准12导联信号,如含有标准12导联信号、那么便执行步骤D11至步骤D18,否则执行步骤DD1至步骤DD3;
D11. X、Y、Z导联逐搏QRS波定位;
D12. 标准12导联逐搏QRS波定位;
D13. X、Y、Z导联逐搏QRS波校正;
D14. X、Y、Z导联平均模板合成;
D15. 标准12导联平均模板合成;
D16. X、Y、Z导联平均模板波形分析;
D17. 标准12导联平均模板波形分析;
D18. X、Y、Z导联模板特征点校正;
DD1. X、Y、Z导联逐搏QRS波定位;
DD2. X、Y、Z导联平均模板合成;
DD3. X、Y、Z导联平均模板波形分析。
作为本发明的进一步改进,所述步骤D包括如下步骤:
DDD1.X、Y、Z导联逐搏QRS波定位;
DDD2. X、Y、Z导联平均模板合成;
DDD3. X、Y、Z导联平均模板波形分析;
DDD4.判断是否含标准12导联信号,如含有标准12导联信号、那么执行步骤DDD5至步骤DDD8,否则结束;
DDD5. 标准12导联逐搏QRS波定位;
DDD6. 标准12导联平均模板合成;
DDD7. 标准12导联平均模板波形分析;
DDD8. X、Y、Z导联模板特征点校正。
作为本发明的进一步改进,该心电向量图检测分析方法还包括步骤:
E. 根据分析导联中X、Y、Z导联的模板及特征点合成心电向量环并分析、计算各环参数;
所述步骤E包括如下步骤:
E1. X、Y、Z导联模板特征点归一化;
E2. 额面、横面 、侧面向量环合成;
E3. 心电向量环各参数计算。
作为本发明的进一步改进,该心电向量图检测分析方法还包括如下步骤:
F.各向量环的参数在心电向量图准则库中进行分析;
G.分析结果输出。
作为本发明的进一步改进,所述步骤F包括如下步骤:
F1. 将各向量环参数与心电向量图准则库进行比较;
F2. 判断是否含有12导联信号,如果有、则执行步骤F3,否则执行步骤G。
F3. 标准12导联心电图准则库联合比较、校正。
本发明还提供了一种心电向量图检测分析系统,包括:
接收模块,用于接收心电信号;
导联体系类型判断模块,与接收模块相连、用于判断导联体系类型。
作为本发明的进一步改进,该心电向量图检测分析系统还包括:
分析导联设定模块,与导联体系类型判断模块相连、用于根据导联体系类型设定相应的分析导联信号;
导联信号自动检测分析模块,与分析导联设定模块相连、用于对各分析导联的信号进行自动检测分析,获取平均模板及特征点。
作为本发明的进一步改进,该心电向量图检测分析系统还包括:
心电向量环自动分析模块,与导联信号自动检测分析模块相连、用于根据分析导联中X、Y、Z导联的模板及特征点合成心电向量环并分析、计算各环参数。
作为本发明的进一步改进,所述心电向量图检测分析系统还包括自动比较模块和分析输出模块,
所述自动比较模块与心电向量环自动分析模块相连、且用于将各向量环的参数与心电向量图准则库进行自动比较;
所述分析输出模块与所述导联信号自动检测分析模块相连,且所述分析输出模块用于将导联信号自动检测分析模块得到的X、Y、Z导联模板进行输出;
所述分析输出模块与心电向量环自动分析模块相连,且所述分析输出模块用于将心电向量环自动分析模块得到的心电向量环及各参数进行输出;
所述分析输出模块与所述自动比较模块相连,且所述分析输出模块用于将自动比较模块得到的心电向量图自动比较结果进行输出。
本发明的有益效果是:本发明公开的心电向量图检测分析方法及系统能够对不同应用场景下采集到的心电信号自动选择合适的导联体系判断方法判断当前所采用的导联体系,判断准确快速,为后续的工作做好准备。
附图说明
图1是本发明心电向量图检测分析方法的流程图。
图2是本发明心电向量图检测分析系统原理框图。
图3是本发明心电向量图检测分析方法一实施例流程图。
图4是本发明心电向量图检测分析方法另一实施例流程图。
图5是本发明心电向量图检测分析方法又一实施例流程图。
图6是本发明心电向量图检测分析方法中步骤S4的一实施例流程图。
图7是本发明心电向量图检测分析方法中步骤S4的另一实施例流程图。
图8是本发明心电向量图检测分析系统的一实施例原理框图。
具体实施方式
如图1所示,本发明公开了一种心电向量图检测分析方法,包括如下步骤:在步骤S1中,接收心电信号。在步骤S2中,判断导联体系类型。其中在步骤S2中具体包括如下步骤:判断是否有预设信息,如有、那么根据预设信息判断导联体系类型,否则根据导联数目判断导联体系类型。
在根据导联数目判断导联体系类型步骤中包括如下两个步骤:第一步:读取采集信号导联数据标识;第二步:根据导联数目判断所采用的导联体系,若导联数目为3、则所采用为Frank导联体系,若导联数目为12、则所采用为标准12导联体系,若导联数为15、则为标准12导联体系和Frank导联体系共存。
本发明所述的心电向量图检测分析方法为计算机或其他具备信息处理能力的设备所执行和处理的方法。
如图2所示,本发明还公开了一种心电向量图检测分析系统,包括:接收模块10,用于接收心电信号;导联体系类型判断模块11,导联体系类型判断模块11与接收模块10相连、导联体系类型判断模块11用于判断导联体系类型。
临床操作者根据需求正确连接导联线,并在心电图机装置(或系统)界面输入病人信息,如年龄、性别等,由硬件采集系统完成心电信号采集和相应处理。硬件采集系统将处理后的心电信号传输至接收模块10。
临床操作者可以通过心电图机装置(或系统)能界面操作预设导联体系信息,亦可以紧急情况下省略该步骤。
其中预设导联体系信息仅方便具体实现时标记所用,可采用多种规则。一个实施例规则如下:以标记取值不同标识不同的导联体系,标记= 0时,未预设导联体系信息,即没有预设信息。标记= 1时,标识仅Frank导联体系采集;标记= 2时,标识仅用标准12导联体系采集;标记= 3时,标识Frank导联体系采集和标准12导联体系两种导联体系并存。导联体系类型判断模块11根据标记的取值即可判定是何种导联体系在工作。
本发明公开的心电向量图检测分析方法及系统能够对不同应用场景下采集到的心电信号自动选择合适的导联体系判断方法判断当前所采用的导联体系,判断准确快速,为后续的工作做好准备。
如图3所示,该心电向量图检测分析方法还包括步骤:S3. 根据导联体系类型设定相应的分析导联信号。
在步骤S3中包括步骤S31至步骤S33,在步骤S31中,判断是否需要推导X、Y、Z导联,如需要推导、那么执行步骤S32,否则执行步骤S33。在步骤S32中,标准12导联推导X、Y、Z导联。在步骤S33中,设定分析导联信号。
在步骤S31中,若导联体系为标准12导联体系,此时采集到的心电信号不含X、Y、Z导联信号,此时需要推导得到,执行步骤S32,否则,直接执行步骤S33。
在步骤S32中,标准12导联推导X、Y、Z导联;以12导联为输入,通过转换矩阵推导X、Y、Z导联,该转换矩阵在相关文献中均有所描述,如张开滋 郭继鸿等主编,湖南科学技术出版社出版的《临床心电信息学》及MS Guillem, AV Sahakian, S Swiryn撰写的《Derivation of Orthogonal Leads from the 12-lead ECG. Accuracy of a Single Transform for the Derivation of Atrial and Ventricular Waves》等。亦可以以文献中所介绍的矩阵为基础进行改进。
在步骤S33中,设定分析导联信号;Frank导联体系时,分析导联为采集得到的X、Y、Z导联;标准12导联体系时,分析导联为采集得到的12导联和推导得到的X、Y、Z导联;Frank导联体系和标准12导联体系共存时,分析导联为采集得到X、Y、Z导联和标准12导联。
通过在不同应用场景下选择合适的导联体系判断方法判断当前所采用的导联体系并选择相应的分析导联,能够适用于只有标准12导联体系、只有Frank导联体系和两个导联体系共存时的情况,通过步骤S32完成数据的转换,从而可以适应多种情况。
该心电向量图检测分析方法还包括步骤S4和步骤S5,在步骤S4中,对各分析导联的信号进行自动检测分析,获取平均模板及特征点。对步骤S3中的分析导联进行单导联或者多导联联合分析,定位逐搏QRS波,通过有效平均叠加,合成各分析导联的平均模板,并对该平均模板进行QRS波及P、T波检测以及各波起点、终点定位。
在步骤S5中,根据分析导联中X、Y、Z导联的模板及特征点合成心电向量环并分析、计算各环参数。根据前述步骤S4得到分析导联中X、Y、Z导联的平均模板及其P波、QRS波和T波起点、终点,合成额面、横面、侧面P向量环、QRS向量环和T向量环。对各环进行检测分析,并计算各特征参数。作为公知的方法,X、Y导联合成额面,X、Z导联合成横面(又称水平面),Y、Z导联合成侧面。
通过执行步骤S3至步骤S5,可以得到更为准确的X、Y、Z导联模板及对应特征点,确保心电向量环合成的准确性,进而提高心电向量环各环参数的准确性。
如图4所示,该心电向量图检测分析方法还包括步骤S6和步骤S7。在步骤S6中,各向量环的参数在心电向量图准则库中进行分析;按照病人年龄和性别选择预设的心电向量图准则库,将步骤S5中得到的各向量环的众多参数与该准则库进行自动分析,并给出分析结果。
在步骤S7中,分析结果输出;将步骤S5得到的各向量环的参数和步骤S6得到的自动比较结果按需输出。
如图5所示,在步骤S5中包括步骤S51至步骤S53。
在步骤S51中,X、Y、Z导联模板特征点归一化;为了保证准确生成向量环,根据前述步骤S4得到的X、Y、Z导联模板各波的特征点进行归一化处理; 对步骤S33设定的分析导联进行分析,分别对P波起点、P波终点、QRS波起点、QRS波终点、T波起点和T波终点进行多导联有效平均计算; 以P波起点为例,按照设定的剔除规则,去掉各导联P波起点误差较大的值,对其余P波起点取均值,作为X、Y、Z导联统一的P波起点,方便说明,下文中以Pb_all代替该统一的P波起点,以区分前述步骤S4中未经归一化处理的P波起点; P波终点Pe_all,QRS波起点QRSb_all、QRS波终点QRSe_all、T波起点Tb_all和T波终点Te_all亦采用同样的方法处理得到。
在步骤S52中,额面、横面 、侧面向量环合成;根据前述步骤S4得到分析导联中X、Y、Z导联的模板及其P波、QRS波和T波起点、终点,合成额面、横面、侧面P向量环、QRS向量环和T向量环;该步骤以公知的方法合成各面的向量环,其中:X、Y导联模板合成额面向量环,X、Z导联模板合成横面(又称水平面)向量环,Y、Z导联合成侧面向量环。其中侧面可以分为右侧面和左侧面,具体选择可由临床操作者通过界面实现。各面向量环组成为:由P波起点到P波终点形成P向量环,由QRS波起点至QRS波终点形成QRS环,由T波起点到T波终点形成T向量环。
为了方便向量环画图及后续的参数计算,在合成向量环时,可以先行对X、Y、Z导联的模板进行基准等电位线(下述以“基线”简称)调整,以保证各环的起点处在坐标零点的位置,然后再根据上述说明合成各面的P、QRS、T向量环。其中,基线调整有多种方法实现可供选择,举例说明如下:
方法1:各导联直接以QRS波起点QRSb_all为基准点,进行基线调整,调整后,QRS波起点处幅值为零。
方法2:根据步骤SS1或步骤QQ1至步骤QQ3中得到的X、Y、Z导联逐搏QRS波信息计算心率,根据心率水平设定不通的持续时间阈值TH,在QRS波起点QRSb_all前设定搜索窗,若持续TH毫秒内,信号幅度低于20uV,则将该时间段内的平均信号幅值作为基准等电位线。若所搜窗内没有满足上述条件的基准等电位线,在直接取PR段(即P波终点Pe_all到QRS波起点QRSb_all)为基准等电位线。
在步骤S53中,心电向量环各特征参数计算;对步骤S52得到的各向量环进行检测、分析。并计算各面P向量环、QRS向量环、T向量环的最大向量、运行方向、各向力、各象限面积占比、各环时限等一百多个参数。
本发明通过执行步骤S1至步骤S5,实现了在不同应用场景下选择合适的导联体系判断方法判断当前所采用的导联体系并选择相应的分析导联,能够适用于只有标准12导联体系、只有Frank导联体系(X、Y、Z导联)和两个导联体系共存时获取准确的心电向量图,即可以适用于12导及以下的心电图机装置(或系统),也适用于15导及以上的心电图机装置或(系统),满足获得心电向量图的不同应用需求。
在步骤S6中包括步骤S61至步骤S63。在步骤S61中,将各向量环参数与心电向量图准则库进行比较;按照病人年龄和性别选择预设的心电向量图准则库,将步骤S53得到的各向量环的众多参数与该准则库进行自动比较,并给出比较结果。
在步骤S62中,判断是否含有12导联信号,如果有,则执行步骤S63,否则,执行步骤S7。
在步骤S63中,标准12导联心电图准则库联合比较、校正;利用12导联心电图与心电图准则库的比较结果对步骤S61得到的心电向量图比较结果进行校正,以得到更准确的心电向量图比较结果。
在步骤S7中,自动分析、比较结果输出;将步骤S53得到的各向量环的参数和步骤S63得到的自动比较结果按需输出。
本发明通过执行步骤S6至步骤S7,实现了将获取的心电向量图进行自动检测分析,并将所得各参数与心电向量图准则库进行自动比较,给出可靠、准确的参数及比较结果。
如图6所示,作为步骤S4的一个实施例,所述步骤S4包括如下步骤:首先判断是否含标准12导联信号,如含有标准12导联信号、那么便执行步骤QQ1至步骤QQ8,否则执行步骤SS1至步骤SS3。
在步骤SS1中,X、Y、Z导联逐搏QRS波定位;对选定的节律导联(操作者可将X、Y、Z导联之一预设为节律导联,未作选择的情况下默认为Y导联)进行逐搏QRS波检测,定位QRS波主峰位置;联合未选做节律导联的两个导联,进行多导联联合QRS波定位,确定X、Y、Z各导联的QRS波主峰位置。
在步骤SS2中,X、Y、Z导联平均模板合成;X、Y、Z各导联合成一个平均模板心搏,包括下述三个步骤:第一步,对X、Y、Z导联逐搏QRS波进行形态分类,对占主导优势的部分心搏利用中值法挑选中值心搏;第二步,依次选择上述步骤SS2中的第一步中的其他优势心搏,与中值心搏进行逐点相关运算,记录所得最大相关系数作为当前心搏的相关系数;第三步,对满足相关系数阈值条件的心搏进行叠加平均,合成平均模板;若无满足相关系数阈值条件的心搏,则取中值心搏作为平均模板。
在步骤SS3中,X、Y、Z导联平均模板波形分析;分别对步骤SS2得到的X、Y、Z导联的平均模板进行检测分析,完成对各模板P、QRS、T波的起、止点定位;首先检测QRS波起点、终点,然后在QRS波起点前一段时间范围内搜索P波,并判断其起点、终点,在QRS波终点之后一段时间范围内搜索T波,并判断其起点、终点。
在步骤QQ1中,X、Y、Z导联逐搏QRS波定位;对选定的节律导联(操作者可将X、Y、Z导联之一预设为节律导联,未作选择的情况下默认为Y导联)进行逐搏QRS波检测,定位QRS波主峰位置;联合未选做节律导联的两个导联,进行多导联联合QRS波定位,确定X、Y、Z各导联的QRS波主峰位置。
在步骤QQ2中,标准12导联逐搏QRS波定位;对选定的节律导联(操作者可将标准12导联的任一导联预设为节律导联,未作预设时,默认为II导联)进行逐搏QRS波检测,定位QRS波主峰位置;对其他未选做节律导联的导联进行多导联联合QRS波定位,确定标准12导联各导联的QRS波主峰位置。
在步骤QQ3中,X、Y、Z导联逐搏QRS波校正;为了提高QRS波定位的准确性,避免干扰引起的定位误差或者QRS波误检,利用步骤QQ2得到的标准12导联的定位结果对X、Y、Z导联的逐搏QRS波进行校正,包含下述步骤:第一步,节律导联的QRS波主峰位置对X、Y、Z导联中节律导联的QRS波定位进行校正,对于主峰位置相差超过阈值的QRS波进行定位调整或者剔除;第二步,利用校正的X、Y、Z导联节律导联对非节律导联进行校正。
在步骤QQ4中,X、Y、Z导联平均模板合成;X、Y、Z各导联合成一个平均模板心搏,包括下述三个步骤:第一步,对X、Y、Z导联逐搏QRS波进行形态分类,对占主导优势的部分心搏利用中值法挑选中值心搏;第二步,依次选择上述步骤QQ4中的第一步中的其他优势心搏,与中值心搏进行逐点相关运算,记录所得最大相关系数作为当前心搏的相关系数;第三步,对满足相关系数阈值条件的心搏进行叠加平均,合成平均模板;若无满足相关系数阈值条件的心搏,则取中值心搏作为平均模板。
在步骤QQ5中,标准12导联平均模板合成;标准12导联各导联合成一个平均模板心搏,包括下述步骤:第一步,根据步骤QQ2中得到的逐搏QRS波定位结果,对各QRS波进行形态分类,对占主导优势的部分心搏利用中值法挑选中值心搏;第二步,依次选择步骤QQ5中的第一步中的其他优势心搏,与中值心搏进行逐点相关运算,记录所得最大相关系数作为当前心搏的相关系数ρ;第三步,对满足相关系数阈值条件的心搏进行叠加平均,合成平均模板;若无满足相关系数阈值条件的心搏,则取中值心搏作为平均模板。
在步骤QQ6中,X、Y、Z导联平均模板波形分析;分别对步骤QQ4中得到的X、Y、Z导联的平均模板进行检测分析,完成对各模板P、QRS、T波的起、止点定位;首先检测QRS波起点、终点,然后在QRS波起点前一段时间范围内搜索P波,并判断其起点、终点,在QRS波终点之后一段时间范围内搜索T波,并判断其起点、终点。
在步骤QQ7中,标准12导联平均模板波形分析;分别对步骤QQ5中得到的标准12导联平均模板进行分析,完成对各模板P、QRS、T波的起、止点定位;最终得到12个导联的P波起点、终点,QRS波起点、终点和T波起点、终点。
在步骤QQ8中,X、Y、Z导联模板特征点校正;联合步骤QQ6和QQ7的结果,对X、Y、Z导联和标准12导联共15个导联模板的各特征点进行校正;以P波起点为例,首先按照设定的剔除规则去掉15个P波起点中误差较大的,用剩下满足误差条件的P波起点计算平均值;然后判断当前导联P波起点与该平均值是否满足误差要求,如果满足,则当前导联P波起点不作校正和调整,否则,用该平均值代替当前导联的P波起点,完成校正。P波终点,QRS波起点、QRS波终点、T波起点和T波终点亦采用同样的方法处理。
如图7所示,作为步骤S4的另一个实施例,所述步骤S4包括步骤W1至步骤W8。
在步骤W1中,X、Y、Z导联逐搏QRS波定位;对选定的节律导联(操作者可将X、Y、Z导联之一预设为节律导联,未作选择的情况下默认为Y导联)进行逐搏QRS波检测,定位QRS波主峰位置;联合未选做节律导联的两个导联,进行多导联联合QRS波定位,确定X、Y、Z各导联的QRS波主峰位置。
在步骤W2中,X、Y、Z导联平均模板合成;X、Y、Z各导联合成一个平均模板心搏,包括下述三个步骤:第一步,对X、Y、Z导联逐搏QRS波进行形态分类,对占主导优势的部分心搏利用中值法挑选中值心搏;第二步,依次选择上述步骤W2中的第一步中的其他优势心搏,与中值心搏进行逐点相关运算,记录所得最大相关系数作为当前心搏的相关系数;第三步,对满足相关系数阈值条件的心搏进行叠加平均,合成平均模板;若无满足相关系数阈值条件的心搏,则取中值心搏作为平均模板。
在步骤W3中,X、Y、Z导联平均模板波形分析;分别对步骤W2得到的X、Y、Z导联的平均模板进行检测分析,完成对各模板P、QRS、T波的起、止点定位;首先检测QRS波起点、终点,然后在QRS波起点前一段时间范围内搜索P波,并判断其起点、终点,在QRS波终点之后一段时间范围内搜索T波,并判断其起点、终点。
在步骤W4中,判断是否含12导联信号,若不含标准12导联信号、则分析结束,否则执行步骤W5至步骤W8;执行步骤W5至步骤W8是对对标准12导联进行检测分析,并将其分析结果用于X、Y、Z导联检测分析的校正,以进一步提高检测分析的准确性。
在步骤W5中,标准12导联逐搏QRS波定位;对选定的节律导联(操作者可将标准12导联的任一导联预设为节律导联,未作预设时,默认为II导联)进行逐搏QRS波检测,定位QRS波主峰位置;对其他未选做节律导联的导联进行多导联联合QRS波定位,确定标准12导联各导联的QRS波主峰位置。
在步骤W6中,标准12导联平均模板合成;标准12导联各导联合成一个平均模板心搏,包括下述步骤:第一步,根据步骤W5中得到的逐搏QRS波定位结果,对各QRS波进行形态分类,对占主导优势的部分心搏利用中值法挑选中值心搏;第二步,依次选择步骤W6中的第一步中的其他优势心搏,与中值心搏进行逐点相关运算,记录所得最大相关系数作为当前心搏的相关系数ρ;第三步,对满足相关系数阈值条件的心搏进行叠加平均,合成平均模板;若无满足相关系数阈值条件的心搏,则取中值心搏作为平均模板。
在步骤W7中,标准12导联平均模板波形分析;分别对步骤W6中得到的标准12导联平均模板进行分析,完成对各模板P、QRS、T波的起、止点定位;最终得到12个导联的P波起点、终点,QRS波起点、终点和T波起点、终点。
在步骤W8中,X、Y、Z导联模板特征点校正;联合步骤W3和W7的结果,对X、Y、Z导联和标准12导联共15个导联模板的各特征点进行校正;以P波起点为例,首先按照设定的剔除规则去掉15个P波起点中误差较大的,用剩下满足误差条件的P波起点计算平均值;然后判断当前导联P波起点与该平均值是否满足误差要求,如果满足,则当前导联P波起点不作校正和调整,否则,用该平均值代替当前导联的P波起点,完成校正。P波终点,QRS波起点、QRS波终点、T波起点和T波终点亦采用同样的方法处理。
如图8所示,本发明公开的心电向量图检测分析系统还包括分析导联设定模块12,与导联体系类型判断模块11相连、用于根据导联体系类型设定相应的分析导联信号;导联信号自动检测分析模块13,与分析导联设定模块12相连、用于对各分析导联的信号进行自动检测分析,获取平均模板及特征点;心电向量环自动分析模块14,与导联信号自动检测分析模块13相连、用于根据分析导联中X、Y、Z导联的模板及特征点合成心电向量环并分析、计算各环参数;自动比较模块15,与心电向量环自动分析模块14相连、且用于将各向量环的参数与心电向量图准则库进行自动比较;
所述分析输出模块16与所述导联信号自动检测分析模块13相连,且所述分析输出模块16用于将导联信号自动检测分析模块13得到的X、Y、Z导联模板进行输出;
所述分析输出模块16与心电向量环自动分析模块14相连,且所述分析输出模块16用于将心电向量环自动分析模块14得到的心电向量环及各参数进行输出;
所述分析输出模块16与所述自动比较模块15相连,且所述分析输出模块16用于将自动比较模块15得到的心电向量图自动比较结果进行输出。
所述输出包括显示和/或打印。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明,不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干简单推演或替换,都应当视为属于本发明的保护范围。
Claims (13)
1.一种心电向量图检测分析方法,其特征在于,包括如下步骤:
A. 接收心电信号;
B. 判断导联体系类型;
在所述步骤B中包括如下步骤:
判断是否有预设信息,如有、那么根据预设信息判断导联体系类型,否则根据导联数目判断导联体系类型。
2.根据权利要求1所述的心电向量图检测分析方法,其特征在于,该心电向量图检测分析方法还包括步骤:
C. 根据导联体系类型设定相应的分析导联信号。
3. 根据权利要求2所述的心电向量图检测分析方法,其特征在于,在所述步骤C中包括如下步骤:
C1. 判断是否需要推导X、Y、Z导联,如需要推导、那么执行步骤C2,否则执行步骤C3;
C2. 标准12导联推导X、Y、Z导联;
C3. 设定分析导联信号。
4.根据权利要求3所述的心电向量图检测分析方法,其特征在于,该心电向量图检测分析方法还包括步骤:
D. 对各分析导联的信号进行自动检测分析,获取平均模板及特征点。
5. 根据权利要求4所述的心电向量图检测分析方法,其特征在于,所述步骤D包括如下步骤:
首先判断是否含标准12导联信号,如含有标准12导联信号、那么便执行步骤D11至步骤D18,否则执行步骤DD1至步骤DD3;
D11. X、Y、Z导联逐搏QRS波定位;
D12. 标准12导联逐搏QRS波定位;
D13. X、Y、Z导联逐搏QRS波校正;
D14. X、Y、Z导联平均模板合成;
D15. 标准12导联平均模板合成;
D16. X、Y、Z导联平均模板波形分析;
D17. 标准12导联平均模板波形分析;
D18. X、Y、Z导联模板特征点校正;
DD1. X、Y、Z导联逐搏QRS波定位;
DD2. X、Y、Z导联平均模板合成;
DD3. X、Y、Z导联平均模板波形分析。
6. 根据权利要求4所述的心电向量图检测分析方法,其特征在于,所述步骤D包括如下步骤:
DDD1.X、Y、Z导联逐搏QRS波定位;
DDD2. X、Y、Z导联平均模板合成;
DDD3. X、Y、Z导联平均模板波形分析;
DDD4.判断是否含标准12导联信号,如含有标准12导联信号、那么执行步骤DDD5至步骤DDD8,否则结束;
DDD5. 标准12导联逐搏QRS波定位;
DDD6. 标准12导联平均模板合成;
DDD7. 标准12导联平均模板波形分析;
DDD8. X、Y、Z导联模板特征点校正。
7. 根据权利要求1至6任一项所述的心电向量图检测分析方法,其特征在于,该心电向量图检测分析方法还包括步骤:
E. 根据分析导联中X、Y、Z导联的模板及特征点合成心电向量环并分析、计算各环参数;
所述步骤E包括如下步骤:
E1. X、Y、Z导联模板特征点归一化;
E2. 额面、横面 、侧面向量环合成;
E3. 心电向量环各参数计算。
8. 根据权利要求7所述的心电向量图检测分析方法,其特征在于,该心电向量图检测分析方法还包括如下步骤:
F.各向量环的参数在心电向量图准则库中进行分析;
G.分析结果输出。
9. 根据权利要求8所述的心电向量图检测分析方法,其特征在于,所述步骤F包括如下步骤:
F1. 将各向量环参数与心电向量图准则库进行比较;
F2. 判断是否含有12导联信号,如果有、则执行步骤F3,否则执行步骤G;
F3. 标准12导联心电图准则库联合比较、校正。
10. 一种心电向量图检测分析系统,其特征在于,包括:
接收模块,用于接收心电信号;
导联体系类型判断模块,与接收模块相连、用于判断导联体系类型。
11. 根据权利要求10所述的心电向量图检测分析系统,其特征在于,该心电向量图检测分析系统还包括:
分析导联设定模块,与导联体系类型判断模块相连、用于根据导联体系类型设定相应的分析导联信号;
导联信号自动检测分析模块,与分析导联设定模块相连、用于对各分析导联的信号进行自动检测分析,获取平均模板及特征点。
12. 根据权利要求11所述的心电向量图检测分析系统,其特征在于,该心电向量图检测分析系统还包括:
心电向量环自动分析模块,与导联信号自动检测分析模块相连、用于根据分析导联中X、Y、Z导联的模板及特征点合成心电向量环并分析、计算各环参数。
13. 根据权利要求12所述的心电向量图检测分析系统,其特征在于:所述心电向量图检测分析系统还包括自动比较模块和分析输出模块,
所述自动比较模块与心电向量环自动分析模块相连、且用于将各向量环的参数与心电向量图准则库进行自动比较;
所述分析输出模块与所述导联信号自动检测分析模块相连,且所述分析输出模块用于将导联信号自动检测分析模块得到的X、Y、Z导联模板进行输出;
所述分析输出模块与心电向量环自动分析模块相连,且所述分析输出模块用于将心电向量环自动分析模块得到的心电向量环及各参数进行输出;
所述分析输出模块与所述自动比较模块相连,且所述分析输出模块用于将自动比较模块得到的心电向量图自动比较结果进行输出。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210317717.XA CN102793539B (zh) | 2012-08-31 | 2012-08-31 | 心电向量图检测分析方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210317717.XA CN102793539B (zh) | 2012-08-31 | 2012-08-31 | 心电向量图检测分析方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102793539A true CN102793539A (zh) | 2012-11-28 |
CN102793539B CN102793539B (zh) | 2014-07-09 |
Family
ID=47193032
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210317717.XA Active CN102793539B (zh) | 2012-08-31 | 2012-08-31 | 心电向量图检测分析方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102793539B (zh) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103810393A (zh) * | 2014-02-27 | 2014-05-21 | 厦门纳龙科技有限公司 | 一种基于心电向量的心电波形特征点定位方法 |
CN105411575B (zh) * | 2015-12-01 | 2016-08-17 | 浙江铭众科技有限公司 | 一种心电电极位置监控方法 |
CN106889984A (zh) * | 2017-01-22 | 2017-06-27 | 河北大学 | 一种心电信号自动降噪方法 |
CN107260161A (zh) * | 2017-07-18 | 2017-10-20 | 上海图灵医疗科技有限公司 | 一种心电动力学数据量化分析方法 |
CN108125678A (zh) * | 2017-12-01 | 2018-06-08 | 北京顺源开华科技有限公司 | 心电信号的方向检测方法、装置及电子设备 |
CN110013247A (zh) * | 2019-05-24 | 2019-07-16 | 东北大学 | 一种心电图p波和t波的检测、区分与定位方法 |
CN110090012A (zh) * | 2019-03-15 | 2019-08-06 | 上海图灵医疗科技有限公司 | 一种基于机器学习的人体疾病检测方法及检测产品 |
CN111603170A (zh) * | 2020-06-01 | 2020-09-01 | 河南华南医电科技有限公司 | 一种基于心电向量图的人体体位检测方法 |
WO2021148045A1 (zh) * | 2020-01-20 | 2021-07-29 | 深圳市理邦精密仪器股份有限公司 | 心电波形的数据测量方法、心电图机和可读存储介质 |
CN113768511A (zh) * | 2020-06-04 | 2021-12-10 | 深圳市理邦精密仪器股份有限公司 | 生理参数检测方法及电子设备 |
CN114869296A (zh) * | 2022-07-08 | 2022-08-09 | 长沙博为软件技术股份有限公司 | 定位qrs-t波参考水平的方法、系统、设备及介质 |
CN116077785A (zh) * | 2023-02-08 | 2023-05-09 | 北京新尖科技有限公司 | 一种基于心电矢量弧的快速心搏识别算法 |
CN116077785B (zh) * | 2023-02-08 | 2024-05-24 | 北京新尖科技有限公司 | 一种基于心电矢量弧的快速心搏识别算法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5713367A (en) * | 1994-01-26 | 1998-02-03 | Cambridge Heart, Inc. | Measuring and assessing cardiac electrical stability |
CN1206580A (zh) * | 1998-08-13 | 1999-02-03 | 赵峰 | 立体心电图的成像方法及成像仪 |
CN2530301Y (zh) * | 2001-03-09 | 2003-01-08 | 北京卡迪欧医疗设备有限责任公司 | 用于对心电信号进行分析的分析器 |
US20040006265A1 (en) * | 2002-04-30 | 2004-01-08 | Karim Alhussiny | Wireless transmission-ST-segment preserved of the standard 12 leads EKG apparatus for the remote administration of thrrombolytic therapy under severe cellular channel impairment |
CN101011242A (zh) * | 2007-02-06 | 2007-08-08 | 赵峰 | 一种向量心电图仪及实现方法 |
JP2011072542A (ja) * | 2009-09-30 | 2011-04-14 | Fukuda Denshi Co Ltd | 心電図データ処理装置、心電図データ処理方法および心電図データ処理プログラム |
CN102028458A (zh) * | 2011-01-13 | 2011-04-27 | 深圳市科曼医疗设备有限公司 | 一种实现同步输出心电图和心电向量图的方法 |
-
2012
- 2012-08-31 CN CN201210317717.XA patent/CN102793539B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5713367A (en) * | 1994-01-26 | 1998-02-03 | Cambridge Heart, Inc. | Measuring and assessing cardiac electrical stability |
CN1206580A (zh) * | 1998-08-13 | 1999-02-03 | 赵峰 | 立体心电图的成像方法及成像仪 |
CN2530301Y (zh) * | 2001-03-09 | 2003-01-08 | 北京卡迪欧医疗设备有限责任公司 | 用于对心电信号进行分析的分析器 |
US20040006265A1 (en) * | 2002-04-30 | 2004-01-08 | Karim Alhussiny | Wireless transmission-ST-segment preserved of the standard 12 leads EKG apparatus for the remote administration of thrrombolytic therapy under severe cellular channel impairment |
CN101011242A (zh) * | 2007-02-06 | 2007-08-08 | 赵峰 | 一种向量心电图仪及实现方法 |
JP2011072542A (ja) * | 2009-09-30 | 2011-04-14 | Fukuda Denshi Co Ltd | 心電図データ処理装置、心電図データ処理方法および心電図データ処理プログラム |
CN102028458A (zh) * | 2011-01-13 | 2011-04-27 | 深圳市科曼医疗设备有限公司 | 一种实现同步输出心电图和心电向量图的方法 |
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103810393A (zh) * | 2014-02-27 | 2014-05-21 | 厦门纳龙科技有限公司 | 一种基于心电向量的心电波形特征点定位方法 |
CN105411575B (zh) * | 2015-12-01 | 2016-08-17 | 浙江铭众科技有限公司 | 一种心电电极位置监控方法 |
CN106889984A (zh) * | 2017-01-22 | 2017-06-27 | 河北大学 | 一种心电信号自动降噪方法 |
CN107260161A (zh) * | 2017-07-18 | 2017-10-20 | 上海图灵医疗科技有限公司 | 一种心电动力学数据量化分析方法 |
WO2019015017A1 (zh) * | 2017-07-18 | 2019-01-24 | 上海图灵医疗科技有限公司 | 一种心电动力学数据量化分析方法 |
CN107260161B (zh) * | 2017-07-18 | 2019-02-19 | 上海图灵医疗科技有限公司 | 一种心电动力学数据量化分析方法 |
CN108125678A (zh) * | 2017-12-01 | 2018-06-08 | 北京顺源开华科技有限公司 | 心电信号的方向检测方法、装置及电子设备 |
CN108125678B (zh) * | 2017-12-01 | 2021-01-26 | 北京顺源开华科技有限公司 | 心电信号的方向检测方法、装置及电子设备 |
CN110090012A (zh) * | 2019-03-15 | 2019-08-06 | 上海图灵医疗科技有限公司 | 一种基于机器学习的人体疾病检测方法及检测产品 |
CN110013247B (zh) * | 2019-05-24 | 2020-06-05 | 东北大学 | 一种心电图p波和t波的检测、区分与定位方法 |
CN110013247A (zh) * | 2019-05-24 | 2019-07-16 | 东北大学 | 一种心电图p波和t波的检测、区分与定位方法 |
WO2021148045A1 (zh) * | 2020-01-20 | 2021-07-29 | 深圳市理邦精密仪器股份有限公司 | 心电波形的数据测量方法、心电图机和可读存储介质 |
CN111603170A (zh) * | 2020-06-01 | 2020-09-01 | 河南华南医电科技有限公司 | 一种基于心电向量图的人体体位检测方法 |
CN113768511A (zh) * | 2020-06-04 | 2021-12-10 | 深圳市理邦精密仪器股份有限公司 | 生理参数检测方法及电子设备 |
CN113768511B (zh) * | 2020-06-04 | 2023-09-22 | 深圳市理邦精密仪器股份有限公司 | 生理参数检测方法及电子设备 |
CN114869296A (zh) * | 2022-07-08 | 2022-08-09 | 长沙博为软件技术股份有限公司 | 定位qrs-t波参考水平的方法、系统、设备及介质 |
CN114869296B (zh) * | 2022-07-08 | 2022-10-04 | 长沙博为软件技术股份有限公司 | 定位qrs-t波参考水平的方法、系统、设备及介质 |
CN116077785A (zh) * | 2023-02-08 | 2023-05-09 | 北京新尖科技有限公司 | 一种基于心电矢量弧的快速心搏识别算法 |
CN116077785B (zh) * | 2023-02-08 | 2024-05-24 | 北京新尖科技有限公司 | 一种基于心电矢量弧的快速心搏识别算法 |
Also Published As
Publication number | Publication date |
---|---|
CN102793539B (zh) | 2014-07-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102793539B (zh) | 心电向量图检测分析方法及系统 | |
CN110013247B (zh) | 一种心电图p波和t波的检测、区分与定位方法 | |
CN101810476B (zh) | 一种动态心电图的心拍模板分类方法 | |
Kuchar et al. | Electrocardiographic localization of the site of origin of ventricular tachycardia in patients with prior myocardial infarction | |
Madeiro et al. | An innovative approach of QRS segmentation based on first-derivative, Hilbert and Wavelet Transforms | |
CN102834050B (zh) | 儿童心电图自动检测分析方法及系统 | |
US5967994A (en) | Method and system for characterizing the quality of signals indicative of heart function | |
US20200113474A1 (en) | Matching and tracking time sequences of heart activation | |
JP6901399B2 (ja) | Ecgデバイス用のecgコントローラ、該ecgコントローラの作動の方法 | |
CN103815897B (zh) | 一种心电图特征提取方法 | |
CN105411567A (zh) | 一种心电图数据分析方法及系统 | |
CN103110417A (zh) | 一种心电图自动识别系统 | |
CN101059793A (zh) | 用于心脏病学和病人监护数据分析的多层系统 | |
CN104856669A (zh) | 一种心电图导联错接的纠正方法 | |
CN102379692A (zh) | 一种面向区端网络的远程心电监护系统及方法 | |
Yazdani et al. | Adaptive mathematical morphology for QRS fiducial points detection in the ECG | |
CN104680186B (zh) | 一种心电信号中st段形态的自动分类方法 | |
CA2567282A1 (en) | Differentiating ischemic from non-ischemic t-wave inversion | |
CN109077720B (zh) | 信号处理方法、装置、设备和存储介质 | |
CN111419212A (zh) | 心电图数据处理的方法、装置、存储介质及计算机设备 | |
CN101773392B (zh) | 动态心电数据自适应高效存储方法 | |
Christov et al. | Ranking of the most reliable beat morphology and heart rate variability features for the detection of atrial fibrillation in short single-lead ECG | |
CN105877742A (zh) | 一种针对心电信号ⅰ导联电极的错放检测方法 | |
CN102349828A (zh) | 一种心电图导联纠错方法 | |
CN112971800B (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 |