CN111631704A

CN111631704A - 基于心电与脑电信息结合的糖尿病前期检测系统和方法

Info

Publication number: CN111631704A
Application number: CN202010475003.6A
Authority: CN
Inventors: 聂泽东; 李景振; 刘宇航
Original assignee: Shenzhen Institute of Advanced Technology of CAS
Current assignee: Shenzhen Institute of Advanced Technology of CAS
Priority date: 2020-05-29
Filing date: 2020-05-29
Publication date: 2020-09-08
Anticipated expiration: 2040-05-29
Also published as: WO2021238092A1; US20220313172A1; CN111631704B

Abstract

本发明公开了一种基于心电与脑电信息结合的糖尿病前期检测系统和方法。该系统包括：信号获取模块，其利用可穿戴设备以无创方式同步获取用户的心电信号和脑电信号；特征提取模块，其利用多种方式对由心电特征和脑电特征构成的组合特征集进行降维处理，获得多个降维的组合特征集，并通过分析多个降维的组合特征集与血糖浓度值的相关性选出满足设定相关性标准的心电特征和脑电特征，构成优化的组合特征集；多模融合模块，用于将优化的组合特征集分别输入至经训练的多种类型的神经网络模型，通过融合多种类型的神经网络的输出结果，获得用户糖尿病前期的检测结果。本发明能够提供无创、无疼痛、便捷舒适、低成本的糖尿病前期检测方案。

Description

基于心电与脑电信息结合的糖尿病前期检测系统和方法

技术领域

本发明涉及医疗健康技术领域，更具体地，涉及一种基于心电与脑电信息结合的糖尿病前期检测系统和方法。

背景技术

在由血糖正常的健康群体，逐渐发展成糖尿病人群的过程中存在一个时期，即糖尿病前期。糖尿病前期是指血糖调节功能受损,包括空腹血糖受损以及糖耐量受损，但是尚未达到糖尿病诊断标准的一段时期。据统计，有约25％的年轻人,以及约20％的青少年都患有糖尿病前期。对处于糖尿病前期的人，如果不进行干预，每年会有10％左右进展为糖尿病；而如果在糖尿病前期能及时采用相应措施进行干预，例如使用药物、控制其饮食、加强运动等，发生糖尿病的危险能够下降30％-75％，恢复为正常血糖状况的几率最高可升至70％左右。因此，对糖尿病前期的检测具有非常重要的意义。

目前糖尿病前期的判定方法包括：1)在空腹情况下通过采血获取血糖浓度值，如果空腹血糖值在5.6mmol/L至7.0mmol/L之间，则可判定为糖尿病前期；2)开展口服葡萄糖耐量试验，口服葡萄糖2小时后，通过采血获取血糖浓度，如果血糖值在7.8mmol/L至11.1mmol/L，则判定为糖尿病前期。然而，现有技术方案均需采集静脉血或者指尖血，会给患者带来较大疼痛和感染的风险，并且检测成本较高。

发明内容

本发明的目的是克服上述现有技术的缺陷，提供一种基于心电和脑电信息相结合的糖尿病前期检测方法，通过开展口服葡萄糖耐量试验，并利用可穿戴设备同步获取心电和脑电信息，进而提取心电和脑电的相关特征，实现糖尿病前期的检测。

根据本发明的第一方面，提供一种基于心电与脑电信息结合的糖尿病前期检测系统。该系统包括：

信号获取模块：利用可穿戴设备以无创方式同步获取用户的心电信号和脑电信号；

特征提取模块：利用多种方式对由心电特征和脑电特征构成的组合特征集进行降维处理，获得多个降维的组合特征集，并通过分析所述多个降维的组合特征集与血糖浓度值的相关性筛选出满足设定相关性标准的心电特征和脑电特征，构成优化的组合特征集；

多模融合模块：用于将所述优化的组合特征集分别输入至经训练的多种类型的神经网络模型，通过融合所述多种类型的神经网络的输出结果，获得用户糖尿病前期的检测结果。

在一个实施例中，所述同步获取用户的心电信号和脑电信号包括：

将用于监测心电信号的六个电极片V1至V6分别放置在用户胸部，其中V1设置在胸骨右缘第四肋间，V2设置在胸骨左缘第四肋间，V3设置在V2与V4之间连线的中点，V4设置在左锁骨中线与第五肋间交叉处，V5水平于腋前线，V6水平于腋中线；

用户在头部佩戴脑电电极帽，该脑电电极帽中分布设置用于监测脑电信号的六个电极，分别对应大脑左半球和右半球的额叶、枕叶、顶叶；

开展葡萄糖耐量试验，并开启心电采集设备和脑电采集设备，以同步获取用户的心电信号和脑电信号。

在一个实施例中，所述特征提取模块执行以下过程：

从所述心电信号中提取多个不同片段的特征信息，从所述脑电信号中分别提取对应大脑不同位置的不同频段的多个脑电特征信息，构成所述组合特征集；

利用主成分分析对所述组合特征集进行降维处理，获得第一组合特征集；

利用独立成分分析对所述组合特征集进行降维处理，获得第二组合特征集；

利用套索回归分析对所述组合特征集进行降维处理，获得第三组合特征集；

分别分析所述第一组合特征集、所述第二组合特征集和所述第三组合特征集与血糖浓度的相关性，进而筛选出满足设定相关性标准的心电特征和脑电特征，构成所述优化的组合特征集。

在一个实施例中，利用Pearson相关分析法分析相关性，所述相关性标准设置为相关性k＞0.2且P≤0.05，P表示对相关性系数进行假设检验的概率。

在一个实施例中，所述利用主成分分析对所述组合特征集进行降维处理包括：

计算所述组合特征集中各特征的特征点的协方差矩阵；

计算所述协方差矩阵的特征向量和对应的特征值：

将特征向量按特征值的大小排序，组成矩阵u＝[u₁,u₂,u₃,...,u_n]，对应的特征值由大到小分别是λ₁,λ₂,λ₃,...,λ_n，从矩阵u中截取排序靠前的设定比例的值作为每个特征的新特征点，以实现数据降维。

在一个实施例中，所述多种类型的神经网络模型包括支持向量机、随机森林、卷积神经网络、长短期记忆网络、递归神经网络中的至少两种类型。

在一个实施例中，所述多模融合模块采用投票法融合所述多种类型的神经网络的输出结果，获得用户是否属于糖尿病前期的检测结果。

根据本发明的第二方面，提供一种基于心电与脑电信息结合的糖尿病前期检测方法。该方法包括以下步骤：

利用可穿戴设备以无创方式同步获取用户的心电信号和脑电信号；

利用多种方式对由心电特征和脑电特征构成的组合特征集进行降维处理，获得多个降维的组合特征集，并通过分析所述多个降维的组合特征集与血糖浓度值的相关性筛选出满足设定相关性标准的心电特征和脑电特征，构成优化的组合特征集；

将所述优化的组合特征集分别输入至经训练的多种类型的神经网络模型，通过融合所述多种类型的神经网络的输出结果，获得用户糖尿病前期的检测结果。

与现有技术相比，本发明的优点在于，具有无创、无疼痛、便捷舒适、低成本、实时监测等优势，能够广泛运用于儿童、青少年、老人等不同群体的糖尿病前期检测。

通过以下参照附图对本发明的示例性实施例的详细描述，本发明的其它特征及其优点将会变得清楚。

附图说明

被结合在说明书中并构成说明书的一部分的附图示出了本发明的实施例，并且连同其说明一起用于解释本发明的原理。

图1是根据本发明一个实施例的基于心电和脑电信息相结合的糖尿病前期检测系统的示意图；

图2是根据本发明一个实施例的从心电信号中获取其不同片段的特征信息的示意图；

图3是根据本发明一个实施例的所选取的心电特征和脑电特征的示意图。

具体实施方式

现在将参照附图来详细描述本发明的各种示例性实施例。应注意到：除非另外具体说明，否则在这些实施例中阐述的部件和步骤的相对布置、数字表达式和数值不限制本发明的范围。

以下对至少一个示例性实施例的描述实际上仅仅是说明性的，决不作为对本发明及其应用或使用的任何限制。

对于相关领域普通技术人员已知的技术、方法和设备可能不作详细讨论，但在适当情况下，所述技术、方法和设备应当被视为说明书的一部分。

在这里示出和讨论的所有例子中，任何具体值应被解释为仅仅是示例性的，而不是作为限制。因此，示例性实施例的其它例子可以具有不同的值。

应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步讨论。

参见图1所示，本发明实施例提供的基于心电与脑电信息结合的糖尿病前期检测系统包括信号获取模块110、特征提取模块120、多模融合模块130和结果输出模块140。

信号获取模块110用于采集用户(或称被测者)的心电信号和脑电信号。例如，利用可穿戴设备实现同步获取用户在开展口服葡萄糖耐量试验时的心电和脑电信号。在一个实施例中，同步获取心电和脑电信号的具体过程是：首先，将六个电极片分别放置在用户胸部，将电极位置(或电极标识)标记为V1-V6，分别是胸骨右缘第四肋间V1；胸骨左缘第四肋间V2；V3的位置是V2和V4之间连线的中点；V4位于左锁骨中线和第五肋间交叉处；V5水平于腋前线；V6水平于腋中线。并且，用户在头部佩戴一个一次性脑电电极帽，电极帽中分布有6个电极用于监测脑电信息，6个电极分别位于大脑左半球及右半球的额叶、枕叶、顶叶。接下来，在用户佩戴好心电电极及脑电电极且检查无误后，静坐10分钟，待气息平稳后，开展口服葡萄糖耐量试验，例如服用75克葡萄糖。与此同时，开启心电采集设备以及脑电采集设备，实现同步获取用户的心电和脑电信号。2小时后，停止采集心电和脑电信号，并将数据导出保存。

特征提取模块120用于实现心电特征和脑电特征的降维，并选择出与血糖浓度变化最为相关的一些心电特征及脑电特征用于后续分析，以提高计算速度。

特征提取模块120可通过多种方式实现降维，例如利用主成分分析、独立成分分析、套索回归分析等算法，进而通过对多种方式降维后的特征进行相关性分析，选出与血糖浓度变化最为相关的一些心电特征及脑电特征。

在一个实施例中，特征提取模块120的具体实现过程包括：

步骤S210，从心电信号中提取心电特征。

例如，对心电信号进行分析，从心电信号中提取其不同片段的特征信息，包括RRI，R-H，P-H，QRS，PRQ，QT，QTC，ST，HR等9个不同片段的特征信息，部分片段信息的标注如图2所示，详细含义如下表1。

表1：心电信号所提取的不同片段的特征信息

心电信号的不同片段特征信息名称	含义
		RRI	心电信号中相邻两个R点的时间长度
R-H	心电信号中R波的高度
		P-H	心电信号中P波的高度
QRS	QRS波群的周期长度
		PRQ	心电信号中R点到P点的时间长度
QT	从Q波开始到T波结束的时间长度
		QTC	经过校准后的QT值
ST	S点到T点中的直线段的时间距离
		HR	一分钟内心跳的次数

步骤S220，从脑电信号中提取脑电特征。

具体地，对脑电信号进行分析，分别提取大脑左半球和右半球的额叶、枕叶、顶叶的不同频段的信号特征，例如，共提取30个脑电特征信息，如下表2。

表2：脑电信号所提取的不同频段下的特征信息

步骤S230，采用多种方式对心电特征和脑电特征集进行降维。

具体地，考虑到所提取的心电特征和脑电特征维度过长，为了实现快速计算，可采用主成分分析、独立成分分析、套索回归分析三种算法分别对心电特征和脑电特征数据进行降维处理。

在一个实施例中，基于主成分分析的降维过程如下：

以获取了共39个特征(包括9个心电特征和30个脑电特征)为例，假设每个特征有n个特征点，则所有的心电、脑电特征集(或称组合特征集)可表示为：

X＝{x⁽¹⁾,x⁽²⁾,x⁽³⁾,...,x⁽³⁹⁾} (1)

其中第i个特征可表示为：

x⁽ⁱ⁾＝(x₁ ⁽ⁱ⁾,x₂ ⁽ⁱ⁾,x₃ ⁽ⁱ⁾,...,x_n ⁽ⁱ⁾)^T (2)

则特征集的协方差矩阵为：

其中，协方差矩阵为n×n大小的方阵，具有n个特征点。

计算协方差矩阵的特征向量及对应的特征值：

Uu＝λu (4)

进一步地，将特征向量按特征值的大小按列排放，组成矩阵u＝[u₁,u₂,u₃,...,u_n]，对应的特征值由大到小分别为：λ₁,λ₂,λ₃,...,λ_n，则特征向量u₁为主特征向量(对应的特征值最大)，u₂为次特征向量，以此类推。通过截取矩阵u中例如前10％的值作为每个特征的新特征点，从而实现数据的降维。

在一个实施例中，基于独立成分分析的降维过程如下：

同样假设所有的心电、脑电特征集可表示为

X＝{x⁽¹⁾,x⁽²⁾,x⁽³⁾,...,x⁽³⁹⁾} (5)

其中第i个特征可表示为：

x⁽ⁱ⁾＝(x₁ ⁽ⁱ⁾,x₂ ⁽ⁱ⁾,x₃ ⁽ⁱ⁾,...,x_n ⁽ⁱ⁾)^T (6)

假设经过独立成分分析降维后的特征集可表示为：

Q＝{q⁽¹⁾,q⁽²⁾,q⁽³⁾,...,q⁽³⁹⁾} (7)

根据独立成分分析，则有：

X＝AQ (8)

令W＝A^-1，则q⁽ⁱ⁾＝A^-1x⁽ⁱ⁾＝Wx⁽ⁱ⁾。

其中

(即将W表示成另一种形式)，α，g为比例系数，其值大小采用极大似然法确定。当确定W矩阵的值后，则可计算出Q矩阵的值，该值即为经过独立成分分析降维后的特征集。

在一个实施例中，基于套索回归分析的降维过程如下：

同样假设所有的心电、脑电特征集可表示为

X＝{x⁽¹⁾,x⁽²⁾,x⁽³⁾,...,x⁽³⁹⁾} (9)

其中第i个特征可表示为：

x⁽ⁱ⁾＝(x₁ ⁽ⁱ⁾,x₂ ⁽ⁱ⁾,x₃ ⁽ⁱ⁾,...,x_n ⁽ⁱ⁾)^T (10)

假设经过套索回归分析降维后的特征集可表示为：

Y＝{y⁽¹⁾,y⁽²⁾,y⁽³⁾,...,y⁽³⁹⁾} (11)

套索回归分析过程可以视为凸优化的求解过程，即：

其中J为代价函数，K为待求解的相关矩阵，通过利用含有范数的拉格朗日函数求解方法，解析出K矩阵的大小，则可以得到J函数的表达式。

步骤S240，对采用多种方式进行降维后的组合特征集进行相关性分析，进而选取满足设定相关性标准的心电特征和脑电特征构成优化的组合特征集。

具体地，对于利用主成分分析得到降维后的心电、脑电特征集u，利用独立成分分析得到降维后的心电、脑电特征集Q，以及利用套索回归分析降维后的心电、脑电特征集Y，利用Pearson相关分析方法分别分析上述降维后的3个心电、脑电特征集与血糖浓度的相关性，即分别分析血糖浓度值与RRI，R-H，P-H，QRS，PRQ，QT，QTC，ST，HR等9个心电特征以及F1_T，F1_D，F1_A，F1_B，F1_G，O1_T，O1_D，O1_A，O1_B，O1_G等30个脑电特征的相关性。例如，如果满足相关系数k＞0.2且P值(双尾)≤0.05(其中k用于表示相关程度，P是对相关系数进行假设检验的概率，P≤0.05表示有显著相关关系)，则认为血糖值与该特征具有较强的相关性，则该特征暂时保留。

相关分析完成后，通过选取在主成分分析、独立成分分析、套索回归分析等三种算法中均满足相关性(k＞0.2，P值(双尾)≤0.05)的心电、脑电特征，将这些特征作为可使用的特征，如图3所示，最终所选取的特征包括RRI，QT，QRS三个心电特征以及O1_A，O1_D，P2_A，P2_B，P2_T五个脑电特征。

在本发明实施例中，通过主成分分析选择特征点投影具有最大方差的方向，能够获得低维特征集，主成分分析适用于样本呈高斯分布的情况，而独立成分分析不要求样本呈高斯分布，套索回归分析能够有效进行数据降维并准确识别更重要的特征。通过采用多种降维方法并结合相关性分析进行特征筛选，在有效进行数据降维的同时，能够准确识别出与血糖浓度具有强相关性的特征，提高了后续处理速度并且不影响检测精度，并且对特征数据的分布不敏感，适用范围更广泛。

多模融合模块130用于基于选出的心电特征和脑电特征集，利用不同方式获得糖尿病前期的多个检测结果。

具体地，可利用支持向量机、随机森林、卷积神经网络、长短期记忆网络、递归神经网络等，并将所提取的RRI，QT，QRS三个心电特征以及O1_A，O1_D，P2_A，P2_B，P2_T五个脑电特征分别作为该五种网络模型的输入，选择合适的参数，经过多次迭代训练后，分别得出每种网络模型的输出结果，即是否属于糖尿病前期。进一步地，通过采用投票法获得融合结果，例如，如果五种网络模型有三种或三种以上的输出结果为糖尿病前期，则判定为该用户为糖尿病前期，否则判定该用户为非糖尿病前期。采用融合结果作为最终的检测结果，提高了检测准确性。

结果输出模块140，用于向用户显示最终的判断结果。例如，可通过语音、文本等显示检测结果。

综上所述，本发明利用可穿戴设备同步获取心电信号和脑电信号，通过对心电特征和脑电特征进行降维处理并选择优化的组合特征集，进而基于多种算法融合的糖尿病前期检测方法，能够为用户提供无创、无疼痛、便捷舒适、低成本的糖尿病前期检测方法，可广泛运用于儿童、青少年、老人等不同群体的糖尿病前期检测。

本发明可以是系统、方法和/或计算机程序产品。计算机程序产品可以包括计算机可读存储介质，其上载有用于使处理器实现本发明的各个方面的计算机可读程序指令。

计算机可读存储介质可以是可以保持和存储由指令执行设备使用的指令的有形设备。计算机可读存储介质例如可以是――但不限于――电存储设备、磁存储设备、光存储设备、电磁存储设备、半导体存储设备或者上述的任意合适的组合。计算机可读存储介质的更具体的例子(非穷举的列表)包括：便携式计算机盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、静态随机存取存储器(SRAM)、便携式压缩盘只读存储器(CD-ROM)、数字多功能盘(DVD)、记忆棒、软盘、机械编码设备、例如其上存储有指令的打孔卡或凹槽内凸起结构、以及上述的任意合适的组合。这里所使用的计算机可读存储介质不被解释为瞬时信号本身，诸如无线电波或者其他自由传播的电磁波、通过波导或其他传输媒介传播的电磁波(例如，通过光纤电缆的光脉冲)、或者通过电线传输的电信号。

这里所描述的计算机可读程序指令可以从计算机可读存储介质下载到各个计算/处理设备，或者通过网络、例如因特网、局域网、广域网和/或无线网下载到外部计算机或外部存储设备。网络可以包括铜传输电缆、光纤传输、无线传输、路由器、防火墙、交换机、网关计算机和/或边缘服务器。每个计算/处理设备中的网络适配卡或者网络接口从网络接收计算机可读程序指令，并转发该计算机可读程序指令，以供存储在各个计算/处理设备中的计算机可读存储介质中。

用于执行本发明操作的计算机程序指令可以是汇编指令、指令集架构(ISA)指令、机器指令、机器相关指令、微代码、固件指令、状态设置数据、或者以一种或多种编程语言的任意组合编写的源代码或目标代码，所述编程语言包括面向对象的编程语言—诸如Smalltalk、C++等，以及常规的过程式编程语言—诸如“C”语言或类似的编程语言。计算机可读程序指令可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。在涉及远程计算机的情形中，远程计算机可以通过任意种类的网络—包括局域网(LAN)或广域网(WAN)—连接到用户计算机，或者，可以连接到外部计算机(例如利用因特网服务提供商来通过因特网连接)。在一些实施例中，通过利用计算机可读程序指令的状态信息来个性化定制电子电路，例如可编程逻辑电路、现场可编程门阵列(FPGA)或可编程逻辑阵列(PLA)，该电子电路可以执行计算机可读程序指令，从而实现本发明的各个方面。

这里参照根据本发明实施例的方法、装置(系统)和计算机程序产品的流程图和/或框图描述了本发明的各个方面。应当理解，流程图和/或框图的每个方框以及流程图和/或框图中各方框的组合，都可以由计算机可读程序指令实现。

这些计算机可读程序指令可以提供给通用计算机、专用计算机或其它可编程数据处理装置的处理器，从而生产出一种机器，使得这些指令在通过计算机或其它可编程数据处理装置的处理器执行时，产生了实现流程图和/或框图中的一个或多个方框中规定的功能/动作的装置。也可以把这些计算机可读程序指令存储在计算机可读存储介质中，这些指令使得计算机、可编程数据处理装置和/或其他设备以特定方式工作，从而，存储有指令的计算机可读介质则包括一个制造品，其包括实现流程图和/或框图中的一个或多个方框中规定的功能/动作的各个方面的指令。

也可以把计算机可读程序指令加载到计算机、其它可编程数据处理装置、或其它设备上，使得在计算机、其它可编程数据处理装置或其它设备上执行一系列操作步骤，以产生计算机实现的过程，从而使得在计算机、其它可编程数据处理装置、或其它设备上执行的指令实现流程图和/或框图中的一个或多个方框中规定的功能/动作。

附图中的流程图和框图显示了根据本发明的多个实施例的系统、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上，流程图或框图中的每个方框可以代表一个模块、程序段或指令的一部分，所述模块、程序段或指令的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。在有些作为替换的实现中，方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如，两个连续的方框实际上可以基本并行地执行，它们有时也可以按相反的顺序执行，这依所涉及的功能而定。也要注意的是，框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合，可以用执行规定的功能或动作的专用的基于硬件的系统来实现，或者可以用专用硬件与计算机指令的组合来实现。对于本领域技术人员来说公知的是，通过硬件方式实现、通过软件方式实现以及通过软件和硬件结合的方式实现都是等价的。

以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择，旨在最好地解释各实施例的原理、实际应用或对市场中的技术改进，或者使本技术领域的其它普通技术人员能理解本文披露的各实施例。本发明的范围由所附权利要求来限定。

Claims

1.一种基于心电与脑电信息结合的糖尿病前期检测系统，包括：

2.根据权利要求1所述的系统，其中，所述同步获取用户的心电信号和脑电信号包括：

3.根据权利要求1所述的系统，其中，所述特征提取模块执行以下过程：

4.根据权利要求3所述的系统，其中，利用Pearson相关分析法分析相关性，所述相关性标准设置为相关性k＞0.2且P≤0.05，P表示对相关性系数进行假设检验的概率。

5.根据权利要求3所述的系统，其中，所述利用主成分分析对所述组合特征集进行降维处理包括：

计算所述组合特征集中各特征的特征点的协方差矩阵；

计算所述协方差矩阵的特征向量和对应的特征值：

6.根据权利要求1所述的系统，其中，所述多种类型的神经网络模型包括支持向量机、随机森林、卷积神经网络、长短期记忆网络、递归神经网络中的至少两种类型。

7.根据权利要求1所述的系统，其中，所述多模融合模块采用投票法融合所述多种类型的神经网络的输出结果，获得用户是否属于糖尿病前期的检测结果。

8.一种基于心电与脑电信息结合的糖尿病前期检测方法，包括以下步骤：

9.一种计算机可读存储介质，其上存储有计算机程序，其中，该程序被处理器执行时实现以下步骤：