CN116491960B - Brain transient monitoring device, electronic device, and storage medium - Google Patents

Abstract

Translated fromChinese

本申请公开了一种脑瞬态监测设备、电子设备及存储介质。在当前时刻采集待监测对象的多个原始脑电信号；基于重组矩阵对多个原始脑电信号进行溯源分析，得到多个第一脑电源信号；从每个第一脑电源信号中截取与当前时刻对应的待分析脑电信号，得到多个待分析脑电信号；对多个待分析脑电信号进行处理，得到第一延时自协方差矩阵；确定第一延时自协方差矩阵与多个预设延时自协方差矩阵之间的多个第一相似度；根据多个第一相似度、上一时刻的脑瞬态、多种脑瞬态状态转移概率矩阵、多种脑瞬态的初始状态概率，确定当前时刻处于各种脑瞬态的概率；根据当前时刻处于各种脑瞬态的概率，确定当前时刻的脑瞬态。

This application discloses a brain transient monitoring device, electronic device and storage medium. Collect multiple original EEG signals of the object to be monitored at the current moment; perform traceability analysis on the multiple original EEG signals based on the reorganization matrix to obtain multiple first brain power signals; intercept and extract from each first brain power signal the current The corresponding EEG signals to be analyzed at each moment are used to obtain multiple EEG signals to be analyzed; the multiple EEG signals to be analyzed are processed to obtain the first delayed autocovariance matrix; the first delayed autocovariance matrix and the multiple EEG signals are determined. Multiple first similarities between preset delay autocovariance matrices; based on multiple first similarities, the brain transient at the previous moment, multiple brain transient state transition probability matrices, multiple brain transients The initial state probability is used to determine the probability of being in various brain transient states at the current moment; based on the probability of being in various brain transient states at the current moment, the brain transient state at the current moment is determined.

Description

Translated fromChinese

脑瞬态监测设备、电子设备及存储介质Brain transient monitoring equipment, electronic equipment and storage media

技术领域Technical field

本申请涉及医疗技术领域，具体涉及一种脑瞬态监测设备、电子设备及存储介质。This application relates to the field of medical technology, specifically to a brain transient monitoring device, electronic device and storage medium.

背景技术Background technique

人类大脑是一个复杂动态变化的非线性系统，由于大脑自发活动可以呈现为一些规律性的状态，并且这些状态与脑疾病的机制和发展密切关联，因此大脑活动状态的实时监测对于大脑疾病的评估和治疗具有极其重要的意义。The human brain is a complex and dynamically changing nonlinear system. Since spontaneous brain activity can present some regular states, and these states are closely related to the mechanism and development of brain diseases, real-time monitoring of brain activity states is important for the assessment of brain diseases. and treatment are extremely important.

目前临床上应用最为广泛的基于脑电的脑状态监测技术是针对癫痫的视频脑电技术以及麻醉监护技术。视频脑电技术使用视频录像与脑电结合，同步观察癫痫发作的行为特征与脑电特征，进而确定癫痫类型。但是该技术只限于数据采集，脑状态的分析需要经验丰富的神经内科医生进行离线状态下的肉眼判断，这种监测受人工经验的影响，导致监测精度较低，并且无法对脑状态进行实时性监测，也无法在线对脑状态进行监测。麻醉监护技术一般通过脑电的计算实时评估麻醉深度并指导麻醉医生调整麻醉药物。但是这种技术只针对前额脑电进行监测和计算，不对全脑状态进行评估，并且该技术只利用一段时间内的脑电数据进行脑状态分析，导致脑状态监测精度较低，而且无法对脑状态进行实时性监测。Currently, the most widely used brain state monitoring technology based on EEG in clinical practice is video EEG technology for epilepsy and anesthesia monitoring technology. Video EEG technology uses video recording and EEG to simultaneously observe the behavioral characteristics and EEG characteristics of epileptic seizures, thereby determining the type of epilepsy. However, this technology is limited to data collection. The analysis of brain states requires experienced neurologists to make offline visual judgments. This kind of monitoring is affected by manual experience, resulting in low monitoring accuracy and the inability to conduct real-time analysis of brain states. Monitoring, and brain status cannot be monitored online. Anesthesia monitoring technology generally evaluates the depth of anesthesia in real time through EEG calculations and guides anesthesiologists to adjust anesthetic drugs. However, this technology only monitors and calculates the forehead EEG and does not evaluate the whole brain state. Moreover, this technology only uses EEG data within a period of time to analyze the brain state, resulting in low accuracy of brain state monitoring and the inability to analyze the brain state. Status is monitored in real time.

因此，如何进行实时在线监测大脑的脑状态，以及提高大脑的脑状态的监测精度是目前亟待解决的技术问题。Therefore, how to monitor the brain's brain state online in real time and improve the monitoring accuracy of the brain's brain state are technical issues that need to be solved urgently.

发明内容Contents of the invention

本申请实施例提供了一种脑瞬态监测设备、电子设备及存储介质，通过离线处理的参数，并结合上一时刻的脑瞬态，提高了脑瞬态监测的精度，以及能够实时在线进行脑状态监测。Embodiments of the present application provide a brain transient monitoring device, an electronic device and a storage medium, which improve the accuracy of brain transient monitoring through offline processing of parameters and combined with the brain transient at the previous moment, and can be carried out online in real time. Brain state monitoring.

第一方面，本申请实施例提供一种脑瞬态监测设备，所述脑瞬态监测设备包括图像获取模块、离线处理模块、脑电采集模块和脑瞬态分析模块；In a first aspect, embodiments of the present application provide a brain transient monitoring device, which includes an image acquisition module, an offline processing module, an EEG acquisition module, and a brain transient analysis module;

所述图像获取模块，用于获取待监测对象的脑部的原始核磁共振图像；The image acquisition module is used to acquire original MRI images of the brain of the subject to be monitored;

所述脑电采集模块，用于在预设时间段内在多个通道上采集所述待监测对象的多个离线脑电信号；The EEG collection module is used to collect multiple offline EEG signals of the object to be monitored on multiple channels within a preset time period;

所述离线处理模块，用于基于所述原始核磁共振图像，得到所述待监测对象的三维脑模型；The offline processing module is used to obtain a three-dimensional brain model of the object to be monitored based on the original MRI image;

获取所述三维脑模型中的多个皮层网格；将所述多个通道上用于脑电信号采集的多个电极在所述待监测对象的脑部的位置与所述多个皮层网格在脑部的位置进行对齐，得到与所述待监测对象对应的Lead-field矩阵；Obtain multiple cortical grids in the three-dimensional brain model; compare the positions of the multiple electrodes used for brain electrical signal collection on the multiple channels in the brain of the subject to be monitored with the multiple cortical grids. Align the position of the brain to obtain a Lead-field matrix corresponding to the object to be monitored;

获取所述多个离线脑电信号的协方差矩阵；Obtain the covariance matrix of the multiple offline EEG signals;

基于所述协方差矩阵以及所述Lead-field矩阵，确定重组矩阵；Based on the covariance matrix and the lead-field matrix, determine the reorganization matrix;

基于所述Lead-field矩阵，对所述多个离线脑电信号进行逆分解，得到所述多个皮层网格对应的多个皮层信号；将所述三维脑模型划分为多个源极子；Based on the lead-field matrix, inversely decompose the multiple offline EEG signals to obtain multiple cortical signals corresponding to the multiple cortical grids; divide the three-dimensional brain model into multiple source sub-elements;

对所述多个皮层信号进行加权映射，得到与所述多个源极子对应的多个第二脑电源信号；Perform weighted mapping on the plurality of cortical signals to obtain a plurality of second brain power signals corresponding to the plurality of source sub-elements;

基于预设时间窗口，对每个第二脑电源信号进行多次加窗，以对每个第二脑电源信号进行分割，得到与每个第二脑电源信号对应的多个子脑电源信号；Based on the preset time window, perform multiple windowing on each second brain power signal to segment each second brain power signal and obtain multiple sub-brain power signals corresponding to each second brain power signal;

基于每个第二脑电源信号对应的多个子脑电源信号，得到多个预设延时自协方差矩阵、多种脑瞬态之间的状态转移概率矩阵以及所述多种脑瞬态对应的初始状态概率，其中，每个预设延时自协方差矩阵用于表征一种脑瞬态；Based on multiple sub-brain power signals corresponding to each second brain power signal, multiple preset delay auto-covariance matrices, state transition probability matrices between multiple brain transient states, and corresponding values of the multiple brain transient states are obtained. Initial state probability, where each preset delay autocovariance matrix is used to characterize a brain transient;

所述脑电采集模块，用于在当前时刻从多个通道上采集待监测对象的多个原始脑电信号；The EEG acquisition module is used to collect multiple original EEG signals of the object to be monitored from multiple channels at the current moment;

所述脑瞬态分析模块，用于基于所述重组矩阵对所述多个原始脑电信号进行溯源分析，得到多个第一脑电源信号；The brain transient analysis module is used to perform traceability analysis on the plurality of original EEG signals based on the reorganization matrix to obtain a plurality of first brain power signals;

从每个第一脑电源信号中截取与当前时刻对应的待分析脑电信号，得到多个待分析脑电信号；Intercept the EEG signal to be analyzed corresponding to the current moment from each first brain power signal to obtain multiple EEG signals to be analyzed;

对所述多个待分析脑电信号进行处理，得到第一延时自协方差矩阵；Process the plurality of EEG signals to be analyzed to obtain a first delayed autocovariance matrix;

确定所述第一延时自协方差矩阵与所述多个预设延时自协方差矩阵之间的相似度，得到多个第一相似度；Determine the similarity between the first delay autocovariance matrix and the plurality of preset delay autocovariance matrices, and obtain a plurality of first similarities;

根据所述多个第一相似度、所述待监测对象在上一时刻的脑瞬态、所述状态转移概率矩阵以及所述初始状态概率，确定所述待监测对象当前时刻处于各种脑瞬态的概率；According to the plurality of first similarities, the brain transient state of the object to be monitored at the previous moment, the state transition probability matrix and the initial state probability, it is determined that the object to be monitored is in various brain transient states at the current moment. probability of state;

根据当前时刻处于各种脑瞬态的概率，确定所述待监测对象在当前时刻的脑瞬态。According to the probability of being in various brain transient states at the current moment, the brain transient state of the object to be monitored at the current moment is determined.

第二方面，本申请实施例提供一种脑瞬态监测方法，所述方法应用于所述脑瞬态监测设备，所述脑瞬态监测设备包括图像获取模块、离线处理模块、脑电采集模块和脑瞬态分析模块；所述方法包括：In the second aspect, embodiments of the present application provide a brain transient monitoring method, which method is applied to the brain transient monitoring device. The brain transient monitoring device includes an image acquisition module, an offline processing module, and an EEG acquisition module. and a brain transient analysis module; the method includes:

获取待监测对象的脑部的原始核磁共振图像；Obtaining raw MRI images of the brain of the subject to be monitored;

在预设时间段内在多个通道上采集所述待监测对象的多个离线脑电信号；Collect multiple offline EEG signals of the object to be monitored on multiple channels within a preset time period;

基于所述原始核磁共振图像，得到所述待监测对象的三维脑模型；Based on the original MRI image, obtain a three-dimensional brain model of the object to be monitored;

获取所述三维脑模型中的多个皮层网格；将所述多个通道上用于脑电信号采集的多个电极在所述待监测对象的脑部的位置与所述多个皮层网格在脑部的位置进行对齐，得到与所述待监测对象对应的Lead-field矩阵；Obtain multiple cortical grids in the three-dimensional brain model; compare the positions of the multiple electrodes used for brain electrical signal collection on the multiple channels in the brain of the subject to be monitored with the multiple cortical grids. Align the position of the brain to obtain a Lead-field matrix corresponding to the object to be monitored;

获取所述多个离线脑电信号的协方差矩阵；Obtain the covariance matrix of the multiple offline EEG signals;

基于所述协方差矩阵以及所述Lead-field矩阵，确定重组矩阵；Based on the covariance matrix and the lead-field matrix, determine the reorganization matrix;

基于所述Lead-field矩阵，对所述多个离线脑电信号进行逆分解，得到所述多个皮层网格对应的多个皮层信号；将所述三维脑模型划分为多个源极子；Based on the lead-field matrix, inversely decompose the multiple offline EEG signals to obtain multiple cortical signals corresponding to the multiple cortical grids; divide the three-dimensional brain model into multiple source sub-elements;

对所述多个皮层信号进行加权映射，得到与所述多个源极子对应的多个第二脑电源信号；Perform weighted mapping on the plurality of cortical signals to obtain a plurality of second brain power signals corresponding to the plurality of source sub-elements;

基于预设时间窗口，对每个第二脑电源信号进行多次加窗，以对每个第二脑电源信号进行分割，得到与每个第二脑电源信号对应的多个子脑电源信号；Based on the preset time window, perform multiple windowing on each second brain power signal to segment each second brain power signal and obtain multiple sub-brain power signals corresponding to each second brain power signal;

基于每个第二脑电源信号对应的多个子脑电源信号，得到多个预设延时自协方差矩阵、多种脑瞬态之间的状态转移概率矩阵以及所述多种脑瞬态对应的初始状态概率，其中，每个预设延时自协方差矩阵用于表征一种脑瞬态；Based on multiple sub-brain power signals corresponding to each second brain power signal, multiple preset delay auto-covariance matrices, state transition probability matrices between multiple brain transient states, and corresponding values of the multiple brain transient states are obtained. Initial state probability, where each preset delay autocovariance matrix is used to characterize a brain transient;

在当前时刻从多个通道上采集待监测对象的多个原始脑电信号；Collect multiple raw EEG signals of the object to be monitored from multiple channels at the current moment;

基于所述重组矩阵对所述多个原始脑电信号进行溯源分析，得到多个第一脑电源信号；Perform traceability analysis on the plurality of original EEG signals based on the recombination matrix to obtain a plurality of first brain power signals;

从每个第一脑电源信号中截取与当前时刻对应的待分析脑电信号，得到多个待分析脑电信号；Intercept the EEG signal to be analyzed corresponding to the current moment from each first brain power signal to obtain multiple EEG signals to be analyzed;

对所述多个待分析脑电信号进行处理，得到第一延时自协方差矩阵；Process the plurality of EEG signals to be analyzed to obtain a first delayed autocovariance matrix;

确定所述第一延时自协方差矩阵与所述多个预设延时自协方差矩阵之间的相似度，得到多个第一相似度；Determine the similarity between the first delay autocovariance matrix and the plurality of preset delay autocovariance matrices, and obtain a plurality of first similarities;

根据所述多个第一相似度、所述待监测对象在上一时刻的脑瞬态、所述状态转移概率矩阵以及所述初始状态概率，确定所述待监测对象当前时刻处于各种脑瞬态的概率；According to the plurality of first similarities, the brain transient state of the object to be monitored at the previous moment, the state transition probability matrix and the initial state probability, it is determined that the object to be monitored is in various brain transient states at the current moment. probability of state;

根据当前时刻处于各种脑瞬态的概率，确定所述待监测对象在当前时刻的脑瞬态。According to the probability of being in various brain transient states at the current moment, the brain transient state of the object to be monitored at the current moment is determined.

第三方面，本申请提供一种电子设备，包括：处理器和存储器，所述处理器与存储器相连，所述存储器用于存储计算机程序，所述处理器用于执行所述存储器中存储的计算机程序，以使得所述电子设备执行如第二方面所述的方法。In a third aspect, the present application provides an electronic device, including: a processor and a memory. The processor is connected to the memory. The memory is used to store a computer program. The processor is used to execute the computer program stored in the memory. , so that the electronic device performs the method described in the second aspect.

第四方面，本申请提供一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序使得计算机执行如第二方面所述的方法。In a fourth aspect, the present application provides a computer-readable storage medium that stores a computer program, and the computer program causes a computer to execute the method as described in the second aspect.

第五方面，本申请提供一种计算机程序产品，所述计算机程序产品包括存储了计算机程序的非瞬时性计算机可读存储介质，所述计算机可操作来使计算机执行如第二方面所述的方法。In a fifth aspect, the present application provides a computer program product. The computer program product includes a non-transitory computer-readable storage medium storing a computer program, and the computer is operable to cause the computer to perform the method as described in the second aspect. .

实施本申请实施例，具有如下有益效果：Implementing the embodiments of this application has the following beneficial effects:

可以看出，在本申请实施例中，首先离线计算出脑电信号到脑电源信号的重组矩阵，以及用于表征脑瞬态的预设延时自协方差矩阵，以及多种脑瞬态之间的状态转移概率矩阵，以及多种脑瞬态的初始状态概率，这样进行实时的在线监测时，可以利用离线计算好的重组矩阵，直接对在线采集到的多个原始脑电信号进行溯源，得到多个第一脑电源信号；以及从每个第一脑电源信号中截取与当前时刻对应的待分析脑电信号，得到多个待分析脑电信号，并获取多个待分析脑电信号对应的第一延时自协方差矩阵。然后，直接使用多个预设延时自协方差矩阵，确定出多个第一相似度。最后，直接基于多个第一相似度、待监测对象在上一时刻的脑瞬态、多种脑瞬态之间的状态转移概率矩阵、多种脑瞬态对应的初始状态概率，确定出待监测对象当前时刻处于各种脑瞬态的概率，进而确定出当前时刻的脑状态。由于在线监测的大部分复杂参数都是离线计算好的，大幅度的减轻了在线计算的复杂度，因此在线监测时可以快速地计算出当前时刻的脑状态，从而实现实时在线地监测大脑的脑状态。并且在确定当前时刻的脑瞬态时，不会只利用当前时刻下的脑电信号，与各种脑瞬态之间的相似度，还会结合上一个时刻的脑瞬态，以及多种脑瞬态之间的转移概率，从而可以精确地确定出当前时刻的脑瞬态，提高了对脑瞬态的监测精度。It can be seen that in the embodiment of the present application, the recombination matrix from the brain electrical signal to the brain power signal is first calculated offline, as well as the preset delay autocovariance matrix used to characterize brain transients, and a variety of brain transients. The state transition probability matrix between states, as well as the initial state probabilities of various brain transients, so that when performing real-time online monitoring, the offline calculated reorganization matrix can be used to directly trace the source of multiple original EEG signals collected online. Obtain a plurality of first brain power signals; and intercept an EEG signal to be analyzed corresponding to the current moment from each first brain power signal, obtain a plurality of EEG signals to be analyzed, and obtain the corresponding EEG signals of the multiple EEG signals to be analyzed. The first delay autocovariance matrix of . Then, multiple preset delay autocovariance matrices are directly used to determine multiple first similarities. Finally, directly based on multiple first similarities, the brain transient of the object to be monitored at the previous moment, the state transition probability matrix between multiple brain transients, and the initial state probabilities corresponding to multiple brain transients, the target to be monitored is determined. Monitor the probability that the subject is in various brain transient states at the current moment, and then determine the brain state at the current moment. Since most of the complex parameters for online monitoring are calculated offline, the complexity of online calculations is greatly reduced. Therefore, the brain state at the current moment can be quickly calculated during online monitoring, thereby achieving real-time online monitoring of the brain. state. And when determining the brain transient at the current moment, we will not only use the similarity between the EEG signal at the current moment and various brain transients, but also combine the brain transient at the previous moment and various brain transients. The transition probability between transient states can accurately determine the brain transient at the current moment, improving the accuracy of monitoring brain transients.

附图说明Description of drawings

为了更清楚地说明本申请实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。In order to explain the technical solutions in the embodiments of the present application more clearly, the drawings needed to be used in the description of the embodiments will be briefly introduced below. Obviously, the drawings in the following description are some embodiments of the present application. Those of ordinary skill in the art can also obtain other drawings based on these drawings without exerting creative efforts.

图1为本申请实施例提供的一种脑瞬态监测的场景示意图；Figure 1 is a schematic diagram of a brain transient monitoring scene provided by an embodiment of the present application;

图2为本申请实施例提供的一种脑瞬态监测设备的示意图；Figure 2 is a schematic diagram of a brain transient monitoring device provided by an embodiment of the present application;

图3为本申请实施例提供的一种基于预设时间窗口获取每个第二脑电源信号对应的多个子脑电源信号的示意图；Figure 3 is a schematic diagram of obtaining multiple sub-brain power signals corresponding to each second brain power signal based on a preset time window according to an embodiment of the present application;

图4为本申请实施例提供的一种获取每个时间窗口的第二延时自协方差矩阵的示意图；Figure 4 is a schematic diagram of obtaining the second delay autocovariance matrix of each time window provided by the embodiment of the present application;

图5为本申请实施例提供的一种脑瞬态监测的示意图；Figure 5 is a schematic diagram of brain transient monitoring provided by an embodiment of the present application;

图6为本申请实施例提供的一种脑瞬态监测方法的流程示意图；Figure 6 is a schematic flow chart of a brain transient monitoring method provided by an embodiment of the present application;

图7为本申请实施例提供的一种电子设备的示意图。FIG. 7 is a schematic diagram of an electronic device provided by an embodiment of the present application.

具体实施方式Detailed ways

下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。The technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present application. Obviously, the described embodiments are part of the embodiments of the present application, rather than all of the embodiments. Based on the embodiments in this application, all other embodiments obtained by those of ordinary skill in the art without creative efforts fall within the scope of protection of this application.

本申请的说明书和权利要求书及所述附图中的术语“第一”、“第二”、“第三”和“第四”等是用于区别不同对象，而不是用于描述特定顺序。此外，术语“包括”和“具有”以及它们任何变形，意图在于覆盖不排他的包含。例如包含了一系列步骤或单元的过程、方法、系统、产品或设备没有限定于已列出的步骤或单元，而是可选地还包括没有列出的步骤或单元，或可选地还包括对于这些过程、方法、产品或设备固有的其它步骤或单元。The terms “first”, “second”, “third” and “fourth” in the description, claims and drawings of this application are used to distinguish different objects, rather than to describe a specific sequence. . Furthermore, the terms "including" and "having" and any variations thereof are intended to cover non-exclusive inclusion. For example, a process, method, system, product or device that includes a series of steps or units is not limited to the listed steps or units, but optionally also includes steps or units that are not listed, or optionally also includes Other steps or units inherent to such processes, methods, products or devices.

在本文中提及“实施例”意味着，结合实施例描述的特定特征、结果或特性可以包含在本申请的至少一个实施例中。在说明书中的各个位置出现该短语并不一定均是指相同的实施例，也不是与其它实施例互斥的独立的或备选的实施例。本领域技术人员显式地和隐式地理解的是，本文所描述的实施例可以与其它实施例相结合。Reference herein to "an embodiment" means that a particular feature, result or characteristic described in connection with the embodiment may be included in at least one embodiment of the present application. The appearances of this phrase in various places in the specification are not necessarily all referring to the same embodiment, nor are separate or alternative embodiments mutually exclusive of other embodiments. Those skilled in the art understand, both explicitly and implicitly, that the embodiments described herein may be combined with other embodiments.

参阅图1，图1为本申请实施例提供的一种脑瞬态监测的场景示意图。Refer to Figure 1, which is a schematic diagram of a brain transient monitoring scenario provided by an embodiment of the present application.

如图1所示，脑瞬态监测设备通过放置在待监测对象的脑皮层的电极（图1中脑皮层的圆圈）在当前时刻从多个通道上采集待监测对象的多个原始脑电信号；然后，基于离线处理得到的重组矩阵对所述多个原始脑电信号进行溯源分析，得到多个第一脑电源信号。然后，从每个第一脑电源信号中截取与当前时刻对应的待分析脑电信号，得到多个待分析脑电信号；对所述多个待分析脑电信号进行处理，得到第一延时自协方差矩阵；确定所述第一延时自协方差矩阵与多个预设延时自协方差矩阵之间的相似度，得到多个第一相似度，其中，每个预设延时自协方差矩阵用于表征一种脑瞬态；根据所述多个第一相似度、所述待监测对象在上一时刻的脑瞬态、多种脑瞬态之间的状态转移概率矩阵、所述多种脑瞬态对应的初始状态概率，确定所述待监测对象当前时刻处于各种脑瞬态的概率；根据当前时刻处于各种脑瞬态的概率，确定所述待监测对象在当前时刻的脑瞬态。最后，可以在脑瞬态监测设备显示界面显示待监测对象在当前时刻的脑瞬态。As shown in Figure 1, the brain transient monitoring equipment collects multiple raw EEG signals of the subject to be monitored from multiple channels at the current moment through electrodes placed on the cerebral cortex of the subject to be monitored (the circle in the cerebral cortex in Figure 1). ; Then, perform traceability analysis on the multiple original EEG signals based on the reorganization matrix obtained through offline processing, and obtain multiple first EEG signals. Then, intercept the EEG signal to be analyzed corresponding to the current moment from each first brain power signal to obtain multiple EEG signals to be analyzed; process the multiple EEG signals to be analyzed to obtain the first delay Autocovariance matrix; determine the similarity between the first delay autocovariance matrix and multiple preset delay autocovariance matrices to obtain multiple first similarities, wherein each preset delay autocovariance matrix The covariance matrix is used to represent a brain transient; according to the multiple first similarities, the brain transient of the object to be monitored at the previous moment, and the state transition probability matrix between multiple brain transients, the According to the initial state probabilities corresponding to the various brain transient states, the probability that the subject to be monitored is in various brain transient states at the current moment is determined; based on the probability that the subject to be monitored is in various brain transient states at the current moment, the probability that the subject to be monitored is in the various brain transient states at the current moment is determined. of brain transients. Finally, the brain transient state of the subject to be monitored at the current moment can be displayed on the display interface of the brain transient monitoring device.

应说明，本申请所提到的脑瞬态也可以称为脑状态，两者在本质上是一致的，可以不用区分。It should be noted that the brain transients mentioned in this application can also be called brain states. The two are essentially the same and do not need to be distinguished.

参阅图2，图2为本申请实施例提供的一种脑瞬态监测设备的示意图。如图2所示，脑瞬态监测设备包括：脑电采集模块、脑瞬态分析模块、图像获取模块和离线处理模块。Refer to Figure 2, which is a schematic diagram of a brain transient monitoring device provided by an embodiment of the present application. As shown in Figure 2, the brain transient monitoring equipment includes: EEG acquisition module, brain transient analysis module, image acquisition module and offline processing module.

为了便于理解本申请的技术方案，首先对本申请的离线处理过程进行解释和说明。In order to facilitate understanding of the technical solution of the present application, the offline processing process of the present application is first explained and described.

本申请的离线处理过程，主要是为了准备在线脑瞬态监测所需的参数。其中，在线脑瞬态所需的参数至少包括：重组矩阵、每种脑瞬态对应的预设延时自协方差矩阵、多种脑瞬态之间的状态转移概率矩阵以及多种脑瞬态对应的初始状态概率。The offline processing process of this application is mainly to prepare the parameters required for online brain transient monitoring. Among them, the parameters required for online brain transients include at least: a reorganization matrix, a preset delay autocovariance matrix corresponding to each brain transient, a state transition probability matrix between multiple brain transients, and multiple brain transients. The corresponding initial state probability.

下面结合附图说明本申请的离线处理过程。The offline processing process of this application will be described below with reference to the accompanying drawings.

示例性的，图像获取模块用于获取待监测对象的脑部的原始核磁共振图像（Magnetic Resonance Imaging，MRI）。For example, the image acquisition module is used to acquire the original magnetic resonance image (Magnetic Resonance Imaging, MRI) of the brain of the subject to be monitored.

脑电采集模块，用于在预设时间段内在多个通道上采集待监测对象的多个离线脑电信号；该预设时间段可以是对待监测对象进行脑瞬态在线监测之前的一个时间段。例如，在待监测对象的脑部安装好电极（即采集脑电信号的脑电电极），且t时刻准备对待监测对象进行脑瞬态监测，则可以在位于t时刻之前的一个预设时间段内采集待监测对象的多个脑电信号，作为多个离线脑电信号进行离线分析。当然，预设时间段可以为进行脑瞬态监测之前的任意一个时间段，比如，在t时刻之前的任意一个预设时间段，本申请不对预设时间段进行具体限定。然后，离线处理模块用于基于核磁共振图像和多个离线脑电信号，确定多个预设延时自协方差矩阵、重组矩阵、状态转移概率矩阵以及初始状态概率。The EEG acquisition module is used to collect multiple offline EEG signals of the subject to be monitored on multiple channels within a preset time period; the preset time period can be a time period before the subject to be monitored undergoes online brain transient monitoring. . For example, if electrodes (i.e., EEG electrodes that collect EEG signals) are installed on the brain of the subject to be monitored, and the subject is prepared to perform brain transient monitoring at time t, then the subject can be monitored in a preset time period before time t. Multiple EEG signals of the object to be monitored are collected within the system and analyzed offline as multiple offline EEG signals. Of course, the preset time period can be any time period before brain transient monitoring is performed, for example, any preset time period before time t. This application does not specifically limit the preset time period. Then, the offline processing module is used to determine multiple preset time-delay autocovariance matrices, reorganization matrices, state transition probability matrices, and initial state probabilities based on the MRI image and multiple offline EEG signals.

具体地，离线处理模块先基于原始核磁共振图像，得到所述待监测对象的三维脑模型。具体地，离线处理模块对原始核磁共振图像进行重切片，得到目标核磁共振图像。然后，离线处理模块对目标核磁共振图像进行分割，得到脑组织、颅骨以及头皮，并基于所述分割出脑组织、颅骨以及头皮的所述目标核磁共振图像进行脑模型建立，即基于该目标核磁共振图像，以及使用各向一致性的传导参数建立信号传导正向脑模型，得到待监测对象的三维脑模型。然后，离线处理模块获取三维脑模型中的所有皮层网格，并对所有皮层网格进行降采样（即减少皮层网格的数量），并在降采样后将皮层网络与标准脑模型的皮层网格进行匹配校准，得到多个皮层网格。应说明，随着算力的不断提升，若离线处理模块能够处理的皮层网格的数量不受限制，则可以不用对皮层网格进行降采样，仅需对所有皮层网格进行校准即可。本申请中主要以降采样以及校准后得到的多个皮层网格为例进行说明。Specifically, the offline processing module first obtains the three-dimensional brain model of the object to be monitored based on the original MRI image. Specifically, the offline processing module reslices the original MRI image to obtain the target MRI image. Then, the offline processing module segments the target MRI image to obtain the brain tissue, skull and scalp, and establishes a brain model based on the target MRI image segmented into the brain tissue, skull and scalp, that is, based on the target MRI Resonance images, and the use of isotropically consistent conduction parameters to establish signal conduction forward brain models, to obtain a three-dimensional brain model of the object to be monitored. Then, the offline processing module obtains all cortical grids in the three-dimensional brain model, downsamples all cortical grids (i.e., reduces the number of cortical grids), and compares the cortical network with that of the standard brain model after downsampling. The grids are matched and calibrated to obtain multiple cortical grids. It should be noted that with the continuous improvement of computing power, if the number of cortical grids that the offline processing module can process is not limited, there is no need to downsample the cortical grids and only need to calibrate all cortical grids. This application mainly takes multiple cortical grids obtained after downsampling and calibration as an example for explanation.

然后，离线处理模块获取多个通道上用于脑电信号采集的多个电极在待监测对象的脑部的位置，以及多个皮层网格在待监测对象的脑部的位置，并将这两个位置进行对齐，得到与待监测对象对应的Lead-field矩阵（也可称为导联场矩阵）。最后，离线处理模块获取该多个离线脑电信号的协方差矩阵。例如，将每个离线脑电信号在各个时刻的幅值作为矩阵中的一行或一列元素，可得到第一幅值矩阵。基于该第一幅值矩阵，得到该多个离线脑电信号的协方差矩阵。最后，离线处理模块基于该Lead-field矩阵和该协方差矩阵，得到重组矩阵。Then, the offline processing module obtains the positions of the multiple electrodes used for EEG signal collection on multiple channels in the brain of the subject to be monitored, and the positions of multiple cortical grids in the brain of the subject to be monitored, and combines these two Align the positions to obtain the Lead-field matrix (also called the lead field matrix) corresponding to the object to be monitored. Finally, the offline processing module obtains the covariance matrices of the multiple offline EEG signals. For example, the first amplitude matrix can be obtained by taking the amplitude of each offline EEG signal at each moment as a row or column element in the matrix. Based on the first amplitude matrix, covariance matrices of the multiple offline EEG signals are obtained. Finally, the offline processing module obtains the reorganization matrix based on the lead-field matrix and the covariance matrix.

示例性的，重组矩阵可以通过公式（1）表示：Illustratively, the reorganization matrix can be expressed by formula (1):

公式（1）； Formula 1);

其中，W为重组矩阵，R上述协方差矩阵，L为上述Lead-field矩阵，T为转置操作。Among them, W is the reorganization matrix, R is the above-mentioned covariance matrix, L is the above-mentioned Lead-field matrix, and T is the transpose operation.

应说明，重组矩阵主要用于对脑电信号进行逆向分解，即将脑电信号分解到多个源极子对应的多个脑电源信号。由于对于一个个体来说，脑电信号的逆向分解主要与所使用的正向脑模型、电极位置、Lead-field矩阵和脑电信号的噪声水平相关。因此针对同一个个体，在相同的脑电环境（即电极位置相同以及使用的脑模型相同）下，则每次采集到的脑电信号逆向分解所使用的重组矩阵是相同的，也就是将每次采集到的脑电信号乘以该重组矩阵即可完成从脑电信号到脑电源信号的逆向分解。It should be noted that the recombination matrix is mainly used to reversely decompose the EEG signal, that is, to decompose the EEG signal into multiple brain power signals corresponding to multiple source sub-elements. Because for an individual, the inverse decomposition of the EEG signal is mainly related to the forward brain model used, the electrode position, the lead-field matrix and the noise level of the EEG signal. Therefore, for the same individual, in the same EEG environment (that is, the electrode positions are the same and the brain model used is the same), the reorganization matrix used for reverse decomposition of the EEG signals collected each time is the same, that is, each EEG signal is The collected EEG signals are multiplied by the reorganization matrix to complete the reverse decomposition from EEG signals to brain power signals.

应说明，本申请中在线脑瞬态监测和离线处理过程，对待监测对象所使用的脑电环境是相同，即所使用的脑模型、电极位置都是相同的。这样也就保证本申请后续在线进行脑瞬态监测时，可以使用离线分析出的重组矩阵直接进行脑电信号到脑电源信号的逆分解，无需再次计算重组矩阵。It should be noted that in the online brain transient monitoring and offline processing processes in this application, the EEG environment used for the object to be monitored is the same, that is, the brain model and electrode positions used are the same. This also ensures that when this application performs subsequent online brain transient monitoring, the reorganization matrix analyzed offline can be used to directly perform the inverse decomposition of the EEG signal into the brain power signal, without the need to calculate the reorganization matrix again.

进一步地，在得到Lead-field矩阵后，离线处理模块还基于该Lead-field矩阵，对多个离线脑电信号进行逆分解，得到多个皮层网格对应的多个皮层信号，即使用该Lead-field矩阵，通过线性约束最小方差算法对该多个离线脑电信号进行逆分解，得到多个皮层信号。Further, after obtaining the Lead-field matrix, the offline processing module also performs inverse decomposition of multiple offline EEG signals based on the Lead-field matrix to obtain multiple cortical signals corresponding to multiple cortical grids, that is, using the Lead-field matrix -field matrix, use the linear constrained minimum variance algorithm to inversely decompose the multiple offline EEG signals to obtain multiple cortical signals.

然后，离线处理模块将三维脑模型划分为多个源极子，即将脑模型中的一块皮层区域划分为一个源极子，可得到多个源极子。然后，对多个皮层信号进行加权映射，得到与多个源极子对应的多个第二脑电源信号。可选地，可以使用主成分分析方法，将多个皮层信号进行加权映射，得到多个第二脑电信号。可选地，在对多个皮层信号进行加权映射后，还需要对加权映射得到的脑电信号的空间泄露和源极子极性进行校正，将校正后的脑电信号作为多个第二脑电信号。Then, the offline processing module divides the three-dimensional brain model into multiple source sub-units, that is, divides a cortical area in the brain model into one source sub-unit, and multiple source sub-units can be obtained. Then, weighted mapping is performed on the plurality of cortical signals to obtain a plurality of second brain power signals corresponding to the plurality of source sub-elements. Alternatively, the principal component analysis method can be used to perform weighted mapping on multiple cortical signals to obtain multiple second EEG signals. Optionally, after performing weighted mapping on multiple cortical signals, it is also necessary to correct the spatial leakage and source sub-polarity of the EEG signals obtained by the weighted mapping, and use the corrected EEG signals as multiple second brain signals. electric signal.

最后，离线处理模块基于预设时间窗口和所述多个第二脑电源信号，得到多个预设延时自协方差矩阵、多种脑瞬态之间的状态转移概率矩阵以及所述多种脑瞬态对应的初始状态概率，其中，每个预设延时自协方差矩阵用于表征一种脑瞬态。Finally, based on the preset time window and the plurality of second brain power signals, the offline processing module obtains a plurality of preset delay autocovariance matrices, a state transition probability matrix between various brain transient states, and the plurality of second brain power signals. The initial state probability corresponding to the brain transient, where each preset delay autocovariance matrix is used to represent a brain transient.

具体地，离线处理模块基于预设时间窗口，对每个第二脑电源信号进行多次加窗，以对每个第二脑电源信号进行分割，得到与每个第二脑电源信号对应的多个子脑电源信号。如图3所示，以预设时间窗口在每个第二脑电源信号上进行滑动，这样在滑动过程中就会产生多个时间窗口，并且每个时间窗口会从第二脑电源信号上分割出一个子脑电源信号。基于每个第二脑电源信号对应的多个子脑电源信号，确定多个预设延时自协方差矩阵、状态转移概率矩阵以及初始状态概率。Specifically, the offline processing module performs multiple windowing on each second brain power signal based on a preset time window to segment each second brain power signal and obtain multiple brain power signals corresponding to each second brain power signal. Subbrain power signal. As shown in Figure 3, each second brain power signal is slid with a preset time window, so that multiple time windows will be generated during the sliding process, and each time window will be segmented from the second brain power signal. Outputs a sub-brain power signal. Based on multiple sub-brain power signals corresponding to each second brain power signal, a plurality of preset delay autocovariance matrices, state transition probability matrices and initial state probabilities are determined.

具体地，离线处理模块获取多个第二脑电源信号在每个时间窗口下的多个子脑电源信号对应的第二延时自协方差矩阵，即获取该多个子脑电源信号的自协方差矩阵，其中，由于该多个子脑电源信号是从不同脑皮层区域的采集到的，并且该多个子脑电源信号是在一个时间窗口下的脑电信号，因此，该第二延时自协方差矩阵同时包含了空间维度和时间维度上的脑节律信息。举例来说，如图4所示，针对第一个时间窗口，获取每个第二脑电源信号在该第一个时间窗口下的子脑电源信号，可得到在第一个时间窗口下的多个子脑电源信号。然后，将每个子脑电源信号在各个时刻的幅值作为矩阵中的一行或者一列的元素，可得到与第一个时间窗口对应的第二幅值矩阵。然后，基于该第二幅值矩阵，得到与第一个时间窗口对应的第二延时自协方差矩阵。最后，针对每个时间窗口，可获取与每个时间窗口下的第二延时自协方差矩阵，进而得到多次加窗对应的多个时间窗口下的多个第二延时自协方差矩阵。Specifically, the offline processing module obtains the second delay auto-covariance matrix corresponding to the multiple sub-brain power signals of the multiple second brain power signals in each time window, that is, obtains the auto-covariance matrix of the multiple sub-brain power signals. , wherein, since the multiple sub-brain power signals are collected from different cerebral cortex areas, and the multiple sub-brain power signals are EEG signals in a time window, therefore, the second delay auto-covariance matrix It also contains brain rhythm information in both spatial and temporal dimensions. For example, as shown in Figure 4, for the first time window, obtain the sub-brain power signal of each second brain power signal under the first time window, and obtain multiple sub-brain power signals under the first time window. Subbrain power signal. Then, using the amplitude of each sub-brain power signal at each moment as an element in a row or column of the matrix, a second amplitude matrix corresponding to the first time window can be obtained. Then, based on the second amplitude matrix, a second delay autocovariance matrix corresponding to the first time window is obtained. Finally, for each time window, the second delay autocovariance matrix under each time window can be obtained, and then multiple second delay autocovariance matrices under multiple time windows corresponding to multiple windowings can be obtained .

进一步地，离线处理模块对多个第二延时自协方差矩阵进行分组，得到多个延时自协方差矩阵组，其中，每个延时自协方差矩阵组包括该多个第二延时自协方差矩阵中的一个或多个。例如，离线处理模块可以对多个第二延时自协方差矩阵进行聚类，得到多个延时自协方差矩阵组。然后，获取每个延时自协方差矩阵组的中心，将每个延时自协方差矩阵组的中心作为一个预设延时自协方差矩阵组，进而得到多个预设延时自协方差矩阵组。例如，将每个延时自协方差矩阵组对应的聚类中心作为该延时自协方差矩阵组的中心，或者，获取每个延时自协方差矩阵组中的第二延时自协方差矩阵的平均值，将该平均值作为每个延时自协方差矩阵组的中心。然后，离线处理模块使用每个预设延时自协方差矩阵表征一种脑瞬态。在本申请中，也可以将每个预设延时自协方差矩阵称作一个高斯观测模型，故每种脑瞬态也可以通过一个高斯观测模型表征，则多种脑瞬态可以通过多个高斯观测模型，即高斯观测模型群进行表征。在本申请的一个实施方式中，在得到每种脑瞬态对应的预设延时自协方差矩阵后，为了使后续的在线脑瞬态监测泛化，可以使用主成分分析对预设延时自协方差矩阵降维，仅保留预设延时自协方差矩阵中的关键时间信息和关键空间信息，得到降维后的延时自协方差矩阵，然后可以将降维后的延时自协方差矩阵作为脑瞬态对应的预设延时自协方差矩阵。本申请中主要以不对预设延时自协方差矩阵进行降维为例进行说明。Further, the offline processing module groups a plurality of second delay autocovariance matrices to obtain a plurality of delay autocovariance matrix groups, wherein each delay autocovariance matrix group includes the plurality of second delay autocovariance matrices. One or more of the autocovariance matrices. For example, the offline processing module can cluster multiple second delay autocovariance matrices to obtain multiple delay autocovariance matrix groups. Then, obtain the center of each delay autocovariance matrix group, and use the center of each delay autocovariance matrix group as a preset delay autocovariance matrix group to obtain multiple preset delay autocovariances. matrix group. For example, use the cluster center corresponding to each delay autocovariance matrix group as the center of the delay autocovariance matrix group, or obtain the second delay autocovariance in each delay autocovariance matrix group. The average of the matrices, which is used as the center of each delay autocovariance matrix group. The offline processing module then uses each preset delay autocovariance matrix to characterize a brain transient. In this application, each preset delay autocovariance matrix can also be called a Gaussian observation model, so each brain transient can also be characterized by a Gaussian observation model, and multiple brain transients can be represented by multiple Gaussian observation model, that is, Gaussian observation model group is represented. In one embodiment of the present application, after obtaining the preset delay autocovariance matrix corresponding to each brain transient, in order to generalize subsequent online brain transient monitoring, principal component analysis can be used to analyze the preset delay The dimensionality of the autocovariance matrix is reduced, and only the key time information and key spatial information in the preset delay autocovariance matrix are retained, and the dimensionally reduced delay autocovariance matrix is obtained. The dimensionally reduced delay autocovariance matrix can then be The variance matrix serves as a preset delay autocovariance matrix corresponding to brain transients. This application mainly explains by taking the example of not reducing the dimensionality of the preset delay autocovariance matrix.

示例性的，本申请中的多种脑瞬态包括但不限于：前额瞬态、感觉运动瞬态、顶叶瞬态和视觉瞬态。进一步地，在得到多个预设延时自协方差矩阵后，为了将每种脑瞬态进行具化表示，可以将每个预设延时自协方差矩阵表征的脑瞬态进行编号，比如，可以分别编号为：脑瞬态1、脑瞬态2、……。为了方便解释，本申请中主要以分组出4种脑瞬态为例进行说明，并且这四种脑瞬态分别编号为脑瞬态1、脑瞬态2、脑瞬态3以及脑瞬态4。By way of example, various brain transients in this application include, but are not limited to: frontal transients, sensorimotor transients, parietal lobe transients and visual transients. Furthermore, after obtaining multiple preset delay autocovariance matrices, in order to specifically represent each brain transient, the brain transients represented by each preset delay autocovariance matrix can be numbered, such as , which can be numbered respectively: brain transient 1, brain transient 2, .... For the convenience of explanation, this application mainly takes the example of grouping four brain transient states, and these four brain transient states are numbered as brain transient 1, brain transient 2, brain transient 3 and brain transient 4 respectively. .

进一步地，离线处理模块获取每个时间窗口下的第二延时自协方差矩阵，与多个预设延时自协方差矩阵之间的多个第二相似度。然后，离线处理模块根据每个时间窗口下的多个第二相似度，确定所述待监测对象在每个时间窗口下的脑瞬态，即将第二相似度最大所对应的预设延时自协方差矩阵表征的脑瞬态作为每个时间窗口下的脑瞬态。然后，离线处理模块按照时间的先后顺序，对多个时间窗口下的多个脑瞬态进行排列，得到脑瞬态序列，即将多个时间窗口下的多个脑瞬态按照时间顺序组合起来，得到该脑瞬态序列。最后，基于脑瞬态序列，确定多种脑瞬态之间的状态转移概率矩阵以及多种脑瞬态对应的初始状态概率。Further, the offline processing module obtains a plurality of second similarities between the second delay autocovariance matrix in each time window and a plurality of preset delay autocovariance matrices. Then, the offline processing module determines the brain transient state of the object to be monitored in each time window based on the plurality of second similarities in each time window, that is, the preset delay corresponding to the maximum second similarity is automatically The covariance matrix represents the brain transient as the brain transient under each time window. Then, the offline processing module arranges multiple brain transients in multiple time windows in order of time to obtain a brain transient sequence, that is, combines multiple brain transients in multiple time windows in chronological order. Obtain this brain transient sequence. Finally, based on the brain transient sequence, the state transition probability matrix between multiple brain transient states and the initial state probabilities corresponding to multiple brain transient states are determined.

具体地，离线处理模块确定该脑瞬态序列中每种脑瞬态出现的次数，基于每种脑瞬态出现的次数，以及所述脑瞬态序列中脑瞬态的数量（也就是上述多个脑瞬态的数量或者可以理解为时间窗口的数量）确定每种脑瞬态的出现概率，即将每种脑瞬态出现的次数与脑瞬态序列中脑瞬态的数量之间的比值作为每种脑瞬态的出现概率。最后，将每种脑瞬态的出现概率，作为多种脑瞬态对应的初始状态概率。Specifically, the offline processing module determines the number of occurrences of each brain transient in the brain transient sequence, based on the number of occurrences of each brain transient and the number of brain transients in the brain transient sequence (that is, the above-mentioned number of brain transients). The number of individual brain transients (or can be understood as the number of time windows) determines the occurrence probability of each brain transient, that is, the ratio between the number of occurrences of each brain transient and the number of brain transients in the brain transient sequence is Probability of occurrence of each brain transient. Finally, the occurrence probability of each brain transient is used as the initial state probability corresponding to multiple brain transients.

具体地，针对每种脑瞬态，离线处理模块确定所述脑瞬态序列中下一时刻与该脑瞬态相邻的脑瞬态，基于下一时刻与该脑瞬态相邻的脑瞬态，确每种脑瞬态下一时刻转移为各种脑瞬态的次数；基于每种脑瞬态下一时刻转移为各种脑瞬态的次数，以及每种脑瞬态出现的次数，确定每种脑瞬态转移为各种脑瞬态的概率。例如，将每种脑瞬态下一时刻转移为各种脑瞬态的次数与每种脑瞬态出现的次数之间的比值，作为每种脑瞬态转移为各种脑瞬态的概率；基于每种脑瞬态转移为各种脑瞬态的概率，确定多种脑瞬态之间的状态转移概率矩阵，即将每种脑瞬态转移为各种脑瞬态的概率分别作为状态转移概率矩阵中的每行元素，得到状态转移概率矩阵。Specifically, for each brain transient, the offline processing module determines the brain transient adjacent to the brain transient at the next moment in the brain transient sequence, based on the brain transient adjacent to the brain transient at the next moment. state, to determine the number of times each brain transient transitions to various brain transients at the next moment; based on the number of times each brain transient transitions to various brain transients at the next moment, and the number of occurrences of each brain transient, Determine the probability that each brain transient transitions into various brain transients. For example, the ratio between the number of times each brain transient transitions to various brain transients at the next moment and the number of occurrences of each brain transient is used as the probability that each brain transient transitions to various brain transients; Based on the probability that each brain transient transitions to various brain transients, the state transition probability matrix between multiple brain transient states is determined, that is, the probability of each brain transient transitioning to various brain transients is regarded as the state transition probability. For each row of elements in the matrix, the state transition probability matrix is obtained.

举例来说，时间窗口的数量为8，则脑瞬态序列中包含了8个脑瞬态，例如，脑瞬态序列为[脑瞬态1、脑瞬态2、脑瞬态1、脑瞬态3、脑瞬态4、脑瞬态2、脑瞬态1、脑瞬态4]。则脑瞬态1的出现次数为3，脑瞬态2出现的次数为2、脑瞬态3出现的次数为1，脑瞬态4出现的次数为2，故这四种脑瞬态对应的初始状态概率分别为：[3/8、1/4、1/8、1/4]。For example, if the number of time windows is 8, then the brain transient sequence contains 8 brain transients. For example, the brain transient sequence is [brain transient 1, brain transient 2, brain transient 1, brain transient state 3, brain transient 4, brain transient 2, brain transient 1, brain transient 4]. Then the number of occurrences of brain transient 1 is 3, the number of occurrences of brain transient 2 is 2, the number of occurrences of brain transient 3 is 1, and the number of occurrences of brain transient 4 is 2, so the corresponding numbers of these four brain transients are The initial state probabilities are: [3/8, 1/4, 1/8, 1/4].

然后，针对脑瞬态1，下一时刻转移为脑瞬态1的次数为0，下一时刻转移为脑瞬态2的次数为1，下一时刻转移为脑瞬态3的次数为1，下一时刻转移为脑瞬态4的次数为1，则脑瞬态1下一时刻转移为四种脑瞬态的概率分别为：0、1/3、1/3、1/3。针对脑瞬态2，下一时刻转移为脑瞬态1的次数为2，下一时刻转移为脑瞬态2的次数为0，下一时刻转移为脑瞬态3的次数为0，下一时刻转移为脑瞬态4的次数为0，则脑瞬态2下一时刻转移为四种脑瞬态的概率分别为：1、0、0、0。针对脑瞬态3，下一时刻转移为脑瞬态1的次数为0，下一时刻转移为脑瞬态2的次数为0，下一时刻转移为脑瞬态3的次数为0，下一时刻转移为脑瞬态4的次数为1，则脑瞬态3下一时刻转移为四种脑瞬态的概率分别为：0、0、0、1。针对脑瞬态4，下一时刻转移为脑瞬态1的次数为0，下一时刻转移为脑瞬态2的次数为1，下一时刻转移为脑瞬态3的次数为0，下一时刻转移为脑瞬态4的次数为0，则脑瞬态4下一时刻转移为四种脑瞬态的概率分别为：0、1、0、0。Then, for brain transient 1, the number of transitions to brain transient 1 at the next moment is 0, the number of transitions to brain transient 2 at the next moment is 1, and the number of transitions to brain transient 3 at the next moment is 1. The number of times that brain transient 1 transitions to brain transient 4 at the next moment is 1, then the probabilities that brain transient 1 transitions to four brain transient states at the next moment are: 0, 1/3, 1/3, and 1/3. For brain transient 2, the number of transitions to brain transient 1 at the next moment is 2, the number of transitions to brain transient 2 at the next moment is 0, the number of transitions to brain transient 3 at the next moment is 0, and the number of transitions to brain transient 3 at the next moment is 0. The number of times a moment transitions to brain transient 4 is 0, then the probabilities of brain transient 2 transitioning to four brain transient states at the next moment are: 1, 0, 0, 0. For brain transient 3, the number of transitions to brain transient 1 at the next moment is 0, the number of transitions to brain transient 2 at the next moment is 0, the number of transitions to brain transient 3 at the next moment is 0, and the number of transitions to brain transient 3 at the next moment is 0. The number of times a moment transitions to brain transient 4 is 1, then the probabilities that brain transient 3 transitions to four brain transient states at the next moment are: 0, 0, 0, and 1 respectively. For brain transient 4, the number of transitions to brain transient 1 at the next moment is 0, the number of transitions to brain transient 2 at the next moment is 1, the number of transitions to brain transient 3 at the next moment is 0, and the number of transitions to brain transient 3 at the next moment is 0. The number of times a moment transitions to brain transient 4 is 0, then the probabilities that brain transient 4 transitions to four brain transient states at the next moment are: 0, 1, 0, 0.

故这四种脑瞬态之间的状态转移概率矩阵为：Therefore, the state transition probability matrix between these four brain transient states is:

上面叙述了离线处理获取在线脑瞬态监测所需的参数的过程，下面结合离线过程获取到的参数，详细说明本申请的在线脑瞬态监测过程。The above describes the process of offline processing to obtain parameters required for online brain transient monitoring. The following describes the online brain transient monitoring process of this application in detail based on the parameters obtained in the offline process.

示例性的，进行在线监测时，脑电采集模块用于在当前时刻从多个通道上采集待监测对象的多个原始脑电信号。具体地，脑瞬态监测设备上可设置有脑瞬态监测按钮，当按下该按钮后，脑电采集模块就会一直采集待监测对象的脑电信号。本申请中主要以分析当前时刻的脑瞬态为例进行说明，故主要以当前时刻采集到的脑电信号为例进行说明，其中，当前时刻可以为进行在线监测后的任意一个时刻。For example, when performing online monitoring, the EEG collection module is used to collect multiple original EEG signals of the object to be monitored from multiple channels at the current moment. Specifically, the brain transient monitoring device can be provided with a brain transient monitoring button. When the button is pressed, the EEG collection module will continue to collect the EEG signals of the subject to be monitored. In this application, the analysis of brain transients at the current moment is mainly used as an example. Therefore, the electroencephalogram signal collected at the current moment is mainly used as an example. The current moment can be any moment after online monitoring.

进一步地，脑瞬态分析模块，用于基于重组矩阵对多个待分析脑电信号进行溯源分析，得到多个第一脑电源信号，即基于重组矩阵对多个原始脑电信号进行逆分解，得到与多个源极子对应的多个第一脑电源信号，即基于重组矩阵对多个原始脑电信号进行映射，可得到多个第一脑电源信号。然后，脑瞬态分析模块从每个第一脑电源信号中截取与当前时刻对应的待分析脑电信号，得到多个待分析脑电信号。即以当前时刻开始，分别从每个第一脑电源信号往前截取脑电信号，得到多个待分析脑电信号，其中，截取的长度为预设窗口所对应的长度。应说明的是，由于本申请是通过预设时间窗口进行信号截取，对各个时刻的脑瞬态进行分析，故当某个时刻之前的脑电信号的长度不够加窗时，则这个时刻可以不进行脑瞬态的监测。例如，预设时间窗口为30ms，则在启动在线监测功能后，前30ms的脑电信号无法加窗，故前30ms内可以不去输出各个时刻的脑瞬态，从30ms开始可以实时监测待监测对象在各个时刻下的脑瞬态，并将各个时刻下的脑瞬态进行输出显示。Further, the brain transient analysis module is used to conduct traceability analysis of multiple EEG signals to be analyzed based on the reorganization matrix, and obtain multiple first brain power signals, that is, to perform inverse decomposition of multiple original EEG signals based on the reorganization matrix, To obtain multiple first brain power signals corresponding to multiple source sub-elements, that is, by mapping multiple original brain power signals based on the recombination matrix, multiple first brain power signals can be obtained. Then, the brain transient analysis module intercepts the EEG signal to be analyzed corresponding to the current moment from each first brain power signal, and obtains multiple EEG signals to be analyzed. That is, starting from the current moment, the EEG signals are intercepted forward from each first EEG power signal to obtain multiple EEG signals to be analyzed, where the intercepted length is the length corresponding to the preset window. It should be noted that since this application intercepts signals through a preset time window and analyzes brain transients at each moment, when the length of the EEG signal before a certain moment is not enough to add a window, this moment does not need to be Perform brain transient monitoring. For example, if the preset time window is 30ms, then after the online monitoring function is started, the EEG signal of the first 30ms cannot be windowed, so the brain transients at each moment do not need to be output in the first 30ms, and the brain transients to be monitored can be monitored in real time starting from 30ms. The brain transient state of the subject at each moment, and the brain transient state at each moment is output and displayed.

进一步地，脑瞬态分析模块对多个待分析脑电信号进行处理，得到第一延时自协方差矩阵，即获取该多个待分析脑电信号对应的自协方差矩阵。例如，可以将每个待分析脑电信号在各个时刻下的幅值作为矩阵中的一行或一列元素，得到第三幅值矩阵。然后，基于幅值矩阵，可以确定出第一延时自协方差矩阵。Further, the brain transient analysis module processes multiple EEG signals to be analyzed to obtain a first delayed autocovariance matrix, that is, the autocovariance matrix corresponding to the multiple EEG signals to be analyzed is obtained. For example, the amplitude of each EEG signal to be analyzed at each moment can be used as a row or column element in the matrix to obtain a third amplitude matrix. Then, based on the amplitude matrix, the first delay autocovariance matrix can be determined.

进一步地，脑瞬态分析模块确定第一延时自协方差矩阵与多个预设延时自协方差矩阵之间的相似度，得到多个第一相似度。然后，根据多个第一相似度、所述待监测对象在上一时刻的脑瞬态、多种脑瞬态之间的状态转移概率矩阵、所述多种脑瞬态对应的初始状态概率，确定所述待监测对象当前时刻处于各种脑瞬态的概率。Further, the brain transient analysis module determines the similarity between the first delay autocovariance matrix and a plurality of preset delay autocovariance matrices, and obtains a plurality of first similarities. Then, based on the plurality of first similarities, the brain transient of the object to be monitored at the previous moment, the state transition probability matrix between the multiple brain transients, and the initial state probabilities corresponding to the multiple brain transients, Determine the probability that the subject to be monitored is in various brain transient states at the current moment.

示例性的，若当前时刻是首次确定所述待监测对象的脑瞬态时，则脑瞬态分析模块基于所述多种脑瞬态对应的初始状态概率与所述多个第一相似度，确定所述待监测对象当前时刻处于各种脑瞬态的概率。可选地，脑瞬态分析模块将所述多个第一相似度输入到隐马尔科夫模型中与该当前时刻对应的节点，得到待监测对象当前时刻处于各种脑瞬态的概率，其中，该隐马尔科夫模型中嵌入有与多种脑瞬态对应的初始状态概率以及状态转移概率矩阵。则将多个第一相似度输入到隐马尔科夫模型对应的节点后，可将多种脑瞬态对应的初始状态概率与多个第一相似度进行点乘，得到确定所述待监测对象当前时刻处于各种脑瞬态的概率，即将每种脑瞬态对应的第一相似度以及初始状态概率进行乘积，得到确定所述待监测对象当前时刻处于每种脑瞬态的概率。For example, if the current moment is when the brain transient of the object to be monitored is determined for the first time, the brain transient analysis module is based on the initial state probabilities corresponding to the multiple brain transients and the multiple first similarities, Determine the probability that the subject to be monitored is in various brain transient states at the current moment. Optionally, the brain transient analysis module inputs the plurality of first similarities into the node corresponding to the current moment in the hidden Markov model to obtain the probability that the object to be monitored is in various brain transient states at the current moment, where , the Cain Markov model is embedded with initial state probabilities and state transition probability matrices corresponding to various brain transient states. After inputting multiple first similarities into the corresponding nodes of the hidden Markov model, the initial state probabilities corresponding to various brain transient states can be dot-multiplied by the multiple first similarities to determine the object to be monitored. The probability of being in various brain transient states at the current moment is determined by multiplying the first similarity corresponding to each brain transient state and the initial state probability to determine the probability that the subject to be monitored is in each brain transient state at the current moment.

示例性的，若当前时刻是首次确定所述待监测对象的脑瞬态时，脑瞬态分析模块基于所述状态转移概率矩阵、所述多个第一相似度以及所述待监测对象在上一时刻的脑瞬态，确定所述待监测对象当前时刻处于各种脑瞬态的概率。可选地，脑瞬态分析模块可将所述多个第一相似度输入到隐马尔科夫模型中与该当前时刻对应的节点，得到待监测对象当前时刻处于各种脑瞬态的概率。具体地，将多个第一相似度输入到隐马尔科夫模型对应的节点后，可从隐马尔科夫模型中获取上一个节点（即上一时刻）的脑瞬态；然后，从状态转移概率矩阵中获取所述待监测对象在上一时刻的脑瞬态转移为各种脑瞬态的概率序列；根据所述概率序列与所述多个第一相似度，确定所述待监测对象当前时刻处于各种脑瞬态的概率。例如，将所述概率序列与所述多个第一相似度进行点乘，得到待监测对象当前时刻处于各种脑瞬态的概率，即将每种脑瞬态对应的第一相似度以及上一时刻的脑瞬态转移为这种脑瞬态的概率进行乘积，得到待监测对象当前时刻处于各种脑瞬态的概率。For example, if the current moment is when the brain transient state of the object to be monitored is determined for the first time, the brain transient analysis module is based on the state transition probability matrix, the plurality of first similarities and the state of the object to be monitored. The brain transient state at a moment is determined to determine the probability that the subject to be monitored is in various brain transient states at the current moment. Optionally, the brain transient analysis module can input the plurality of first similarities into the node corresponding to the current moment in the hidden Markov model to obtain the probability that the subject to be monitored is in various brain transient states at the current moment. Specifically, after inputting multiple first similarities to the corresponding node of the hidden Markov model, the brain transient state of the previous node (i.e., the previous moment) can be obtained from the hidden Markov model; then, from the state transition Obtain the probability sequence of the brain transient transition of the subject to be monitored at the previous moment into various brain transients from the probability matrix; determine the current state of the subject to be monitored based on the probability sequence and the plurality of first similarities. The probability of being in various brain transients at any time. For example, dot multiplication of the probability sequence and the plurality of first similarities is performed to obtain the probability that the subject to be monitored is in various brain transient states at the current moment, that is, the first similarity corresponding to each brain transient state and the previous similarity are obtained. The brain transient transition at a moment is multiplied by the probability of such a brain transient to obtain the probability that the subject to be monitored is in various brain transient states at the current moment.

最后，脑瞬态分析模块基于该待监测对象当前时刻处于各种脑瞬态的概率，确定待检测对象在当前时刻的脑瞬态。例如，脑瞬态分析模块将概率最大所对应的脑瞬态作为待检测对象在当前时刻的脑瞬态。Finally, the brain transient analysis module determines the brain transient state of the subject to be detected at the current moment based on the probability that the subject to be monitored is in various brain transient states at the current moment. For example, the brain transient analysis module uses the brain transient corresponding to the maximum probability as the brain transient of the object to be detected at the current moment.

在本申请的一个实施方式中，脑瞬态监测设备还包括显示模块。则在获取待检测对象在当前时刻的脑瞬态后，可通过显示模块显示待检测对象在当前时刻的脑瞬态，以便医生及时观看到待监测对象在当前时刻的脑瞬态。In one embodiment of the present application, the brain transient monitoring device further includes a display module. After obtaining the brain transient state of the object to be detected at the current moment, the display module can display the brain transient state of the object to be detected at the current moment, so that the doctor can view the brain transient state of the object to be monitored at the current moment in time.

在本申请的一个实施方式中，脑瞬态监测设备还可以包括报警模块。则在获取待检测对象在当前时刻的脑瞬态后，可将脑瞬态与预设脑瞬态进行比对，以确定当前时刻的脑瞬态是否存在危险；若是，则通过报警模块进行报警，比如，向远端联系人或者医生发送报警信息，或者播放报警信号。In one embodiment of the present application, the brain transient monitoring device may further include an alarm module. After obtaining the brain transient of the object to be detected at the current moment, the brain transient can be compared with the preset brain transient to determine whether the brain transient at the current moment is dangerous; if so, an alarm is issued through the alarm module , for example, sending alarm information to remote contacts or doctors, or playing alarm signals.

在本申请的一个实施方式中，脑瞬态监测设备还可将待监测对象在各个时刻下获取到的脑瞬态输出给其他设备，从而便于其他设备基于待监测对象在各个时刻的脑瞬态进行疾病分析、健康监测、个体化神经调控治疗策略等等。或者，脑瞬态监测设备还可将待监测对象在各个时刻下的脑瞬态进行缓存以及存档，以便后续进行分析。In one embodiment of the present application, the brain transient monitoring device can also output the brain transients obtained by the subject to be monitored at each moment to other devices, thereby facilitating other devices to base on the brain transients of the subject to be monitored at each moment. Perform disease analysis, health monitoring, personalized neuromodulation treatment strategies, and more. Alternatively, the brain transient monitoring equipment can also cache and archive the brain transients of the subject to be monitored at various moments for subsequent analysis.

可以看出，在本申请实施例中，首先离线计算出脑电信号到脑电源信号的重组矩阵，以及用于表征脑瞬态的预设延时自协方差矩阵，以及多种脑瞬态之间的状态转移概率矩阵，以及多种脑瞬态的初始状态概率，这样进行实时的在线监测时，可以利用离线计算好的重组矩阵，直接对在线采集到的多个原始脑电信号进行溯源，得到多个第一脑电源信号；以及从每个第一脑电源信号中截取与当前时刻对应的待分析脑电信号，得到多个待分析脑电信号，并获取多个待分析脑电信号对应的第一延时自协方差矩阵。然后，直接使用多个预设延时自协方差矩阵，确定出多个第一相似度。最后，直接基于多个第一相似度、待监测对象在上一时刻的脑瞬态、多种脑瞬态之间的状态转移概率矩阵、多种脑瞬态对应的初始状态概率，确定出待监测对象当前时刻处于各种脑瞬态的概率，进而确定出当前时刻的脑状态。由于在线监测的大部分复杂参数都是离线计算好的，大幅度的减轻了在线计算的复杂度，因此在线监测时可以快速地计算出当前时刻的脑状态，从而实现实时在线地监测大脑的脑状态。并且在确定当前时刻的脑瞬态时，不会只利用当前时刻下的脑电信号，与各种脑瞬态之间的相似度，还会结合上一个时刻的脑瞬态，以及多种脑瞬态之间的转移概率，从而可以精确地确定出当前时刻的脑瞬态，提高了对脑瞬态的监测精度。It can be seen that in the embodiment of the present application, the recombination matrix from the brain electrical signal to the brain power signal is first calculated offline, as well as the preset delay autocovariance matrix used to characterize brain transients, and a variety of brain transients. The state transition probability matrix between states, as well as the initial state probabilities of various brain transients, so that when performing real-time online monitoring, the offline calculated reorganization matrix can be used to directly trace the source of multiple original EEG signals collected online. Obtain a plurality of first brain power signals; and intercept an EEG signal to be analyzed corresponding to the current moment from each first brain power signal, obtain a plurality of EEG signals to be analyzed, and obtain the corresponding EEG signals of the multiple EEG signals to be analyzed. The first delay autocovariance matrix of . Then, multiple preset delay autocovariance matrices are directly used to determine multiple first similarities. Finally, directly based on multiple first similarities, the brain transient of the object to be monitored at the previous moment, the state transition probability matrix between multiple brain transients, and the initial state probabilities corresponding to multiple brain transients, the target to be monitored is determined. Monitor the probability that the subject is in various brain transient states at the current moment, and then determine the brain state at the current moment. Since most of the complex parameters for online monitoring are calculated offline, the complexity of online calculations is greatly reduced. Therefore, the brain state at the current moment can be quickly calculated during online monitoring, thereby achieving real-time online monitoring of the brain. state. And when determining the brain transient at the current moment, we will not only use the similarity between the EEG signal at the current moment and various brain transients, but also combine the brain transient at the previous moment and various brain transients. The transition probability between transient states can accurately determine the brain transient at the current moment, improving the accuracy of monitoring brain transients.

上面分别从离线和在线两个角度叙述了对脑瞬态进行监测的过程。下面结合附图，整体叙述本申请的脑瞬态监测过程。The process of monitoring brain transients is described above from two perspectives: offline and online. The brain transient monitoring process of the present application will be described as a whole in conjunction with the accompanying drawings.

如图5所示，首先进行离线处理，获取待监测对象的头皮脑电，得到离线脑电信号，以及获取待监测对象的头部的核磁共振图像；然后，对核磁共振图像进行脑皮层重构，得到三维脑模型，以及三维脑模型中的皮层网格，从而构建出Lead-field矩阵，以及重组矩阵，基于Lead-field矩阵对离线脑电信号进行逆分解，得到多个源极子对应的多个脑电源信号；基于多个脑电源信号，可构建出多个高斯观测模型，其中，每个高斯观测模型用于表征一种脑瞬态。同样，基于多个脑电源信号，可构造出离线处理时的脑瞬态序列；基于该脑瞬态序列，确定出多种脑瞬态对应的初始状态概率，以及多种脑瞬态之间的状态转移概率矩阵。As shown in Figure 5, offline processing is first performed to obtain the scalp EEG of the subject to be monitored, obtain the offline EEG signal, and obtain the MRI image of the head of the subject to be monitored; then, perform cerebral cortex reconstruction on the MRI image , obtain a three-dimensional brain model, and the cortical grid in the three-dimensional brain model, thereby constructing a lead-field matrix and a reorganization matrix. Based on the lead-field matrix, the offline EEG signals are inversely decomposed to obtain multiple source sub-elements. Multiple brain power signals; based on multiple brain power signals, multiple Gaussian observation models can be constructed, where each Gaussian observation model is used to represent a brain transient state. Similarly, based on multiple brain power signals, a brain transient sequence during offline processing can be constructed; based on this brain transient sequence, the initial state probabilities corresponding to multiple brain transients are determined, as well as the relationships between multiple brain transients. State transition probability matrix.

如图5所示，进行在线监测时，可获取当前时刻的脑电信号，基于重组矩阵将脑电信号逆分解（即脑电源信号重构）为多个脑电源信号；使用预设时间窗口从脑电源信号中截取当前时刻的待分析脑电信号，并计算待分析脑电信号的第一延时自协方差矩阵；计算第一延时自协方差矩阵与每个高斯观测模型之间的第一相似度，得到多个第一相似度；然后将多个第一相似度输入到隐马尔科夫模型，基于隐马尔科夫模型中的多种脑瞬态对应的初始状态概率，以及多种脑瞬态之间的状态转移概率矩阵、以及历史脑瞬态序列，得到待监测对象在当前时刻下的脑瞬态。As shown in Figure 5, during online monitoring, the EEG signal at the current moment can be obtained, and the EEG signal is inversely decomposed (i.e., the brain power signal is reconstructed) into multiple brain power signals based on the reorganization matrix; a preset time window is used from Intercept the EEG signal to be analyzed at the current moment from the brain power signal, and calculate the first delayed auto-covariance matrix of the EEG signal to be analyzed; calculate the first delayed auto-covariance matrix and each Gaussian observation model. One similarity, obtain multiple first similarities; then input the multiple first similarities into the hidden Markov model, based on the initial state probabilities corresponding to various brain transients in the hidden Markov model, and multiple The state transition probability matrix between brain transients and the historical brain transient sequence are used to obtain the brain transient state of the object to be monitored at the current moment.

参阅图6，图6为本申请实施例提供的一种脑瞬态监测方法的流程示意图。该方法应用于上述脑瞬态监测设备。该方法包括但不限于以下步骤内容：Refer to Figure 6, which is a schematic flow chart of a brain transient monitoring method provided by an embodiment of the present application. This method is applied to the brain transient monitoring device described above. The method includes but is not limited to the following steps:

S601：所述脑电采集模块在当前时刻从多个通道上采集待监测对象的多个原始脑电信号。S601: The EEG collection module collects multiple original EEG signals of the object to be monitored from multiple channels at the current moment.

S602：所述脑瞬态分析模块基于重组矩阵对所述多个原始脑电信号进行溯源分析，得到多个第一脑电源信号。S602: The brain transient analysis module performs traceability analysis on the multiple original EEG signals based on the reorganization matrix, and obtains multiple first brain power signals.

S603：所述脑瞬态分析模块从每个第一脑电源信号中截取与当前时刻对应的待分析脑电信号，得到多个待分析脑电信号。S603: The brain transient analysis module intercepts the EEG signal to be analyzed corresponding to the current moment from each first brain power signal, and obtains multiple EEG signals to be analyzed.

S604：所述脑瞬态分析模块对所述多个待分析脑电信号进行处理，得到第一延时自协方差矩阵。S604: The brain transient analysis module processes the multiple EEG signals to be analyzed to obtain a first delayed autocovariance matrix.

S605：所述脑瞬态分析模块确定所述第一延时自协方差矩阵与多个预设延时自协方差矩阵之间的相似度，得到多个第一相似度，其中，每个预设延时自协方差矩阵用于表征一种脑瞬态。S605: The brain transient analysis module determines the similarity between the first delayed autocovariance matrix and multiple preset delayed autocovariance matrices, and obtains multiple first similarities, where each preset Let the delayed autocovariance matrix be used to represent a brain transient state.

S606：所述脑瞬态分析模块根据所述多个第一相似度、所述待监测对象在上一时刻的脑瞬态、多种脑瞬态之间的状态转移概率矩阵、所述多种脑瞬态对应的初始状态概率，确定所述待监测对象当前时刻处于各种脑瞬态的概率。S606: The brain transient analysis module analyzes the brain transient state according to the plurality of first similarities, the brain transient state of the object to be monitored at the previous moment, the state transition probability matrix between the multiple brain transient states, and the multiple brain transient states. The initial state probability corresponding to the brain transient determines the probability that the subject to be monitored is in various brain transient states at the current moment.

S607：所述脑瞬态分析模块根据当前时刻处于各种脑瞬态的概率，确定所述待监测对象在当前时刻的脑瞬态。S607: The brain transient analysis module determines the brain transient state of the object to be monitored at the current moment based on the probability of being in various brain transient states at the current moment.

其中，步骤S601~步骤S607的具体实现过程，可参照上述离线分析模块、图像获取模块、脑电采集模块以及脑瞬态分析模块的具体功能，不再赘述。For the specific implementation process of steps S601 to S607, reference can be made to the specific functions of the above-mentioned offline analysis module, image acquisition module, EEG acquisition module and brain transient analysis module, and will not be described again.

参阅图7，图7是本申请实施例提供的一种电子设备的示意图。图7所示的电子设备700包括存储器701、处理器702、通信接口703以及总线704。其中，存储器701、处理器702、通信接口703通过总线704实现彼此之间的通信连接。电子设备700可以为上述脑瞬态监测装置。处理器702可以实现上述脑瞬态分析模块以及离线处理模块的功能。通信接口703可以实现上述图像获取模块和脑电采集模块的功能。Refer to FIG. 7 , which is a schematic diagram of an electronic device provided by an embodiment of the present application. The electronic device 700 shown in FIG. 7 includes a memory 701, a processor 702, a communication interface 703, and a bus 704. Among them, the memory 701, the processor 702, and the communication interface 703 realize communication connections between each other through the bus 704. The electronic device 700 may be the above-mentioned brain transient monitoring device. The processor 702 can implement the functions of the above-mentioned brain transient analysis module and offline processing module. The communication interface 703 can realize the functions of the above-mentioned image acquisition module and EEG acquisition module.

其中，处理器702可以集成有上述脑瞬态分析模块以及离线处理模块的功能，用于离线计算（即离线计算重组矩阵、预设延时自协方差矩阵、种脑瞬态之间的状态转移概率矩阵、多种脑瞬态对应的初始状态概率）和在线计算（即计算当前时刻的脑瞬态）。通信接口703可以集成有上述图像获取模块和脑电采集模块的功能，例如，从待监测对象采集脑电信号，或者获取待监测对象的核磁共振图像。Among them, the processor 702 can integrate the functions of the above-mentioned brain transient analysis module and offline processing module for offline calculation (ie, offline calculation of reorganization matrix, preset delay autocovariance matrix, and state transfer between brain transient states). Probability matrix, initial state probabilities corresponding to various brain transients) and online calculation (that is, calculating the brain transient at the current moment). The communication interface 703 can integrate the functions of the above-mentioned image acquisition module and EEG acquisition module, for example, to acquire EEG signals from the object to be monitored, or to acquire MRI images of the object to be monitored.

存储器701可以是只读存储器（Read Only Memory，ROM），静态存储设备，动态存储设备或者随机存取存储器（Random Access Memory，RAM）。存储器701可以存储程序，当存储器701中存储的程序被处理器702执行时，处理器702和通信接口703用于执行本申请实施例的脑瞬态监测方法的各个步骤。The memory 701 may be a read-only memory (Read Only Memory, ROM), a static storage device, a dynamic storage device or a random access memory (Random Access Memory, RAM). The memory 701 can store programs. When the program stored in the memory 701 is executed by the processor 702, the processor 702 and the communication interface 703 are used to execute various steps of the brain transient monitoring method according to the embodiment of the present application.

处理器702可以采用通用的中央处理器（Central Processing Unit，CPU），微处理器，应用专用集成电路（Application Specific Integrated Circuit，ASIC），图形处理器（graphics processing unit，GPU）或者一个或多个集成电路，用于执行相关程序，以实现本申请实施例的音频特征补偿装置或音频识别装置中的单元所需执行的功能，或者执行本申请方法实施例的脑瞬态监测方法。The processor 702 may be a general central processing unit (CPU), a microprocessor, an application specific integrated circuit (Application Specific Integrated Circuit, ASIC), a graphics processing unit (GPU), or one or more The integrated circuit is used to execute relevant programs to realize the functions required to be performed by the units in the audio characteristic compensation device or the audio recognition device according to the embodiment of the present application, or to perform the brain transient monitoring method according to the method embodiment of the present application.

处理器702还可以是一种集成电路芯片，具有信号的处理能力。在实现过程中，本申请的刺激强度设定方法中的各个步骤可以通过处理器702中的硬件的集成逻辑电路或者软件形式的指令完成。上述的处理器702还可以是通用处理器、数字信号处理器（DigitalSignal Processing，DSP）、专用集成电路（ASIC）、现成可编程门阵列（Field ProgrammableGate Array，FPGA）或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。可以实现或者执行本申请实施例中的公开的各方法、步骤及逻辑框图。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。结合本申请实施例所公开的方法的步骤可以直接体现为硬件译码处理器执行完成，或者用译码处理器中的硬件及软件模块组合执行完成。软件模块可以位于随机存储器，闪存、只读存储器，可编程只读存储器或者电可擦写可编程存储器、寄存器等本领域成熟的存储介质中。该存储介质位于存储器701，处理器702读取存储器701中的信息，结合其硬件完成本申请实施例的音频特征补偿装置或音频识别装置中包括的单元所需执行的功能，或者执行本申请方法实施例的脑瞬态监测方法中的各个步骤。The processor 702 may also be an integrated circuit chip with signal processing capabilities. During the implementation process, each step in the stimulation intensity setting method of the present application can be completed by instructions in the form of hardware integrated logic circuits or software in the processor 702 . The above-mentioned processor 702 can also be a general-purpose processor, a digital signal processor (Digital Signal Processing, DSP), an application specific integrated circuit (ASIC), an off-the-shelf programmable gate array (Field ProgrammableGate Array, FPGA) or other programmable logic devices, discrete Gate or transistor logic devices, discrete hardware components. Each method, step and logical block diagram disclosed in the embodiment of this application can be implemented or executed. A general-purpose processor may be a microprocessor or the processor may be any conventional processor, etc. The steps of the method disclosed in conjunction with the embodiments of the present application can be directly implemented by a hardware decoding processor, or executed by a combination of hardware and software modules in the decoding processor. The software module can be located in random access memory, flash memory, read-only memory, programmable read-only memory or electrically erasable programmable memory, registers and other mature storage media in this field. The storage medium is located in the memory 701. The processor 702 reads the information in the memory 701, and combines its hardware to complete the functions required to be performed by the units included in the audio feature compensation device or the audio recognition device in the embodiment of the present application, or to perform the method of the present application. Various steps in the brain transient monitoring method of the embodiment.

通信接口703使用例如但不限于收发器、输入-输出设备一类的收发装置，来实现电子设备700与其他设备或通信网络之间的通信。例如，可以通过通信接口703采集脑电信号。The communication interface 703 uses transceiver devices such as but not limited to transceivers and input-output devices to implement communication between the electronic device 700 and other devices or communication networks. For example, EEG signals can be collected through the communication interface 703.

总线704可包括在电子设备700各个部件（例如，存储器701、处理器702、通信接口703）之间传送信息的通路。Bus 704 may include a path that carries information between various components of electronic device 700 (eg, memory 701, processor 702, communication interface 703).

应注意，尽管图7所示电子设备700仅仅示出了存储器、处理器、通信接口，但是在具体实现过程中，本领域的技术人员应当理解，电子设备700还包括实现正常运行所必须的其他器件。同时，根据具体需要，本领域的技术人员应当理解，电子设备700还可包括实现其他附加功能的硬件器件。此外，本领域的技术人员应当理解，电子设备700也可仅仅包括实现本申请实施例所必须的器件，而不必包括图7中所示的全部器件。It should be noted that although the electronic device 700 shown in FIG. 7 only shows a memory, a processor, and a communication interface, during specific implementation, those skilled in the art will understand that the electronic device 700 also includes other components necessary for normal operation. device. At the same time, according to specific needs, those skilled in the art should understand that the electronic device 700 may also include hardware devices that implement other additional functions. In addition, those skilled in the art should understand that the electronic device 700 may only include components necessary to implement the embodiments of the present application, and does not necessarily include all components shown in FIG. 7 .

所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，上述描述的系统、装置和单元的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。Those skilled in the art can clearly understand that for the convenience and simplicity of description, the specific working processes of the systems, devices and units described above can be referred to the corresponding processes in the foregoing method embodiments, and will not be described again here.

在本申请所提供的几个实施例中，应该理解到，所揭露的系统、装置和方法，可以通过其它的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如，单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。In the several embodiments provided in this application, it should be understood that the disclosed systems, devices and methods can be implemented in other ways. For example, the device embodiments described above are only illustrative. For example, the division of units is only a logical function division. In actual implementation, there may be other division methods. For example, multiple units or components may be combined or integrated. to another system, or some features can be ignored, or not implemented. On the other hand, the coupling or direct coupling or communication connection between each other shown or discussed may be through some interfaces, and the indirect coupling or communication connection of the devices or units may be in electrical, mechanical or other forms.

作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。A unit described as a separate component may or may not be physically separate. A component shown as a unit may or may not be a physical unit, that is, it may be located in one place, or it may be distributed to multiple network units. Some or all of the units can be selected according to actual needs to achieve the purpose of the solution of this embodiment.

另外，在本申请各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。In addition, each functional unit in each embodiment of the present application can be integrated into one processing unit, each unit can exist physically alone, or two or more units can be integrated into one unit.

所述功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备（可以是个人计算机，服务器，或者网络设备等）执行本申请各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器（Read-Only Memory，ROM）、随机存取存储器（Random Access Memory，RAM）、磁碟或者光盘等各种可以存储程序代码的介质。If the functions are implemented in the form of software functional units and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present application is essentially or the part that contributes to the existing technology or the part of the technical solution can be embodied in the form of a software product. The computer software product is stored in a storage medium, including Several instructions are used to cause a computer device (which may be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the methods described in various embodiments of this application. The aforementioned storage media include: U disk, mobile hard disk, read-only memory (ROM), random access memory (RAM), magnetic disk or optical disk and other media that can store program code. .

以上所述，仅为本申请的具体实施方式，但本申请的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本申请揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本申请的保护范围之内。因此，本申请的保护范围应以所述权利要求的保护范围为准。The above are only specific embodiments of the present application, but the protection scope of the present application is not limited thereto. Any person familiar with the technical field can easily think of changes or substitutions within the technical scope disclosed in the present application. should be covered by the protection scope of this application. Therefore, the protection scope of this application should be subject to the protection scope of the claims.

本申请实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行以实现如上述方法实施例中记载的任何一种脑瞬态监测方法的部分或全部步骤。Embodiments of the present application also provide a computer-readable storage medium that stores a computer program, and the computer program is executed by a processor to implement any brain transient state as described in the above method embodiments. Monitor some or all steps of a method.

本申请实施例还提供一种计算机程序产品，所述计算机程序产品包括存储了计算机程序的非瞬时性计算机可读存储介质，所述计算机程序可操作来使计算机执行如上述方法实施例中记载的任何一种脑瞬态监测方法的部分或全部步骤。Embodiments of the present application also provide a computer program product. The computer program product includes a non-transitory computer-readable storage medium storing a computer program. The computer program is operable to cause the computer to execute the steps described in the above method embodiments. Some or all steps of any brain transient monitoring method.

需要说明的是，对于前述的各方法实施例，为了简单描述，故将其都表述为一系列的动作组合，但是本领域技术人员应该知悉，本申请并不受所描述的动作顺序的限制，因为依据本申请，某些步骤可以采用其他顺序或者同时进行。其次，本领域技术人员也应该知悉，说明书中所描述的实施例均属于可选实施例，所涉及的动作和模块并不一定是本申请所必须的。It should be noted that for the sake of simple description, the foregoing method embodiments are expressed as a series of action combinations. However, those skilled in the art should know that the present application is not limited by the described action sequence. Because in accordance with this application, certain steps may be performed in other orders or simultaneously. Secondly, those skilled in the art should also know that the embodiments described in the specification are all optional embodiments, and the actions and modules involved are not necessarily necessary for this application.

在上述实施例中，对各个实施例的描述都各有侧重，某个实施例中没有详述的部分，可以参见其他实施例的相关描述。In the above embodiments, each embodiment is described with its own emphasis. For parts that are not described in detail in a certain embodiment, please refer to the relevant descriptions of other embodiments.

在本申请所提供的几个实施例中，应该理解到，所揭露的装置，可通过其它的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性或其它的形式。In the several embodiments provided in this application, it should be understood that the disclosed device can be implemented in other ways. For example, the device embodiments described above are only illustrative. For example, the division of the units is only a logical function division. In actual implementation, there may be other division methods. For example, multiple units or components may be combined or may be Integrated into another system, or some features can be ignored, or not implemented. On the other hand, the coupling or direct coupling or communication connection between each other shown or discussed may be through some interfaces, and the indirect coupling or communication connection of the devices or units may be in electrical or other forms.

所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。The units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in one place, or they may be distributed to multiple network units. Some or all of the units can be selected according to actual needs to achieve the purpose of the solution of this embodiment.

另外，在本申请各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件程序模块的形式实现。In addition, each functional unit in each embodiment of the present application can be integrated into one processing unit, each unit can exist physically alone, or two or more units can be integrated into one unit. The above integrated units can be implemented in the form of hardware or software program modules.

所述集成的单元如果以软件程序模块的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储器中。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储器中，包括若干指令用以使得一台计算机设备（可为个人计算机、服务器或者网络设备等）执行本申请各个实施例所述方法的全部或部分步骤。而前述的存储器包括：U盘、只读存储器（ROM，Read-Only Memory）、随机存取存储器（RAM，Random Access Memory）、移动硬盘、磁碟或者光盘等各种可以存储程序代码的介质。The integrated unit, if implemented in the form of a software program module and sold or used as an independent product, may be stored in a computer-readable memory. Based on this understanding, the technical solution of the present application is essentially or contributes to the existing technology, or all or part of the technical solution can be embodied in the form of a software product, and the computer software product is stored in a memory, It includes several instructions to cause a computer device (which can be a personal computer, a server or a network device, etc.) to execute all or part of the steps of the methods described in various embodiments of this application. The aforementioned memory includes: U disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), mobile hard disk, magnetic disk or optical disk and other media that can store program code.

本领域普通技术人员可以理解上述实施例的各种方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成，该程序可以存储于一计算机可读存储器中，存储器可以包括：闪存盘、只读存储器（英文：Read-Only Memory，简称：ROM）、随机存取器（英文：Random Access Memory，简称：RAM）、磁盘或光盘等。Those of ordinary skill in the art can understand that all or part of the steps in the various methods of the above embodiments can be completed by instructing relevant hardware through a program. The program can be stored in a computer-readable memory. The memory can include: a flash disk. , read-only memory (English: Read-Only Memory, abbreviation: ROM), random access device (English: Random Access Memory, abbreviation: RAM), magnetic disk or optical disk, etc.

以上对本申请实施例进行了详细介绍，本文中应用了具体个例对本申请的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本申请的方法及其核心思想；同时，对于本领域的一般技术人员，依据本申请的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本申请的限制。The embodiments of the present application have been introduced in detail above. Specific examples are used in this article to illustrate the principles and implementation methods of the present application. The description of the above embodiments is only used to help understand the method and the core idea of the present application; at the same time, for Those of ordinary skill in the art will have changes in the specific implementation and application scope based on the ideas of the present application. In summary, the content of this description should not be understood as a limitation of the present application.

Claims (6)

Translated fromChinese

1.一种脑瞬态监测设备，其特征在于，1. A brain transient monitoring device, characterized by:所述脑瞬态监测设备包括图像获取模块、离线处理模块、脑电采集模块和脑瞬态分析模块；The brain transient monitoring equipment includes an image acquisition module, an offline processing module, an EEG acquisition module and a brain transient analysis module;所述图像获取模块，用于获取待监测对象的脑部的原始核磁共振图像；The image acquisition module is used to acquire original MRI images of the brain of the subject to be monitored;所述脑电采集模块，用于在预设时间段内在多个通道上采集所述待监测对象的多个离线脑电信号；The EEG collection module is used to collect multiple offline EEG signals of the object to be monitored on multiple channels within a preset time period;所述离线处理模块，用于基于所述原始核磁共振图像，得到所述待监测对象的三维脑模型；The offline processing module is used to obtain a three-dimensional brain model of the object to be monitored based on the original MRI image;获取所述三维脑模型中的多个皮层网格；Obtaining a plurality of cortical grids in the three-dimensional brain model;将所述多个通道上用于脑电信号采集的多个电极在所述待监测对象的脑部的位置与所述多个皮层网格在脑部的位置进行对齐，得到与所述待监测对象对应的Lead-field矩阵；Align the positions of the multiple electrodes on the multiple channels used for EEG signal collection on the brain of the subject to be monitored with the positions of the multiple cortical grids on the brain to obtain the Lead-field matrix corresponding to the object;获取所述多个离线脑电信号的协方差矩阵；Obtain the covariance matrix of the multiple offline EEG signals;基于所述协方差矩阵以及所述Lead-field矩阵，确定重组矩阵；Based on the covariance matrix and the lead-field matrix, determine the reorganization matrix;基于所述Lead-field矩阵，对所述多个离线脑电信号进行逆分解，得到所述多个皮层网格对应的多个皮层信号；将所述三维脑模型划分为多个源极子；Based on the lead-field matrix, inversely decompose the multiple offline EEG signals to obtain multiple cortical signals corresponding to the multiple cortical grids; divide the three-dimensional brain model into multiple source sub-elements;对所述多个皮层信号进行加权映射，得到与所述多个源极子对应的多个第二脑电源信号；Perform weighted mapping on the plurality of cortical signals to obtain a plurality of second brain power signals corresponding to the plurality of source sub-elements;基于预设时间窗口，对每个第二脑电源信号进行多次加窗，以对每个第二脑电源信号进行分割，得到与每个第二脑电源信号对应的多个子脑电源信号；Based on the preset time window, perform multiple windowing on each second brain power signal to segment each second brain power signal and obtain multiple sub-brain power signals corresponding to each second brain power signal;基于每个第二脑电源信号对应的多个子脑电源信号，得到多个预设延时自协方差矩阵、多种脑瞬态之间的状态转移概率矩阵以及所述多种脑瞬态对应的初始状态概率，其中，每个预设延时自协方差矩阵用于表征一种脑瞬态；Based on multiple sub-brain power signals corresponding to each second brain power signal, multiple preset delay auto-covariance matrices, state transition probability matrices between multiple brain transient states, and corresponding values of the multiple brain transient states are obtained. Initial state probability, where each preset delay autocovariance matrix is used to characterize a brain transient;所述脑电采集模块，用于在当前时刻从所述多个通道上采集所述待监测对象的多个原始脑电信号；The EEG collection module is used to collect multiple original EEG signals of the object to be monitored from the multiple channels at the current moment;所述脑瞬态分析模块，用于基于所述重组矩阵对所述多个原始脑电信号进行溯源分析，得到多个第一脑电源信号；The brain transient analysis module is used to perform traceability analysis on the plurality of original EEG signals based on the reorganization matrix to obtain a plurality of first brain power signals;从每个第一脑电源信号中截取与当前时刻对应的待分析脑电信号，得到多个待分析脑电信号；Intercept the EEG signal to be analyzed corresponding to the current moment from each first brain power signal to obtain multiple EEG signals to be analyzed;对所述多个待分析脑电信号进行处理，得到第一延时自协方差矩阵；Process the plurality of EEG signals to be analyzed to obtain a first delayed autocovariance matrix;确定所述第一延时自协方差矩阵与所述多个预设延时自协方差矩阵之间的相似度，得到多个第一相似度；Determine the similarity between the first delay autocovariance matrix and the plurality of preset delay autocovariance matrices, and obtain a plurality of first similarities;根据所述多个第一相似度、所述待监测对象在上一时刻的脑瞬态、所述状态转移概率矩阵以及所述初始状态概率，确定所述待监测对象当前时刻处于各种脑瞬态的概率；具体用于：According to the plurality of first similarities, the brain transient state of the object to be monitored at the previous moment, the state transition probability matrix and the initial state probability, it is determined that the object to be monitored is in various brain transient states at the current moment. The probability of state; specifically used for:若当前时刻为首次确定所述待监测对象的脑瞬态时，基于所述多种脑瞬态对应的初始状态概率与所述多个第一相似度，确定所述待监测对象当前时刻处于各种脑瞬态的概率；If the current moment is the first time the brain transient state of the subject to be monitored is determined, based on the initial state probabilities corresponding to the multiple brain transient states and the plurality of first similarities, it is determined that the subject to be monitored is in each state at the current moment. the probability of brain transients;若当前时刻不是首次确定所述待监测对象的脑瞬态时，基于所述状态转移概率矩阵，确定所述待监测对象在上一时刻的脑瞬态转移为各种脑瞬态的概率序列；根据所述概率序列与所述多个第一相似度，确定所述待监测对象当前时刻处于各种脑瞬态的概率；If the current moment is not the first time that the brain transient state of the object to be monitored is determined, based on the state transition probability matrix, it is determined that the brain transient state of the object to be monitored at the previous moment is transformed into a probability sequence of various brain transient states; Determine the probability that the subject to be monitored is in various brain transient states at the current moment according to the probability sequence and the plurality of first similarities;根据当前时刻处于各种脑瞬态的概率，确定所述待监测对象在当前时刻的脑瞬态。According to the probability of being in various brain transient states at the current moment, the brain transient state of the object to be monitored at the current moment is determined.2.根据权利要求1所述的设备，其特征在于，2. The device according to claim 1, characterized in that,在基于每个第二脑电源信号对应的多个子脑电源信号，得到多个预设延时自协方差矩阵、多种脑瞬态之间的状态转移概率矩阵以及所述多种脑瞬态对应的初始状态概率方面，所述离线处理模块，具体用于：Based on multiple sub-brain power signals corresponding to each second brain power signal, multiple preset delay auto-covariance matrices, state transition probability matrices between multiple brain transient states, and correspondences between the multiple brain transient states are obtained In terms of the initial state probability, the offline processing module is specifically used for:获取所述多个第二脑电源信号在每个时间窗口下的多个子脑电源信号；Obtain a plurality of sub-brain power signals of the plurality of second brain power signals in each time window;基于每个时间窗口下的多个子脑电源信号，确定每个时间窗口下的第二延时自协方差矩阵；Based on the multiple sub-brain power signals under each time window, determine the second delay autocovariance matrix under each time window;基于每个时间窗口的第二延时自协方差矩阵，得到多次加窗对应的多个时间窗口下的多个第二延时自协方差矩阵；Based on the second delay autocovariance matrix of each time window, obtain multiple second delay autocovariance matrices under multiple time windows corresponding to multiple windowings;对所述多个第二延时自协方差矩阵进行分组，得到多个延时自协方差矩阵组；Group the plurality of second delay autocovariance matrices to obtain a plurality of delay autocovariance matrix groups;将每个延时自协方差矩阵组的中心，作为一个预设延时自协方差矩阵，得到所述多个预设延时自协方差矩阵；Use the center of each delay autocovariance matrix group as a preset delay autocovariance matrix to obtain the plurality of preset delay autocovariance matrices;确定每个时间窗口下的第二延时自协方差矩阵，与所述多个预设延时自协方差矩阵之间的多个第二相似度；Determine a plurality of second similarities between the second delay autocovariance matrix in each time window and the plurality of preset delay autocovariance matrices;基于每个时间窗口下的多个第二相似度，确定所述待监测对象在每个时间窗口下的脑瞬态；Based on the plurality of second similarities under each time window, determine the brain transient state of the object to be monitored under each time window;按照时间的先后顺序，对所述多个时间窗口下的多个脑瞬态进行排列，得到脑瞬态序列；Arranging multiple brain transients under the multiple time windows in order of time to obtain a brain transient sequence;基于所述脑瞬态序列，确定所述多种脑瞬态之间的状态转移概率矩阵，以及所述多种脑瞬态对应的初始状态概率。Based on the brain transient sequence, a state transition probability matrix between the multiple brain transient states is determined, as well as an initial state probability corresponding to the multiple brain transient states.3.根据权利要求2所述的设备，其特征在于，3. The device according to claim 2, characterized in that,在基于所述脑瞬态序列，确定所述多种脑瞬态之间的状态转移概率矩阵，以及所述多种脑瞬态对应的初始状态概率方面，所述离线处理模块，具体用于：In terms of determining the state transition probability matrix between the multiple brain transients based on the brain transient sequence, and the initial state probabilities corresponding to the multiple brain transients, the offline processing module is specifically used to:确定所述脑瞬态序列中每种脑瞬态出现的次数；determining the number of occurrences of each brain transient in the sequence of brain transients;基于每种脑瞬态出现的次数，以及所述脑瞬态序列中脑瞬态的数量，确定每种脑瞬态的出现概率；determining the probability of occurrence of each brain transient based on the number of occurrences of each brain transient and the number of brain transients in the sequence of brain transients;将每种脑瞬态的出现概率，作为所述多种脑瞬态对应的初始状态概率；The occurrence probability of each brain transient is used as the initial state probability corresponding to the multiple brain transients;针对每种脑瞬态，确定所述脑瞬态序列中下一时刻与该脑瞬态相邻的脑瞬态，基于下一时刻与该脑瞬态相邻的脑瞬态，确定该种脑瞬态下一时刻转移为各种脑瞬态的次数；For each brain transient, determine the brain transient adjacent to the brain transient at the next moment in the brain transient sequence, and determine the brain transient based on the brain transient adjacent to the brain transient at the next moment. The number of times a transient transitions into various brain transients at the next moment;基于该种脑瞬态下一时刻转移为各种脑瞬态的次数，以及每种脑瞬态出现的次数，确定每种脑瞬态转移为各种脑瞬态的概率；Based on the number of times this brain transient transitions to various brain transients at the next moment and the number of occurrences of each brain transient, determine the probability that each brain transient transitions to various brain transients;基于每种脑瞬态转移为各种脑瞬态的概率，确定所述多种脑瞬态之间的状态转移概率矩阵。Based on the probability that each brain transient transitions into various brain transients, a state transition probability matrix between the various brain transients is determined.4.根据权利要求1所述的设备，其特征在于，4. The device according to claim 1, characterized in that,在基于所述原始核磁共振图像，得到所述待监测对象的三维脑模型方面，所述离线处理模块，具体用于：In terms of obtaining a three-dimensional brain model of the object to be monitored based on the original MRI image, the offline processing module is specifically used to:对所述原始核磁共振图像进行重切片，得到目标核磁共振图像；Re-slice the original nuclear magnetic resonance image to obtain a target nuclear magnetic resonance image;对所述目标核磁共振图像进行分割，得到脑组织、颅骨以及头皮；Segment the target MRI image to obtain brain tissue, skull and scalp;基于分割出脑组织、颅骨以及头皮的所述目标核磁共振图像进行脑模型建立，得到所述待监测对象的三维脑模型。A brain model is established based on the target MRI image segmented into brain tissue, skull and scalp to obtain a three-dimensional brain model of the object to be monitored.5.一种电子设备，其特征在于，包括：处理器和存储器，所述处理器与所述存储器相连，所述存储器用于存储计算机程序，所述处理器用于执行所述存储器中存储的计算机程序，以使得所述电子设备执行如下步骤：5. An electronic device, characterized in that it includes: a processor and a memory, the processor is connected to the memory, the memory is used to store computer programs, and the processor is used to execute the computer program stored in the memory. Program, so that the electronic device performs the following steps:获取待监测对象的脑部的原始核磁共振图像；Obtaining raw MRI images of the brain of the subject to be monitored;在预设时间段内在多个通道上采集所述待监测对象的多个离线脑电信号；Collect multiple offline EEG signals of the object to be monitored on multiple channels within a preset time period;基于所述原始核磁共振图像，得到所述待监测对象的三维脑模型；Based on the original MRI image, obtain a three-dimensional brain model of the object to be monitored;获取所述三维脑模型中的多个皮层网格；将所述多个通道上用于脑电信号采集的多个电极在所述待监测对象的脑部的位置与所述多个皮层网格在脑部的位置进行对齐，得到与所述待监测对象对应的Lead-field矩阵；Obtain multiple cortical grids in the three-dimensional brain model; compare the positions of the multiple electrodes used for brain electrical signal collection on the multiple channels in the brain of the subject to be monitored with the multiple cortical grids. Align the position of the brain to obtain a Lead-field matrix corresponding to the object to be monitored;获取所述多个离线脑电信号的协方差矩阵；Obtain the covariance matrix of the multiple offline EEG signals;基于所述协方差矩阵以及所述Lead-field矩阵，确定重组矩阵；Based on the covariance matrix and the lead-field matrix, determine the reorganization matrix;基于所述Lead-field矩阵，对所述多个离线脑电信号进行逆分解，得到所述多个皮层网格对应的多个皮层信号；将所述三维脑模型划分为多个源极子；Based on the lead-field matrix, inversely decompose the multiple offline EEG signals to obtain multiple cortical signals corresponding to the multiple cortical grids; divide the three-dimensional brain model into multiple source sub-elements;对所述多个皮层信号进行加权映射，得到与所述多个源极子对应的多个第二脑电源信号；Perform weighted mapping on the plurality of cortical signals to obtain a plurality of second brain power signals corresponding to the plurality of source sub-elements;基于预设时间窗口，对每个第二脑电源信号进行多次加窗，以对每个第二脑电源信号进行分割，得到与每个第二脑电源信号对应的多个子脑电源信号；Based on the preset time window, perform multiple windowing on each second brain power signal to segment each second brain power signal and obtain multiple sub-brain power signals corresponding to each second brain power signal;基于每个第二脑电源信号对应的多个子脑电源信号，得到多个预设延时自协方差矩阵、多种脑瞬态之间的状态转移概率矩阵以及所述多种脑瞬态对应的初始状态概率，其中，每个预设延时自协方差矩阵用于表征一种脑瞬态；Based on multiple sub-brain power signals corresponding to each second brain power signal, multiple preset delay auto-covariance matrices, state transition probability matrices between multiple brain transient states, and corresponding values of the multiple brain transient states are obtained. Initial state probability, where each preset delay autocovariance matrix is used to characterize a brain transient;在当前时刻从所述多个通道上采集所述待监测对象的多个原始脑电信号；Collect multiple raw EEG signals of the object to be monitored from the multiple channels at the current moment;基于所述重组矩阵对所述多个原始脑电信号进行溯源分析，得到多个第一脑电源信号；Perform traceability analysis on the plurality of original EEG signals based on the recombination matrix to obtain a plurality of first brain power signals;从每个第一脑电源信号中截取与当前时刻对应的待分析脑电信号，得到多个待分析脑电信号；Intercept the EEG signal to be analyzed corresponding to the current moment from each first brain power signal to obtain multiple EEG signals to be analyzed;对所述多个待分析脑电信号进行处理，得到第一延时自协方差矩阵；Process the plurality of EEG signals to be analyzed to obtain a first delayed autocovariance matrix;确定所述第一延时自协方差矩阵与所述多个预设延时自协方差矩阵之间的相似度，得到多个第一相似度；Determine the similarity between the first delay autocovariance matrix and the plurality of preset delay autocovariance matrices, and obtain a plurality of first similarities;根据所述多个第一相似度、所述待监测对象在上一时刻的脑瞬态、所述状态转移概率矩阵以及所述初始状态概率，确定所述待监测对象当前时刻处于各种脑瞬态的概率；包括：According to the plurality of first similarities, the brain transient state of the object to be monitored at the previous moment, the state transition probability matrix and the initial state probability, it is determined that the object to be monitored is in various brain transient states at the current moment. The probability of the state; including:若当前时刻为首次确定所述待监测对象的脑瞬态时，基于所述多种脑瞬态对应的初始状态概率与所述多个第一相似度，确定所述待监测对象当前时刻处于各种脑瞬态的概率；If the current moment is the first time the brain transient state of the subject to be monitored is determined, based on the initial state probabilities corresponding to the multiple brain transient states and the plurality of first similarities, it is determined that the subject to be monitored is in each state at the current moment. the probability of brain transients;若当前时刻不是首次确定所述待监测对象的脑瞬态时，基于所述状态转移概率矩阵，确定所述待监测对象在上一时刻的脑瞬态转移为各种脑瞬态的概率序列；根据所述概率序列与所述多个第一相似度，确定所述待监测对象当前时刻处于各种脑瞬态的概率；If the current moment is not the first time that the brain transient state of the object to be monitored is determined, based on the state transition probability matrix, it is determined that the brain transient state of the object to be monitored at the previous moment is transformed into a probability sequence of various brain transient states; Determine the probability that the subject to be monitored is in various brain transient states at the current moment according to the probability sequence and the plurality of first similarities;根据当前时刻处于各种脑瞬态的概率，确定所述待监测对象在当前时刻的脑瞬态。According to the probability of being in various brain transient states at the current moment, the brain transient state of the object to be monitored at the current moment is determined.6.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行以实现如下步骤：6. A computer-readable storage medium, characterized in that the computer-readable storage medium stores a computer program, and the computer program is executed by a processor to implement the following steps:获取待监测对象的脑部的原始核磁共振图像；Obtaining raw MRI images of the brain of the subject to be monitored;在预设时间段内在多个通道上采集所述待监测对象的多个离线脑电信号；Collect multiple offline EEG signals of the object to be monitored on multiple channels within a preset time period;基于所述原始核磁共振图像，得到所述待监测对象的三维脑模型；Based on the original MRI image, obtain a three-dimensional brain model of the object to be monitored;获取所述三维脑模型中的多个皮层网格；将所述多个通道上用于脑电信号采集的多个电极在所述待监测对象的脑部的位置与所述多个皮层网格在脑部的位置进行对齐，得到与所述待监测对象对应的Lead-field矩阵；Obtain multiple cortical grids in the three-dimensional brain model; compare the positions of the multiple electrodes used for brain electrical signal collection on the multiple channels in the brain of the subject to be monitored with the multiple cortical grids. Align the position of the brain to obtain a Lead-field matrix corresponding to the object to be monitored;获取所述多个离线脑电信号的协方差矩阵；Obtain the covariance matrix of the multiple offline EEG signals;基于所述协方差矩阵以及所述Lead-field矩阵，确定重组矩阵；Based on the covariance matrix and the lead-field matrix, determine the reorganization matrix;基于所述Lead-field矩阵，对所述多个离线脑电信号进行逆分解，得到所述多个皮层网格对应的多个皮层信号；将所述三维脑模型划分为多个源极子；Based on the lead-field matrix, inversely decompose the multiple offline EEG signals to obtain multiple cortical signals corresponding to the multiple cortical grids; divide the three-dimensional brain model into multiple source sub-elements;对所述多个皮层信号进行加权映射，得到与所述多个源极子对应的多个第二脑电源信号；Perform weighted mapping on the plurality of cortical signals to obtain a plurality of second brain power signals corresponding to the plurality of source sub-elements;基于预设时间窗口，对每个第二脑电源信号进行多次加窗，以对每个第二脑电源信号进行分割，得到与每个第二脑电源信号对应的多个子脑电源信号；Based on the preset time window, perform multiple windowing on each second brain power signal to segment each second brain power signal and obtain multiple sub-brain power signals corresponding to each second brain power signal;基于每个第二脑电源信号对应的多个子脑电源信号，得到多个预设延时自协方差矩阵、多种脑瞬态之间的状态转移概率矩阵以及所述多种脑瞬态对应的初始状态概率，其中，每个预设延时自协方差矩阵用于表征一种脑瞬态；Based on multiple sub-brain power signals corresponding to each second brain power signal, multiple preset delay auto-covariance matrices, state transition probability matrices between multiple brain transient states, and corresponding values of the multiple brain transient states are obtained. Initial state probability, where each preset delay autocovariance matrix is used to characterize a brain transient;在当前时刻从所述多个通道上采集所述待监测对象的多个原始脑电信号；Collect multiple raw EEG signals of the object to be monitored from the multiple channels at the current moment;基于所述重组矩阵对所述多个原始脑电信号进行溯源分析，得到多个第一脑电源信号；Perform traceability analysis on the plurality of original EEG signals based on the recombination matrix to obtain a plurality of first brain power signals;从每个第一脑电源信号中截取与当前时刻对应的待分析脑电信号，得到多个待分析脑电信号；Intercept the EEG signal to be analyzed corresponding to the current moment from each first brain power signal to obtain multiple EEG signals to be analyzed;对所述多个待分析脑电信号进行处理，得到第一延时自协方差矩阵；Process the plurality of EEG signals to be analyzed to obtain a first delayed autocovariance matrix;确定所述第一延时自协方差矩阵与所述多个预设延时自协方差矩阵之间的相似度，得到多个第一相似度；Determine the similarity between the first delay autocovariance matrix and the plurality of preset delay autocovariance matrices, and obtain a plurality of first similarities;根据所述多个第一相似度、所述待监测对象在上一时刻的脑瞬态、所述状态转移概率矩阵以及所述初始状态概率，确定所述待监测对象当前时刻处于各种脑瞬态的概率；包括：According to the plurality of first similarities, the brain transient state of the object to be monitored at the previous moment, the state transition probability matrix and the initial state probability, it is determined that the object to be monitored is in various brain transient states at the current moment. The probability of the state; including:若当前时刻为首次确定所述待监测对象的脑瞬态时，基于所述多种脑瞬态对应的初始状态概率与所述多个第一相似度，确定所述待监测对象当前时刻处于各种脑瞬态的概率；If the current moment is the first time the brain transient state of the subject to be monitored is determined, based on the initial state probabilities corresponding to the multiple brain transient states and the plurality of first similarities, it is determined that the subject to be monitored is in each state at the current moment. the probability of brain transients;若当前时刻不是首次确定所述待监测对象的脑瞬态时，基于所述状态转移概率矩阵，确定所述待监测对象在上一时刻的脑瞬态转移为各种脑瞬态的概率序列；根据所述概率序列与所述多个第一相似度，确定所述待监测对象当前时刻处于各种脑瞬态的概率；If the current moment is not the first time that the brain transient state of the object to be monitored is determined, based on the state transition probability matrix, it is determined that the brain transient state of the object to be monitored at the previous moment is transformed into a probability sequence of various brain transient states; Determine the probability that the subject to be monitored is in various brain transient states at the current moment according to the probability sequence and the plurality of first similarities;根据当前时刻处于各种脑瞬态的概率，确定所述待监测对象在当前时刻的脑瞬态。According to the probability of being in various brain transient states at the current moment, the brain transient state of the object to be monitored at the current moment is determined.