一种通过计算机实现的心电波形的检测方法及装置制造方法
【专利摘要】本发明提供了一种通过计算机实现的心电波形的检测方法及装置，该方法包括：获取M个心拍波形数据，M为正整数；根据聚类方法将所述M个心拍波形数据聚类成至少一个第一簇；从所述至少一个第一簇当中选取心拍波形数据的个数最多的第一簇；基于选取的第一簇生成标准心拍模板。该装置包括：获取单元，用于获取M个心拍波形数据，M为正整数；第一聚类单元，用于根据聚类方法将所述M个心拍波形数据聚类成至少一个第一簇；第一选取单元，用于从所述至少一个第一簇当中选取心拍波形数据的个数最多的第一簇；生成单元，用于基于选取的第一簇生成标准心拍模板。本发明通过计算机自动地生成标准心拍模板，并且具有较好的准确性。
【专利说明】一种通过计算机实现的心电波形的检测方法及装置
【技术领域】
[0001 ] 本发明涉及心电波形处理，尤其涉及一种心电波形的检测方法及装置。
【背景技术】
[0002]心律失常是一种常见的心脏病，动态心电信可以帮助诊断心律失常的早期诊断。由于动态心电图记录的数据量大，不可能由人工来完成分析，因而心电图的自动分析显得尤为重要，已经成为了当前信号处理领域的研究热点。
[0003]心电波形是一个时域周期信号，每个周期由一个心跳的节拍波形组成。图1为典型的心跳节拍波形的示意图。
[0004]如图1所示，心电图的波形周期划分为P波、P-R段、QRS复合波、S-T段和T波，每个部分都代表着相应的生理意义。如果心脏产生病理性变化，则心电信号将在波形周期和形态上发生畸变。
[0005]现有技术中选用心电波形的前N拍生成心拍模板。在模板没有生成前，通过第一个QRS复合波中的R点和第二个QRS复合波的R点之间的间隔判断心拍是否正常。
[0006]通过上述方法判断正常的心拍会导致以下问题:即使第一个QRS复合波中的R点和第二个QRS复合波的R点之间的间隔正常，由于噪音或其他因素，以此波形作为模板与其他的波形对比，会产生误差。
【发明内容】
[0007]本发明的实施例提供一种心电波形的检测方法及装置，通过采用聚类的方法得到标准心拍模板，从而通过计算机处理在波形分类中自动地、较准确地生成标准心拍模板。
[0008]为达到上述目的，本发明的实施例采用如下技术方案:
[0009]一种通过计算机实现的心电波形的检测方法，所述方法包括:
[0010]获取M个心拍波形数据，M为正整数；
[0011]根据聚类方法将所述M个心拍波形数据聚类成至少一个第一簇；
[0012]从所述至少一个第一簇当中选取心拍波形数据的个数最多的第一簇；
[0013]基于选取的第一簇生成标准心拍模板。
[0014]一种心电波形的检测装置，所述装置包括:
[0015]获取单元，用于获取M个心拍波形数据，M为正整数；
[0016]第一聚类单元，用于根据聚类方法将所述M个心拍波形数据聚类成至少一个第一簇；
[0017]第一选取单元，用于从所述至少一个第一簇当中选取心拍波形数据的个数最多的第一簇；
[0018]生成单元，用于基于选取的第一簇生成标准心拍模板。
[0019]本发明实施例提供的心电波形的检测方法和装置，根据聚类方法将趋于正常的心拍波形数据聚类到一起，从而根据聚类得到的包括心拍个数最多的簇来生成标准心拍模板，并且通过计算机自动地生成标准心拍模板，具有较好的准确性。
【专利附图】
【附图说明】
[0020]图1为典型的心跳节拍波形的示意图；
[0021]图2为本发明实施例提供的通过计算机实现的心电波形的检测方法流程图；
[0022]图3为本发明实施例提供的心电波形的检测装置示意图。
【具体实施方式】
[0023]心跳节拍的正常波形是相似的，而异常波形是千差万别的。同时，由于心电图因人而异，很难取得一个固定的模板。因此，在分析病人的心电图之前，通过对心电波形的分析确定正常的心拍模板，可有效去除心电图分析中的个体差异。本发明提出一种通过聚类分析算法对心电波形进行聚类，以将趋于正常的波形汇聚到一起，从而生成标准心拍模板的技术方案。
[0024]聚类分析是指将物理或抽象对象的集合分组成为由类似的对象组成的多个类的分析过程，所以同一个簇中的对象有很大的相似性，而不同簇间的对象有很大的相异性。聚类分析的目标就是在相似的基础上收集数据来分类。
[0025]传统的聚类算法可以被分为五类:划分方法、层次方法、基于密度方法、基于网格方法和基于模型方法。
[0026]在本申请中，仅以基于密度的聚类方法为例对本发明的心电波形的检测方法和装置进行描述。但是，本领域的普通技术人员可以理解，本发明提出的心电波形的检测技术不限于通过基于密度的聚类方法对心电波形进行分析和聚类，而是可以通过任何基于心电波形的特征进行聚类的算法对心电波形进行聚类。
[0027]下面结合附图对本发明实施例提供的心电波形的检测方法和装置进行详细描述。
[0028]图2为本发明实施例提供的心电波形的检测方法流程图。如图2所示，该方法包括以下步骤:
[0029]步骤S201，获取M个心拍波形数据，M为正整数。
[0030]这里，从心电波形图中获取连续的M个心拍波形数据。
[0031]步骤S202，根据聚类方法将所述M个心拍波形数据聚类成至少一个第一簇。稍后将具体说明步骤S202的处理。
[0032]步骤S203，从所述至少一个第一簇当中选取心拍波形数据的个数最多的第一簇。
[0033]步骤S204，基于选取的第一簇生成标准心拍模板。
[0034]以下将详细描述根据本发明的示例性实施例在步骤S202的处理。
[0035]基于密度的聚类算法(Density-BasedSpatial Clustering of Applicat1nwith Noise, DBSCAN)它将簇定义为密度相连的点的最大集合，能够把具有足够高密度的区域划分为簇，可在噪声的空间数据库中发现任意形状的簇。
[0036]根据本发明的示例性实施例，步骤S202包括两轮聚类处理:
[0037]首先，遍历所述M个心拍波形数据，根据心拍间的距离将所述M个心拍波形数据聚类成至少一个第二簇，以此将波形相近(密度可达)的心拍波形数据初步聚类到一起；
[0038]其次，再次遍历所述M个心拍波形数据，根据所述M个心拍波形数据当中的各心拍与所述至少一个第二簇中的各心拍之间的距离将所述至少一个第二簇进一步聚类成至少一个第一簇，从而由已被聚类的心拍波形数据通过心拍间的距离进行扩展，将密度相连的心拍波形数据也聚类到一起。
[0039]具体地，前述遍历所述M个心拍波形数据，根据心拍间的距离将所述M个心拍波形数据聚类成至少一个第二簇的处理包括，重复执行以下处理，直至全部所述M个心拍波形数据都已被标记为止:
[0040]从所述M个心拍波形数据当中选取任一尚未被标记的心拍波形数据作为当前的心拍波形数据，并且标记当前的心拍波形数据；
[0041]从尚未被标记的心拍波形数据当中选取与所述当前的心拍波形数据的距离在预定的距离阈值内(即密度可达)的心拍波形数据，并且标记所述选取的心拍波形数据；
[0042]如果所述选取的心拍波形数据的个数大于预定的阈值，则将所述当前的心拍波形数据和选取的心拍波形数据聚类成一个第二簇。
[0043]这里，通过计算的距离是否满足预定的条件来判定心拍波形数据间是否满足密度可达的条件。
[0044]此外，如果通过前述处理后，在从任一心拍波形数据聚类得到的心拍波形数据的个数不大于预定的阈值，也就是说，初步聚类得到的心拍波形数据的个数相对少，则将所述聚类得到的心拍波形数据作为噪音处理，丢弃所述心拍波形数据。
[0045]为了再次对所述M个心拍波形数据进行遍历，将所述M个心拍波形数据去除标记。然后，执行前述再次遍历的处理，包括，重复执行以下处理，直至全部所述M个心拍波形数据都已被标记为止:
[0046]从所述至少一个第二簇选取任一尚未被标记的心拍波形数据作为当前的心拍波形数据，并且标记当前的心拍波形数据；
[0047]从尚未被标记的心拍波形数据当中选取与所述当前的心拍波形数据的距离在预定的距离阈值内的心拍波形数据，将选取的心拍波形数据聚类到当前的心拍波形数据所在的第二簇内，并且标记所述选取的心拍波形数据。
[0048]通过所述再次遍历的处理，可从初步聚类得到的至少一个第二簇当中的所有心拍波形数据通过心拍间的距离条件进行扩展，从而将自第二簇密度相连的心拍波形数据进一步聚类到一起。
[0049]这里所述的心拍间的距离可以是，但不限于，心拍波形数据之间的欧几里得距离、马氏距离、岗可夫斯基距离、余弦相似度、皮尔逊相关系数等。
[0050]可例如，通过以下方式计算所述心拍间的距离。以图1中QRS波形中的特征点(如R点)为中心，左边按照等距离取预定个数的点(如取12个点)，右边按照等距离取预定个数的点(如取23个点)，再加上所述特征点(如R点)，由此为每个心拍波形数据建立多维向量(如前述的36维向量)。为了计算两个心拍间的距离，以X= (X1, X2,…Xn)T，Y =(Y1,Y2,…Υη)τ分别表示两个心拍的高维向量，然后可根据任何公知的距离计算公式来计算所述两个心拍间的距离。
[0051]本领域的普通技术人员可以理解，以上计算心拍间的距离的方式仅是示例性的，计算使用的变量不仅限于以R点为中心的36个点，而可以是其他可用于比较心拍波形数据相似性的指标。
[0052]在通过步骤S201?S203的处理后，获得了可代表正常心拍的簇。根据基于密度的聚类方法，在聚类的簇中处于密度核心的心拍波形数据可认为是最标准的心拍波形数据。因此，根据本发明的优选实施例，步骤S204包括:从所述选取的第一簇选取从密度上处于核心的心拍波形数据，并且将所述处于核心的心拍波形数据作为所述标准心拍模板。
[0053]根据本发明的可选实施例，在步骤S204，对所述选取的第一簇当中的心拍波形数据的各个指标求平均，并且基于求得的平均值生成所述标准心拍模板。
[0054]本发明实施例提供的心电波形的检测方法根据聚类方法将趋于正常的心拍波形数据聚类到一起，并且根据聚类得到的包括心拍个数最多的簇来生成标准心拍模板，从而通过计算机自动地生成标准心拍模板，并且具有较好的准确性。
[0055]图3为本发明实施例提供的心电波形的检测装置示意图。如图3所示，该装置包括:获取单元100、第一聚类单元200、第一选取单元300和生成单元400。
[0056]获取单元100，用于获取M个心拍波形数据，M为正整数；
[0057]第一聚类单元200，用于根据聚类方法将所述M个心拍波形数据聚类成至少一个第一簇；
[0058]第一选取单元300，用于从所述至少一个第一簇当中选取心拍波形数据的个数最多的第一簇；
[0059]生成单元400，用于基于选取的第一簇生成标准心拍模板。
[0060]具体地，第一聚类单元200可进一步包括第二聚类单元210和第二聚类单元220。
[0061]第二聚类单元210用于遍历所述M个心拍波形数据，根据心拍间的距离将所述M个心拍波形数据聚类成至少一个第二簇；
[0062]第三聚类单元220用于遍历所述M个心拍波形数据，根据所述M个心拍波形数据当中的各心拍与所述至少一个第二簇中的各心拍之间的距离将所述至少一个第二簇进一步聚类成至少一个第一簇。
[0063]根据本发明的示例性实施例，第二聚类单元210用于重复执行以下处理，直至全部所述M个心拍波形数据都已被标记为止:
[0064]从所述M个心拍波形数据当中选取任一尚未被标记的心拍波形数据作为当前的心拍波形数据，并且标记当前的心拍波形数据；
[0065]从尚未被标记的心拍波形数据当中选取与所述当前的心拍波形数据的距离在预定的距离阈值内的心拍波形数据，并且标记所述选取的心拍波形数据；
[0066]如果所述选取的心拍波形数据的个数大于预定的阈值，则将所述当前的心拍波形数据和选取的心拍波形数据聚类成一个第二簇。
[0067]根据本发明的示例性实施例，第三聚类单元220用于将所述M个心拍波形数据去除标记，然后重复执行以下处理，直至全部所述M个心拍波形数据都已被标记为止:
[0068]从所述至少一个第二簇选取任一尚未被标记的心拍波形数据作为当前的心拍波形数据，并且标记当前的心拍波形数据。
[0069]从尚未被标记的心拍波形数据当中选取与所述当前的心拍波形数据的距离在预定的距离阈值内的心拍波形数据，将选取的心拍波形数据聚类到当前的心拍波形数据所在的第二簇内，并且标记所述选取的心拍波形数据。
[0070]根据本发明的示例性实施例，生成单元400进一步包括:
[0071]第二选取单元(未示出)，用于从所述选取的第一簇选取从密度上处于核心的心拍波形数据。
[0072]第一确定单元(未示出)，将所述处于核心的心拍波形数据作为所述标准心拍模板。
[0073]根据本发明的另一示例性实施例，生成单元400包括:
[0074]平均单元(未示出)，用于对所述选取的第一簇当中的心拍波形数据的各个指标求平均。
[0075]第二确定单元(未示出)，用于基于求得的平均值生成所述标准心拍模板。
[0076]本发明实施例提供的心电波形的检测方法和装置，根据聚类方法将趋于正常的心拍波形数据聚类到一起，并且根据聚类得到的包括心拍个数最多的簇来生成标准心拍模板，从而通过计算机自动地生成标准心拍模板，并且具有较好的准确性。
[0077]需要指出，根据实施的需要，可将本申请中描述的各个步骤拆分为更多步骤，也可将两个或多个步骤或者步骤的部分操作组合成新的步骤，以实现本发明的目的。
[0078]上述根据本发明的方法可在硬件、固件中实现，或者被实现为可存储在记录介质(诸如CD ROM、RAM、软盘、硬盘或磁光盘)中的软件或计算机代码，或者被实现通过网络下载的原始存储在远程记录介质或非暂时机器可读介质中并将被存储在本地记录介质中的计算机代码，从而在此描述的方法可被存储在使用通用计算机、专用处理器或者可编程或专用硬件(诸如ASIC或FPGA)的记录介质上的这样的软件处理。可以理解，计算机、处理器、微处理器控制器或可编程硬件包括可存储或接收软件或计算机代码的存储组件(例如，RAM、R0M、闪存等)，当所述软件或计算机代码被计算机、处理器或硬件访问且执行时，实现在此描述的处理方法。此外，当通用计算机访问用于实现在此示出的处理的代码时，代码的执行将通用计算机转换为用于执行在此示出的处理的专用计算机。
[0079]以上所述，仅为本发明的【具体实施方式】，但本发明的保护范围并不局限于此，任何熟悉本【技术领域】的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以所述权利要求的保护范围为准。
【权利要求】
1.一种通过计算机实现的心电波形的检测方法，其特征在于，所述方法包括: 获取M个心拍波形数据，M为正整数； 根据聚类方法将所述M个心拍波形数据聚类成至少一个第一簇； 从所述至少一个第一簇当中选取心拍波形数据的个数最多的第一簇； 基于选取的第一簇生成标准心拍模板。
2.根据权利要求1所述的方法，其特征在于，所述根据聚类方法将所述M个心拍波形数据聚类成至少一个第一簇的处理包括: 遍历所述M个心拍波形数据，根据心拍间的距离将所述M个心拍波形数据聚类成至少一个第二簇； 遍历所述M个心拍波形数据，根据所述M个心拍波形数据当中的各心拍与所述至少一个第二簇中的各心拍之间的距离将所述至少一个第二簇进一步聚类成至少一个第一簇。
3.根据权利要求2所述的方法，其特征在于，所述遍历所述M个心拍波形数据，根据心拍间的距离将所述M个心拍波形数据聚类成至少一个第二簇的处理包括，重复执行以下处理，直至全部所述M个心拍波形数据都已被标记为止: 从所述M个心拍波形数据当中选取任一尚未被标记的心拍波形数据作为当前的心拍波形数据，并且标记当前的心拍波形数据； 从尚未被标记的心拍波形数据当中选取与所述当前的心拍波形数据的距离在预定的距离阈值内的心拍波形数据，并且标记所述选取的心拍波形数据； 如果所述选取的心拍波形数据的个数大于预定的阈值，则将所述当前的心拍波形数据和选取的心拍波形数据聚类成一个第二簇。
4.根据权利要求3所述的方法，其特征在于，所述遍历所述M个心拍波形数据，根据所述M个心拍波形数据当中的各心拍与所述至少一个第二簇中的各心拍之间的距离将所述至少一个第二簇进一步聚类成至少一个第一簇的处理包括，将所述M个心拍波形数据去除标记，然后重复执行以下处理，直至全部所述M个心拍波形数据都已被标记为止: 从所述至少一个第二簇选取任一尚未被标记的心拍波形数据作为当前的心拍波形数据，并且标记当前的心拍波形数据； 从尚未被标记的心拍波形数据当中选取与所述当前的心拍波形数据的距离在预定的距离阈值内的心拍波形数据，将选取的心拍波形数据聚类到当前的心拍波形数据所在的第二簇内，并且标记所述选取的心拍波形数据。
5.根据权利要求1?4中任一项所述的方法，其特征在于，所述基于选取的第一簇生成标准心拍模板的处理包括: 从所述选取的第一簇选取从密度上处于核心的心拍波形数据； 将所述处于核心的心拍波形数据作为所述标准心拍模板。
6.根据权利要求1?4中任一项所述的方法,其特征在于，所述基于选取的第一簇生成标准心拍模板的处理包括: 对所述选取的第一簇当中的心拍波形数据的各个指标求平均； 基于求得的平均值生成所述标准心拍模板。
7.—种心电波形的检测装置，其特征在于，所述装置包括: 获取单元，用于获取M个心拍波形数据，M为正整数； 第一聚类单元，用于根据聚类方法将所述M个心拍波形数据聚类成至少一个第一簇； 第一选取单元，用于从所述至少一个第一簇当中选取心拍波形数据的个数最多的第一簇； 生成单兀，用于基于选取的第一簇生成标准心拍模板。
8.根据权利要求7所述的装置，其特征在于，所述第一聚类单元包括: 第二聚类单元，用于遍历所述M个心拍波形数据，根据心拍间的距离将所述M个心拍波形数据聚类成至少一个第二簇； 第三聚类单元，用于遍历所述M个心拍波形数据，根据所述M个心拍波形数据当中的各心拍与所述至少一个第二簇中的各心拍之间的距离将所述至少一个第二簇进一步聚类成至少一个第一簇。
9.根据权利要求8所述的装置，其特征在于，所述第二聚类单元用于重复执行以下处理，直至全部所述M个心拍波形数据都已被标记为止: 从所述M个心拍波形数据当中选取任一尚未被标记的心拍波形数据作为当前的心拍波形数据，并且标记当前的心拍波形数据； 从尚未被标记的心拍波形数据当中选取与所述当前的心拍波形数据的距离在预定的距离阈值内的心拍波形数据，并且标记所述选取的心拍波形数据； 如果所述选取的心拍波形数据的个数大于预定的阈值，则将所述当前的心拍波形数据和选取的心拍波形数据聚类成一个第二簇。
10.根据权利要求9所述的装置，其特征在于，所述第三聚类单元用于将所述M个心拍波形数据去除标记，然后重复执行以下处理，直至全部所述M个心拍波形数据都已被标记为止: 从所述至少一个第二簇选取任一尚未被标记的心拍波形数据作为当前的心拍波形数据，并且标记当前的心拍波形数据； 从尚未被标记的心拍波形数据当中选取与所述当前的心拍波形数据的距离在预定的距离阈值内的心拍波形数据，将选取的心拍波形数据聚类到当前的心拍波形数据所在的第二簇内，并且标记所述选取的心拍波形数据。
【文档编号】A61B5/0452GK104398256SQ201410640514
【公开日】2015年3月11日 申请日期:2014年11月13日 优先权日:2014年11月13日
【发明者】侯桂文 申请人:北京海思敏医疗技术有限公司