CN104644166B - 一种基于格子复杂性算法的无线动态麻醉深度检测方法 - Google Patents

一种基于格子复杂性算法的无线动态麻醉深度检测方法 Download PDF

Info

Publication number
CN104644166B
CN104644166B CN201510085532.4A CN201510085532A CN104644166B CN 104644166 B CN104644166 B CN 104644166B CN 201510085532 A CN201510085532 A CN 201510085532A CN 104644166 B CN104644166 B CN 104644166B
Authority
CN
China
Prior art keywords
module
symbol
processing module
sequence
frequency
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
CN201510085532.4A
Other languages
English (en)
Other versions
CN104644166A (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.)
JIANGSU APON MEDICAL TECHNOLOGY Co.,Ltd.
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN201510085532.4A priority Critical patent/CN104644166B/zh
Publication of CN104644166A publication Critical patent/CN104644166A/zh
Application granted granted Critical
Publication of CN104644166B publication Critical patent/CN104644166B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
    • 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/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
    • A61B5/7217—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal of noise originating from a therapeutic or surgical apparatus, e.g. from a pacemaker
    • 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/7225—Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • Surgery (AREA)
  • Medical Informatics (AREA)
  • Psychiatry (AREA)
  • Physics & Mathematics (AREA)
  • Veterinary Medicine (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Public Health (AREA)
  • Molecular Biology (AREA)
  • General Health & Medical Sciences (AREA)
  • Animal Behavior & Ethology (AREA)
  • Artificial Intelligence (AREA)
  • Physiology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Power Engineering (AREA)
  • Psychology (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

本发明公开一种基于格子复杂性算法的无线动态麻醉深度检测方法。该发明方法是基于以下麻醉深度监测装置，该装置包括采集模块、无线传输模块、数据处理模块、中央处理模块、显示模块；采集模块采集脑电信号对该模拟脑电信号进行放大，并将该将其转换为数字脑电信号。无线传输模块通过WiFi无线网络，将采集模块输出的数字脑电信号数据上传到数据处理模块进行预处理；预处理后的数据输入中央处理模块进行分析处理，分别计算格子复杂度、边缘频率、爆发抑制比，并利用决策树算法拟合得到麻醉深度指数并显示。由于使用了无线传输模块和格子复杂性算法，使得该发明方法省去传统检测装置的导线所带来的不便，且算法简单，实时性强。

Description

一种基于格子复杂性算法的无线动态麻醉深度检测方法
技术领域
本发明涉及一种麻醉深度检测装置，尤其涉及一种包含格子复杂性算法及无线传输的用于重症监护病房的麻醉深度检测装置。
背景技术
麻醉是指借助于药物等方法而产生的全身或局部感觉的消失及记忆遗忘状态，它可以确保手术的顺利进行，麻醉过深或过浅都会对患者造成危害。因此，麻醉深度的监测尤为重要。麻醉是通过引起可逆的中枢神经系统的抑制和兴奋，从而达到意识消失和止痛的目的。而脑电可以直接反映出中枢神经系统的活动。因此脑电技术成为确定麻醉深度的最佳手段之一。脑电的许多时域、频域参数都曾被用来检测麻醉深度，但效果都不是很理想。
时域分析是最原始的脑电分析方法,主要是对脑电波波形和幅度的检测分析。时域分析大部分只停留在定性分析上,且系统复杂、计算量大、效果也多不理想，因此限制了其在临床上的进一步应用。
脑电图的频域分析主要包括脑电功率谱、各频率段的脑电功率比、中心频率(MF)、边缘频率(SEF)脚等。MF与SEF的单参数定量化特性也很适合麻醉深度的定量监测研究,大量研究表明在某些药物麻醉中MF和SEF随着麻醉深度的加深作相应的变化。但同时很多研究也证明MF和SEF的个体差异和药物差异性较大。另外其分析基础一傅立叶变换是基于信号的平稳性假设的,这与脑电信号的特性并不符合,因此其可靠性已越来越受到许多学者的质疑。
涉及时域、频域及双谱域的复合指数一一双谱指数(BiPsecrtalnIdexScale,BIS)是目前麻醉深度监测中近年来最受欢迎的一个的EEG参数。它除了表达EEG信号时域和频域信息外,更重要的是增加了相位信息,具有非线性特点。BIS是唯一被美国食物药品管理局认可的麻醉药对大脑作用的监测仪,是目前商业化麻醉深度监测仪中敏感度和特异度量较好的监测仪之一26[,27]。尽管这样,BIS的药物差异性仍然较大,缺乏明确的麻醉安全边界指标,还不能临床推广。
复杂度是近几年出现的一种非线性动力学分析方法，而脑电正是一种非平稳信号，所以复杂度的分析犯法非常适合用于脑电的处理。脑电序列的复杂度表现了EEG序列的随机程度，即大脑神经元处理信息活动的有序程度，反应了决定这段EEG序列的信息量的大小。复杂度算法简单，易于实现而且计算速度很快；只需较短时间的数据就能得到稳定值。已有的近似熵等复杂性算法由于计算所需序列长度长或计算所需时间长的缺点无法实现实时监测。
因此本技术领域人员需要一个问题是：如何找到一种新的复杂性参数来提高麻醉深度监测的实时性和准确性。
同时考虑到设备用于重症监护病房，以往的有线监测显得不便且经济成本较高。这也是本技术领域人员需要考虑的。
发明内容
本发明的目的就是为了解决上述问题，提供一种基于格子复杂性算法的无线动态麻醉深度检测方法，它具有监测方便，准确的优点。
本发明方法是基于以下麻醉深度监测装置，该装置包括采集模块、无线传输模块、数据处理模块、中央处理模块、显示模块；采集模块的输出端与无线传输模块的接收端信号连接，无线传输模块的输出端通过无线WIFI与数据处理模块的接收端连接，数据处理模块的输出端与中央处理模块的接收端信号连接，中央处理模块的输出端与显示模块的输入端信号连接；
所述的采集模块，包括脑电信号采集电极和依次连接的前置放大模块和模数转换模块；所述前置放大模块用于接收脑电信号采集电极所采集的模拟脑电信号，并对该模拟脑电信号进行放大；所述模数转换模块用于接收前置放大模块放大后的模拟脑电信号，将模拟脑电信号转换为数字脑电信号；
所述的无线传输模块，通过WiFi无线网络将采集模块采集到的数字脑电信号数据上传到数据处理模块；
所述的数据处理模块用于对无线传输模块传输的脑电信号数据进行存储并预处理，所述的预处理即为对数字脑电信号数据进行去噪处理；
由于未预处理的脑电信号数据包含工频、眼电、电刀等各类噪声，不可直接用于脑电信号分析，故需要去除这些噪声；
所述的中央处理模块包括参数计算模块、决策树分类器模块；
其中所述的参数计算模块，包括格子复杂性、边缘频率、爆发抑制比三个参数的计算；
所述的格子复杂性用于评价麻醉深度，尤其可以区分病人清醒与麻醉时期；但是若病人处于麻醉时期，本发明采用边缘频率判断其麻醉深度，边缘频率在中度麻醉时期变化最为剧烈，可以较好的区分浅麻、中麻和深麻(分别代表浅度麻醉、中度麻醉和深度麻醉)；在深度麻醉时期，由于爆发抑制比(即幅度变化特征)最为明显，可作为此时麻醉深度量化的指标；
本发明方法包括以下步骤：
步骤(1).利用采集模块脑电信号采集电极采集模拟脑电信号，经前置放大模块对该模拟脑电信号放大后，模数转换模块将其转换为数字脑电信号；
步骤(2).无线传输模块通过WiFi无线网络，将采集模块采集到的数字脑电信号数据上传到数据处理模块；数据处理模块对数据进行存储并预处理，所述的预处理即为对数字脑电信号数据进行去噪处理；
步骤(3).中央处理模块对于步骤(2)数据处理模块预处理后的数据进行分析处理：
3.1格子复杂度的计算：
3.1.1对步骤(2)数据处理模块预处理后的数据进行符号化处理，得到符号序列；其中符号化处理方法为均值分划或排列分划；
3.1.2对符号序列按照窗口大小，结合窗口移动长度进行分段，得到多个子序列；
3.1.3对每个子序列P进行格子复杂性计算：
本发明设子序列P＝s1s2s3...sn，n表示子序列P的长度；lc(n)表示格子复杂度即格子数；S序列，Q序列都是P序列的子序列，SQ是S序列和Q序列的合并序列，SQπ序列是SQ合并序列去掉最后一个符号后剩余部分；
假设P序列中S前r个符号组成的子序列S已经进行格子化操作，其中S＝s1s2...sr，r<n,则下一个格子划分如下：
(a)首先要找到一个序列Q满足Q序列的最后一个符号等于Qπ序列中的任意一个符号(Qπ序列表示将Q序列符号串最后一个符号去掉剩余部分)：令Q＝sr+1，因为Q只有一个符号，不做任何操作；继续令Q＝sr+1sr+2，观察Q的最后一个符号是否等于Qπ中的任意一个符号(Qπ表示将Q符号串最后一个符号去掉剩下的部分，这里有Qπ＝sr+1).若否，重置Q，令Q＝sr+ 1sr+2sr+3，继续观察Q最后一个符号是否等于Qπ中的任意一个符号(此时Qπ＝sr+1sr+2)，若否，则继续重置Q，直到Q的最后一个符号可以由Qπ中的一个符号复制为止，接着进行下一步步骤(b)；
(b)若满足步骤(a)后Q＝sr+1sr+2...sr+j，Qπ中的相等符号为sr+i(0<i<j)，继续观察sr+j+1是否等于sr+i+1，若是则继续观察下一个符号是否相等(即sr+j+2是否等于sr+i+2)，直到两者不相等为止，进行下一步步骤(c)。
(c)若满足步骤(b)后不相等的符号为sr+j+k,0<k<n-r-j，令Q＝sr+1sr+2...sr+j+k，观察Q能否由SQπ的某个子串复制得到，其中SQπ＝s1s2s3...sr+j+k-1；若能则继续令Q＝sr+1sr+ 2sr+3...sr+j+ksr+j+k+1，观察Q能否由SQπ中的子串复制得到，若能则再令Q＝sr+1sr+2sr+ 3...sr+j+ksr+j+k+1sr+j+k+2,直到Q最后一个符号为sn时则格子数lc(n)＝lc(n)+1,同时结束操作；若Q不能由SQπ的某个子串复制得到，则格子数lc(n)＝lc(n)+1，继续执行步骤(a)。
3.2边缘频率的计算：
所述的边缘频率是指一段信号的功率谱中，功率积分达到总功率的95％时的频率；
3.2.1对步骤(2)数据处理模块预处理后的数据xk,k＝1,2,...,N进行离散傅里叶变换，见公式(1)：
3.2.2由于X[m]中前N/2点包含了步骤(2)数据处理模块预处理后脑电信号数据的频率信息，利用此脑电信号的采样率fs与频率fi之间的关系式(2)，将X[m]转换为关于频率fi的表达式(3)：
m＝N*fi/fs,m＝1,2,...,N/2 (2)
X(fi)＝X[m]＝X(N*fi/fs) (3)
在X(fi)中寻找某一频率值fsef，使其满足以下关系式：
此时频率值fsef即边缘频率；
3.3爆发抑制比(BSR)的计算：
所述的爆发抑制比(BSR)为深麻时期内抑制状态占一段EEG信号的比例；抑制状态是指当EEG电位不大于±5uV，时间超过0.5s；
3.4决策树分类器的训练：
所述的决策树分类器是将格子复杂度、边缘频率和爆发抑制比作为决策树算法的输入特征，将麻醉专家给出的麻醉状态分类作为输出，训练出一棵决策树；
所述的决策树的训练采用以下三个终止条件，满足任一终止条件则训练结束：
①信息增益为0时；
②输入特征和麻醉状态分类的集合D中某一类数据的比例大于阈值时；其中阈值根据本领域技术人员经验设定；
③集合D中的数据小于阈值时，这时用D中出现最多的类来统一判断该节点的数据。
训练结束后根据实际情况对决策树进行剪枝，剪枝的目的是：1.在保证精度的同时尽量使决策树简单，以提高推广能力；2.排除过度学习的可能，尽可能的使分类器与生理知识保持一致。
3.5获得决策树输出的麻醉深度类别后，利用格子复杂度、边缘频率和爆发抑制比三个参数，拟合现有的麻醉量化指标，得到此时病人的麻醉深度指数。
本发明的有益效果：
本发明采用无线传输模块的应用使得该方法省去传统检测装置的导线所带来的不便；采用格子复杂性算法使得整个算法的复杂度降低，提高了实时监测性。
附图说明
图1为本发明方法流程图；
图2为电极的贴放位置；
图3为决策树结构；
图4为一例典型患者麻醉过程的本发明监测指数示意图；
图5为一例典型患者麻醉过程中的BIS指数变化。
具体实施方式
下面结合附图与实施例对本发明作进一步说明。
本发明方法是基于以下麻醉深度监测装置，该装置包括采集模块、无线传输模块、数据处理模块、中央处理模块、显示模块；采集模块的输出端与无线传输模块的接收端信号连接，无线传输模块的输出端通过无线WIFI与数据处理模块的接收端连接，数据处理模块的输出端与中央处理模块的接收端信号连接，中央处理模块的输出端与显示模块的输入端信号连接；
所述的采集模块，包括脑电信号采集电极和依次连接的前置放大模块和模数转换模块；脑电信号采集电极贴于病人头部特定的四个位置(如图2所示)，以128HZ采样频率采集数据；所述前置放大模块用于接收脑电信号采集电极所采集的模拟脑电信号，并对该模拟脑电信号进行放大；所述模数转换模块用于接收前置放大模块放大后的模拟脑电信号，将模拟脑电信号转换为数字脑电信号；
所述的无线传输模块，通过WiFi无线网络将采集模块采集到的数字脑电信号数据上传到数据处理模块；
所述的数据处理模块用于对无线传输模块传输的脑电信号数据进行存储并预处理，所述的预处理即为对数字脑电信号数据进行去噪处理；由于未预处理的脑电信号数据包含工频、眼电、电刀等各类噪声，不可直接用于脑电信号分析，故需要去除这些噪声；
本发明去噪处理是数据处理模块利用截止频率为1～47Hz的数字滤波器对脑电信号滤波，然后人为去除滤波无法去除的眼电伪迹和电刀干扰；
所述的中央处理模块包括参数计算模块、决策树分类器模块；
其中所述的参数计算模块，包括格子复杂性、边缘频率、爆发抑制比三个参数的计算；
所述的格子复杂性用于评价麻醉深度，尤其可以区分病人清醒与麻醉时期；但是若病人处于麻醉时期，本发明采用边缘频率判断其麻醉深度，边缘频率在中度麻醉时期变化最为剧烈，可以较好的区分浅麻、中麻和深麻(分别代表浅度麻醉、中度麻醉和深度麻醉)；在深度麻醉时期，由于爆发抑制比(即幅度变化特征)最为明显，可作为此时麻醉深度量化的指标；
如图1所示，本发明方法包括以下步骤：
步骤(1).利用采集模块脑电信号采集电极采集模拟脑电信号，经前置放大模块对该模拟脑电信号放大后，模数转换模块将其转换为数字脑电信号；
步骤(2).无线传输模块通过WiFi无线网络，将采集模块采集到的数字脑电信号数据上传到数据处理模块；数据处理模块对数据进行存储并预处理，所述的预处理即为对数字脑电信号数据进行去噪处理；
步骤(3).中央处理模块对于步骤(2)数据处理模块预处理后的数据进行分析处理：
3.1格子复杂度的计算：
3.1.1对步骤(2)数据处理模块预处理后的数据进行符号化处理，得到符号序列；其中符号化处理方法为均值分划或排列分划；
3.1.2对符号序列按照窗口大小，结合窗口移动长度进行分段，得到多个子序列；
3.1.3对每个子序列P进行格子复杂性计算：
本发明设子序列P＝s1s2s3...sn，n表示子序列P的长度；lc(n)表示格子复杂度即格子数；S序列，Q序列都是P序列的子序列，SQ是S序列和Q序列的合并序列，SQπ序列是SQ合并序列去掉最后一个符号后剩余部分；
假设P序列中S前r个符号组成的子序列S已经进行格子化操作，其中S＝s1s2...sr，r<n,则下一个格子划分如下：
(a)首先要找到一个序列Q满足Q序列的最后一个符号等于Qπ序列中的任意一个符号(Qπ序列表示将Q序列符号串最后一个符号去掉剩余部分)：令Q＝sr+1，因为Q只有一个符号，不做任何操作；继续令Q＝sr+1sr+2，观察Q的最后一个符号是否等于Qπ中的任意一个符号(Qπ表示将Q符号串最后一个符号去掉剩下的部分，这里有Qπ＝sr+1).若否，重置Q，令Q＝sr+ 1sr+2sr+3，继续观察Q最后一个符号是否等于Qπ中的任意一个符号(此时Qπ＝sr+1sr+2)，若否，则继续重置Q，直到Q的最后一个符号可以由Qπ中的一个符号复制为止，接着进行下一步步骤(b)；
(b)若满足步骤(a)后Q＝sr+1sr+2...sr+j，Qπ中的相等符号为sr+i(0<i<j)，继续观察sr+j+1是否等于sr+i+1，若是则继续观察下一个符号是否相等(即sr+j+2是否等于sr+i+2)，直到两者不相等为止，进行下一步步骤(c)。
(c)若满足步骤(b)后不相等的符号为sr+j+k,0<k<n-r-j，令Q＝sr+1sr+2...sr+j+k，观察Q能否由SQπ的某个子串复制得到，其中SQπ＝s1s2s3...sr+j+k-1；若能则继续令Q＝sr+1sr+ 2sr+3...sr+j+ksr+j+k+1，观察Q能否由SQπ中的子串复制得到，若能则再令Q＝sr+1sr+2sr+ 3...sr+j+ksr+j+k+1sr+j+k+2,直到Q最后一个符号为sn时则格子数lc(n)＝lc(n)+1,同时结束操作；若Q不能由SQπ的某个子串复制得到，则格子数lc(n)＝lc(n)+1，继续执行步骤(a)。
3.2边缘频率的计算：
所述的边缘频率是指一段信号的功率谱中，功率积分达到总功率的95％时的频率；
3.2.1对步骤(2)数据处理模块预处理后的数据xk,k＝1,2,...,N进行离散傅里叶变换，见公式(1)：
3.2.2由于X[m]中前N/2点包含了步骤(2)数据处理模块预处理后脑电信号数据的频率信息，利用此脑电信号的采样率fs与频率fi之间的关系式(2)，将X[m]转换为关于频率fi的表达式(3)：
m＝N*fi/fs,m＝1,2,...,N/2； (2)
X(fi)＝X[m]＝X(N*fi/fs) (3)
在X(fi)中寻找某一频率值fsef，使其满足以下关系式：
此时频率值fsef即边缘频率；
3.3爆发抑制比(BSR)的计算：
所述的爆发抑制比(BSR)为深麻时期内抑制状态占一段EEG信号的比例；抑制状态是指当EEG电位不大于±5uV，时间超过0.5s；
例如一段深麻脑电信号长度为N，电位处于±5uV且时间超过0.5s的数据段累计长度为m,m小于等于N；则爆发抑制比BSR＝m/N。
3.4决策树分类器的训练：
所述的决策树分类器是将格子复杂度、边缘频率和爆发抑制比作为决策树算法的输入特征，将麻醉专家给出的麻醉状态分类作为输出，训练出一棵决策树；
所述的决策树的训练采用以下三个终止条件：
⑴信息增益为0时；
⑵D(输入特征和麻醉状态分类的集合)中某一类数据的比例大于一定量时；
⑶D中的数据少于一定量时，这时用D中出现最多的类来统一判断该节点的数据；
训练结束后根据实际情况对决策树进行剪枝，剪枝的目的是：1.在保证精度的同时尽量使决策树简单，以提高推广能力；2.排除过度学习的可能，尽可能的使分类器与生理知识保持一致。最终得到的决策树结构如图3所示，图中LC、Sef、BSR分别表示格子复杂度、边缘频率与爆发抑制比，Th为阈值，1、2、3、4分别代表清醒、浅麻、中麻和深麻。
3.5获得决策树输出的麻醉深度类别后，利用格子复杂度、边缘频率和爆发抑制比三个参数，拟合现有的麻醉量化指标，得到此时病人的麻醉深度指数。
由于在各麻醉深度类别中麻醉深度与格子复杂度、边缘频率和爆发抑制比等参数基本呈线性关系，故采用最小二乘法进行线性拟合。这种方法虽没有一些复杂的拟合方法精度高，但是它的简单性会在一定程度上保证算法的推广能力。
临床试验实施例：
将上述装置以及检测方法应用于临床，观察复合麻醉下经典麻醉深度指标BIS和本发明麻醉深度指标的异同。
选取全身麻醉下手术患者30例，年龄24～73岁，身高155～175cm，体重52～80kg，性别不限，ASA体格状况分级为Ⅰ～Ⅱ级，排除颅面部手术者、精神疾病病史者及无民事行为能力者。
患者入手术室后确定传感器的粘贴位置，并将传感器固定于被测者头部，在另一侧按照说明书安放BIS监护仪专用脑电电极，启动本实施例和BIS监护仪采集分析脑电信号。麻醉诱导期依次注入镇静类药物咪唑安定0.1mg*kg-1，丙泊酚1.5～2mg*kg-1，肌松类药物维库溴铵0.1mg*kg-1，镇痛类药物芬太尼5ug*kg-1。麻醉维持期用药如下：每小时追加2mg维库溴铵，0.1mg芬太尼，以30～50ml/h速率恒速泵入10mg/ml浓度的丙泊酚，10～25ml/h速率恒速泵入20mg/50ml浓度的瑞芬太尼。记录整个麻醉过程中的脑电信号变化。
临床试验结果
本实施例和BIS监护仪显示的麻醉指数分别如图4、图5所示，从图中可以看出实施例的灵敏度较高，能较好的识别麻醉和清醒的转换。这是由于格子复杂性算法对混沌边缘比较敏感，而诱导期和恢复期正对应此状态。在麻醉期间，实施例相较于BIS整体波动较小，局部波动较大，这和格子复杂性均值分化的符号化方法有关，这种方法使得其对细节的刻画较差。
本发明的主要结论
本发明将格子复杂性这一参数用于麻醉深度监测，并与BIS监护仪监测结果进行比较，证明了格子复杂性可作为识别麻醉与清醒的主要参数。本发明的临床试验表明，格子复杂性能够得到比BIS更高的灵敏度，本发明可作为临床麻醉检测的一种手段。

Claims (2)

1.一种基于格子复杂性算法的无线动态麻醉深度检测方法，是基于以下麻醉深度监测装置，该装置包括采集模块、无线传输模块、数据处理模块、中央处理模块、显示模块；采集模块的输出端与无线传输模块的接收端信号连接，无线传输模块的输出端通过无线WIFI与数据处理模块的接收端连接，数据处理模块的输出端与中央处理模块的接收端信号连接，中央处理模块的输出端与显示模块的输入端信号连接，其特征在于该方法包括以下步骤：
步骤(1).利用采集模块脑电信号采集电极采集模拟脑电信号，经前置放大模块对该模拟脑电信号放大后，模数转换模块将其转换为数字脑电信号；
步骤(2).将采集模块采集到的数字脑电信号数据通过无线传输模块，上传到数据处理模块；数据处理模块对数据进行存储并预处理，所述的预处理即为对数字脑电信号数据进行去噪处理；
步骤(3).中央处理模块对于步骤(2)数据处理模块预处理后的数据进行分析处理：
3.1格子复杂度的计算：
3.1.1对步骤(2)数据处理模块预处理后的数据进行符号化处理，得到符号序列；
3.1.2对符号序列按照窗口大小，结合窗口移动长度进行分段，得到多个子序列；
3.1.3对每个子序列P进行格子复杂性计算：
假设子序列P＝s1s2s3...sn，n表示子序列P的长度；lc(n)表示格子复杂度即格子数；S序列，Q序列都是P序列的子序列，SQ是S序列和Q序列的合并序列，SQπ序列是SQ合并序列去掉最后一个符号后剩余部分；
假设P序列中S前r个符号组成的子序列S已经进行格子化操作，其中S＝s1s2...sr，r<n,则下一个格子划分如下：
(a)首先要找到一个序列Q,满足Q序列的最后一个符号等于Qπ序列中的任意一个符号,其中Qπ序列表示将Q序列符号串最后一个符号去掉剩余部分：令Q＝sr+1，因为Q只有一个符号，不做任何操作；继续令Q＝sr+1sr+2，观察Q的最后一个符号是否等于Qπ中的任意一个符号，若否，重置Q，令Q＝sr+1sr+2sr+3，继续观察Q最后一个符号是否等于Qπ中的任意一个符号，若否，则继续重置Q，直到Q的最后一个符号可以由Qπ中的一个符号复制为止，接着进行下一步步骤(b)；
(b)若满足步骤(a)后Q＝sr+1sr+2...sr+j，Qπ中的相等符号为sr+i(0<i<j)，继续观察sr+j+1是否等于sr+i+1，若是则继续观察下一个符号是否相等，直到两者不相等为止，进行下一步步骤(c)；
(c)若满足步骤(b)后不相等的符号为sr+j+k，0<k<n-r-j，令Q＝sr+1sr+2...sr+j+k，观察Q能否由SQp的某个子串复制得到，其中SQp＝s1s2s3...sr+j+k-1；若能则继续令Q＝sr+1sr+2sr+ 3...sr+j+ksr+j+k+1，观察Q能否由SQp中的子串复制得到，若能则再令Q＝sr+1sr+2sr+3...sr+j+ ksr+j+k+1sr+j+k+2，直到Q最后一个符号为sn时，则格子数lc(n)＝lc(n)+1，同时结束操作；若Q不能由SQp的某个子串复制得到，则格子数lc(n)＝lc(n)+1，继续执行步骤(a)；
3.2边缘频率的计算：
所述的边缘频率是指一段信号的功率谱中，功率积分达到总功率的95％时的频率；
3.2.1对步骤(2)数据处理模块预处理后的数据xk，k＝1,2,...,N进行离散傅里叶变换，见公式(1)：
3.2.2由于X[m]中前N/2点包含了步骤(2)数据处理模块预处理后脑电信号数据的频率信息，利用此脑电信号的采样率fs与频率fi之间的关系式(2)，将X[m]转换为关于频率fi的表达式(3)：
m＝N*fi/fs,m＝1,2,...,N/2 (2)
X(fi)＝X[m]＝X(N*fi/fs) (3)
在X(fi)中寻找某一频率值fsef，使其满足以下关系式：
此时频率值fsef即边缘频率；
3.3爆发抑制比的计算：
所述的爆发抑制比为深麻时期内抑制状态占一段脑电信号的比例；抑制状态是指当脑电信号电位不大于±5uV，时间超过0.5s；
3.4决策树分类器的训练：
所述的决策树分类器是将格子复杂度、边缘频率和爆发抑制比作为决策树算法的输入特征，将麻醉专家给出的麻醉状态分类作为输出，训练出一棵决策树；
3.5获得决策树输出的麻醉深度类别后，利用格子复杂度、边缘频率和爆发抑制比三个参数，拟合现有的麻醉量化指标，得到此时病人的麻醉深度指数。
2.如权利要求1所述的一种基于格子复杂性算法的无线动态麻醉深度检测方法，其特征在于步骤3.1.1所述的符号化处理采用均值分划或排列分划方法。
CN201510085532.4A 2015-02-16 2015-02-16 一种基于格子复杂性算法的无线动态麻醉深度检测方法 Active CN104644166B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510085532.4A CN104644166B (zh) 2015-02-16 2015-02-16 一种基于格子复杂性算法的无线动态麻醉深度检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510085532.4A CN104644166B (zh) 2015-02-16 2015-02-16 一种基于格子复杂性算法的无线动态麻醉深度检测方法

Publications (2)

Publication Number Publication Date
CN104644166A CN104644166A (zh) 2015-05-27
CN104644166B true CN104644166B (zh) 2017-03-08

Family

ID=53236177

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510085532.4A Active CN104644166B (zh) 2015-02-16 2015-02-16 一种基于格子复杂性算法的无线动态麻醉深度检测方法

Country Status (1)

Country Link
CN (1) CN104644166B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110234272B (zh) 2017-12-29 2022-06-14 深圳迈瑞生物医疗电子股份有限公司 麻醉阶段识别与麻醉深度计算方法及装置
WO2019127558A1 (zh) * 2017-12-29 2019-07-04 深圳迈瑞生物医疗电子股份有限公司 基于脑电的麻醉深度监测方法和装置
DE102018110275A1 (de) * 2018-04-27 2019-10-31 Susanne Koch Verfahren und Vorrichtung zur Bereitstellung eines Parameters, der auf einen Bewusstseinsverlust eines Patienten unter Narkose hinweist
CN109394210A (zh) * 2018-12-14 2019-03-01 张维维 一种麻醉科用麻醉深度监测装置
CN112438701A (zh) * 2019-08-16 2021-03-05 浙江象立医疗科技有限公司 一种基于光电融合的血液参数和电生理参数的检测方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5252383B2 (ja) * 2007-02-06 2013-07-31 公立大学法人広島市立大学 生体内電流双極子の推定方法
US20110015503A1 (en) * 2009-07-17 2011-01-20 WAVi Medical apparatus for collecting patient electroencephalogram (eeg) data
CN103876734B (zh) * 2014-03-24 2015-09-02 北京工业大学 一种基于决策树的脑电信号特征选择方法

Also Published As

Publication number Publication date
CN104644166A (zh) 2015-05-27

Similar Documents

Publication Publication Date Title
CN104644166B (zh) 一种基于格子复杂性算法的无线动态麻醉深度检测方法
Almazaydeh et al. Obstructive sleep apnea detection using SVM-based classification of ECG signal features
Lu et al. A comparison of photoplethysmography and ECG recording to analyse heart rate variability in healthy subjects
CN104994782B (zh) 基于脑波信号进行分析的方法和设备
CN105147248A (zh) 基于生理信息的抑郁症评估系统及其评估方法
DE102008002933A1 (de) Datenaufzeichnung zur Patientenstatusanalyse
CN108392211A (zh) 一种基于多信息融合的疲劳检测方法
WO2013063053A1 (en) Method of building classifiers for real-time classification of neurological states
US11504056B2 (en) Method for classifying anesthetic depth in operations with total intravenous anesthesia
US20200327995A1 (en) Method for identifying anesthetic drug, and method and device for processing anesthesia electroencephalogram signal
CN112493995B (zh) 一种适用不同年龄患者的麻醉状态评估系统及方法
CN109497997A (zh) 基于多数据采集的癫痫发作检测设备及预警系统
CN112137616B (zh) 一种多感觉脑-体联合刺激的意识检测装置
WO2019127558A1 (zh) 基于脑电的麻醉深度监测方法和装置
Wang et al. Eeg-based real-time drowsiness detection using hilbert-huang transform
Tiwari et al. Detection of insomnia using advanced complexity and entropy features of sleep stage data of EEG recordings
CN113974557A (zh) 基于脑电奇异谱分析的深度神经网络麻醉深度分析方法
Sloboda et al. A simple sleep stage identification technique for incorporation in inexpensive electronic sleep screening devices
Dilber et al. EEG based detection of epilepsy by a mixed design approach
Park et al. A study on the development of a day-to-day mental stress monitoring system using personal physiological data
Raghavendra et al. Multiscale fractal dimension technique for characterization and analysis of biomedical signals
CN115349861A (zh) 一种基于单通道脑电信号的精神压力检测系统及方法
Benzy et al. Approximate entropy and wavelet entropy based depth of anesthesia monitoring
EP3831287A1 (en) Method and device for capacitive touch panel based biosignal measurement
CN114533088A (zh) 多模态大脑信号分类方法、装置、电子设备及存储介质

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20210301

Address after: No. 716, 7th floor, building 4, No. 99, Guangfu Road, Wuhou District, Chengdu, Sichuan 610041

Patentee after: Chengdu yishenrui Technology Co.,Ltd.

Address before: 310027 No. 38, Zhejiang Road, Hangzhou, Zhejiang, Xihu District

Patentee before: ZHEJIANG University

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20210308

Address after: 226000 east side of Yongtong Avenue, Rudong Economic Development Zone, Nantong City, Jiangsu Province

Patentee after: JIANGSU APON MEDICAL TECHNOLOGY Co.,Ltd.

Address before: No. 716, 7th floor, building 4, No. 99, Guangfu Road, Wuhou District, Chengdu, Sichuan 610041

Patentee before: Chengdu yishenrui Technology Co.,Ltd.