一种一维序列升维聚类方法及系统

技术领域

本发明属于人工智能图像处理技术领域，具体涉及一种一维序列升维聚类方法及系统。

背景技术

心电信号(Electrocardiogram,ECG)是对心脏活动的直观记录，通过心电图展现的心率、S-T段、P波、QRS波群的形态及出现位置等特征能够对于窦性心动过速、窦性心律不齐、室性早搏和心房颤动等心律失常进行诊断，是进行心脏疾病判别的主要依据之一。过去，心电信号的判别主要靠医生人为判断，随着医疗标准的不断提高，人工判别的速度与准确度已不足以满足医疗需求。本技术能够极大程度的展现数据隐性特征，提高心电信号识别速度与识别准确率，解决由于数据不足而导致的模型识别能力差的问题，同时由于使用无监督学习，极大地减少了标注成本。同时拥有较强的泛化能力，可以适应不同病人的病理分类。

现有技术存在的问题或缺陷：1、提取特征单一，无法准确描述信号所反映的信息，导致信号分类准确率较低；2、样本较少导致准确率较低，深度学习模型已被证明在数据分类方面具有良好的效果，但深度学习模型的训练对数据规模要求较高，而大多序列信号标注数据较少，会导致模型训练不完全，模型欠拟合，识别准确率较低。目前的一维序列分类方法存在序列分类准确率低、分类模泛化能力差、算力要求过高等问题。

发明内容

针对上述传统心电信号的处理提取特征单一、识别准确率较低的技术问题，本发明提供了一种分类准确率高、特征提取容易、成本低、自动化程度高的一维序列升维聚类方法及系统。

为了解决上述技术问题，本发明采用的技术方案为：

一种一维序列升维聚类方法，包括下列步骤：

S1、一维心电信号收集；

S2、一维心电信号去噪；

S3、对去噪后的一维心电信号进行格罗姆变换；

S4、导入无监督K-means聚类算法；

S5、输出结果，并储存。

所述S1中一维心电信号收集的方法为：导入一维心电信号，将其分割为长度为n的心电信号段，并储存为NPY文件供模型读取。

所述S2中一维心电信号去噪的方法为：首先对一维心电信号按照一维时间序列排序X＝x

所述S2中一维心电信号去噪采用滤波器，所述滤波器的截止频率为0.08，所述滤波器采用阶数为4的低通滤波器。

所述S3中对去噪后的一维心电信号进行格罗姆变换的方法为：将去噪后的一维心电信号导入格罗姆变换模块，格罗姆变换模块的格罗姆角场通过对时间序列进行内积来衡量数据内部的相似性，将一维心电信号转换为格罗姆矩阵，实现一维心电信号转变为二维心电图像。

所述格罗姆变换转换过程为：一维心电信号进行缩放和去噪后，计算每一对的点积并放入格罗姆矩阵G中，

所述S4中导入无监督K-means聚类算法的方法为：首先随机指定类中心，根据所述二维心电图像距离类中心的远近进行类簇的划分，划分完成之后，重新对类中心进行计算，经过循环迭代，直到类中心收敛，完成类簇划分。

所述对类中心进行计算的方法为；通过EM算法的两步走策略计算出类中心，求取最小化平方误差函数E，求取最小化平方误差函数E的公式为：

一种一维序列升维聚类系统，包括心电信号收集模块、预处理模块、去噪模块、格罗姆变换模块、无监督K-means聚类模块、输出结果模块，所述心电信号收集模块与预处理模块，所述预处理模块连接有去噪模块，所述去噪模块连接有格罗姆变换模块，所述格罗姆变换模块连接有无监督K-means聚类模块，所述无监督K-means聚类模块连接有输出结果模块。

本发明与现有技术相比，具有的有益效果是：

本发明将心电信号的识别从时域问题转换图像分类问题，通过对时间序列进行內积来衡量数据内部的相似性，将其转换为格罗姆矩阵，完成从一维到二维的升维过程。可以在保留数据特征的同时，展现出数据更多的隐藏特征，使得后续使用深度神经网络进行特征提取更加容易。本专利采用无监督K-means聚类算法，对数据量的要求较低。由于采用无监督算法，并不需要数据的标注，这在一定程度上节约了经济成本，且此方法只需在确立类簇中心时耗费一定算力，而在中心确立后，只需较少的计算成本即可完成数据的分类工作。同时真正做到了无人参与、智能识别。

附图说明

图1为本发明的工作流程图；

图2为本发明的算法流程总图；

图3为本发明的信号转化流程图；

图4为本发明的聚类流程图；

图5为本发明的结构原理图；

其中：1为心电信号收集模块，2为预处理模块，3为去噪模块，4为格罗姆变换模块，5为无监督K-means聚类模块，6为输出结果模块。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

一种一维序列升维聚类方法，如图1、图2、图3所示，包括下列步骤：

步骤一、心电信号收集；导入一维心电信号，将其分割为长度为n的心电信号段，并储存为NPY文件供模型读取；这里对心电信号进行分割的目的在于，统一输入模型数据的大小，方便模型进行识别。

步骤二、一维心电信号去噪；将所述一维心电信号的时间序列X＝x

步骤三、对去噪后的一维心电信号进行格罗姆变换；将去噪后的心电信号导入格罗姆变换模块，格罗姆角场通过对时间序列进行內积来衡量数据内部的相似性，将其转换为格罗姆矩阵，实现一维心电信号转变为二维心电图像，将心电信号的识别从时域问题转换图像分类问题，通过对时间序列进行內积来衡量数据内部的相似性，将其转换为格罗姆矩阵，完成从一维到二维的升维过程。可以在保留数据特征的同时，展现出数据更多的隐藏特征，使得后续使用深度神经网络进行特征提取更加容易。

步骤四、导入无监督K-means聚类算法；首先随机指定类中心，根据二维心电图像距离类中心的远近进行类簇的划分，划分完成之后，重新对类中心进行计算，经过循环迭代，直到类中心收敛，完成类簇划分；采用无监督K-means聚类算法，对数据量的要求较低。由于采用无监督算法，并不需要数据的标注，这在一定程度上节约了经济成本，且此方法只需在确立类簇中心时耗费一定算力，而在中心确立后，只需较少的计算成本即可完成数据的分类工作。同时真正做到了无人参与、智能识别。

步骤五、输出结果，并储存。

在一种可能的实施方式中，步骤步骤二中去噪时，首先对一维心电信号按照一维时间序列排序X＝x

在一种可能的实施方式中，S2中一维心电信号去噪采用滤波器，所述滤波器的截止频率为0.08，所述滤波器采用阶数为4的低通滤波器。

在一种可能的实施方式中，步骤三中格罗姆变换的具体转换过程如下：一维心电信号进行缩放和去噪后，计算每一对的点积并放入格罗姆矩阵G中，

在一种可能的实施方式中，步骤四中采用的类中心计算方法为；通过EM算法的两步走策略而计算出类中心，求取最小化平方误差函数E，最小化平方误差函数E公式为：

其中具体的聚类过程，参照图4，首先导入二维心电图像，其导入的二维心电图像数据量为k，先通过人为设定类中心μ

一种一维序列升维聚类系统，如图5所示，包括心电信号收集模块1、预处理模块2、去噪模块3、格罗姆变换模块4、无监督K-means聚类模块5、输出结果模块6，所述心电信号收集模块1与预处理模块2，所述预处理模块2连接有去噪模块3，所述去噪模块3连接有格罗姆变换模块4，所述格罗姆变换模块4连接有无监督K-means聚类模块5，所述无监督K-means聚类模块5连接有输出结果模块6。心电信号收集模块1导入一维心电信号，并储存于数据库中，并对一维心电信号依照导入数据库先后时间顺序从X

上面仅对本发明的较佳实施例作了详细说明，但是本发明并不限于上述实施例，在本领域普通技术人员所具备的知识范围内，还可以在不脱离本发明宗旨的前提下作出各种变化，各种变化均应包含在本发明的保护范围之内。