CN110141226B - 自动睡眠分期方法、装置、计算机设备及计算机存储介质 - Google Patents

自动睡眠分期方法、装置、计算机设备及计算机存储介质 Download PDF

Info

Publication number
CN110141226B
CN110141226B CN201910455076.6A CN201910455076A CN110141226B CN 110141226 B CN110141226 B CN 110141226B CN 201910455076 A CN201910455076 A CN 201910455076A CN 110141226 B CN110141226 B CN 110141226B
Authority
CN
China
Prior art keywords
sleep
level representation
automatic
samples
staging
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201910455076.6A
Other languages
English (en)
Other versions
CN110141226A (zh
Inventor
王兴军
徐子尧
邹庆言
冯发润
黄海平
李庆
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Dongguan Jianda Information Technology Co ltd
Shenzhen International Graduate School of Tsinghua University
Original Assignee
Dongguan Jianda Information Technology Co ltd
Shenzhen Graduate School Tsinghua University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Dongguan Jianda Information Technology Co ltd, Shenzhen Graduate School Tsinghua University filed Critical Dongguan Jianda Information Technology Co ltd
Priority to CN201910455076.6A priority Critical patent/CN110141226B/zh
Publication of CN110141226A publication Critical patent/CN110141226A/zh
Application granted granted Critical
Publication of CN110141226B publication Critical patent/CN110141226B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/369Electroencephalography [EEG]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/389Electromyography [EMG]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/398Electrooculography [EOG], e.g. detecting nystagmus; Electroretinography [ERG]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/48Other medical applications
    • A61B5/4806Sleep evaluation
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/48Other medical applications
    • A61B5/4806Sleep evaluation
    • A61B5/4812Detecting sleep stages or cycles
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7264Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
    • A61B5/7267Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems involving training the classification device

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Surgery (AREA)
  • Public Health (AREA)
  • Pathology (AREA)
  • Veterinary Medicine (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Psychiatry (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physiology (AREA)
  • Signal Processing (AREA)
  • Psychology (AREA)
  • Ophthalmology & Optometry (AREA)
  • Evolutionary Computation (AREA)
  • Fuzzy Systems (AREA)
  • Mathematical Physics (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)

Abstract

本发明实施例提供了一种自动睡眠分期方法、装置、计算机设备及计算机可读存储介质,其中,该方法包括:获取待处理的睡眠数据,其中,所述睡眠数据包括脑电信号、眼电信号以及下颌肌电信号;将所述睡眠数据输入自动睡眠分期模型,输出所述睡眠数据的睡眠分期,其中,所述自动睡眠分期模型是以已知的睡眠数据和所述已知的睡眠数据的睡眠分期作为样本训练分类器得到的。该方案的自动睡眠分期过程由于采用了脑电信号、眼电信号以及下颌肌电信号,与现有技术中只根据单通道脑电进行自动睡眠分期的方式相比,可以有利于提高自动睡眠分期的准确率。

Description

自动睡眠分期方法、装置、计算机设备及计算机存储介质
技术领域
本发明涉及数据处理技术领域,特别涉及一种自动睡眠分期方法、装置、计算机设备及计算机可读存储介质。
背景技术
睡眠分期是按照一定的标准对人在睡眠过程中所经历的阶段进行划分。研究表明健康者在睡眠过程中会经历若干个睡眠阶段周期性循环的过程,而存在睡眠问题的患者这种周期性的睡眠结构不够明显或者比较混乱。例如,失眠症患者的睡眠障碍主要体现在入睡过程较为困难,即从清醒阶段到睡眠阶段的时间过长。相反,嗜睡症患者则极易入睡,从清醒阶段到睡眠阶段的时间极短。又如,睡眠呼吸障碍疾病是重大源头疾病,与多种疾病具有高共患率,我国未诊人群及高危人群数量极大。准确的睡眠分期监测能够为睡眠问题的诊断提供有效的辅助工具,进而采取相应的措施进行调理和治疗。
目前,睡眠监测的常见方式是先对患者实施整晚的PSG(多导睡眠图)睡眠参数监测,次日由技师对整晚数据曲线图进行人工分图和事件分析,再由医生根据事件进行分析诊断。除监控视频数据外,每个患者的整晚PSG监测一般在200MB以上,有经验的医师每天完成3例左右患者数据的完整分析和诊断报告。虽然基于规则的计算机辅助分图有了长足进展,但受准确性限制,手工分析仍是目前国内外普遍认可的方式。而受技师差异的影响,手工分图还存在高达5%~20%的变异率。随着睡眠监测的患者数量的增多,有限的睡眠分图技师越来越难以满足日益增长的睡眠分图需求。因此,自动睡眠分期越来越引起研究人员的兴趣。
睡眠分期是PSG分图的基础。对患者整晚的睡眠监测记录进行睡眠分期,不仅可以为分析患者的睡眠结构提供依据,还是进一步分析患者是否发生微觉醒事件的参考依据,而睡眠分期和微觉醒事件,又是进一步分析患者是否发生呼吸暂停、低通气等呼吸事件的参考依据。自动睡眠分期需要从每帧PSG信号中提取特征,再根据特征进行分类。传统地,提取的特征包括时频特征和非线性动力学特征,时频特征的提取方法包括小波变换、经验模式分解等,非线性动力学特征包括多尺度熵、Lyapunov指数等。根据手工设计的特征,使用传统的机器学习模型,例如线性判别分析、随机森林、支持向量机、AdaBoosting,进行自动睡眠分期。基于特征工程的自动睡眠分期模型依赖于人工设计的特征,然而,每种人工设计的特征都存在其局限性,我们很难知道应该提取哪些特征。
发明内容
本发明实施例提供了一种自动睡眠分期方法,以解决现有技术中由于只根据单通道脑电进行自动睡眠分期而存在的睡眠分期准确率低的技术问题。该方法包括:
获取待处理的睡眠数据,其中,所述睡眠数据包括脑电信号、眼电信号以及下颌肌电信号;
将所述睡眠数据输入自动睡眠分期模型,输出所述睡眠数据的睡眠分期,其中,所述自动睡眠分期模型是以已知的睡眠数据和所述已知的睡眠数据的睡眠分期作为样本训练分类器得到的。
本发明实施例还提供了一种自动睡眠分期装置,以解决现有技术中由于只根据单通道脑电进行自动睡眠分期而存在的睡眠分期准确率低的技术问题。该装置包括:
数据获取模块,用于获取待处理的睡眠数据,其中,所述睡眠数据包括脑电信号、眼电信号以及下颌肌电信号;
睡眠分期模块,用于将所述睡眠数据输入自动睡眠分期模型,输出所述睡眠数据的睡眠分期,其中,所述自动睡眠分期模型是以已知的睡眠数据和所述已知的睡眠数据的睡眠分期作为样本训练分类器得到的。
本发明实施例还提供了一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述任意的自动睡眠分期方法,以解决现有技术中由于只根据单通道脑电进行自动睡眠分期而存在的睡眠分期准确率低的技术问题。
本发明实施例还提供了一种计算机可读存储介质,所述计算机可读存储介质存储有执行上述任意的自动睡眠分期方法的计算机程序,以解决现有技术中由于只根据单通道脑电进行自动睡眠分期而存在的睡眠分期准确率低的技术问题。
在本发明实施例中,提出了采用训练的自动睡眠分期模型来实现基于综合多个通道的信息(例如,包括脑电信号、眼电信号以及下颌肌电信号等多个通道的信息)的自动睡眠分期,该自动睡眠分期过程由于采用了脑电信号、眼电信号以及下颌肌电信号,与现有技术中只根据单通道脑电进行自动睡眠分期的方式相比,可以有利于提高自动睡眠分期的准确率。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,并不构成对本发明的限定。在附图中:
图1是本发明实施例提供的一种自动睡眠分期方法的流程图;
图2是本发明实施例提供的一种具体实施上述自动睡眠分期方法的流程图;
图3是本发明实施例提供的一种计算机设备的结构框图;
图4是本发明实施例提供的一种自动睡眠分期装置的结构框图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,下面结合实施方式和附图,对本发明做进一步详细说明。在此,本发明的示意性实施方式及其说明用于解释本发明,但并不作为对本发明的限定。
在本发明实施例中,提供了一种自动睡眠分期方法,如图1所示,该方法包括:
步骤102:获取待处理的睡眠数据,其中,所述睡眠数据包括脑电信号、眼电信号以及下颌肌电信号;
步骤104:将所述睡眠数据输入自动睡眠分期模型,输出所述睡眠数据的睡眠分期,其中,所述自动睡眠分期模型是以已知的睡眠数据和所述已知的睡眠数据的睡眠分期作为样本训练分类器得到的。
由图1所示的流程可知,在本发明实施例中,提出了采用训练的自动睡眠分期模型来实现基于综合多个通道的信息(例如,包括脑电信号、眼电信号以及下颌肌电信号等多个通道的信息)的自动睡眠分期,该自动睡眠分期过程由于采用了脑电信号、眼电信号以及下颌肌电信号,与现有技术中只根据单通道脑电进行自动睡眠分期的方式相比,可以有利于自动提高睡眠分期的准确率。
具体实施时,上述睡眠数据包括的脑电信号、眼电信号以及下颌肌电信号可以从PSG(多导睡眠图)信号中提取得到。
具体实施时,采用已知的睡眠数据和已知的睡眠数据的睡眠分期作为样本训练分类器,其中,已知的睡眠数据也包括脑电信号、眼电信号以及下颌肌电信号,已知的睡眠数据也可以从已知的PSG信号中提取,已知的睡眠数据的睡眠分期可以是由医师或其他分析方式确定的,本申请对此不做限定。
具体实施时,上述分类器可以包括但不限于朴素贝叶斯、决策树、逻辑回归、支持向量机、神经网络(例如,任意形式的卷积神经网络和循环神经网络)等。
具体实施时,为了提高自动睡眠分期模型的准确率,在本实施例中,通过以下方式练分类器得到所述自动睡眠分期模型:
采用所有样本训练分类器,得到初步的自动睡眠分期模型;具体的,在训练分类器的过程中,可以将所有样本划分为训练集、验证集以及测试集,通过训练一个神经网络,得到初步的自动睡眠分期模型。
针对每个样本,将该样本的睡眠数据输入所述初步的自动睡眠分期模型,将所述初步的自动睡眠分期模型中间层的输出作为该样本的高层表示样本,进而得到每个样本的高层表示样本;具体的,所述初步的自动睡眠分期模型中间层的输出包括所述初步的自动睡眠分期模型的一个或多个隐藏层的输出,即可以将初步的自动睡眠分期模型的任意一个隐藏层的输出或任意多个隐藏层组合的输出作为高层表示样本。
针对每类睡眠分期,对属于该类睡眠分期的样本的高层表示样本进行密度聚类,得到属于该类睡眠分期的多个高层表示子集;得到属于每类睡眠分期的高层表示子集;具体的,为了克服综合多个通道的信息存在特征波形不够明显的问题,通过对高层表示样本进行密度聚类,使得属于每类睡眠分期的多个高层表示子集分别具有不同的特征波形明显度。
采用属于每类睡眠分期的高层表示子集训练分类器,得到最终的所述自动睡眠分期模型。具体的,在采用属于每类睡眠分期的高层表示子集训练分类器的过程中,也可以将高层表示子集划分为训练集、验证集、测试集,进而训练并微调分类器得到最终的上述自动睡眠分期模型。
具体实施时,在本实施例中,可以通过以下方式对属于该类睡眠分期的样本的高层表示样本进行密度聚类,得到属于该类睡眠分期的多个高层表示子集:
针对属于该类睡眠分期的各样本的高层表示样本,对维度大于预设维度阈值的高层表示样本进行降维,降维至维度等于所述预设维度阈值;具体的,该维度可以是高层表示样本的数据点的个数,降维方法可以包括但不限于主成分分析、线性判别分析、等距映射、局部线性嵌入、拉普拉斯特征映射、局部保留投影等方法。
针对属于该类睡眠分期的所有高层表示样本,计算每两个高层表示样本之间的距离;
针对每个高层表示样本,计算与该高层表示样本的距离小于预设距离阈值的高层表示样本的个数,将该个数作为该高层表示样本的局部密度,将局部密度最大的高层表示样本作为该类睡眠分期的聚类中心;
按照与聚类中心的距离大小顺序,将属于该类睡眠分期的所有高层表示样本依次分成多个高层表示子集。
具体的,在将属于该类睡眠分期的所有高层表示样本依次分成多个高层表示子集的过程中,针对每类睡眠分期,按照与聚类中心的距离由小到大的顺序形成多个高层表示子集,包含高层表示样本的与聚类中心的距离越大的高层表示子集,特征波形越明显。
例如,将属于每个睡眠分期类别的高层表示样本集合记作
Di={(x,y)|y=i}
其中,x表示高层表示样本,y表示对应的睡眠分期类别,i=1,2,...,num_sleep_stages,num_sleep_stages表示睡眠分期的类别数。
对Di中的每个样本
Figure BDA0002076348650000051
其中xi(i=1,2,...,num_sleep_stages)表示属于第i个睡眠分期类别的高层表示样本,j表示样本序号,j=1,2,...,|Di|,|Di|表示Di集合的元素个数。
对Di集合中任意两个样本dij与dik(dik中的k也表示样本序号,只是k的数值不同于j),计算样本之间的距离disijk
Figure BDA0002076348650000052
对Di集合中每个样本dij,定义局部密度ρij
Figure BDA0002076348650000053
其中,I(·)为指示函数,dT为所有disijk中按从小到大排列的位于T%分位点的数值,T∈[0,100]为预先设定的参数。
Figure BDA0002076348650000061
为第i类睡眠分期的聚类中心,其中
Figure BDA0002076348650000062
对disick按从小到大的顺序划分为N个集合,用Din(n=1,2,...,N)表示Di的第n个高层表示子集,用Dim(m=1,2,...,N)表示Di的第m个高层表示子集,Din,Dim满足
Figure BDA0002076348650000063
Din∩Dim=Φ(1≤m,n≤N,m≠n)
Figure BDA0002076348650000064
表示包含有各个睡眠分期类别的第n个高层表示子集。
具体实施时,本申请发明人发现,如果不进行数据增强,在类别不平衡的数据集上训练,很容易使分类结果更偏向于样本数目更多的类别,造成分类精度下降,为了进一步提高睡眠分期模型的准确率,在本实施例中,对各睡眠分期的高层表示子集进行数据增强,例如,在采用属于每类睡眠分期的高层表示子集训练分类器之前,针对各类睡眠分期的相同序数的高层表示子集,在各类睡眠分期的相同序数的高层表示子集包括的高层表示样本数目不平衡的情况下,进行平衡处理,直至各类睡眠分期的相同序数的高层表示子集包括的高层表示样本数目相同。
具体的,针对各类睡眠分期的第n个高层表示子集En,针对每类睡眠分期的第n个高层表示子集Din,统计该第n个高层表示子集En包括多少个该类睡眠分期的高层表示样本,如各类睡眠分期的第n个高层表示子集En所包含的高层表示样本数目|Din|不同,则进行平衡处理,直至各类睡眠分期的第n个高层表示子集所包含的高层表示样本数目相同,记对Din进行平衡处理之后得到的集合为D'in,即使|D'1n|=|D'2n|=…=|D'Sn|(S=num_sleep_stages,n=1,2,...,N),得到N个各类别样本数目平衡后的高层表示子集
Figure BDA0002076348650000065
平衡方法可以包括但不限于简单上采样、简单下采样、SMOTE算法、训练生成式对抗网络生成样本等方法。
具体实施时,由于采集到的睡眠监测信号有时会受到较强噪声的干扰,有时会出现特征波不明显的情况,这些情况可能导致数据质量的不一致性,为了避免数据质量的不一致性对分类器训练产生影响,提升分类精度,在本实施例中,采用属于每类睡眠分期的高层表示子集训练并微调分类器,包括:
针对属于每类睡眠分期的高层表示子集,按照与聚类中心的距离从小到大的顺序,采用属于该类睡眠分期的高层表示子集训练分类器,即先采用与聚类中心的距离相对小的高层表示子集训练分类器,再采用与聚类中心的距离相对大的高层表示子集训练分类器,以得到最终的自动睡眠分期模型。
具体实施时,为了进一步避免数据受噪声影响,在本实施例中,在以已知的睡眠数据和已知的睡眠数据的睡眠分期作为样本训练分类器之前,将样本中已知的睡眠数据变换为同一固定频率;即将脑电信号EEG_raw、眼电信号EOG_raw以及下颌肌电信号EMG_raw的采样频率变换到某一固定频率,例如变换采样频率到100Hz,得到EEG_100hz,EOG_100hz,EMG_100hz。还可以对变换频率后的已知的睡眠数据进行归一化处理,对采样频率变换后的脑电信号EEG_100hz、眼电信号EOG_100hz、下颌肌电信号EMG_100hz进行数据归一化(数据归一化的方法可以包括但不限于线性函数归一化、零均值归一化等),分别得到EEG_norm,EOG_norm,EMG_norm,将每一帧数据归一化后的脑电信号EEG_norm、眼电信号EOG_norm、下颌肌电信号EMG_norm与前后若干帧信号拼接起来,作为一个样本,从而得到一个脑电信号、眼电信号、下颌肌电信号联合的数据集ExG,进而采用数据集ExG来训练分类器。
具体的,以下详细描述实施上述自动睡眠分期方法的过程,如图2所示,该过程包括以下步骤:
步骤S1、从睡眠监测得到的PSG信号中,提取脑电信号、眼电信号、下颌肌电信号及医师标注的睡眠分期,以此作为样本,并对样本进行预处理,得到预处理后的脑电信号、眼电信号、下颌肌电信号;
具体的,步骤S1可以具体包括S11、S12、S13以及S14:
步骤S11、从睡眠监测记录文件中,根据文件格式协议,将脑电信号、眼电信号、下颌肌电信号提取出来,分别得到EEG_raw,EOG_raw,EMG_raw;
步骤S12、将提取出的脑电信号EEG_raw、眼电信号EOG_raw、下颌肌电信号EMG_raw的采样频率变换到某一固定频率,例如变换采样频率到100Hz,分别得到EEG_100hz,EOG_100hz,EMG_100hz;
步骤S13、对采样频率变换后的脑电信号EEG_100hz、眼电信号EOG_100hz、下颌肌电信号EMG_100hz进行数据归一化,得到EEG_norm,EOG_norm,EMG_norm;
步骤S14、将每一帧数据归一化后的脑电信号EEG_norm、眼电信号EOG_norm、下颌肌电信号EMG_norm与前后若干帧信号拼接起来,作为一个样本,从而得到一个脑电信号、眼电信号、下颌肌电信号联合的样本数据集ExG。
步骤S2、将步骤S1得到的预处理后的脑电信号、眼电信号、下颌肌电信号组成的样本数据集划分为训练集、验证集、测试集,并训练一个神经网络,得到初步的自动睡眠分期模型;
具体的,步骤S2具体可以包括S21、S22:
步骤S21、将得到的脑电信号、眼电信号、下颌肌电信号联合的样本数据集,划分为训练集ExG_train、验证集ExG_validation、测试集ExG_test;
步骤S22、根据划分的训练集、验证集、测试集,训练分类网络,得到初步的自动睡眠分期模型model01。所述的分类网络,包括任意形式的卷积神经网络和循环神经网络。
步骤S3、对所有样本ExGi(i=1,2,…,n_samples),n_samples为样本数,以样本作为步骤S2得到的初步的自动睡眠分期模型model01的输入,将初步的自动睡眠分期模型model01中间层的输出fi=feature(ExGi),作为样本的高层表示样本,进而得到所有样本的高层表示样本;
具体的,步骤S3具体可以包括S31、S32:
步骤S31、将脑电信号、眼电信号、下颌肌电信号联合的样本数据集ExG,输入初步的自动睡眠分期模型model01,得到模型的中间层输出fi=feature(ExGi),作为数据集的高层表示。所述的中间层输出,可以是分类模型任意一个隐藏层的输出或任意几个隐藏层输出的组合;
步骤S32、将样本的高层表示fi=feature(ExGi)和睡眠分期标签yi组合在一起,作为高层表示样本数据集D={(fi,yi),i=1,2,...,n_samples}。
步骤S4、将属于每个睡眠分期类别的高层表示样本集合记作
Di={(x,y)|y=i}
其中,x表示高层表示样本,y表示对应的睡眠分期类别,i=1,2,...,num_sleep_stages,num_sleep_stages表示睡眠分期的类别数。
对Di中的每个样本
Figure BDA0002076348650000081
其中xi(i=1,2,...,num_sleep_stages)表示属于第i个睡眠分期类别的高层表示样本,j表示样本序号,j=1,2,...,|Di|,|Di|表示Di集合的元素个数。
对Di集合中任意两个样本dij与dik(dik中的k也表示样本序号,只是k的数值不同于j),计算样本之间的距离disijk
Figure BDA0002076348650000091
(y(j)=y(k)=i,1≤j,k≤|Di|)
对Di集合中每个样本dij,定义局部密度ρij
Figure BDA0002076348650000092
其中,I(·)为指示函数,dT为所有disijk中按从小到大排列的位于T%分位点的数值,T∈[0,100]为预先设定的参数。
Figure BDA0002076348650000093
为第i类睡眠分期的聚类中心,其中
Figure BDA0002076348650000094
对disick按从小到大的顺序划分为N个集合,用Din(n=1,2,...,N)表示Di的第n个高层表示子集,用
Figure BDA0002076348650000095
表示包含有各个睡眠分期类别的第n个高层表示子集。
步骤S5、对步骤S4得到的各类睡眠分期的第j个高层表示子集Ej,统计其中属于该类睡眠分期的样本数量,进行各类别样本数目的平衡,直至各类睡眠分期的第j个高层表示子集Ej所包含的高层表示样本数目相同,即使|D'1j|=|D'2j|=…=|D'num_sleep_stages,j|(j=1,2,…,N),得到N个睡眠分期类别的样本数目平衡后的高层表示子集
Figure BDA0002076348650000096
具体的,步骤S5具体可以包括S51、S52:
步骤S51、对每类睡眠分期的第j个高层表示子集,统计其中该睡眠分期类别的样本数目;
步骤S52、若各个睡眠分期类别的样本数目不平衡,则进行类别平衡,得到N个类别平衡后的高层表示子集。所述的类别平衡可以包括但不限于简单上采样、简单下采样、SMOTE算法、训练生成式对抗网络生成样本等方法。
步骤S6、对步骤S5得到的N个各类别样本数目平衡后的高层表示子集,对每个子集划分训练集、验证集、测试集,并设计学习策略,根据设计的学习策略分阶段地训练分类器,得到最终的自动睡眠分期模型。
具体的,步骤S6具体可以包括S61、S62:
步骤S61、对已得到的N个类别平衡后的高层表示子集中的每个子集
Figure BDA0002076348650000101
划分训练集
Figure BDA0002076348650000102
验证集
Figure BDA0002076348650000103
测试集
Figure BDA0002076348650000104
步骤S62、设计学习策略,根据设计的学习策略分阶段地训练分类器,先在局部密度大的高层表示子集上训练、再在局部密度小的高层表示子集上训练,得到最终的自动睡眠分期模型。
具体实施时,在本实施例中,还提供了计算机设备,如图3所示,包括存储器302、处理器304及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述任意的自动睡眠分期方法。
具体的,该计算机设备可以是计算机终端、服务器或者类似的运算装置。
具体实施时,在本实施例中,还提供了计算机可读存储介质,所述计算机可读存储介质存储有执行上述任意的自动睡眠分期方法的计算机程序。
具体的,计算机可读存储介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机可读存储介质的例子包括,但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带,磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质,可用于存储可以被计算设备访问的信息。按照本文中的界定,计算机可读存储介质不包括暂存电脑可读媒体(transitory media),如调制的数据信号和载波。
基于同一发明构思,本发明实施例中还提供了一种自动睡眠分期装置,如下面的实施例所述。由于自动睡眠分期装置解决问题的原理与自动睡眠分期方法相似,因此自动睡眠分期装置的实施可以参见自动睡眠分期方法的实施,重复之处不再赘述。以下所使用的,术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现,但是硬件,或者软件和硬件的组合的实现也是可能并被构想的。
图4是本发明实施例的自动睡眠分期装置的一种结构框图,如图4所示,该装置包括:
数据获取模块402,用于获取待处理的睡眠数据,其中,所述睡眠数据包括脑电信号、眼电信号以及下颌肌电信号;
自动睡眠分期模块404,用于将所述睡眠数据输入自动睡眠分期模型,输出所述睡眠数据的睡眠分期,其中,所述自动睡眠分期模型是以已知的睡眠数据和所述已知的睡眠数据的睡眠分期作为样本训练分类器得到的。
在一个实施例中,所述自动睡眠分期模块,包括:
第一训练单元,用于采用所有样本训练分类器,得到初步的自动睡眠分期模型;
高层表示获取单元,用于针对每个样本,将该样本的睡眠数据输入所述初步的自动睡眠分期模型,将所述初步的自动睡眠分期模型中间层的输出作为该样本的高层表示样本,得到每个样本的高层表示样本;
密度聚类单元,用于针对每类睡眠分期,对属于该类睡眠分期的样本的高层表示样本进行密度聚类,得到属于该类睡眠分期的多个高层表示子集;得到属于每类睡眠分期的高层表示子集;
第二训练单元,用于采用属于每类睡眠分期的高层表示子集训练分类器,得到最终的所述自动睡眠分期模型。
在一个实施例中,所述密度聚类单元,包括:
降维子单元,用于针对属于该类睡眠分期的各样本的高层表示样本,对维度大于预设维度阈值的高层表示样本进行降维,降维至维度等于所述预设维度阈值;
距离计算子单元,用于针对属于该类睡眠分期的所有高层表示样本,计算每两个高层表示样本之间的距离;
局部密度确定子单元,用于针对每个高层表示样本,计算与该高层表示样本的距离小于预设距离阈值的高层表示样本的个数,将该个数作为该高层表示样本的局部密度,将局部密度最大的高层表示样本作为该类睡眠分期的聚类中心;
密度聚类子单元,用于按照与聚类中心的距离大小顺序,将属于该类睡眠分期的所有高层表示样本依次分成多个高层表示子集。
在一个实施例中,上述装置还包括:
平衡模块,用于在采用属于每类睡眠分期的高层表示子集训练分类器之前,针对各类睡眠分期的相同序数的高层表示子集,在各类睡眠分期的相同序数的高层表示子集包括的高层表示样本数目不平衡的情况下,进行平衡处理,直至各类睡眠分期的相同序数的高层表示子集包括的高层表示样本数目相同。
在一个实施例中,所述初步的自动睡眠分期模型中间层的输出包括所述初步的自动睡眠分期模型的一个或多个隐藏层的输出。
在一个实施例中,所述第二训练单元具体用于,针对属于每类睡眠分期的高层表示子集,按照与聚类中心的距离由小到大的顺序,采用属于各类睡眠分期的各高层表示子集依次训练并微调分类器。
在一个实施例中,上述装置还包括:
预处理模块,用于在以已知的睡眠数据和所述已知的睡眠数据的睡眠分期作为样本训练分类器之前,将所述样本中所述已知的睡眠数据变换为同一固定频率;对变换频率后的所述已知的睡眠数据进行归一化处理。
本发明实施例实现了如下技术效果:提出了采用训练的自动睡眠分期模型来实现基于综合多个通道的信息(例如,包括脑电信号、眼电信号以及下颌肌电信号等多个通道的信息)的自动睡眠分期,该自动睡眠分期过程由于采用了脑电信号、眼电信号以及下颌肌电信号,与现有技术中只根据单通道脑电进行自动睡眠分期的方式相比,可以有利于提高自动睡眠分期的准确率。
显然,本领域的技术人员应该明白,上述的本发明实施例的各模块或各步骤可以用通用的计算装置来实现,它们可以集中在单个的计算装置上,或者分布在多个计算装置所组成的网络上,可选地,它们可以用计算装置可执行的程序代码来实现,从而,可以将它们存储在存储装置中由计算装置来执行,并且在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤,或者将它们分别制作成各个集成电路模块,或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样,本发明实施例不限制于任何特定的硬件和软件结合。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (14)

1.一种自动睡眠分期方法,其特征在于,包括:
获取待处理的睡眠数据,其中,所述睡眠数据包括脑电信号、眼电信号以及下颌肌电信号;
将所述睡眠数据输入自动睡眠分期模型,输出所述睡眠数据的睡眠分期,其中,所述自动睡眠分期模型是以已知的睡眠数据和所述已知的睡眠数据的睡眠分期作为样本训练分类器得到的;
训练分类器得到所述自动睡眠分期模型,包括:
采用所有样本训练分类器,得到初步的自动睡眠分期模型;
针对每个样本,将该样本的睡眠数据输入所述初步的自动睡眠分期模型,将所述初步的自动睡眠分期模型中间层的输出作为该样本的高层表示样本,得到每个样本的高层表示样本;
针对每类睡眠分期,对属于该类睡眠分期的样本的高层表示样本进行密度聚类,得到属于该类睡眠分期的多个高层表示子集;得到属于每类睡眠分期的高层表示子集;
采用属于每类睡眠分期的高层表示子集训练分类器,得到最终的所述自动睡眠分期模型。
2.如权利要求1所述的自动睡眠分期方法,其特征在于,对属于该类睡眠分期的样本的高层表示样本进行密度聚类,得到属于该类睡眠分期的多个高层表示子集,包括:
针对属于该类睡眠分期的各样本的高层表示样本,对维度大于预设维度阈值的高层表示样本进行降维,降维至维度等于所述预设维度阈值;
针对属于该类睡眠分期的所有高层表示样本,计算每两个高层表示样本之间的距离;
针对每个高层表示样本,计算与该高层表示样本的距离小于预设距离阈值的高层表示样本的个数,将该个数作为该高层表示样本的局部密度,将局部密度最大的高层表示样本作为该类睡眠分期的聚类中心;
针对属于该类睡眠分期的所有高层表示样本,按照与聚类中心的距离大小顺序,将属于该类睡眠分期的所有高层表示样本依次分成多个高层表示子集。
3.如权利要求1所述的自动睡眠分期方法,其特征在于,还包括:
在采用属于每类睡眠分期的高层表示子集训练分类器之前,针对各类睡眠分期的相同序数的高层表示子集,在各类睡眠分期的相同序数的高层表示子集包括的高层表示样本数目不平衡的情况下,进行平衡处理,直至各类睡眠分期的相同序数的高层表示子集包括的高层表示样本数目相同。
4.如权利要求1所述的自动睡眠分期方法,其特征在于,所述初步的自动睡眠分期模型中间层的输出包括所述初步的自动睡眠分期模型的一个或多个隐藏层的输出。
5.如权利要求2所述的自动睡眠分期方法,其特征在于,采用属于每类睡眠分期的高层表示子集训练分类器,包括:
针对属于每类睡眠分期的高层表示子集,按照与聚类中心的距离由小到大的顺序,采用属于各类睡眠分期的各高层表示子集依次训练并微调分类器。
6.如权利要求1至5中任一项所述的自动睡眠分期方法,其特征在于,还包括:
在以已知的睡眠数据和所述已知的睡眠数据的睡眠分期作为样本训练分类器之前,将所述样本中所述已知的睡眠数据去除噪声、去除伪迹以及变换数据采样频率为同一固定频率;
对去除噪声、去除伪迹以及变换数据采样频率后的所述已知的睡眠数据进行归一化处理。
7.一种自动睡眠分期装置,其特征在于,包括:
数据获取模块,用于获取待处理的睡眠数据,其中,所述睡眠数据包括脑电信号、眼电信号以及下颌肌电信号;
自动睡眠分期模块,用于将所述睡眠数据输入自动睡眠分期模型,输出所述睡眠数据的睡眠分期,其中,所述自动睡眠分期模型是以已知的睡眠数据和所述已知的睡眠数据的睡眠分期作为样本训练分类器得到的;
所述自动睡眠分期模块,包括:
第一训练单元,用于采用所有样本训练分类器,得到初步的自动睡眠分期模型;
高层表示获取单元,用于针对每个样本,将该样本的睡眠数据输入所述初步的自动睡眠分期模型,将所述初步的自动睡眠分期模型中间层的输出作为该样本的高层表示样本,得到每个样本的高层表示样本;
密度聚类单元,用于针对每类睡眠分期,对属于该类睡眠分期的样本的高层表示样本进行密度聚类,得到属于该类睡眠分期的多个高层表示子集;得到属于每类睡眠分期的高层表示子集;
第二训练单元,用于采用属于每类睡眠分期的高层表示子集训练分类器,得到最终的所述自动睡眠分期模型。
8.如权利要求7所述的自动睡眠分期装置,其特征在于,所述密度聚类单元,包括:
降维子单元,用于针对属于该类睡眠分期的各样本的高层表示样本,对维度大于预设维度阈值的高层表示样本进行降维,降维至维度等于所述预设维度阈值;
距离计算子单元,用于针对属于该类睡眠分期的所有高层表示样本,计算每两个高层表示样本之间的距离;
局部密度确定子单元,用于针对每个高层表示样本,计算与该高层表示样本的距离小于预设距离阈值的高层表示样本的个数,将该个数作为该高层表示样本的局部密度,将局部密度最大的高层表示样本作为该类睡眠分期的聚类中心;
密度聚类子单元,用于按照与聚类中心的距离大小顺序,将属于该类睡眠分期的所有高层表示样本依次分成多个高层表示子集。
9.如权利要求7所述的自动睡眠分期装置,其特征在于,还包括:
平衡模块,用于在采用属于每类睡眠分期的高层表示子集训练分类器之前,针对各类睡眠分期的相同序数的高层表示子集,在各类睡眠分期的相同序数的高层表示子集包括的高层表示样本数目不平衡的情况下,进行平衡处理,直至各类睡眠分期的相同序数的高层表示子集包括的高层表示样本数目相同。
10.如权利要求7所述的自动睡眠分期装置,其特征在于,所述初步的自动睡眠分期模型中间层的输出包括所述初步的自动睡眠分期模型的一个或多个隐藏层的输出。
11.如权利要求8所述的自动睡眠分期装置,其特征在于,所述第二训练单元具体用于,针对属于每类睡眠分期的高层表示子集,按照与聚类中心的距离由小到大的顺序,采用属于各类睡眠分期的各高层表示子集依次训练并微调分类器。
12.如权利要求7至11中任一项所述的自动睡眠分期装置,其特征在于,还包括:
预处理模块,用于在以已知的睡眠数据和所述已知的睡眠数据的睡眠分期作为样本训练分类器之前,将所述样本中所述已知的睡眠数据将所述样本中所述已知的睡眠数据去除噪声、去除伪迹以及变换数据采样频率为同一固定频率;对去除噪声、去除伪迹以及变换数据采样频率后的所述已知的睡眠数据进行归一化处理。
13.一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至6任一项所述的自动睡眠分期方法。
14.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有执行权利要求1至6任一项所述的自动睡眠分期方法的计算机程序。
CN201910455076.6A 2019-05-29 2019-05-29 自动睡眠分期方法、装置、计算机设备及计算机存储介质 Active CN110141226B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910455076.6A CN110141226B (zh) 2019-05-29 2019-05-29 自动睡眠分期方法、装置、计算机设备及计算机存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910455076.6A CN110141226B (zh) 2019-05-29 2019-05-29 自动睡眠分期方法、装置、计算机设备及计算机存储介质

Publications (2)

Publication Number Publication Date
CN110141226A CN110141226A (zh) 2019-08-20
CN110141226B true CN110141226B (zh) 2022-03-15

Family

ID=67593652

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910455076.6A Active CN110141226B (zh) 2019-05-29 2019-05-29 自动睡眠分期方法、装置、计算机设备及计算机存储介质

Country Status (1)

Country Link
CN (1) CN110141226B (zh)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110772267B (zh) * 2019-11-07 2022-04-19 中国人民解放军63850部队 一种人体生理疲劳数据标记方法及疲劳识别模型
CN111214211A (zh) * 2020-01-16 2020-06-02 珠海格力电器股份有限公司 睡眠监测方法、装置和智能床
CN111358435B (zh) * 2020-03-13 2023-02-28 珠海向量科技有限公司 一种提高深度神经网络精度的数据增强方法
CN111657860B (zh) * 2020-05-14 2021-07-20 清华大学 一种基于对抗学习的睡眠阶段识别方法及系统
CN112057068A (zh) * 2020-08-27 2020-12-11 武汉工程大学 一种癫痫病理数据分类方法、装置及存储介质
CN112426131A (zh) * 2020-11-22 2021-03-02 华东理工大学 一种基于分步聚类模型的睡眠分期自动判读方法
CN112568873B (zh) * 2021-02-25 2021-06-08 南京畅享医疗科技有限公司 一种实时睡眠监测记录与分析方法
CN114041753B (zh) * 2021-11-16 2024-04-26 上海市第六人民医院 睡眠分期方法、装置、计算机设备和存储介质
CN115374882B (zh) * 2022-10-24 2023-01-31 中国医学科学院药用植物研究所 睡眠分期方法、装置、电子设备及存储介质
CN115590477B (zh) * 2022-11-16 2023-03-28 中国医学科学院药用植物研究所 基于自监督的睡眠分期方法、装置、电子设备及存储介质

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1440380A1 (en) * 2001-11-01 2004-07-28 Oxford Biosignals Limited Improvements in or relating to physiological monitoring
US9754471B2 (en) * 2012-11-02 2017-09-05 Koninklijke Philips N.V. Electronic switch for controlling a device in dependency on a sleep stage
CN107495962A (zh) * 2017-09-18 2017-12-22 北京大学 一种单导联脑电的睡眠自动分期方法
CN108451505A (zh) * 2018-04-19 2018-08-28 广西欣歌拉科技有限公司 轻量入耳式睡眠分期系统
CN108595585A (zh) * 2018-04-18 2018-09-28 平安科技(深圳)有限公司 样本数据分类方法、模型训练方法、电子设备及存储介质
CN109044280A (zh) * 2018-08-20 2018-12-21 深圳和而泰数据资源与云技术有限公司 一种睡眠分期方法及相关设备
CN109645983A (zh) * 2019-01-09 2019-04-19 南京航空航天大学 一种基于多模块神经网络的不平衡心拍分类方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1440380A1 (en) * 2001-11-01 2004-07-28 Oxford Biosignals Limited Improvements in or relating to physiological monitoring
US9754471B2 (en) * 2012-11-02 2017-09-05 Koninklijke Philips N.V. Electronic switch for controlling a device in dependency on a sleep stage
CN107495962A (zh) * 2017-09-18 2017-12-22 北京大学 一种单导联脑电的睡眠自动分期方法
CN108595585A (zh) * 2018-04-18 2018-09-28 平安科技(深圳)有限公司 样本数据分类方法、模型训练方法、电子设备及存储介质
CN108451505A (zh) * 2018-04-19 2018-08-28 广西欣歌拉科技有限公司 轻量入耳式睡眠分期系统
CN109044280A (zh) * 2018-08-20 2018-12-21 深圳和而泰数据资源与云技术有限公司 一种睡眠分期方法及相关设备
CN109645983A (zh) * 2019-01-09 2019-04-19 南京航空航天大学 一种基于多模块神经网络的不平衡心拍分类方法

Also Published As

Publication number Publication date
CN110141226A (zh) 2019-08-20

Similar Documents

Publication Publication Date Title
CN110141226B (zh) 自动睡眠分期方法、装置、计算机设备及计算机存储介质
Cui et al. Automatic Sleep Stage Classification Based on Convolutional Neural Network and Fine‐Grained Segments
Andreotti et al. Multichannel sleep stage classification and transfer learning using convolutional neural networks
Yuan et al. A novel channel-aware attention framework for multi-channel eeg seizure detection via multi-view deep learning
CN110801221B (zh) 基于无监督特征学习的睡眠呼吸暂停片段检测设备
CH716863B1 (de) Verfahren zum Erstellen eines Depressionserkennungs-Klassifikationsmodells basierend auf einer Kanalauswahl von Mehrkanal-Elektroenzephalographie.
DE202019005960U1 (de) Verwalten von Atembeschwerden auf Basis von Geräuschen des Atemsystems
Boubchir et al. A review of feature extraction for EEG epileptic seizure detection and classification
Aljalal et al. Detection of Parkinson’s disease from EEG signals using discrete wavelet transform, different entropy measures, and machine learning techniques
CN111493822A (zh) 一种基于睡眠脑电的快速眼动期睡眠行为障碍分类方法
Jacob et al. Diagnosis of encephalopathy based on energies of EEG subbands using discrete wavelet transform and support vector machine
CN107066514A (zh) 老年人的情绪识别方法和系统
Yue et al. Deep learning for diagnosis and classification of obstructive sleep apnea: A nasal airflow-based multi-resolution residual network
CN112842279B (zh) 一种基于多维度特征参数的睡眠质量评估方法及装置
Waser et al. Quantifying synchrony patterns in the EEG of Alzheimer’s patients with linear and non-linear connectivity markers
Rachman et al. Alcoholism classification based on eeg data using independent component analysis (ica), wavelet de-noising and probabilistic neural network (pnn)
Back et al. Intra-and inter-epoch temporal context network (IITNet) for automatic sleep stage scoring
Berezutskaya et al. Open multimodal iEEG-fMRI dataset from naturalistic stimulation with a short audiovisual film
CN112263218A (zh) 睡眠分期方法及装置
Gudmundsson et al. Test–retest reliability and feature selection in physiological time series classification
Narmada et al. A novel adaptive artifacts wavelet Denoising for EEG artifacts removal using deep learning with Meta-heuristic approach
Siuly et al. Exploring Rhythms and Channels-Based EEG Biomarkers for Early Detection of Alzheimer's Disease
Pandey et al. Nocturnal sleep sounds classification with artificial neural network for sleep monitoring
CN112842363A (zh) 一种睡眠脑电检测方法及系统
Lee et al. Improving sleep stage classification performance by single-channel eeg data augmentation via spectral band blending

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CP01 Change in the name or title of a patent holder
CP01 Change in the name or title of a patent holder

Address after: Second floor, building a, Tsinghua campus, Shenzhen University Town, Xili street, Nanshan District, Shenzhen City, Guangdong Province

Patentee after: Shenzhen International Graduate School of Tsinghua University

Patentee after: Dongguan JIANDA INFORMATION Technology Co.,Ltd.

Address before: Second floor, building a, Tsinghua campus, Shenzhen University Town, Xili street, Nanshan District, Shenzhen City, Guangdong Province

Patentee before: GRADUATE SCHOOL AT SHENZHEN, TSINGHUA University

Patentee before: DONGGUAN JIANDA INFORMATION TECHNOLOGY Co.,Ltd.

PP01 Preservation of patent right
PP01 Preservation of patent right

Effective date of registration: 20240328

Granted publication date: 20220315