【技术实现步骤摘要】
一种改善N1期类别混淆的多模态多尺度睡眠分期方法
[0001]本专利技术涉及信号处理和模式识别等
,涉及一种改善N1期类别混淆的多模态多尺度睡眠分期方法。
技术介绍
[0002]根据美国睡眠医学会(American Academy of Sleep Medicine,AASM)标准和睡眠手册(Rechtschaffen and Kales,R&K),将睡眠分期分为清醒(W)期、快速眼动(REM)期、非快速眼动1(N1)期、非快速眼动2(N2)期和非快速眼动3(N3)期。AASM标准将整晚的多导睡眠图(Polysomnogram,PSG)数据划分为30秒一段的不重叠片段,逐一对这些数据片段人工判别为不同的睡眠分期,这导致睡眠分期会耗费大量的时间和精力。为解决这一问题,众多自动化睡眠分期方法被研究出来进行睡眠分期。但是,现有自动化睡眠分期方法都存在N1期准确率低的问题。
[0003]对于睡眠数据,每个睡眠分期的持续时间是不相等的。具体来说,N2期占大多数,约占总睡眠时间的45%
‑
55%,而N1期仅占的2%
‑
5%。这个问题存在于所有可用的睡眠数据集中,如蒙特利尔睡眠研究(MontrealArchiveofSleepStudies,MASS)档案数据库和PhysioNet的Sleep
‑
EDF数据库。有一些研究采用过采样、类不平衡抽样器以及GAN等方法弥补少数类,但是这些方法只扩充了样本数量,未保证样本的真实性。
[0004]现有 ...
【技术保护点】
【技术特征摘要】
1.一种改善N1期类别混淆的多模态多尺度睡眠分期方法,其特征在于,包括如下步骤:步骤一,原始信号预处理操作首先,选用脑电EEG信号中的Fpz
‑
Cz通道,眼电EOG信号和肌电EMG信号作为实验数据,即通道数Q取3,将采样频率为T的睡眠数据以30秒划分片段,删除MOVEMENT和UNKNOWN标签以及这些标签对应的数据,合并N3期和N4期数据统称为N3期;在训练集中只保留睡眠时间之前和之后的30分钟的W期;经过预处理得到N
old
个样本,每个样本数据x
i
均为长度为30
×
T的数据,其中i∈[1,2,
…
,N
old
]表示当前样本的序号,假定得到的样本集为X
old
,标签集为Y
old
;步骤二,数据增强对步骤一得到的训练集样本,使用叠取策略,对相邻的两个N1期样本,取两个样本的中间部分组成新的样本;将新生成的样本放入原来的数据集X
old
中,得到新的样本集X,相应地更新标签集Y
old
为新标签集Y;步骤三,多模态多尺度特征提取本发明设计了多模态多尺度特征提取方法充分利用睡眠数据并提取更具代表性的特征,;对步骤二得到的新的数据集,进行特征提取;将步骤二得到的数据按照模态分别进行不同处理,并对EEG信号分频率进行多尺度特征提取;拼接各个模态提取到的特征,送入SE注意力模块,选出具有代表性的特征并去除冗余特征;最后使用双层Bi
‑
LSTM学习睡眠分期的过渡规则,得到最终特征;步骤四,对比及分类经过步骤三可以获得睡眠分期的代表性特征,借助对比学习扩大睡眠分期边界,提升同期睡眠数据特征之间的相似性,降低非同期睡眠数据特征之间的相似性,从而使得各个分期之间更加容易区分;同时,使用分类器得到预测结果;本发明按照国际标准AASM睡眠分期准则进行5类睡眠分期,分别为清醒期(W),快速眼动期(REM),非快速眼动1期(N1),非快速眼动2期(N2),非快速眼动3期(N3)。2.根据权利要求1所述的一种改善N1期类别混淆的多模态多尺度睡眠分期方法,其特征在于,步骤二所述的数据增强,具体包括:使用叠取策略,选取X
old
中连续的两个同类样本和将的后v秒数据和的前(30
‑
v)秒数据拼接在一起,组成一个完整的标签为的新睡眠数据在本发明中v取15,表示生成的第i个新睡眠数据,表示第i个新睡眠数据对应的标签;以此法选取所有相邻的N1期样本,生成新的N1期样本,并将新生成的M个样本放入原来的数据集X
old
中,得到新的样本集X,相应地更新其标签集Y以及样本总数N;基于叠取策略的生成算法公式如下:式如下:N=N
old
+M其中,表示为c类生成的新样本,本发明中c取1;O(
·
)表示上述过程的叠取函数;M为针对类别c生成的新样本数量;N为原始片段数N
old
经过数据增强生成M个新N1期样本后的
训练集X中样本数量,X的维度为(N,Q,(30
×
T))。3.根据权利要求1所述的一种改善N1期类别混淆的多模态多尺度睡眠分期方法,其特征在于,步骤三所述的多模态多尺度特征提取方法,具体包括:将经过数据增强后的新样本集X作为输入数据,送入特征提取模块进行多模态多尺度特征提取,该模块结构图如图4所示;特征提取过程分为4个步骤:(1)多模态特征提取;对步骤2得到的睡眠数据X分模态提取睡眠特征;人类专家对于不同模态的数据有着不同的分期依据,具体地,根据EMG信号的有无、EEG信号中的各种特征波、EOG信号变化快慢来判断睡眠分期;对于输入数据的不同模态数据送入对应模态的特征提取分支中去:将EMG数据直接送入卷积模块进行特征提取,将EEG数据送入不同的卷积核模块中进行多尺度特征提取,对EOG数据先进行短时傅里叶变换,再送入卷积模块进行特征提取;短时傅里叶采用窗口大小为100的汉宁窗进行变换,提取EOG信号的时频特征;EMG和EOG信号特征提取公式如下所示:f
EMG
=E
EMG
(X
EMG
,s
EMG
)其中,f
EMG
、f
EOG
分别为EMG信号和EOG信号提取出的特征;X
EEG
和X
EMG
分别是数据集X的EEG信号和EMG信号数据,是数据集X经过短时傅里叶变换后的EOG信号;s
EMG
、s
EOG
分别是提取EMG信号、EOG信号特征选取的卷积核,其中s
EMG
=10、s
EOG
=(8,8);E
u
(
·
)是使用卷积模块ConvBlock进行特征提取的函数,其中u∈[EMG,EOG,EEG1,EEG2,EEG3],ConvBlock包括一个卷积层,一个批处理归一化层,一个GELU激活函数层,一个最大池化层和为防止过拟合增加参数为0.5的Dropout层;E
EMG
为EMG卷积模块,其中包括三个一维ConvBlock,每个卷积模块中的卷积核大小分别为10、8、8,卷积步长均为3,最大池化层的核分别为6、6、4,池化步长为3、3、4;E
EOG
为EOG卷积模块,其中包括三个二维ConvBlock,每个卷积模块中的卷积核大小分别为8
×
8、7
×
7、7
×
7,卷积步长分别为3
×
3、1
×
1、1
×
1,最大池化层的核分别为6
×
6、4
×
4、2
×
2,池化步长为3
×
3、2
×
2、2
×
2;(2)多尺度特征提取;使用(1)中的EEG信号数据X
EEG
作为输入,根据特征波选取不同尺度的卷积核大小从而提取EEG信号的多尺度睡眠特征,输出EEG信号特征f
EEG
,其中v∈[1,2,3];E...
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。