CN109645989B - 一种麻醉深度估计系统 - Google Patents

一种麻醉深度估计系统 Download PDF

Info

Publication number
CN109645989B
CN109645989B CN201811502030.7A CN201811502030A CN109645989B CN 109645989 B CN109645989 B CN 109645989B CN 201811502030 A CN201811502030 A CN 201811502030A CN 109645989 B CN109645989 B CN 109645989B
Authority
CN
China
Prior art keywords
neural network
convolutional neural
model
electroencephalogram signal
spectrum characteristic
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
CN201811502030.7A
Other languages
English (en)
Other versions
CN109645989A (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.)
Newrise Suzhou Medical Technology Co Ltd
Original Assignee
Yanshan 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 Yanshan University filed Critical Yanshan University
Priority to CN201811502030.7A priority Critical patent/CN109645989B/zh
Publication of CN109645989A publication Critical patent/CN109645989A/zh
Application granted granted Critical
Publication of CN109645989B publication Critical patent/CN109645989B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316—Modalities, i.e. specific diagnostic methods
    • A61B5/369—Electroencephalography [EEG]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/48—Other medical applications
    • A61B5/4821—Determining level or depth of anaesthesia
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235—Details of waveform analysis
    • A61B5/7264—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
    • A61B5/7267—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems involving training the classification device
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/02—Neural networks
    • G06N3/04—Architecture, e.g. interconnection topology
    • G06N3/044—Recurrent networks, e.g. Hopfield networks
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/02—Neural networks
    • G06N3/04—Architecture, e.g. interconnection topology
    • G06N3/045—Combinations of networks
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/02—Neural networks
    • G06N3/08—Learning methods
    • G06N3/084—Backpropagation, e.g. using gradient descent

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Artificial Intelligence (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Theoretical Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Mathematical Physics (AREA)
  • Evolutionary Computation (AREA)
  • Software Systems (AREA)
  • Public Health (AREA)
  • General Physics & Mathematics (AREA)
  • Computing Systems (AREA)
  • Data Mining & Analysis (AREA)
  • Computational Linguistics (AREA)
  • Veterinary Medicine (AREA)
  • General Engineering & Computer Science (AREA)
  • Animal Behavior & Ethology (AREA)
  • Surgery (AREA)
  • Medical Informatics (AREA)
  • Pathology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Psychiatry (AREA)
  • Signal Processing (AREA)
  • Physiology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Fuzzy Systems (AREA)
  • Psychology (AREA)
  • Anesthesiology (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
  • Image Analysis (AREA)

Abstract

本发明公开了一种麻醉深度估计方法及系统。所述方法包括：获取脑电信号训练数据；提取脑电信号训练数据的频谱特征图；将频谱特征图输入至卷积神经网络模型中，采用遗传算法对卷积神经网络模型中的初始权重进行优化，得到卷积神经网络优化模型；将频谱特征图作为卷积神经网络优化模型的输入，采用反向传播算法对卷积神经网络优化模型中的权重进行调节，得到训练好的卷积神经网络优化模型；获取待检测者当前的脑电信号数据；提取当前的脑电信号数据的频谱特征图，得到待检测频谱特征图；将待检测频谱特征图输入至训练好的卷积神经网络优化模型中，得到待检测者的当前状态。本发明能够提高麻醉深度预测的准确性，且计算复杂度低。

Description

一种麻醉深度估计系统
技术领域
本发明涉及麻醉监测技术领域，特别是涉及一种麻醉深度估计方法及系统。
背景技术
麻醉是临床手术过程中必不可少的关键环节，麻醉不当可能会造成患者术中知晓、术后疼痛、术后呕心等不良后果，严重时可致死。因而对麻醉深度的精确估计对于提高麻醉质量，保障手术安全有着重要意义。
自Gibbs等人首次提出用脑电信号(Electroencephalogram，EEG)监测麻醉深度的可能性，已引起越来越多的人对应用EEG监测麻醉深度的研究。以往对EEG信号的研究大多需要从预处理后的信号中计算量化指标，由于人脑思维模式和大脑信号非常复杂，很多采集的信号没有找到明确的意义，手动提取特征会造成信息的损失。近年来，基于EEG的图像识别开始在麻醉学领域被关注。
随着深度学习的兴起，卷积神经网络(Convolutional Neural Network，CNN)受到研究者的广泛关注，并在图像、语音、视频等诸多领域都取得了一定的应用成果。对于生物电信号的分类，目前很少有人使用CNN去替代传统分类方法。CNN的权重共享网络结构使之更类似于生物神经网络，降低了网络模型的复杂度，减少了权重的数量，能够避免人工设计的特征抽取器，依靠反向传播学习算法将不同层变为合适的特征抽取器。CNN是比较适合于像EEG一样随着时间变化很大的信号，可用于分析识别EEG频谱信息。然而，传统的CNN对初始权重的依赖性很大，若初始值设置不当，CNN的训练过程容易陷入局部最优，预测准确性会受到很大影响。
因此，目前麻醉深度估计方法的预测准确度有待提高。
发明内容
基于此，有必要提供一种麻醉深度估计方法及系统，以提高麻醉深度预测的准确性。
为实现上述目的，本发明提供了如下方案：
一种麻醉深度估计方法，所述方法包括：
获取脑电信号训练数据；所述脑电信号训练数据为检测者在整个麻醉时期的脑电信号；
提取所述脑电信号训练数据的频谱特征图；所述频谱特征图包括清醒期频谱图、麻醉期频谱图和恢复期频谱图；
将所述频谱特征图输入至卷积神经网络模型中，采用遗传算法对所述卷积神经网络模型中的初始权重进行优化，得到卷积神经网络优化模型；
将所述频谱特征图作为所述卷积神经网络优化模型的输入，采用反向传播算法对所述卷积神经网络优化模型中的权重进行调节，得到训练好的卷积神经网络优化模型；
获取待检测者当前的脑电信号数据；
提取所述当前的脑电信号数据的频谱特征图，得到待检测频谱特征图；
将所述待检测频谱特征图输入至所述训练好的卷积神经网络优化模型中，得到待检测者的当前状态；所述待检测者的当前状态为清醒期、麻醉期或恢复期。
可选的，所述将所述频谱特征图输入至卷积神经网络模型中，采用遗传算法对所述卷积神经网络模型中的初始权重进行优化，得到卷积神经网络优化模型，具体包括：
构建卷积神经网络模型；所述神经网络模型包括卷积层、池化层和全连接层；
确定所述卷积神经网络模型中卷积层和全连接层的初始权重；
将所述初始权重进行二进制编码，得到初始化种群；所述初始化种群中每个个体对应一个卷积神经网络模型；
将所述频谱特征图作为所述初始化种群中各卷积神经网络模型的输入，计算所述初始化种群中每个卷积神经网络模型的适应度值；
判断是否存在适应度值小于预设分类误差值或当前迭代次数是否大于或等于预设迭代次数；
若是，则将最大的适应度值对应的卷积神经网络模型确定为卷积神经网络优化模型；
若否，则对所述初始化种群中的个体进行选择、交叉和变异操作，更新初始化种群以及当前迭代次数，并返回所述将所述频谱特征图作为所述初始化种群中各卷积神经网络模型的输入，计算所述初始化种群个体中每个卷积神经网络模型的适应度值。
可选的，所述将所述频谱特征图作为所述卷积神经网络优化模型的输入，采用反向传播算法对所述卷积神经网络优化模型中的权重进行调节，得到训练好的卷积神经网络优化模型，具体包括：
将所述频谱特征图作为所述卷积神经网络优化模型的输入，前向计算所述卷积神经网络优化模型中每个神经元的输出值；第j个神经元的输出值为aj；
反向计算所述卷积神经网络优化模型中每个神经元的误差项；第j个神经元的误差项为
Figure GDA0002789199250000031
其中，Ed为卷积神经网络优化模型的损失函数，netj为第j个神经元的加权输入；
依据各神经元的输出值以及各神经元的误差项计算各个神经元连接权重的梯度值；从神经元i连接到神经元j的连接权重wji的梯度值
Figure GDA0002789199250000032
其中，ai为第i个神经元的输出值；
依据所述梯度值，采用梯度下降法对所述卷积神经网络优化模型中的权重进行调节，得到训练好的卷积神经网络优化模型。
可选的，在所述获取脑电信号训练数据之后，还包括：
对所述脑电信号训练数据进行预处理。
可选的，所述对所述脑电信号训练数据进行预处理，具体包括：
采用第一高通滤波器去除所述脑电信号训练数据中的头动噪声，得到第一处理脑电信号；
采用自适应滤波器去除所述第一处理脑电信号中的工频噪声，得到第二处理脑电信号；
采用低通滤波器去除所述第二处理脑电信号中大于第一预设频率的噪声信号，得到第三处理脑电信号；
采用第二高通滤波器去除所述第三处理脑电信号中小于第二预设频率的噪声信号，得到第四处理脑电信号。
可选的，所述脑电信号训练数据以及所述待检测者当前的脑电信号数据均是采用头戴式脑电信号采集系统采集得到的，且采样频率为128Hz。
本发明还提供了一种麻醉深度估计系统，所述系统包括：
第一数据获取模块，用于获取脑电信号训练数据；所述脑电信号训练数据为检测者在整个麻醉时期的脑电信号；
第一提取模块，用于提取所述脑电信号训练数据的频谱特征图；所述频谱特征图包括清醒期频谱图、麻醉期频谱图和恢复期频谱图；
权重优化模块，用于将所述频谱特征图输入至卷积神经网络模型中，采用遗传算法对所述卷积神经网络模型中的初始权重进行优化，得到卷积神经网络优化模型；
权重调节模块，用于将所述频谱特征图作为所述卷积神经网络优化模型的输入，采用反向传播算法对所述卷积神经网络优化模型中的权重进行调节，得到训练好的卷积神经网络优化模型；
第二数据获取模块，用于获取待检测者当前的脑电信号数据；
第二提取模块，用于提取所述当前的脑电信号数据的频谱特征图，得到待检测频谱特征图；
预测模块，用于将所述待检测频谱特征图输入至所述训练好的卷积神经网络优化模型中，得到待检测者的当前状态；所述待检测者的当前状态为清醒期、麻醉期或恢复期。
可选的，所权重优化模块，具体包括：
模型构建单元，用于构建卷积神经网络模型；所述神经网络模型包括卷积层、池化层和全连接层；
初始权重确定单元，用于确定所述卷积神经网络模型中卷积层和全连接层的初始权重；
编码单元，用于将所述初始权重进行二进制编码，得到初始化种群；所述初始化种群中每个个体对应一个卷积神经网络模型；
第一计算单元，用于将所述频谱特征图作为所述初始化种群中各卷积神经网络模型的输入，计算所述初始化种群中每个卷积神经网络模型的适应度值；
判断单元，用于判断是否存在适应度值小于预设分类误差值或当前迭代次数是否大于或等于预设迭代次数；
迭代单元，用于若不存在适应度值小于预设分类误差值且当前迭代次数小于预设迭代次数，则对所述初始化种群中的个体进行选择、交叉和变异操作，更新初始化种群以及当前迭代次数，并返回所述将所述频谱特征图作为所述初始化种群中各卷积神经网络模型的输入，计算所述初始化种群个体中每个卷积神经网络模型的适应度值；
优化模型确定单元，用于若存在适应度值小于预设分类误差值或当前迭代次数大于或等于预设迭代次数，则将最大的适应度值对应的卷积神经网络模型确定为卷积神经网络优化模型。
可选的，所述权重调节模块，具体包括：
第二计算单元，用于将所述频谱特征图作为所述卷积神经网络优化模型的输入，前向计算所述卷积神经网络优化模型中每个神经元的输出值；第j个神经元的输出值为aj；
第三计算单元，用于反向计算所述卷积神经网络优化模型中每个神经元的误差项；第j个神经元的误差项为
Figure GDA0002789199250000051
其中，Ed为卷积神经网络优化模型的损失函数，netj为第j个神经元的加权输入；
第四计算单元，用于依据各神经元的输出值以及各神经元的误差项计算各个神经元连接权重的梯度值；从神经元i连接到神经元j的连接权重wji的梯度值
Figure GDA0002789199250000052
其中，ai为第i个神经元的输出值；
调节单元，用于依据所述梯度值，采用梯度下降法对所述卷积神经网络优化模型中的权重进行调节，得到训练好的卷积神经网络优化模型。
可选的，所述系统还包括：预处理模块，用于对所述脑电信号训练数据进行预处理；
所述预处理模块，具体包括：
第一处理单元，用于采用第一高通滤波器去除所述脑电信号训练数据中的头动噪声，得到第一处理脑电信号；
第二处理单元，用于采用自适应滤波器去除所述第一处理脑电信号中的工频噪声，得到第二处理脑电信号；
第三处理单元，用于采用低通滤波器去除所述第二处理脑电信号中大于第一预设频率的噪声信号，得到第三处理脑电信号；
第四处理单元，用于采用第二高通滤波器去除所述第三处理脑电信号中小于第二预设频率的噪声信号，得到第四处理脑电信号。
与现有技术相比，本发明的有益效果是：
本发明提出了一种麻醉深度估计方法及系统，所述方法通过遗传算法优化卷积神经网络中各层的权重，避免训练过程陷入局部最优，提高了预测精度；采用反向传播算法调整权重，保证了预测模型的可靠性，能够更加精确地识别待检测者的麻醉深度状态，为术中麻醉深度的监测提供了更加方便可靠的途径。并且本发明还能够降低传统的特征指标计算的复杂程度，多个角度无监督学习EEG的图像特征，避免了信息丢失。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例一种麻醉深度估计方法的流程图；
图2为本发明实施例脑电信号训练数据示意图；
图3为本发明实施例脑电信号训练数据的整体频谱特征图；
图4为本发明实施例脑电信号训练数据不同时期的频谱特征图；
图5为本发明实施例卷积神经网络模型的结构示意图；
图6为本发明实施例一种麻醉深度估计系统的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1为本发明实施例一种麻醉深度估计方法的流程图。
参见图1，实施例的麻醉深度估计方法，包括：
步骤S1：获取脑电信号训练数据。所述脑电信号训练数据为检测者在整个麻醉时期的脑电信号，如图2所示。
步骤S2：提取所述脑电信号训练数据的频谱特征图。
该步骤首先依据所述脑电信号训练数据得到整个麻醉时期对应的整体频谱特征图，如图3所示；然后，将整体频谱特征图划分为不同时期的频谱特征图，所述频谱特征图包括清醒期频谱图、麻醉期频谱图和恢复期频谱图，如图4所示，其中图4(a)部分为清醒期频谱图，图4(b)部分为麻醉期频谱图，图4(c)部分为恢复期频谱图。
步骤S3：将所述频谱特征图输入至卷积神经网络模型中，采用遗传算法对所述卷积神经网络模型中的初始权重进行优化，得到卷积神经网络优化模型。
所述步骤S3，具体包括：
301：构建卷积神经网络模型。参见图5，所述神经网络模型由输入、卷积层、池化层、全连接层和输出构成，所述卷积层和所述全连接层各设置三层，利用卷积核提取输入图片的每个特征平面的特征，然后在池化层进行降采样处理，重复叠加三层卷积层加池化层的结构进行学习，最后叠加三个全连接层，输出分类结果。所述卷积层的激活函数采用Relu函数，所述全连接层利用softma x函数输出分类结果。
其中，所述Relu函数为：f(x)＝max(0,x)，
所述softmax函数为：
Figure GDA0002789199250000081
x表示对应函数的输入，即上一层的输出与本层权重的积再加偏置项。
302：确定所述卷积神经网络模型中卷积层和全连接层的初始权重。该步骤随机选择均值为零的较小区间中的数作为初始权重，例如选取区间[-0.3,+0.3]中的数作为初始权重。
303：将所述初始权重进行二进制编码，得到初始化种群。所述初始化种群中每个个体对应一个卷积神经网络模型，选择初始化种群中个体数目为50。
304：将所述频谱特征图作为所述初始化种群中各卷积神经网络模型的输入，计算所述初始化种群中每个卷积神经网络模型的适应度值。本实施例中将分类正确率作为适应度值，所述分类正确率为正确分类训练样本数占训练集总样本数的百分比。
305：判断是否存在适应度值小于预设分类误差值或当前迭代次数是否大于或等于预设迭代次数；若是，则将最大的适应度值对应的卷积神经网络模型确定为卷积神经网络优化模型；若否，则对所述初始化种群中的个体进行选择、交叉和变异操作，更新初始化种群以及当前迭代次数，并返回所述将所述频谱特征图作为所述初始化种群中各卷积神经网络模型的输入，计算所述初始化种群个体中每个卷积神经网络模型的适应度值。本实施例中，分类误差值表示错误分类的样本占输入训练样本总数的百分比，所述预设迭代次数的取值为[100,500]。
步骤S4：将所述频谱特征图作为所述卷积神经网络优化模型的输入，采用反向传播算法对所述卷积神经网络优化模型中的权重进行调节，得到训练好的卷积神经网络优化模型。
所述步骤S4，具体包括：
401：将所述频谱特征图作为所述卷积神经网络优化模型的输入，前向计算所述卷积神经网络优化模型中每个神经元的输出值；第j个神经元的输出值为aj。
402：反向计算所述卷积神经网络优化模型中每个神经元的误差项；第j个神经元的误差项为
Figure GDA0002789199250000091
其中，Ed为卷积神经网络优化模型的损失函数，netj为第j个神经元的加权输入。
403：依据各神经元的输出值以及各神经元的误差项计算各个神经元连接权重的梯度值；从神经元i连接到神经元j的连接权重wji的梯度值
Figure GDA0002789199250000092
其中，ai为第i个神经元的输出值。
404：依据所述梯度值，采用梯度下降法对所述卷积神经网络优化模型中的权重进行调节，得到训练好的卷积神经网络优化模型。具体的：
Figure GDA0002789199250000093
其中，wji′表示调节后的连接权重，η为学习速率。
步骤S5：获取待检测者当前的脑电信号数据。
步骤S6：提取所述当前的脑电信号数据的频谱特征图，得到待检测频谱特征图。
步骤S7：将所述待检测频谱特征图输入至所述训练好的卷积神经网络优化模型中，得到待检测者的当前状态；所述待检测者的当前状态为清醒期、麻醉期或恢复期。
作为一种可选的实施方式，由于原始的麻醉EEG信号中存在肌电、眼电、头动噪声、基线漂移等噪声，因此，在所述获取脑电信号训练数据之后，还要对所述脑电信号训练数据进行预处理。预处理过程具体为：采用第一高通滤波器去除所述脑电信号训练数据中的头动噪声，得到第一处理脑电信号；采用自适应滤波器去除所述第一处理脑电信号中的工频噪声，得到第二处理脑电信号；采用低通滤波器去除所述第二处理脑电信号中大于第一预设频率的噪声信号，得到第三处理脑电信号；采用第二高通滤波器去除所述第三处理脑电信号中小于第二预设频率的噪声信号，得到第四处理脑电信号。
同样的，在所述当前的脑电信号数据之后，还包括：对所述当前的脑电信号数据进行预处理。具体预处理过程与上述相同，即依次经过第一高通滤波器、自适应滤波器、低通滤波器和第二高通滤波器去除噪声。
作为一种可选的实施方式，所述脑电信号训练数据以及所述待检测者当前的脑电信号数据均是采用头戴式脑电信号采集系统采集得到的，且采样频率为128Hz。该实施方式中，第一高通滤波器采用eegfilt函数去除所述脑电信号训练数据中的头动噪声，去掉幅值较大的噪声信号，采用自适应滤波器去除的是第一处理脑电信号中的50Hz工频噪声，第一预设频率为45Hz，第二预设频率为0.1Hz。
本实施例的麻醉深度估计方法，通过遗传算法优化卷积神经网络中各层的权重，利用麻醉不同时期频谱特征作为麻醉深度估计指标，能够实现基于图像的麻醉状态识别，摆脱卷积神经网络初始权重对分类准确率的影响，使麻醉深度的监测更加准确避；采用反向传播算法调整权重，保证了预测模型的可靠性；能够降低传统的特征指标计算的复杂程度，多个角度无监督学习EEG的图像特征，避免了信息丢失。
本发明还提供了一种麻醉深度估计系统，图6为本发明实施例一种麻醉深度估计系统的结构示意图。参见图6，所述系统包括：
第一数据获取模块601，用于获取脑电信号训练数据；所述脑电信号训练数据为检测者在整个麻醉时期的脑电信号。
第一提取模块602，用于提取所述脑电信号训练数据的频谱特征图；所述频谱特征图包括清醒期频谱图、麻醉期频谱图和恢复期频谱图。
权重优化模块603，用于将所述频谱特征图输入至卷积神经网络模型中，采用遗传算法对所述卷积神经网络模型中的初始权重进行优化，得到卷积神经网络优化模型。
所权重优化模块603，具体包括：
模型构建单元，用于构建卷积神经网络模型；所述神经网络模型包括卷积层、池化层和全连接层；
初始权重确定单元，用于确定所述卷积神经网络模型中卷积层和全连接层的初始权重；
编码单元，用于将所述初始权重进行二进制编码，得到初始化种群；所述初始化种群中每个个体对应一个卷积神经网络模型；
第一计算单元，用于将所述频谱特征图作为所述初始化种群中各卷积神经网络模型的输入，计算所述初始化种群中每个卷积神经网络模型的适应度值；
判断单元，用于判断是否存在适应度值小于预设分类误差值或当前迭代次数是否大于或等于预设迭代次数；
迭代单元，用于若不存在适应度值小于预设分类误差值且当前迭代次数小于预设迭代次数，则对所述初始化种群中的个体进行选择、交叉和变异操作，更新初始化种群以及当前迭代次数，并返回所述将所述频谱特征图作为所述初始化种群中各卷积神经网络模型的输入，计算所述初始化种群个体中每个卷积神经网络模型的适应度值；
优化模型确定单元，用于若存在适应度值小于预设分类误差值或当前迭代次数大于或等于预设迭代次数，则将最大的适应度值对应的卷积神经网络模型确定为卷积神经网络优化模型。
权重调节模块604，用于将所述频谱特征图作为所述卷积神经网络优化模型的输入，采用反向传播算法对所述卷积神经网络优化模型中的权重进行调节，得到训练好的卷积神经网络优化模型。
所述权重调节模块604，具体包括：
第二计算单元，用于将所述频谱特征图作为所述卷积神经网络优化模型的输入，前向计算所述卷积神经网络优化模型中每个神经元的输出值；第j个神经元的输出值为aj；
第三计算单元，用于反向计算所述卷积神经网络优化模型中每个神经元的误差项；第j个神经元的误差项为
Figure GDA0002789199250000111
其中，Ed为卷积神经网络优化模型的损失函数，netj为第j个神经元的加权输入；
第四计算单元，用于依据各神经元的输出值以及各神经元的误差项计算各个神经元连接权重的梯度值；从神经元i连接到神经元j的连接权重wji的梯度值
Figure GDA0002789199250000121
其中，ai为第i个神经元的输出值；
调节单元，用于依据所述梯度值，采用梯度下降法对所述卷积神经网络优化模型中的权重进行调节，得到训练好的卷积神经网络优化模型。
第二数据获取模块605，用于获取待检测者当前的脑电信号数据。
第二提取模块606，用于提取所述当前的脑电信号数据的频谱特征图，得到待检测频谱特征图。
预测模块607，用于将所述待检测频谱特征图输入至所述训练好的卷积神经网络优化模型中，得到待检测者的当前状态；所述待检测者的当前状态为清醒期、麻醉期或恢复期。
作为一种可选的实时方式，所述系统还包括：预处理模块，用于对所述脑电信号训练数据进行预处理。所述预处理模块，具体包括：
第一处理单元，用于采用第一高通滤波器去除所述脑电信号训练数据中的头动噪声，得到第一处理脑电信号；
第二处理单元，用于采用自适应滤波器去除所述第一处理脑电信号中的工频噪声，得到第二处理脑电信号；
第三处理单元，用于采用低通滤波器去除所述第二处理脑电信号中大于第一预设频率的噪声信号，得到第三处理脑电信号；
第四处理单元，用于采用第二高通滤波器去除所述第三处理脑电信号中小于第二预设频率的噪声信号，得到第四处理脑电信号。
本实施例的麻醉深度估计系统，通过遗传算法优化卷积神经网络中各层的权重，避免训练过程陷入局部最优，提高了预测精度；采用反向传播算法调整权重，保证了预测模型的可靠性；能够降低传统的特征指标计算的复杂程度，多个角度无监督学习EEG的图像特征，避免了信息丢失。
对于实施例公开的系统而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。

Claims (3)

1.一种麻醉深度估计系统，其特征在于，所述系统包括：
第一数据获取模块，用于获取脑电信号训练数据；所述脑电信号训练数据为检测者在整个麻醉时期的脑电信号；
第一提取模块，用于提取所述脑电信号训练数据的频谱特征图；所述频谱特征图包括清醒期频谱图、麻醉期频谱图和恢复期频谱图；
权重优化模块，用于将所述频谱特征图输入至卷积神经网络模型中，采用遗传算法对所述卷积神经网络模型中的初始权重进行优化，得到卷积神经网络优化模型；
权重调节模块，用于将所述频谱特征图作为所述卷积神经网络优化模型的输入，采用反向传播算法对所述卷积神经网络优化模型中的权重进行调节，得到训练好的卷积神经网络优化模型；
第二数据获取模块，用于获取待检测者当前的脑电信号数据；
第二提取模块，用于提取所述当前的脑电信号数据的频谱特征图，得到待检测频谱特征图；
预测模块，用于将所述待检测频谱特征图输入至所述训练好的卷积神经网络优化模型中，得到待检测者的当前状态；所述待检测者的当前状态为清醒期、麻醉期或恢复期；
所述权重优化模块，具体包括：
模型构建单元，用于构建卷积神经网络模型；所述神经网络模型包括卷积层、池化层和全连接层；
初始权重确定单元，用于确定所述卷积神经网络模型中卷积层和全连接层的初始权重；
编码单元，用于将所述初始权重进行二进制编码，得到初始化种群；所述初始化种群中每个个体对应一个卷积神经网络模型；
第一计算单元，用于将所述频谱特征图作为所述初始化种群中各卷积神经网络模型的输入，计算所述初始化种群中每个卷积神经网络模型的适应度值；
判断单元，用于判断是否存在适应度值小于预设分类误差值或当前迭代次数是否大于或等于预设迭代次数；
迭代单元，用于若不存在适应度值小于预设分类误差值且当前迭代次数小于预设迭代次数，则对所述初始化种群中的个体进行选择、交叉和变异操作，更新初始化种群以及当前迭代次数，并返回所述将所述频谱特征图作为所述初始化种群中各卷积神经网络模型的输入，计算所述初始化种群个体中每个卷积神经网络模型的适应度值；
优化模型确定单元，用于若存在适应度值小于预设分类误差值或当前迭代次数大于或等于预设迭代次数，则将最大的适应度值对应的卷积神经网络模型确定为卷积神经网络优化模型。
2.根据权利要求1所述的一种麻醉深度估计系统，其特征在于，所述权重调节模块，具体包括：
第二计算单元，用于将所述频谱特征图作为所述卷积神经网络优化模型的输入，前向计算所述卷积神经网络优化模型中每个神经元的输出值；第j个神经元的输出值为aj；
第三计算单元，用于反向计算所述卷积神经网络优化模型中每个神经元的误差项；第j个神经元的误差项为
Figure FDA0002766633020000021
其中，Ed为卷积神经网络优化模型的损失函数，netj为第j个神经元的加权输入；
第四计算单元，用于依据各神经元的输出值以及各神经元的误差项计算各个神经元连接权重的梯度值；从神经元i连接到神经元j的连接权重wji的梯度值
Figure FDA0002766633020000022
其中，ai为第i个神经元的输出值；
调节单元，用于依据所述梯度值，采用梯度下降法对所述卷积神经网络优化模型中的权重进行调节，得到训练好的卷积神经网络优化模型。
3.根据权利要求1所述的一种麻醉深度估计系统，其特征在于，所述系统还包括：预处理模块，用于对所述脑电信号训练数据进行预处理；
所述预处理模块，具体包括：
第一处理单元，用于采用第一高通滤波器去除所述脑电信号训练数据中的头动噪声，得到第一处理脑电信号；
第二处理单元，用于采用自适应滤波器去除所述第一处理脑电信号中的工频噪声，得到第二处理脑电信号；
第三处理单元，用于采用低通滤波器去除所述第二处理脑电信号中大于第一预设频率的噪声信号，得到第三处理脑电信号；
第四处理单元，用于采用第二高通滤波器去除所述第三处理脑电信号中小于第二预设频率的噪声信号，得到第四处理脑电信号。
CN201811502030.7A 2018-12-10 2018-12-10 一种麻醉深度估计系统 Active CN109645989B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811502030.7A CN109645989B (zh) 2018-12-10 2018-12-10 一种麻醉深度估计系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811502030.7A CN109645989B (zh) 2018-12-10 2018-12-10 一种麻醉深度估计系统

Publications (2)

Publication Number Publication Date
CN109645989A CN109645989A (zh) 2019-04-19
CN109645989B true CN109645989B (zh) 2021-01-08

Family

ID=66112997

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811502030.7A Active CN109645989B (zh) 2018-12-10 2018-12-10 一种麻醉深度估计系统

Country Status (1)

Country Link
CN (1) CN109645989B (zh)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110232961A (zh) * 2019-06-10 2019-09-13 上海电机学院 一种基于大数据的语音识别智能麻醉系统
CN110811556A (zh) * 2019-11-15 2020-02-21 西安交通大学 一种基于脑电微状态分析的麻醉深度监测系统及方法
CN110811557A (zh) * 2019-11-15 2020-02-21 西安交通大学 一种基于微状态功率谱分析的麻醉深度监测系统及方法
CN110840411B (zh) * 2019-12-06 2022-03-11 深圳市德力凯医疗设备股份有限公司 一种麻醉深度的测量装置、存储介质及电子设备
CN112085174A (zh) * 2020-10-21 2020-12-15 广东工业大学 一种卷积神经网络优化方法、装置和设备
CN112657028A (zh) * 2020-12-03 2021-04-16 遵义医科大学 一种基于物联网的实验室麻醉操作模拟控制系统及方法
CN113208624A (zh) * 2021-04-07 2021-08-06 北京脑陆科技有限公司 一种基于卷积神经网络的疲劳检测方法、系统
CN113491506B (zh) * 2021-08-06 2024-05-17 中科搏锐（北京）科技有限公司 一种基于深度卷积网络的脑电信号处理方法及其设备
CN115105681B (zh) * 2022-06-24 2024-05-24 纽锐思(苏州)医疗科技有限公司 基于多模态生理指标的麻醉闭环输注系统
WO2024196928A1 (en) * 2023-03-21 2024-09-26 The Board Of Regents, The University Of Texas System Method and system of monitoring anesthesia depth
CN118452838B (zh) * 2024-07-09 2024-10-15 上海岩思类脑人工智能研究院有限公司 一种基于神经尖峰编码的麻醉状态检测方法及装置
CN118468166A (zh) * 2024-07-12 2024-08-09 福建医科大学附属第一医院 一种基于人工智能的麻醉深度分类方法及系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2002032305A1 (en) * 2000-10-16 2002-04-25 Instrumentarium Corporation Method and apparatus for determining the cerebral state of a patient with fast response
CN103153178A (zh) * 2010-07-23 2013-06-12 昆腾医疗公司 用于结合麻醉药和止痛药的药物效应相互作用与评估麻醉期间意识水平的脑电图的装置
CN103637798A (zh) * 2013-12-17 2014-03-19 山东大学齐鲁医院 一种麻醉深度监测装置

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101111498B1 (ko) * 2010-07-19 2012-02-22 주식회사 멕 아이씨에스 마취 심도 모니터링 시스템 및 방법
KR101400362B1 (ko) * 2013-02-18 2014-05-30 한국과학기술원 마취 심도 측정 방법 및 장치
BR102013006865B1 (pt) * 2013-03-25 2022-03-03 Universidade Estadual De Campinas - Unicamp Estimulador de nervos periféricos automatizado e método de automatização
CN104887225B (zh) * 2015-06-04 2017-10-10 卞汉道 麻醉精度监护仪器及方法
CN105303252A (zh) * 2015-10-12 2016-02-03 国家计算机网络与信息安全管理中心 基于遗传算法的多阶段神经网络模型训练方法
US11266780B2 (en) * 2015-11-23 2022-03-08 University Of British Columbia Methods and systems for closed-loop control of drug administration
CN106902421B (zh) * 2017-01-17 2020-03-31 燕山大学 一种个人化麻醉闭环控制系统
CN106821337B (zh) * 2017-04-13 2020-01-03 南京理工大学 一种有监督的鼾声来源识别方法
CN108378845A (zh) * 2018-01-29 2018-08-10 湖北民族学院附属民大医院 一种多功能麻醉科用麻醉深度监测装置
CN108415560B (zh) * 2018-02-11 2020-12-04 Oppo广东移动通信有限公司 电子装置、操作控制方法及相关产品

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2002032305A1 (en) * 2000-10-16 2002-04-25 Instrumentarium Corporation Method and apparatus for determining the cerebral state of a patient with fast response
CN103153178A (zh) * 2010-07-23 2013-06-12 昆腾医疗公司 用于结合麻醉药和止痛药的药物效应相互作用与评估麻醉期间意识水平的脑电图的装置
CN103637798A (zh) * 2013-12-17 2014-03-19 山东大学齐鲁医院 一种麻醉深度监测装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
丙泊酚对癫痫患者皮层及皮层下脑电活动的影响;顾娟娟;《中国优秀硕士论文全文数据库医药卫生科技辑》;20180315(第3期);E070-26 *
基于卷积神经网络的运动想象脑电信号特征提取与分类;范明莉;《中国优秀硕士论文全文数据库医药卫生科技辑》;20171115(第11期);E080-6 *

Also Published As

Publication number Publication date
CN109645989A (zh) 2019-04-19

Similar Documents

Publication Publication Date Title
CN109645989B (zh) 一种麻醉深度估计系统
CN109389059B (zh) 一种基于cnn-lstm网络的p300检测方法
KR102221264B1 (ko) 인간 감정 인식을 위한 딥 생리적 정서 네트워크를 이용한 인간 감정 추정 방법 및 그 시스템
CN106503799B (zh) 基于多尺度网络的深度学习模型及在脑状态监测中的应用
CN110353702A (zh) 一种基于浅层卷积神经网络的情感识别方法及系统
CN110353673B (zh) 一种基于标准互信息的脑电通道选择方法
Übeyli Analysis of EEG signals by implementing eigenvector methods/recurrent neural networks
CN114533086B (zh) 一种基于空域特征时频变换的运动想象脑电解码方法
CN111407243B (zh) 一种基于深度学习的脉搏信号压力识别方法
CN113705398B (zh) 基于卷积-长短期记忆网络音乐脑电时空特征分类方法
CN110135244B (zh) 一种基于脑-机协同智能的表情识别方法
Thenmozhi et al. Feature selection using extreme gradient boosting Bayesian optimization to upgrade the classification performance of motor imagery signals for BCI
CN111091074A (zh) 最优区域共空间模式的运动想象脑电信号分类方法
CN113017645B (zh) 一种基于空洞卷积神经网络的p300信号检测方法
CN115969329B (zh) 一种睡眠分期方法、系统、装置及介质
CN111950441A (zh) 上肢运动意图的fnirs实时解码方法及系统
CN116522106A (zh) 一种基于迁移学习并行多尺度滤波器组时域卷积的运动想象脑电信号分类方法
CN117883082A (zh) 一种异常情绪识别方法、系统、设备及介质
Al-dabag et al. EEG motor movement classification based on cross-correlation with effective channel
Khalkhali et al. Low latency real-time seizure detection using transfer deep learning
CN117725367A (zh) 一种源域可迁移性自适应学习的言语想象脑电解码方法
Wang et al. A shallow convolutional neural network for classifying MI-EEG
CN115905819B (zh) 基于生成对抗网络的rPPG信号生成方法及装置
Veerabhadrappa et al. A deep convolutional neural network model for classification of emotions from electroencephalography data
He et al. HMT: An EEG Signal Classification Method Based on CNN Architecture

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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20220714

Address after: Room 101, building 13, No. 158, Jinfeng Road, science and Technology City, high tech Zone, Suzhou, Jiangsu 215163

Patentee after: Jiangsu University of technology collaborative medical robot Co.,Ltd.

Address before: 066000 No. 438, Hebei Avenue, Qinhuangdao, Hebei

Patentee before: Yanshan University

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20220818

Address after: Room 301-6, Building 13, No. 158, Jinfeng Road, Science and Technology City, High-tech Zone, Suzhou City, Jiangsu Province 215151

Patentee after: Newrise (Suzhou) Medical Technology Co., Ltd.

Address before: Room 101, building 13, No. 158, Jinfeng Road, science and Technology City, high tech Zone, Suzhou, Jiangsu 215163

Patentee before: Jiangsu University of technology collaborative medical robot Co.,Ltd.