CN113796830A - 一种睡眠信号分期可信度自动评价方法 - Google Patents

Abstract

一种睡眠信号分期可信度自动评价方法，步骤一：选择主分类器和辅分类器，采用监督学习对主、辅分类器进行训练，获得训练后的主、辅分类器；步骤二：输入睡眠数据，利用步骤一中训练后的主、辅分类器，对睡眠数据进行多粒度分类视图预测，并输出主分类器的预测结果；步骤三：根据步骤二的主、辅分类器预测结果，计算并输出逻辑一致性系数即可信度R；本发明能够有效自动评价睡眠分期结果，对海量无标记睡眠数据自动标记来更新训练集规模来训练健壮性更高、结果可信度高的睡眠分期算法，同时本发明技术可以移植至可穿戴便携式睡眠监测设备中，辅助及促进移动医疗的发展。

Description

一种睡眠信号分期可信度自动评价方法

技术领域

本发明属于生物医学工程技术领域，特别涉及睡眠监测领域中的一种睡眠信号分期可信度自动评价方法。

背景技术

随着现代生活节奏的加快和生活方式的改变，睡眠日益成为一个突出医疗及公共卫生问题。中国睡眠研究会2021年发布睡眠调查报告显示，我国有超3亿人存在睡眠障碍，成年人失眠发生率高达38.2％，其中老年人睡眠障碍发生率占56.7％。临床上通过睡眠分期研究被试的睡眠结构和睡眠质量，对干预改善人们睡眠和辅助医生诊断治疗睡眠疾病至关重要。

现阶段睡眠医学领域中，多导睡眠图(polysomnography，PSG)是睡眠分期的“金标准”，该方法同步采集受试者睡眠过程中的脑电(electroencephalogram，EEG)、心电(electrocardiogram，EEG)、肌电(electromyogram，EMG)、眼电(electrooculogram，EOG)、呼吸等被试的整晚多通道生理信号，然后专业医生按照美国睡眠医学学会(americanacademy of sleep medicine,AASM)制定的睡眠分期标准将夜间睡眠以30s时间间隔为单位分为清醒期(wake，W)、非快速眼动睡眠期(non-rapid eye movement，NREM)和快速眼动睡眠期(rapid eye movement，REM)，其中NREM期又细分为Ⅰ期(N1)、Ⅱ期(N2)、Ⅲ期(N3)。由于PSG价格昂贵、受试者睡眠时需要佩戴大量传感器、需要搭建专业睡眠实验室及专业医生判读结果等缺陷，严重限制了其普及和推广，以及利用人工智能技术在睡眠领域开展海量数据规模的新应用。近年来，很多研究人员开始将注意力集中在通过可穿戴睡眠监测装置，利用人工智能技术进行自动睡眠分期算法研究。

自动睡眠分期算法主要包括生理信号预处理、特征提取和睡眠状态分类三个步骤，对易于采集的ECG、呼吸与光电容积脉搏波(photoplethysmography，PPG)等生理信号预处理后，提取与睡眠分期相关的时域、频域及非线性特征，然后对分类器进行模型训练。现有训练分类器的方法往往采用有监督技术，分类精度与训练集规模成正相关。但由于医疗数据的特殊性和人工标记的复杂性，导致可用的标记数据非常稀缺。利用数量有限的标记数据训练出的分期模型，其健壮性和结果可信度很难保证。对睡眠分期结果可信度进行有效评价具有重要意义和应用价值。

近年来，睡眠健康监测产品的市场化和家庭化，快速产生了海量的无标记睡眠数据。人工标记这些海量的无标记数据需要花费大量时间和人力，该类方案缺乏现实可行性。(1)如何对海量数据进行自动分析，形成可用于模型算法训练的有标记数据是一个巨大的挑战；(2)如何充分有效利用这些海量无标记数据，提高睡眠分期算法的健壮性和结果可信度是另一个需要解决的问题。

发明内容

为了克服上述现有技术的缺陷，本发明的目的在于提供一种睡眠信号分期可信度自动评价方法，利用睡眠多视图概念，对于获得的睡眠数据，按不同分类粒度对同一睡眠时刻分别进行分类，采用五分类方式，把睡眠状态分为W、REM、N1、N2、N3五个期，由于不同的分类视图是用不同的分类器得到的同一组数据的不同视图，使用训练好的多个不同粒度分类器，对睡眠数据进行多分类视图预测，然后计算多个分类视图的逻辑一致性系数即可信度R，作为睡眠分期可信度评价指标。

为了达到上述目的，本发明采用的技术方案为：

一种睡眠信号分期可信度自动评价方法，包括以下步骤：

步骤一：选择主分类器和辅分类器，采用监督学习对主、辅分类器进行训练，获得训练后的主、辅分类器。

步骤二：输入睡眠数据，利用步骤一中训练后的主、辅分类器，对睡眠数据进行多粒度分类视图预测，并输出主分类器的预测结果。

步骤三：根据步骤二的主、辅分类器预测结果，计算并输出逻辑一致性系数即可信度R。

所述的步骤一具体为：

首先根据需要的分类粒度，选择该粒度的分类器作为主分类器，其他粒度的分类器为辅分类器；

用C_m表示主分类器；C₁、...、C_p表示p个辅分类器，然后用标记数据库对主、辅分类器分别进行训练，得到训练后的主、辅分类器。

所述的步骤二具体为：

使用步骤一中训练好的主、辅分类器，对需要分析的睡眠数据分别进行预测，得到不同分类粒度的预测结果，并输出主分类器的预测结果，具体过程：

睡眠数据用X(i)表示，其中i＝1，2，3，......，M，M代表数据长度；对X(i)按固定时间间隔划分为N个睡眠时段，得到Y(j)，其中j＝1，2，3，......，N；把Y(j)分别输入到主、辅分类器C_m和C₁、...、C_p，分别得到所有睡眠时段的预测结果Z_mj、Z_1j、…、Z_pj，并输出主分类器Z_mj的预测结果，其中j＝1，2，3，......，N。

所述的步骤三具体为：

利用步骤二得到的不同分类粒度预测结果Z_mj、Z_1j、...、Z_pj，根据逻辑一致性关系图，计算所有睡眠时段的逻辑一致性值；

所述的逻辑一致性关系图为：把多个睡眠视图，按粒度从细到粗排列，如果细视图的分类结果等效为相邻粗分类视图的结果，与该相邻粗分类视图的时间分类结果相等，就认为这两个视图在逻辑上一致。

第j睡眠时段的逻辑一致性值定义为：

根据L_j睡眠分期值，用公式(2)计算并输出逻辑一致性系数：

式中，j＝1，2，3，......，N，R的取值范围是0～1，R值越大，代表可信度越高。

本发明的优点是：

(1)、多视图是指对于同一事物，可以从多种不同途径或不同角度对其进行描述，以便更全面和深入地揭示问题本质。本发明利用多视图概念，提出一种基于多分类视图包含一致性的睡眠信号分期可信度自动评价方法，适用于对基于各种生理信号和分类器的睡眠分期结果进行可信度评价。同时可利用输出的可信度值，对海量无标记睡眠数据自动标记来更新训练集规模；以及利用海量无标记数据直接训练分类器训练过程，提高分类器精度。主要可用于以下三种应用场景：1)对睡眠分期结果给出可信度指标；2)对无标记睡眠数据进行自动标记，迭代更新标记数据库规模；3)直接利用海量无标记数据，训练健壮性更好、结果可信度更高的睡眠分期算法。

(2)、本发明由睡眠数据X(i)得到的睡眠分段数据Y(j)，可以是对X(i)直接进行分段获得，也可以是对X(i)分段后提取的特征向量。

本发明中的辅分类器可以根据具体应用灵活选择。例如应用于家庭睡眠评价时，选择四分类器为主分类器，即可以选择二、三分类器为辅分类器，也可选择二、三、五分类器为辅分类器。

附图说明

图1为本发明睡眠信号分期可信度自动评价方法流程图。

图2为本发明实施例的分类器结构示意图。

图3为本发明实施例的逻辑一致性系数R计算过程。

图4为专业医师睡眠分期结果图和本发明实施例的自动分期结果图。

图5为逻辑一致性关系图。

具体实施方式

参照图1，本发明结合附图提供一种睡眠信号分期可信度自动评价方法的具体实施例，即将上述可信度自动评价方法，用于对一种利用单导ECG信号、通过双向长短时记忆(bi-directional long short-term memory，Bi-LSTM)深度神经网络构建的睡眠分期模型得到的睡眠分期结果进行评价，包括以下步骤：

一种睡眠信号分期可信度自动评价方法，包括以下步骤：

步骤一：选择主分类器和辅分类器，采用监督学习对主、辅分类器进行训练，获得训练后的主、辅分类器。

所述的步骤一具体为：

结合具体的应用场景，首先根据需要的分类粒度，选择该粒度的分类器作为主分类器，其他粒度的分类器为辅分类器。

例如，需要获得与临床分期类型数一致的结果，则选择五分类器作为主分类器，二、三、四分类器作为辅分类器；若应用于家庭睡眠评价，则选择四分类器为主分类器，二、三分类器为辅分类器，也可选择二、三、五分类器为辅分类器。

用C_m表示主分类器；C₁、...、C_p表示p个辅分类器，然后用标记数据库对主、辅分类器分别进行训练，得到训练后的主、辅分类器。

结合临床睡眠分期的需求，参照图2，首先基于Bi-LSTM网络分别构建二、三、四和五分类器器C₁、C₂、C₃、C_m，其中主分类器C_m为五分类器，其他三个为辅分类器，然后对标记数据库中的ECG信号以30s为间隔划分为N个睡眠时段，每个睡眠时段提取与睡眠分期相关的25个特征，最后利用提取的特征向量和标签对主、辅分类器进行训练。

步骤二：输入睡眠数据，利用步骤一中训练后的主、辅分类器，对睡眠数据进行多粒度分类视图预测，并输出主分类器的预测结果。

使用步骤一中训练好的主、辅分类器，对需要分析的睡眠数据分别进行预测，得到不同分类粒度的预测结果，具体过程：

睡眠数据用X(i)表示，其中i＝1，2，3，......，M，M代表数据长度；对X(i)按固定时间间隔划分为N个睡眠时段，得到Y(j)，其中j＝1，2，3，......，N；把Y(j)分别输入到主、辅分类器C_m和C₁、...、C_p，分别得到所有睡眠时段的预测结果Z_mj、Z_1j、...、Z_pj，并且输出主分类器Z_mj的结果，其中j＝1，2，3，......，N。

具体为：对一段8小时13分钟30s的ECG信号数据X(i)按照30s间隔划分为987个睡眠时段，得到Y(j)，其中i＝1，2，3，......，29610s，M代表数据长度，j＝1，2，3，......，987。将Y(j)输入步骤一中训练好的主、辅分类器C_m和C₁、C₂、C₃，分别得到所有睡眠时段的预测结果Z_mj、Z_1j、Z_2j、Z_3j，，如图3中(a)、(b)、(c)、(d)所示，并且输出主分类器Z_mj的结果。如图4所示，为专业医师睡眠分期结果图和本发明自动分期结果图，说明本发明的方法基本准确，可信度越高，与专业医生做出的判断基本一致。

步骤三：根据步骤二的主、辅分类器预测结果，计算并输出逻辑一致性系数即可信度R。

利用步骤二得到的不同分类粒度预测结果Z_mj、Z_1j、...、Z_pj，根据逻辑一致性关系图，计算所有睡眠时段的逻辑一致性值；

所述的逻辑一致性关系图为：把多个睡眠视图，按粒度从细到粗排列，如果细视图的分类结果等效为相邻粗分类视图的结果，与该相邻粗分类视图的时间分类结果相等，就认为这两个视图在逻辑上一致。

如图5所示，这就是五分类视图，如果把N1、N2合并为浅睡(light sleep，LS)期，把N3称为慢波睡眠(slow wave sleep，SWS)期，其他分期状态不变就得到四分类视图；进一步粗粒化，把四分类的LS期和SWS期合并为NREM期，其他分期状态不变就得到三分类视图；更进一步粗粒化，把三分类的LS期和SWS期合并为睡眠(sleep，S)期，W期状态不变就得到二分类视图。

例如五分类视图的结果为N2，四分类视图的结果为LS，由于五分类视图的N2可以等效为四分类视图LS，与四分类视图的结果LS相等，所以认为这两个视图逻辑上一致；否则，例如五分类视图的结果为W，四分类视图的结果为REM，由于五分类视图的W等效为四分类视图W，与四分类视图的结果REM不等，所以认为这两个视图逻辑上不一致。按这种方式，验证所有相邻视图的一致性，如果所有结果都一致，就认为该时刻分期结果逻辑上一致；否则有任何一个相邻视图逻辑上不一致，就认为该时刻分期结果逻辑上不一致。

按公式(1)计算所有睡眠时段的逻辑一致性值L_j，结果如图3(e)所示。

第j睡眠时段的逻辑一致性值定义为：

最后按公式(2)，

根据L_j睡眠分期值，用公式(2)计算并输出逻辑一致性系数：

式中，j＝1，2，3，......，N，R的取值范围是0～1，R值越大，代表可信度越高。

计算得到R＝0.97，代表可信度很高。

Claims (4)

1.一种睡眠信号分期可信度自动评价方法，其特征在于，包括以下步骤：

步骤一：选择主分类器和辅分类器，采用监督学习对主、辅分类器进行训练，获得训练后的主、辅分类器；

步骤二：输入睡眠数据，利用步骤一中训练后的主、辅分类器，对睡眠数据进行多粒度分类视图预测，并输出主分类器的预测结果；

步骤三：根据步骤二的主、辅分类器预测结果，计算并输出逻辑一致性系数即可信度R。

2.根据权利要求1所述的一种睡眠信号分期可信度自动评价方法，其特征在于，所述的步骤一具体为：

首先根据需要的分类粒度，选择该粒度的分类器作为主分类器，其他粒度的分类器为辅分类器；

用C_m表示主分类器；C₁、...、C_p表示p个辅分类器，然后用标记数据库对主、辅分类器分别进行训练，得到训练后的主、辅分类器。

3.根据权利要求1所述的一种睡眠信号分期可信度自动评价方法，其特征在于，所述的步骤二具体为：

使用步骤一中训练好的主、辅分类器，对需要分析的睡眠数据分别进行预测，得到不同分类粒度的预测结果，并输出主分类器的预测结果，具体过程：

睡眠数据用X(i)表示,其中i＝1，2，3，......，M，M代表数据长度；对X(i)按固定时间间隔划分为N个睡眠时段，得到Y(j)，其中j＝1，2，3，......，N；把Y(j)分别输入到主、辅分类器C_m和C₁、...、C_p，分别得到所有睡眠时段的预测结果Z_mj、Z_1j、...、Z_pj，并输出主分类器Z_mj的预测结果，其中j＝1，2，3，......，N。

4.根据权利要求1所述的一种睡眠信号分期可信度自动评价方法，其特征在于，所述的步骤三具体为：

利用步骤二得到的不同分类粒度预测结果Z_mj、Z_1j、...、Z_pj，根据逻辑一致性关系图，计算所有睡眠时段的逻辑一致性值；

所述的逻辑一致性关系图为：把多个睡眠视图，按粒度从细到粗排列，如果细视图的分类结果等效为相邻粗分类视图的结果，与该相邻粗分类视图的时间分类结果相等，就认为这两个视图在逻辑上一致；

第j睡眠时段的逻辑一致性值定义为：

根据L_j睡眠分期值，用公式(2)计算并输出逻辑一致性系数：

式中，j＝1，2，3，......，N，R的取值范围是0～1，R值越大，代表可信度越高。

Images