CN114742116A - 一种模拟采集脑电信号的生成方法和系统 - Google Patents

一种模拟采集脑电信号的生成方法和系统 Download PDF

Info

Publication number
CN114742116A
CN114742116A CN202210659471.8A CN202210659471A CN114742116A CN 114742116 A CN114742116 A CN 114742116A CN 202210659471 A CN202210659471 A CN 202210659471A CN 114742116 A CN114742116 A CN 114742116A
Authority
CN
China
Prior art keywords
signal
electroencephalogram
noise signal
interest
noise
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202210659471.8A
Other languages
English (en)
Other versions
CN114742116B (zh
Inventor
张家伟
陈超
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sichuan Neosource Biotektronics Ltd
Original Assignee
Sichuan Neosource Biotektronics Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Sichuan Neosource Biotektronics Ltd filed Critical Sichuan Neosource Biotektronics Ltd
Priority to CN202210910794.XA priority Critical patent/CN115238745A/zh
Priority to CN202210659471.8A priority patent/CN114742116B/zh
Priority to CN202210910733.3A priority patent/CN115186718A/zh
Priority to CN202210666384.5A priority patent/CN115034266A/zh
Publication of CN114742116A publication Critical patent/CN114742116A/zh
Application granted granted Critical
Publication of CN114742116B publication Critical patent/CN114742116B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/30Input circuits therefor
    • A61B5/307Input circuits therefor specially adapted for particular uses
    • A61B5/31Input circuits therefor specially adapted for particular uses for electroencephalography [EEG]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/369Electroencephalography [EEG]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7203Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal

Landscapes

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

Abstract

本说明书实施例提供一种模拟采集脑电信号的生成方法和系统,该方法包括:基于基础信号和第一参数组,获得无噪音的模拟脑电信号,该模拟脑电信号用于作为人体中无噪音的脑电信号参与脑电信号的分析与研究;基于初始噪音信号和第二参数组,获得目标噪音信号,该目标噪音信号用于作为人体的采集脑电信号中的噪音信号参与该脑电信号的分析与研究;获取模拟脑电信号的第一权重以及目标噪音信号的第二权重;基于模拟脑电信号、目标噪音信号、第一权重以及第二权重,获得目标模拟采集脑电信号,该目标模拟采集脑电信号用于作为该采集脑电信号中参与该脑电信号的分析与研究以及减少真实脑电信号的采集。

Description

一种模拟采集脑电信号的生成方法和系统
技术领域
本说明书涉及模拟脑电信号领域,特别涉及一种模拟采集脑电信号的生成方法和系统。
背景技术
脑电信号是大脑皮层大量神经元的突触后电位总和的结果。脑电信号同步节律的形成与皮层丘脑非特异性投射系统的活动有关。由于脑电信号的信号幅值极其低(约1uV~100uV),在脑电信号的采集过程中容易引入各种干扰噪声,从而对脑电信号的波形质量造成影响,降低脑电信号的波形分析的准确性。当前可以利用机器学习模型对所采集的脑电信号进行去噪处理,但是由于采集时长和采集手段限制,较难获取模型训练所需的大量训练数据。
因此,有必要提出一种模拟采集脑电信号的生成方法,使得模拟采集脑电信号尽可能接近真实的采集脑电信号,从而可以获得大量模型训练所需的训练数据,从而使训练好的去噪模型有更好的泛化效果,也可以在研究脑电信号时降低采集人体的脑电信号所需的成本。
发明内容
本说明书一个或多个实施例提供一种模拟采集脑电信号的生成方法,所述方法包括基于基础信号和第一参数组,获得无噪音的模拟脑电信号,所述模拟脑电信号用于作为人体中无噪音的脑电信号参与脑电信号的分析与研究;基于初始噪音信号和第二参数组,获得目标噪音信号,所述目标噪音信号用于作为人体的采集脑电信号中的噪音信号参与所述脑电信号的分析与研究;获取所述模拟脑电信号的第一权重以及所述目标噪音信号的第二权重;基于所述模拟脑电信号、所述目标噪音信号、所述第一权重以及所述第二权重,获得目标模拟采集脑电信号,所述目标模拟采集脑电信号用于减少真实脑电信号的采集,所述目标模拟采集脑电信号用于作为所述采集脑电信号参与所述脑电信号的分析与研究以及减少真实脑电信号的采集。
本说明书一个或多个实施例提供一种模拟采集脑电信号的生成系统,所述系统包括:第一获取模块,用于基于基础信号和第一参数组,获得无噪音的模拟脑电信号,所述模拟脑电信号用于作为人体中无噪音的脑电信号参与脑电信号的分析与研究;第二获取模块,用于基于初始噪音信号和第二参数组,获得目标噪音信号,所述目标噪音信号用于作为人体的采集脑电信号中的噪音信号参与所述脑电信号的分析与研究;第三获取模块,用于获取所述模拟脑电信号的第一权重以及所述目标噪音信号的第二权重;第四获取模块,用于基于所述模拟脑电信号、所述目标噪音信号、所述第一权重以及所述第二权重,获得目标模拟采集脑电信号,所述目标模拟采集脑电信号用于作为所述采集脑电信号参与所述脑电信号的分析与研究以及减少真实脑电信号的采集。
本说明书一个或多个实施例提供一种模拟采集脑电信号的生成装置,包括处理器,所述处理器用于执行以实现如上述实施例中任一项所述的模拟采集脑电信号的生成方法。
本说明书一个或多个实施例提供一种计算机可读存储介质,所述存储介质存储计算机指令,当计算机读取存储介质中的计算机指令后,计算机执行以实现如上述实施例中任一项所述的模拟采集脑电信号的生成方法。
附图说明
本说明书将以示例性实施例的方式进一步说明,这些示例性实施例将通过附图进行详细描述。这些实施例并非限制性的,在这些实施例中,相同的编号表示相同的结构,其中:
图1是根据本说明书一些实施例所示的模拟采集脑电信号的生成系统的应用场景示意图;
图2是根据本说明书一些实施例所示的模拟采集脑电信号的生成系统的示例性模块图;
图3是根据本说明书一些实施例所示的生成模拟采集脑电信号的示意图;
图4是根据本说明书一些实施例所示的获得无噪音的模拟脑电信号的示例性流程图;
图5是根据本说明书一些实施例所示的获得目标噪音信号的示例性流程图;
图6是根据本说明书一些实施例所示的获得第一参数组的示例性流程图;
图7是根据本说明书一些实施例所示的采集脑电信号去噪方法的示例性流程图;
图8是根据本说明书一些实施例所示的获取信号转移关系的示例性流程图;
图9是根据本说明书一些实施例所示的信号转移模型的示意图;
图10是根据本说明书一些实施例所示的又一信号转移模型的示意图;
图11是根据本说明书一些实施例所示的确定目标噪音信号的示例性流程图。
具体实施方式
为了更清楚地说明本说明书实施例的技术方案,下面将对实施例描述中所需要使用的附图作简单的介绍。显而易见地,下面描述中的附图仅仅是本说明书的一些示例或实施例,对于本领域的普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图将本说明书应用于其它类似情景。除非从语言环境中显而易见或另做说明,图中相同标号代表相同结构或操作。
应当理解,本文使用的“系统”、“装置”、“单元”和/或“模块”是用于区分不同级别的不同组件、元件、部件、部分或装配的一种方法。然而,如果其他词语可实现相同的目的,则可通过其他表达来替换所述词语。
如本说明书和权利要求书中所示,除非上下文明确提示例外情形,“一”、“一个”、“一种”和/或“该”等词并非特指单数,也可包括复数。一般说来,术语“包括”与“包含”仅提示包括已明确标识的步骤和元素,而这些步骤和元素不构成一个排它性的罗列,方法或者设备也可能包含其它的步骤或元素。
本说明书中使用了流程图用来说明根据本说明书的实施例的系统所执行的操作。应当理解的是,前面或后面操作不一定按照顺序来精确地执行。相反,可以按照倒序或同时处理各个步骤。同时,也可以将其他操作添加到这些过程中,或从这些过程移除某一步或数步操作。
图1是根据本说明书一些实施例所示的模拟采集脑电信号的生成系统的应用场景示意图。
如图1所示,模拟采集脑电信号的生成系统的应用场景100可以包括服务器110、基础信号120、初始噪音信号130、网络140、存储设备150以及终端设备160。模拟采集脑电信号的生成系统的应用场景可以通过实施本说明书中披露的方法和/或过程获取模拟采集脑电信号。
服务器110可以与基础信号120、初始噪音信号130、网络140、存储设备150以及终端设备160进行通信,以实现模拟采集脑电信号的生成系统的各种功能。在一些实施例中,服务器110可以经由,例如,网络140,接收基础信号120和初始噪音信号130,并可以对其进行处理。在一些实施例中,服务器110可以经由,例如,网络140,输出相关数据至存储设备150以及终端设备160。在一些实施例中,服务器110可以是单个服务器,也可以是服务器组。在一些实施例中,服务器110可以本地连接到网络140或者与网络140远程连接。在一些实施例中,服务器110可以在云平台上实施。
基础信号120可以为简单信号。例如,基础信号120可以正弦波。基础信号120可以由波形发生器、模拟信号发生器、波形生成器等任意一种或多种生成。在一些实施例中,基础信号120可以通过网络140发送至服务器110、存储设备150以及终端设备160。
初始噪音信号130可以为采集脑电信号时未经过传播的噪音信号。初始噪音信号130可以通过噪音发生器、模拟信号发生器、噪音生成器等任意一种或多种生成。在一些实施例中,初始噪音信号130可以通过网络140发送至服务器110、存储设备150以及终端设备160。
网络140可以用于信息和/或数据的传输。在一些实施例中,应用场景100中的一个或以上组件(例如,服务器110和/或存储设备150等)可以经由网络140将信息和/或数据发送到应用场景100中的另一个组件。
存储设备150可以储存数据和/或指令。数据可以包括与服务器110、终端设备160、基础信号120和初始噪音信号130等有关的数据。例如,存储设备150可以存储基础信号120和初始噪音信号130。在一些实施例中,存储设备150可以储存服务器110用来执行或使用以完成本说明书中描述的示例性方法的数据和/或指令。在一些实施例中,存储设备150可以连接到网络140以与应用场景100的一个或以上组件(例如,服务器110和/或终端设备160)通信。在一些实施例中,存储设备150可以是服务器110的一部分。在一些实施例中,存储设备150可包括大容量存储器、可移动存储器、易失性读写内存、只读内存(ROM)等或其任意组合。在一些实施例中,存储设备150可在云平台上实现。
终端设备160可以指用户所使用的一个或多个终端或软件。其中,用户可以指研究员。用户可以通过终端设备160查看应用场景100中各个组件的相关数据。在一些实施例中,终端设备160可以包括移动设备160-1、平板计算机160-2、膝上型计算机160-3等或其任意组合。在一些实施例中,终端设备160可以是固定的和/或移动的。例如,终端设备160可以直接安装在服务器110上,成为服务器110的一部分。再例如,终端设备160可以是可移动的设备,用户可以携带终端设备160位于相对于服务器110、基础信号120和初始噪音信号130距离较远的位置,终端设备160可以通过网络140与服务器110、存储设备150连接和/或通信。
应当注意的是,应用场景100仅仅是为了说明的目的而提供,并不意图限制本说明书的范围。对于本领域的普通技术人员来说,可以根据本说明书的描述,做出多种修改或变化。例如,应用场景100还可以包括数据库。然而,这些变化和修改不会背离本说明书的范围。
图2是根据本说明书一些实施例所示的模拟采集脑电信号的生成系统的示例性模块图。
如图2所示,模拟采集脑电信号的生成系统200可以包括第一获取模块210、第二获取模块220、第三获取模块230及第四获取模块240。
第一获取模块210可以用于基于基础信号和第一参数组,获得无噪音的模拟脑电信号,模拟脑电信号用于作为人体中无噪音的脑电信号参与脑电信号的分析与研究。关于基础信号、第一参数组以及模拟脑电信号的更多内容参见图3及其相关描述。在一些实施例中,第一获取模块210还可以用于获取初始第一参数组;基于初始第一参数组与基础信号,获得中间模拟脑电信号;获取参考脑电信号;基于中间模拟脑电信号以及参考脑电信号构建损失函数;基于损失函数迭代更新初始第一参数组的参数,直至满足预设条件,获得第一参数组。关于初始第一参数组、中间模拟脑电信号以及参考脑电信号的更多内容参见图6及其相关描述。在一些实施例中,第一参数组包括感兴趣区域转移矩阵以及信号源转移矩阵,感兴趣区域转移矩阵中元素表征信号在不同感兴趣区域之间进行传播时的转移关系,信号源转移矩阵中的元素表征同一感兴趣区域中的信号源与电极点之间进行传播时的转移关系,其中,不同感兴趣区域表示人体大脑中的不同脑功能区,第一获取模块210还可以用于针对虚拟脑部模型中的每一信号源,基于感兴趣区域转移矩阵以及该信号源发出的基础信号,确定该信号源的目标基础信号;针对虚拟脑部模型中每一感兴趣区域,基于该感兴趣区域中每一信号源的目标基础信号以及信号源转移矩阵,确定该感兴趣区域对应的初始子模拟脑电信号;基于各个感兴趣区域对应的初始子模拟脑电信号,获得模拟脑电信号。关于感兴趣区域转移矩阵、信号源转移矩阵、目标基础信号以及初始子模拟脑电信号的更多内容参见图4及其相关描述。
第二获取模块220可以用于基于初始噪音信号和第二参数组,获得目标噪音信号,目标噪音信号用于作为人体的采集脑电信号中的噪音信号参与脑电信号的分析与研究。关于初始噪音信号、第二参数组以及目标噪音信号的更多内容参见图3及其相关描述。在一些实施例中,初始噪音信号包括模拟脑部噪音信号以及电极点噪音信号,第二参数组包括感兴趣区域转移矩阵以及电极点转移矩阵,感兴趣区域转移矩阵中元素表征信号在不同感兴趣区域之间进行传播时的转移关系,电极点转移矩阵中的元素表征信号在不同电极点之间进行传播时的转移关系,第二获取模块220还可以用于基于感兴趣区域转移矩阵以及模拟脑部噪音信号,确定第一子噪音信号;基于电极点转移矩阵以及电极点噪音信号,确定第二子噪音信号;基于第一子噪音信号以及第二子噪音信号,获得目标噪音信号。关于模拟脑部噪音信号、电极点噪音信号、电极点转移矩阵、第一子噪音信号以及第二子噪音信号的更多内容参见图5及其相关描述。
在一些实施例中,第二获取模块220可以包括第一确定模块、第二确定模块以及信号处理模块等子模块。第一确定模块用于确定被测对象的信号转移关系。其中,信号转移关系表征初始噪音信号在被测对象的各个位置之间进行传播时的信号变化关系。关于初始噪音信号、被测对象以及信号转移关系的更多内容可以参见图7及其相关描述。在一些实施例中,第一确定模块还可以用于针对被测对象的每一位置,从该位置输入初始测试信号,其中,初始测试信号的数量级大于被测对象的脑电信号的数量级;从被测对象的其他位置采集信号,获得采集测试信号;基于初始测试信号以及采集测试信号,确定信号转移关系。关于初始测试信号以及采集测试信号的更多内容可以参见图8及其相关描述。第二确定模块用于基于初始噪音信号以及信号转移关系,确定目标噪音信号。关于目标噪音信号更多内容可以参见图7及其相关描述。在一些实施例中,信号转移关系包括多个简单信号转移关系,简单信号转移关系表征第一简单噪声信号在被测对象的各个位置之间进行传播时的信号变化关系,其中,第一简单噪声信号是对初始噪音信号进行变换而获得的,第二确定模块还可以用于将初始噪音信号变换为多个不同频率和强度的第一简单噪声信号;针对每个所述第一简单噪声信号,基于该第一简单噪声信号对应的简单信号转移关系对该第一简单噪声信号进行处理,确定该第一简单噪声信号对应的第二简单噪音信号;合成各个第二简单噪音信号,获得目标噪音信号。关于简单信号转移关系、第一简单噪声信号以及第二简单噪音信号的更多内容可以参见图10及其相关描述。信号处理模块用于从被测对象中采集到的采集脑电信号去除目标噪音信号,获得无噪音的目标脑电信号。关于采集脑电信号和目标脑电信号更多内容可以参见图7及其相关描述。
第三获取模块230可以用于获取模拟脑电信号的第一权重以及目标噪音信号的第二权重。关于第一权重、第二权重的更多内容参见图3及其相关描述。
第四获取模块240可以用于基于模拟脑电信号、目标噪音信号、第一权重以及第二权重,获得目标模拟采集脑电信号,目标模拟采集脑电信号用于作为采集脑电信号参与脑电信号的分析与研究以及减少真实脑电信号的采集。关于目标模拟采集脑电信号的更多内容参见图3及其相关描述。
需要注意的是,以上对于模拟采集脑电信号的生成系统200及其模块的描述,仅为描述方便,并不能把本说明书限制在所举实施例范围之内。可以理解,对于本领域的技术人员来说,在了解该系统的原理后,可能在不背离这一原理的情况下,对各个模块进行任意组合,或者构成子系统与其他模块连接。在一些实施例中,图2中披露的第一获取模块210、第二获取模块220、第三获取模块230和第四获取模块240可以是一个系统中的不同模块,也可以是一个模块实现上述的两个或两个以上模块的功能。例如,各个模块可以共用一个存储模块,各个模块也可以分别具有各自的存储模块。诸如此类的变形,均在本说明书的保护范围之内。
图3是根据本说明书一些实施例所示的生成模拟采集脑电信号的示意图。在一些实施例中,流程300可以由服务器110执行。如图3所示,流程300可以包括以下步骤:
步骤310,基于基础信号和第一参数组,获得无噪音的模拟脑电信号。在一些实施例中,步骤310可以由第一获取模块210执行。
基础信号可以指信号变化规律的简单信号。模拟脑电信号可以用于作为人体中无噪音的脑电信号参与脑电信号的分析与研究。基础信号可以包括但不限于正弦信号、余弦信号以及方波信号等。第一获取模块210可以将基础信号的特征输入模拟信号发生器以获取基础信号。基础信号的特征可以包括但不限于基础信号的频率、幅值等,其大小可以通过预先设定,以使得基于基础信息获取的模拟脑电信号符合真实的脑电信号的基础特征。例如,可以预先设定基础信号的频率小于100Hz,幅值的取值范围为-200uV~+200 uV。
第一参数组可以指表征基础信号在传播时,虚拟脑部模型中的感兴趣区域、信号源以及安装在虚拟脑部模型上的电极点等因素对基础信号的产生的影响的参数组。在一些实施例中,可以通过多个不同的第一参数组对基础信号进行处理以获得多个不同的模拟脑电信号。其中,虚拟脑部模型中的不同感兴趣区域可以分别表示人体大脑中的各个脑功能区。例如,虚拟脑部模型中的某一感兴趣区域可以表示人体大脑前部负责专注力、短时记忆任务等相关任务的额叶。
应当理解的是,信号在人体大脑中进行传播时,不同感兴趣区域的大脑结构、不同信号源以及用于采集脑电信号的电极点均会对信号产生影响,使得信号发生变化。本说明书的一些实施例中可以通过第一参数组确定感兴趣区域、信号源以及电极点等对信号的影响,以使得获得的模拟脑电信号更接近于真实的无噪音的脑电信号。
在一些实施例中,第一参数组可以包括信号在虚拟脑部模型的不同感兴趣区域之间进行传播时的变化关系。第一参数组可以包括感兴趣区域转移矩阵,该感兴趣区域转移矩阵可以表示基础信号在虚拟脑部模型的不同感兴趣区域之间进行传播时的变化关系。关于感兴趣区域转移矩阵的更多内容参见图4及其相关描述。
在一些实施例中,第一参数组还可以包括虚拟脑部模型中某一感兴趣区域中的某一信号源发出的信号在该感兴趣区域中的电极点间进行传播时的变化关系。第一参数组还可以包括信号源转移矩阵,该信号源转移矩阵可以表示虚拟脑部模型中某一感兴趣区域中的某一信号源发出的基础信号在该感兴趣区域中的电极点之间进行传播时的变化关系。关于信号源转移矩阵的更多内容参见图4及其相关描述。
在一些实施例中,第一参数组中的元素可以是确定值,也可以为区间范围值。当第一参数组为区间范围值时,在基于第一参数组对基础信号进行处理时,可以从该区间范围值内随机选择某一具体值以对基础信号进行处理,从而获取模拟脑电信号。不同人体或者不同虚拟脑部模型对应的第一参数组不同。
在一些实施例中,可以通过多种方式获取第一参数组。例如,可以通过网络获取第一参数组。再例如,还可以通过随机初始化获取第一参数组。
在一些实施例中,第一参数组基于第一预设参数区间获得。可以通过设置第一参数组的第一预设参数区间,并迭代更新第一预设参数区间的范围从而获得第一参数组。其中,第一预设参数区间可以指第一参数组中各元素的预设区间范围。第一预设参数区间可以通过多种方式获取,例如,可以通过随机选取初始范围确定。再例如,可以基于现有第一参数组的范围,并对该范围进行扩展确定。第一预设参数区间对应参数组可以为初始第一参数组。在一些实施例中,可以基于初始第一参数组与基础信号,获得中间模拟脑电信号;基于中间模拟脑电信号以及参考脑电信号构建损失函数;基于损失函数迭代更新初始第一参数组,直至满足预设条件,获得第一参数组。关于上述获得第一参数组的实施例的更多内容参见图6及其相关描述。
本说明书的一些实施例中通过设置第一参数组的第一预设参数区间可以确定第一参数组的大致范围,减小确定第一参数组的计算量,更快确定第一参数组。
模拟脑电信号可以指模拟的人体大脑的脑电信号的信号。在一些实施例中,可以基于第一参数组对基础信号进行处理,确定采集该基础信号时,虚拟脑部模型中的感兴趣区域、信号源以及安装在虚拟脑部模型上的电极点等对该基础信号产生的影响,以获得模拟脑电信号。获取无噪音的模拟脑电信号。在一些实施例中,针对虚拟脑部模型中的每一感兴趣区域,基于信号源转移矩阵以及基础信号,确定该感兴趣区域对应的初始子模拟脑电信号;基于感兴趣区域转移矩阵以及各个感兴趣区域对应的初始子模拟脑电信号,获得模拟脑电信号。关于上述实施例的更多内容参见图4及其相关描述。
在一些实施例中,通过在虚拟脑部模型的多个信号源发出基础信号,多种基础信号进行叠加,可以使得获取的模拟脑电信号具有脑电信号的基础特征。其中,脑电信号的基础特征可以包括但不限于脑电信号的幅值特征(如幅值范围)、频率特征(如频率范围)等。例如,脑电信号的基础特征可以包括幅值范围为-200uV~+200 uV,频率小于100Hz。通过虚拟脑部模型中的3个信号源可以发出幅值分别为10uV、30uV、100uV,频率分别为5Hz、20Hz、40Hz的3个正弦信号,前述3个正弦信号进行叠加后可以得到一个更加复杂的模拟脑电信号,该模拟脑电信号的幅值小于100uV、频率有3个成分且均小于100 Hz,故该模拟脑电信号可以前述脑电信号的基础特征。
应当理解的是,模拟脑电信号是通过第一参数组对基础信号进行处理而直接获取的,因此,该模拟脑电信号中不含有对人体脑电信号进行采集时,采集设备或人体其他部位的活动所带来的噪音信号。
步骤320,基于初始噪音信号和第二参数组,获得目标噪音信号。在一些实施例中,步骤320可以由第二获取模块220执行。
初始噪音信号可以指采集脑电信号时未经过传播的噪音信号。目标噪音信号可以用于作为人体的采集脑电信号中的噪音信号参与脑电信号的分析与研究。
在采集脑电信号时,存在多种因素对采集产生影响,使得得到的采集脑电信号中包含有多种噪音信号。对应的,初始噪音信号可以包括多种噪音信号。在一些实施例中,初始噪音信号可以通过采集或模拟确定。关于确定初始噪音信号的更多内容参见图5及其相关描述。
第二参数组可以指表征初始噪音信号在传播时,虚拟脑部模型中的感兴趣区域、信号源以及安装在虚拟脑部模型上的电极点等因素对初始噪音信号的产生的影响的参数组。
与第一参数组类似的,第二参数组也可以包括感兴趣区域转移矩阵。在一些实施例中,第二参数组还可以包括电极点转移矩阵,该电极点转移矩阵可以表示信号在虚拟脑部模型中电极点之间进行传播时的变化关系。关于电极点转移矩阵的更多内容参见图4及其相关描述。
在一些实施例中,可以通过多种方式获取第二参数组。例如,可以通过网络获取感兴趣区域转移矩阵。
在一些实施例中,第二参数组基于第二预设参数区间获得。可以通过设置第二参数组的第二预设参数区间,并迭代更新第二预设参数区间的区间范围从而获得第二参数组。其中,第二预设参数区间可以指第二参数值中元素的预设区间范围,第二预设参数区间可以通过多种方式获取,例如,可以通过随机选取初始范围确定。第二预设参数区间对应参数组可以为初始第二参数组。在一些实施例中,可以基于初始第二参数组确定第二参数组。关于基于初始第二参数组确定第二参数组的更多内容可以参见图6及其相关说明。
在一些实施例中,还可以通过针对被测对象的每一电极点,从该电极点处输入初始测试信号,其中,初始测试信号的数量级大于被测对象的脑电信号的数量级;从被测对象的其他电极点处采集信号,获得采集测试信号;基于初始测试信号以及采集测试信号,确定该电极点与其他电极点之间的信号转移关系;基于各个电极点与其他电极点之间的信号转移关系,确定电极点转移矩阵。关于上述确定电极点转移矩阵的实施例的更多内容参见图8及其相关描述。
目标噪音信号可以指模拟的采集脑电信号中的噪音信号。在一些实施例中,可以基于第二参数组对初始噪音信号进行处理,确定采集该初始噪音信号时,虚拟脑部模型中的感兴趣区域、信号源以及安装在虚拟脑部模型上的电极点等对该初始噪音信号产生的影响,以获得目标噪音信号。在一些实施例中,基于感兴趣区域转移矩阵以及模拟脑部噪音信号,确定第一子噪音信号;基于电极点转移矩阵以及电极点噪音信号,确定第二子噪音信号;基于第一子噪音信号以及第二子噪音信号,获得目标噪音信号。关于上述实施例的更多内容参见图5及其相关描述。
步骤330,获取模拟脑电信号的第一权重以及目标噪音信号的第二权重。在一些实施例中,步骤330可以由第三获取模块230执行。
第一权重可以指模拟脑电信号在目标模拟采集脑电信号中的比重。第二权重可以指目标噪音信号在目标模拟采集脑电信号中的比重。其中,第一权重以及第二权重可以通过分析真实采集脑电信号中脑电信号与噪音信号的比值而预先确定。
步骤340,基于模拟脑电信号、目标噪音信号、第一权重以及第二权重,获得目标模拟采集脑电信号。
目标模拟采集脑电信号可以指模拟的采集脑电信号。目标模拟采集脑电信号可以用于作为采集脑电信号参与脑电信号的分析与研究以及减少真实脑电信号的采集。采集脑电信号中包含有脑电信号以及噪音信号,对应的,目标模拟采集脑电信号中也可以包含有模拟脑电信号以及目标噪音信号。
在一些实施例中,可以通过公式(1)确定目标模拟采集脑电信号:
Figure 469631DEST_PATH_IMAGE001
其中,x(t)表示目标模拟采集脑电信号,
Figure 779914DEST_PATH_IMAGE002
表示模拟脑电信号,
Figure 488546DEST_PATH_IMAGE003
表示目标噪音信号,
Figure 824368DEST_PATH_IMAGE004
分别表示第一权重与第二权重。
本说明书的一些实施例通过设置第一参数组与第二参数组,使得可以确定虚拟脑部模型中的感兴趣区域、信号源以及安装在虚拟脑部模型上的电极点等因素对信号的影响,进一步可以获得更加真实的目标模拟采集脑电信号。本说明书的一些实施例获得模拟脑电信号、目标噪音信号以及目标采集脑电信号均可以用于脑电信号的分析与研究,避免对脑电信号的分析与研究时缺乏研究数据的问题。同时,目标采集脑电信号还可以用于减少真实脑电信号的采集,降低关于脑电信号的分析与研究中真实脑电信号的采集成本。
图4是根据本说明书一些实施例所示的获得无噪音的模拟脑电信号的示例性流程图。在一些实施例中,流程400可以由第一获取模块210执行。如图4所示,流程400可以包括以下步骤:
步骤410,针对虚拟脑部模型中的每一信号源,基于感兴趣区域转移矩阵以及该信号源发出的基础信号,确定该信号源的目标基础信号。
目标基础信号可以为某一信号源发出的基础信号受到其他感兴趣区域的影响后的信号。
在一些实施例中,第一参数组可以包括感兴趣区域转移矩阵。感兴趣区域转移矩阵中的元素可以表征信号在不同感兴趣区域之间进行传播时的转移关系。感兴趣区域转移矩阵可以为n*n的矩阵,其中,n表示感兴趣区域的数量。在一些实施例中,感兴趣区域转移矩阵中的元素可以为实数,该实数为对应的感兴趣区域之间的调整参数,可以用于确定信号在对应感兴趣区域之间进行传播时的转移关系。例如,感兴趣区域转移矩阵中的某一元素可以为实数
Figure 591072DEST_PATH_IMAGE005
,目标基础信号可以基于如下公式获取:
Figure 213201DEST_PATH_IMAGE006
其中,
Figure 612391DEST_PATH_IMAGE007
表示第j个感兴趣区域中的第i个信号源的目标基础信号,
Figure 610959DEST_PATH_IMAGE007
表示该信号源发出的基础信号,
Figure 497748DEST_PATH_IMAGE005
为信号从第
Figure 946571DEST_PATH_IMAGE008
个感兴趣区域传播至第j个感兴趣区域的调整参数,
Figure 629795DEST_PATH_IMAGE007
表示第
Figure 759949DEST_PATH_IMAGE008
个感兴趣区域中的第i个信号源发出的基础信号。
在一些实施例中,感兴趣区域转移矩阵中的元素还可以为传导函数,该传导函数表征对应的感兴趣区域对信号的影响。前述传导函数可以为线性函数,也可以为非线性函数。关于感兴趣区域转移矩阵的获取方式可以参见图6及其相关描述。
在一些实施例中,目标基础信号可以基于如下公式获取:
Figure 879740DEST_PATH_IMAGE009
其中,
Figure 871153DEST_PATH_IMAGE007
表示第j个感兴趣区域中的第i个信号源的目标基础信号,
Figure 833818DEST_PATH_IMAGE007
表示该信号源发出的基础信号,
Figure 428048DEST_PATH_IMAGE010
表示其他感兴趣区域对该信号源的影响,T表示感兴趣区域转移矩阵,
Figure 526979DEST_PATH_IMAGE011
表示感兴趣区域转移矩阵中的第
Figure 519204DEST_PATH_IMAGE008
行,
Figure 396593DEST_PATH_IMAGE012
可以为第
Figure 275075DEST_PATH_IMAGE008
个感兴趣区域的m个信号源发出的信号的均值,
Figure 584439DEST_PATH_IMAGE013
表示第
Figure 890174DEST_PATH_IMAGE008
个感兴趣区域中的第i个信号源发出的基础信号。
步骤420,针对虚拟脑部模型中每一感兴趣区域,基于该感兴趣区域中每一信号源的目标基础信号以及信号源转移矩阵,确定该感兴趣区域对应的初始子模拟脑电信号。
初始子模拟脑电信号可以为虚拟脑部模型中某一感兴趣区域发出的信号。初始子模拟脑电信号可以为该感兴趣区域中多个信号源的目标基础信号在受到电极点的影响下进行叠加形成的信号。
在一些实施例中,第一参数组还可以包括信号源转移矩阵。信号源转移矩阵中的元素可以表征同一感兴趣区域中的信号源与电极点之间进行传播时的转移关系。信号源转移矩阵可以为m*k的矩阵,其中,m表示该感兴趣区域中信号源的数量,k表示虚拟脑部模型中电极点的数量。与感兴趣区域转移矩阵类似的,信号源转移矩阵中的元素可以为实数,该实数为对应的信号源与电极点之间的调整参数,可以用于确定信号在对应的信号源与电极点之间进行传播时的转移关系。信号源转移矩阵中的元素也可以为传导函数,该传导函数可以表征对应的电极点对对应的信号源的目标基础信号的影响。关于信号源转移矩阵的获取方式可以参见图6及其相关描述。
在一些实施例中,初始子模拟脑电信号可以基于如下公式获取:
Figure 176231DEST_PATH_IMAGE014
其中,
Figure 937032DEST_PATH_IMAGE015
表示第j个感兴趣区域的初始子模拟脑电信号,
Figure 773006DEST_PATH_IMAGE016
表示第j个感兴趣区域中的第i个信号源的目标基础信号,A表示信号源转移矩阵,
Figure 328271DEST_PATH_IMAGE017
表示信号源转移矩阵中的第i行,m表示该感兴趣区域中信号源的数量。
步骤430,基于各个感兴趣区域对应的初始子模拟脑电信号,获得模拟脑电信号。
在一些实施例中,可以基于各个感兴趣区域对应的初始子模拟脑电信号进行叠加处理,获取模拟脑电信号。模拟脑电信号可以基于以下公式获取:
Figure 350892DEST_PATH_IMAGE018
其中,
Figure 446543DEST_PATH_IMAGE019
表示模拟脑电信号,
Figure 917449DEST_PATH_IMAGE015
表示第j个感兴趣区域的初始子模拟脑电信号,n表示感兴趣区域的数量。
图5是根据本说明书一些实施例所示的获得目标噪音信号的示例性流程图。在一些实施例中,流程500可以由第二获取模块220执行。如图5所示,流程500可以包括以下步骤:
步骤510,基于感兴趣区域转移矩阵以及模拟脑部噪音信号,确定第一子噪音信号。
初始噪音信号可以包括模拟脑部噪音信号以及电极点噪音信号。其中,模拟脑部噪音信号可以指模拟的采集到的人体脑电信号中人体大脑中的噪音信号,电极点噪音信号可以指采集到的人体脑电信号中电极点带来的噪音信号。模拟脑部噪音信号可以包括背景噪音信号以及粉红噪音信号。背景噪声信号以及粉红噪音信号可以通过历史数据或网络获取。
在一些实施例中,第二参数组也可以包括感兴趣区域转移矩阵。
第一子噪音信号可以为基于背景噪音信号以及粉红噪音信号确定的噪音信号。第一子噪音信号可以基于以下公式确定:
Figure 283096DEST_PATH_IMAGE020
其中,
Figure 178633DEST_PATH_IMAGE021
表示第一子噪音信号中背景噪音信号带来的噪音信号,A表示信号源转移矩阵,
Figure 66256DEST_PATH_IMAGE022
表示信号源转移矩阵中的第i行,m表示该感兴趣区域中信号源的数量,n表示虚拟脑部模型中感兴趣区域的数量,
Figure 740469DEST_PATH_IMAGE023
表示第j个感兴趣区域中的第i个信号源发出的背景噪音信号;
Figure 757578DEST_PATH_IMAGE024
表示第一子噪音信号中粉红噪音信号带来的噪音信号,f表示该粉红噪音信号的频率,
Figure 459518DEST_PATH_IMAGE025
表示第j个感兴趣区域中的第i个信号源发出的粉红噪音信号;
Figure 37916DEST_PATH_IMAGE026
表示第一子噪音信号,
Figure 889940DEST_PATH_IMAGE027
Figure 692329DEST_PATH_IMAGE028
分别表示第一子噪音信号中背景噪音信号以及粉红噪音信号带来的噪音信号对应的权重参数,可以根据预先设置确定。
步骤520,基于电极点转移矩阵以及电极点噪音信号,确定第二子噪音信号。
在一些实施例中,第二参数组可以包括电极点转移矩阵。电极点转移矩阵中的元素可以表征信号在不同电极点之间进行传播时的转移关系。电极点转移矩阵可以为k*k的矩阵,其中,k为虚拟脑部模型中电极点的数量。与感兴趣区域转移矩阵类似的,电极点转移矩阵中的元素可以为实数,该实数为对应的电极点之间的调整参数,可以用于确定信号在对应电极点之间进行传播时的转移关系。电极点转移矩阵中的元素也可以为传导函数,该传导函数可以表征对应的电极点对电极点噪音信号的影响。
矩阵的中每一个元素可以为一个传导函数,该传导函数可以表征对应电机点对电极点产生的电极点噪音信号的影响。应当理解的是,位于电极点转移矩阵对角线位置的元素可以为1,表示该电极点对自身不产生影响。关于获取电极点转移矩阵的更多内容可以参见图8及其相关描述。
第二子噪音信号可以为基于电极点噪音信号确定的噪音信号。在一些实施例中,第二子噪音信号可以根据以下公式确定:
Figure 976506DEST_PATH_IMAGE029
其中,
Figure 652205DEST_PATH_IMAGE030
表示第二子噪音信号中电极点噪音信号带来的噪音信号,
Figure 691224DEST_PATH_IMAGE031
表示第p个电极点的电极点噪音信号,H表示电极点转移矩阵,
Figure 839747DEST_PATH_IMAGE032
表示电极点转移矩阵中的第p行,k表示电极点的数量;
Figure 433058DEST_PATH_IMAGE033
表示第二子噪音信号,
Figure 754641DEST_PATH_IMAGE034
表示第二子噪音信号中第二子噪音信号中电极点噪音信号带来的噪音信号的权重参数,可以通过预先设置确定。
步骤530,基于第一子噪音信号以及第二子噪音信号,获得目标噪音信号。
在一些实施例中,目标噪音信号可以根据以下公式确定:
Figure 382587DEST_PATH_IMAGE035
其中,
Figure 515061DEST_PATH_IMAGE036
表示目标噪音信号,
Figure 861247DEST_PATH_IMAGE026
表示第一子噪音信号,
Figure 846126DEST_PATH_IMAGE033
表示第二子噪音信号。
本说明书的一些实施例通过设置感兴趣区域转移矩阵、信号源转移矩阵以及电极点转移矩阵,可以明确信号源、感兴趣区域以及电极点等因素对信号的影响,使得获取到的模拟脑电信号以及目标噪音信号更加真实,更符合实际采集到的采集脑电信号。
值得说明的是,在现代技术中,可以基于机器学习模型,即去噪模型,对采集脑电信号进行去噪处理,得到无噪音的脑电信号。此时,需要大量的训练数据对去噪模型进行训练时,训练数据包括采集脑电信号以及去噪后的采集脑电信号。但实际上存在采集脑电信号的现有数据少、采集成本高,且由于人体中本来存在的脑电信号是未知的,故采集脑电信号中的噪音信号也是未知的,使得无法对采集脑电信号进行去噪处理,因此,在实际中使得难以获取上述训练数据。本说明书的一些实施例可以生成大量的目标模拟采集脑电信号以及无噪音的模拟脑电信号以分别对采集脑电信号以及去噪后的采集脑电信号进行模拟,从而可以作为训练数据对去噪模型进行训练,使得训练后的去噪模型有更好的泛化效果,去噪模型的输出更加准确。
图6是根据本说明书一些实施例所示的获得第一参数组的示例性流程图。在一些实施例中,流程600可以由服务器110或第一获取模块210执行。如图6所示,流程600可以包括以下步骤:
步骤610,获取初始第一参数组。
初始第一参数组可以指第一预设参数区间对应的参数组。当第一预设参数区间确定后,对应的初始第一参数组可以被确定。在一些实施例中,第一参数组可以包括感兴趣区域转移矩阵以及信号源转移矩阵,对应的,初始第一参数组可以包括初始感兴趣区域转移矩阵以及初始信号源转移矩阵。
步骤620,基于初始第一参数组与基础信号,获得中间模拟脑电信号。
中间脑电信号可以指通过初始第一参数组对基础信号进行处理后的得到的模拟的人体脑电信号。在一些实施例中,可以基于初始第一参数组对基础信号进行处理,得到中间脑电信号。基于初始第一参数组对基础信号进行处理的方法可以参考本说明书中基于第一参数组对基础信号进行处理得到模拟脑电信号的部分。例如,基于初始第一参数组对基础信号进行处理的方法可以参考图4中的公式(2)~公式(5),将公式中的感兴趣区域转移矩阵替换为初始感兴趣区域转移矩阵,将信号源转移矩阵替换为初始信号源转移矩阵,最终得到的信号可以为中间模拟脑电信号。
步骤630,获取参考脑电信号。
参考脑电信号可以指人体中无噪音的脑电信号。参考脑电信号可以通过历史采集到的人体脑电信号获取,或通过网络获取。
步骤640,基于中间模拟脑电信号以及参考脑电信号构建损失函数。
步骤650,基于损失函数迭代更新初始第一参数组的参数,直至满足预设条件,获得第一参数组。
在一些实施例中,可以基于损失函数迭代更新初始第一参数组的参数。当初始第一参数组与参考脑电信号的损失函数满足预设条件时,获得第一参数组。其中,预设条件可以是中间模拟信号与参考脑电信号的差异小于阈值。预设条件还可以是其他条件,例如,损失函数收敛、迭代的次数达到阈值等。更新初始第一参数组的参数的方法可以为梯度下降法。在一些实施例中,还可以为损失函数添加正则化项,以提高迭代效率,加速损失函数收敛。
在一些实施例中,可以基于初始感兴趣区域转移矩阵、初始信号源转移矩阵以及基础信号获得中间模拟脑电信号,并基于该中间模拟脑电信号以及参考脑电信号构建损失函数,并基于损失函数迭代更新初始感兴趣区域转移矩阵以及初始信号源转移矩阵的参数,当该损失函数满足预设条件时,感兴趣区域转移矩阵以及信号源转移矩阵中具体的参数可以被确定。
在一些实施例中,损失函数可以包括第一损失项和第二损失项。第一损失项可以反映初始感兴趣区域转移矩阵与信号源发出的基础信号的关系,第二损失项可以反映初始信号源转移矩阵与感兴趣区域中每一信号源的目标基础信号的关系。第一损失项与第二损失项可以以各种方式结合。例如,加权求和,可以为第一损失项和第二损失项设置不同的权重以体现不同因素对脑电信号的影响不同。
本说明书的一些实施例可以为初始参数组预先设定一个较大的参数区间范围,即第一预设参数区间,可以避免更多计算量,基于该参数区间进行迭代更新,可以更快得到迭代计算结果,从而可以得到第一参数组。
与基于第一预设参数区间获取第一参数组类似的,可以基于第二预设参数区间获取第二参数组。即,可以获取与第二预设参数区间对应的初始第二参数组,并基于初始第二参数组以及初始噪音信号,获取中间噪音信号;获取参考噪音信号;基于中间噪音信号参与参考噪音信号构建损失函数;基于损失函数迭代更新所述初始第二参数组,直至满足预设条件,获得第二参数组。其中,中间噪音信号可以指通过初始第二参数组对初始噪音信号进行处理后的得到的模拟的采集脑电信号中的噪音信号,参考噪音信号可以指实际的采集脑电信号中的噪音信号。关于上述实施例的更多内容可以参考图6相关描述。
信号转移关系可以用于表示信号在不同电极点之间进行传播时的转移关系。在一些实施例中,信号转移关系可以包括电极点转移矩阵,前述电极点转移矩阵又被称为信号转移矩阵。信号转移关系还可以包括信号转移模型。在一些实施例中,可以基于信号转移关系对采集脑电信号进行去噪处理。
图7是根据本说明书一些实施例所示的采集脑电信号去噪方法的示例性流程图。在一些实施例中,流程700可以由服务器110执行。如图7所示,流程700可以包括以下步骤:
步骤710,确定被测对象的信号转移关系,其中,信号转移关系表征初始噪音信号在被测对象的各个位置之间进行传播时的信号变化关系。在一些实施例中,步骤710可以由第一确定模块执行。
被测对象可以指需要确定信号转移关系的对象。被测对象可以是生物体,例如,人体的头部区域。被测对象也可以是非生物体,例如,虚拟头部模型。
初始噪音信号可以指采集脑电信号中的噪音信号未经过传播前的信号。采集脑电信号可以指通过电极点采集到的被测对象的脑电信号。
应当理解的是,在对被测对象的脑电信号进行采集时,存在多种因素(例如,脑电采集设备的运作或被测对象的生理活动)对采集过程产生影响,从而使得采集脑电信号中包含有多种噪音信号。采集脑电信号中常见的噪音信号可以包括非生物伪迹信号与生物伪迹信号。其中,非生物伪迹信号可以为外部环境、设备等带来的噪音信号,可以包括但不限于脑电采集设备所带来的噪音信号、市电干扰所带来的噪音信号等。非生物伪迹信号可以通过查询相关资料获取。例如,可以通过网络确定该型号的脑电采集设备带来的噪音信号。生物伪迹信号可以为被测对象的生理活动所带来的噪音信号,可以包括但不限于心脏跳动所带来的心电信号、肌肉活动所带来的肌电信号以及眼球转动所带来的眼电信号等。生物伪迹信号可以通过相关采集设备对被测对象进行采集获取。因生物个体的差异性,不同被测对象对应的生物伪迹信号应不同。
信号转移关系可以用于确定初始噪音信号在被测对象的各个位置之间进行传播时的变化。在一些实施例中,基于信号转移关系可以确定出被测对象的某一位置的信号传播到另一位置时所发生的变化。例如,基于信号转移关系,可以确定从被测对象的位置
Figure 203944DEST_PATH_IMAGE037
传入的某一信号
Figure 507319DEST_PATH_IMAGE038
传播到被测对象的位置
Figure 290203DEST_PATH_IMAGE039
后,变化为了信号
Figure 37360DEST_PATH_IMAGE040
。在一些实施例中,由于被测对象的生理情况存在不同、脑部构造形状存在不同,故不同的被测对象对应的信号转移关系不同。例如,老人和儿童的信号转移关系不同。再例如,26岁的女士与26岁的男士的信号转移关系也不同。
在一些实施例中,信号转移关系可以表征为N*N的信号转移矩阵,其中,N为被测对象的电极数量,每个电极可以表示被测对象的不同位置。信号转移矩阵中的各个元素可以为被测对象的对应两个位置之间的传导函数。应当理解的是,该信号转移矩阵中某一位置与该位置自身之间的传导函数可以为1,以表示信号不进行传播时其本身不会发生变化。在一些实施例中,关于获取信号转移矩阵的更多内容可以参见图4及其相关描述。
在一些实施例中,信号转移关系还可以表征为信号转移模型。可以基于信号转移模型对初始噪音信号进行处理,获得目标噪音信号。关于信号转移模型的更多内容可以参见图5及其相关描述。
在一些实施例中,当所述被测对象为虚拟头部模型时,可以基于该虚拟头部模型的材质与形状,确定信号该虚拟头部模型中各个位置的衰减关系,从而确定出该虚拟头部模型的信号转移关系。
在一些实施例中,还可以从被测对象的第一位置处输入初始测试信号,其中,初始测试信号的数量级大于被测对象的脑电信号的数量级;从被测对象的第二位置处采集信号,获得采集测试信号;基于初始测试信号以及采集测试信号,确定信号转移关系。关于上述实施例的更多内容可参见图4及其相关描述。
步骤720,基于初始噪音信号以及信号转移关系,确定目标噪音信号。在一些实施例中,步骤720可以由第二确定模块执行。
目标噪音信号可以指采集脑电信号中的噪音信号。在一些实施例中,目标噪音信号可以为噪音集合,该噪音集合中被测对象的各个位置采集到的噪音信号。
在一些实施例中,在确定初始噪音信号后,可以确定出各个初始噪音信号对应的基础特征。其中,初始噪音信号对应的基础特征可以包括但不限于信号的幅度、方差等。针对多个电极点的每一个电极点,在采集脑电信号时,可以对该电极点采集到的信号进行动态评估,通过比较采集到的各个信号的基础特征与各个初始噪音信号对应的基础特征,判断该电极点采集到的信号中是否存在初始噪音信号。当该电极点采集到的信号中存在初始噪音信号时,可以将该初始噪音信号确定为该电极点对应采集脑电信号中的目标噪音信号。同时,还可以通过信号转移关系确定其他电极点采集到的传播后该初始噪音信号,将传播后的初始噪音信号确定为其他电极点采集到的采集脑电信号中的目标噪音信号。
步骤730,从被测对象中采集到的采集脑电信号去除目标噪音信号,获得无噪音的目标脑电信号。在一些实施例中,步骤730可以由信号处理模块执行。
目标脑电信号可以指被测对象无噪音的脑电信号。目标脑电信号可以用于对被测对象的脑部疾病进行分析、诊断。
在一些实施例中,当某一电极点采集到的目标噪音信号为该电极点处未经传播的初始噪音信号时,可以直接从该电极点采集到的采集脑电信号中去除该初始噪音信号,获得该电极点的目标脑电信号。
在一些实施例中,当某一电极点采集到的目标噪音信号为传播后的初始噪音信号时,可以基于以下公式确定该电极点的目标脑电信号:
Figure 365530DEST_PATH_IMAGE041
其中,
Figure 256387DEST_PATH_IMAGE042
为位置为j的电极点的目标脑电信号,
Figure 690080DEST_PATH_IMAGE043
为该电极点的采集脑电信号,h为信号转移矩阵,
Figure 5043DEST_PATH_IMAGE044
表示信号转移矩阵中信号由位置i传播至位置j的传导函数,
Figure 427801DEST_PATH_IMAGE045
为源于位置为i的电极点的初始噪音信号,k为初始噪音信号的数量。
在一些实施例中,可以基于信号转移模型对某一位置的初始噪音信号进行处理,获取其他位置的采集脑电信号中的目标噪音信号。再去除其他位置的采集脑电信号中的目标噪音信号,以获得对应位置的目标脑电信号。
本说明书的一些实施例中通过确定信号在不同位置间进行传播时的信号转移关系,可以基于初始噪音信号和信号转移关系,可以准确而快速的确定传播后的目标噪音信号,从而可以对采集脑电信号中的目标噪音信号进行去除,获得无噪音的目标脑电信号。相较与现有的对采集脑电信号进行去噪的方法,成本大幅降低。同时,基于每一被测对象确定对应的信号转移关系,可以使得确定出的信号转移关系更加准确,符合该被测对象的实际情况,在准确、完整地去除采集脑电信号中的目标噪音信号的基础上,不会将目标脑电信号去除。
图8是根据本说明书一些实施例所示的获取信号转移关系的示例性流程图。在一些实施例中,流程800可以由第一确定模块执行。如图8所示,流程800可以包括以下步骤:
步骤810,针对被测对象的每一位置,从该位置输入初始测试信号,其中,初始测试信号的数量级大于被测对象的脑电信号的数量级。
初始测试信号可以指输入被测对象的测试信号,初始测试信号可以用于测试信号在不同电极点之间进行传播时的变化关系。在一些实施例中,初始测试信号的数量级为大于人体的脑电信号的数量级。例如,人体的脑电信号的数量级为
Figure 661861DEST_PATH_IMAGE046
,初始测试信号的数量级可以比人体的脑电信号大3个数量级,即选取的初始测试信号的数量级可以为
Figure 848430DEST_PATH_IMAGE047
。在一些实施例中,初始测试信号可以包括多个不同频率、不同强度的信号。在一些实施例中,初始测试信号可以通过预先设置确定。例如,初始测试信号可以为预先设置的脉冲信号,该脉冲信号的数量级大于被测对象的脑电信号的数量级。被测对象的每一位置可以指被测对象的每一电极点,可以在被测对象的一电极点处输入初始测试信号。
应当理解的是,由于初始测试信号的数量级大于被测对象的脑电信号的数量级,故当对信号进行采集时,采集到的信号受到脑电信号的影响较小,可以保证测试结果的精确度。同时,初始测试信号也应该不高于人体能够承受的信号的数量级,以保证被测对象的安全,例如,初始测试信号的幅值范围可以为20mV~200mV。
步骤820,从被测对象的其他位置采集信号,获得采集测试信号。
采集测试信号可以指在其他位置处采集到上述位置输入的初始测试信号经传播后的信号。其他位置可以为其他电极点。对应的,采集测试信号可以包括在其他电极点处基于该位置的电极点采集到的信号。例如,可以在被测对象上设置有电极点A~D,其中,输入初始测试信号的位置为电极点A,其他位置可以指电机点B、C、D。对应的,可以从电机点B、C、D处采集信号以获得采集测试信号。
应当理解的是,初始测试信号在电极点之间进行传播时会受到其他电极点的影响,所以采集测试信号可以不同于初始测试信号。
步骤830,基于初始测试信号以及采集测试信号,确定信号转移关系。
在一些实施例中,可以进行建模或采用各种数据分析算法,例如回归分析法、判别分析法等,对初始测试信号以及采集测试信号进行分析处理,确定信号转移关系。
在一些实施例中,当信号转移矩阵中的元素可以实数时,可以通过基于如下公式确定矩阵中的各个元素:
Figure 29400DEST_PATH_IMAGE048
其中,
Figure 149407DEST_PATH_IMAGE049
表示电极点j采集到的采集测试信号,
Figure 882264DEST_PATH_IMAGE050
为信号从电极点i传播至电极点j的调整参数,
Figure 161324DEST_PATH_IMAGE051
表示从电极点i输入的初始测试信号。
在一些实施例中,信号转移关系可以包括多个简单信号转移关系。简单信号转移关系可以指信号变化较为规律的简单信号在被测对象的多个位置之间进行传播时的变化关系。在一些实施例中,可以将输入的初始测试信号变换为多个不同频率和强度的初始简单测试信号,将其他各个电极点处采集到的采集测试信号变换为多个不同频率和强度的采集简单测试信号;针对每个初始简单测试信号,基于该初始简单测试信号及其对应的采集简单测试信号,确定该初始简单测试信号对应的简单信号转移关系;合成各个简单信号转移关系,获得包含有多个简单信号转移关系的信号转移关系。初始简单测试信号可以为对初始测试信号进行变换后获得的简单信号,采集简单测试信号可以为对采集测试信号进行变换后获得的简单信号,其中,信号的变换方法可以为傅里叶变换。例如,初始简单测试信号和采集简单测试信号可以均为正弦波信号。
在一些实施例中,可以对位置i处输入的初始测试信号进行傅里叶变换,获得M个不同相位和频率的初始简单测试信号集合S,其中,S为N*M的矩阵,N为电极点的数量。在被测对象其他位置未输入信号且初始测试信号未进行传播的情况下,S中的第i行中的各个元素分别各个初始简单测试信号,S中的其他位置可以为0,表征其他位置不包含有信号。在初始测试信号进行传播后,对其他位置处采集到的采集测试信号进行傅里叶变换,可以获得包含多个采集简单测试信号集合S',其中,S'也为N*M的矩阵,S'中的各元素表示在对应电极点处采集到的采集简单测试信号,该采集简单测试信号为对应的初始简单测试信号经过传播后的信号。然后可以基于S与S'进行计算,确定出信号转移矩阵H以使得下列公式成立:
Figure 166889DEST_PATH_IMAGE052
其中,H的元素均为传导函数,表示信号在对应位置进行传播时的变化关系。例如,
Figure 236344DEST_PATH_IMAGE053
可以为H中的元素,表示信号从位置i传播至位置j的传导函数。在一些实施例中,传导函数可以表现为多种形式。如下列公式所示,传导函数
Figure 272258DEST_PATH_IMAGE053
可以等价于一个FIR(FiniteImpulse Response,有限长单位冲激响应)滤波器:
Figure 320504DEST_PATH_IMAGE054
其中,y[n]为采集简单测试信号,x[n]为初始简单测试信号,n为信号长度,L为参数数量,可以通过预先设置确定,当
Figure 854996DEST_PATH_IMAGE055
时, 可以通过解多元一次方程组的方式确定出前述传导函数
Figure 563976DEST_PATH_IMAGE053
中的所有参数。
在一些实施例中,还可以基于初始测试信号以及采集测试信号对信号转移模型进行训练,获得训练好的信号转移模型。关于信号转移模型的更多内容可以参见图9及其相关描述。
应当理解的是,被测对象的脑部结构几乎是固定的,故即使两个信号的强度不同,其在被测对象的位置之间进行传播的信号变化关系应该相同。因此,本说明书的一些实施例可以通过对初始测试信号以及采集测试信号进行分析处理,确定被测对象的信号转移关系,前述信号转移关系也可以表征初始噪音信号在被测对象的各个位置之间进行传播时的信号变化关系。
本说明书的一些实施例针对每一个被测对象,通过初始测试信号以及采集测试信号,可以快速确定出该被测对象的信号转移关系,保证了信号转移关系的准确性,从而避免了在去噪时因被测对象个体的差异所带来的误差。
图9是根据本说明书一些实施例所示的信号转移模型的示意图。
在一些实施例中,可以基于信号转移模型对初始噪音信号进行处理,获得目标噪音信号。如图9所示,信号转移模型930的输入可以包括初始噪音信号910以及初始噪音信号的位置920,输出可以包括目标噪音信号940。在一些实施例中,信号转移模型可以包括但不限于深度神经网络模型、支持向量机模型等。
在一些实施例中,信号转移模型可以是利用训练样本对机器学习模型进行训练而获得的。如图9所述,信号转移模型可以是利用训练样本960以及标签970对初始信号转移模型950进行训练获得的,其中,初始信号转移模型可以为未设置参数的信号转移模型。训练样本960可以包括初始测试信号以及初始测试信号的位置,标签970可以包括采集测试信号。训练样本及标签的获取方式可以参见图8及其相关描述。将多组带有标签970的训练样本960输入初始信号转移模型950,基于初始信号转移模型950的输出和标签970构建损失函数,基于损失函数迭代更新初始信号转移模型950的参数,直至满足预设条件,训练结束,获取训练好的信号转移模型930。预设条件可以包括但不限于损失函数小于阈值、收敛或训练周期达到阈值。
在一些实施例中,信号转移模型可以包括对象信息嵌入层与信号确定层,对象信息嵌入层用于对被测对象的形状信息以及头发信息进行特征提取,获得头部特征,信号确定层用于对头部特征、初始噪音信号以及初始噪音信号的位置信息处理,获得目标噪音信号。其中,形状信息可以包括但不限于被测对象的形状、尺寸等,头发信息可以包括但不限于被测对象头发长度、头发硬度、头发弯曲程度以及头发油脂情况等。被测对象的形状信息以及头发信息可以通过预先设置确定。
如图10所示,信号转移模型930还可以包括依次连接的对象信息嵌入层930-1以及信号确定层930-2,对象信息嵌入层930-1的输入可以包括形状信息1010以及头发信息1020,输出为头部特征931,信号确定层的输入可以包括头部特征931、初始噪音信号910以及初始噪音信号的位置920,输出可以为目标噪音信号940。其中,对象信息嵌入层可以为朴素贝叶斯模型,信号确定层可以为深度神经网络模型。
在一些实施例中,可以通过训练获取对象信息嵌入层:训练样本可以包括测试对象的历史形状信息以及历史头发状况信息,标签可以包括前述测试对象的历史头部特征,其中,测试对象可以指用于获取其数据以训练对象信息嵌入层的对象。与被测对象类似,测试对象可以是生物体或非生物体。前述训练样本和标签可以通过对人工表述测试对象的相关信息获取。可以将多组训练样本输入初始对象信息嵌入层,基于初始对象信息嵌入层的输出与标签构建损失函数,基于损失函数迭代更新初始对象信息嵌入层的参数,直至预设条件被满足,获得训练好的对象信息嵌入层。预设条件可以包括但不限于损失函数小于阈值、收敛或训练周期达到阈值。
在一些实施例中,可以通过训练获取信号确定层:训练样本可以包括历史头部特征、被测对象的初始测试信号以及初始测试信号的位置,标签可以包括采集测试信号,前述训练样本及其标签的获取方式可以参考本说明书前文。可以将多组训练样本输入初始信号确定层,基于初始信号确定层的输出与标签构建损失函数,基于损失函数迭代更新初始信号确定层的参数,直至预设条件被满足,获得训练好的信号确定层。预设条件可以包括但不限于损失函数小于阈值、收敛或训练周期达到阈值。
应当理解的是,用于训练对象信息嵌入层与信号确定层的训练数据对应的对象存在差异性。训练好的对象信息嵌入层可以对所有对象(包括测试对象以及被测对象)的头部特征进行提取,对象信息嵌入层为各个被测对象通用的模型,故基于测试对象的相关数据进行训练获取的对象信息嵌入层可以提取出不同被测对象的头部特征。而训练好的信号确定层用于确定某一被测对象的信号转移关系。如前文所述,当被测对象不同时,其对应的信号转移关系也不同。故基于某一被测对象进行训练获取的信号确定层仅可以用于确定该被测对象的信号转移关系。当需要确定另一被测对象的信号转移关系时,可以在训练好的对象信息嵌入层的基础上,基于该被测对象的初始测试信号以及初始测试信号的位置以及采集测试信号进行训练。
本说明书的一些实施例通过信号转移模型可以更加快速、便捷的确定出被测对象的信号转移关系,以便于对采集脑电信号进行去噪处理。另外,本说明书的一些实施例还可以对被测对象的头部的相关特征进行提取,确定出被测对象头部相关信息的差异性对被测对象的信号转移关系的影响,使得对采集脑电信号的去噪处理更加准确。
图11是根据本说明书一些实施例所示的确定目标噪音信号的示例性流程图。在一些实施例中,流程1100可以由第二确定模块执行。
在一些实施例中,信号转移关系可以包括多个简单信号转移关系。简单信号转移关系表征第一简单噪音信号在所述被测对象的各个位置之间进行传播时的信号变化关系。关于简单信号转移关系的更多内容可以参见图8及其相关描述。第一简单噪声信号是对初始噪音信号进行变换而获得的。
当信号转移关系包括多个简单信号转移关系时,可以执行流程1100以确定目标噪音信号。如图11所示,流程1100可以包括以下步骤:
步骤1110,将初始噪音信号变换为多个不同频率和强度的第一简单噪声信号。
第一简单噪声信号可以是通过对初始噪音信号进行变换的简单信号。例如,第一简单噪声信号可以包括正弦波信号、余弦波信号和方波信号等中的任意一种或组合。可以利用傅里叶变换、拉普拉斯变换以及离散余弦变换等对初始噪音信号进行处理,以获取多个不同频率和强度的第一简单噪声信号。
步骤1120,针对每个第一简单噪声信号,基于该第一简单噪声信号对应的简单信号转移关系对该第一简单噪声信号进行处理,确定该第一简单噪声信号对应的第二简单噪音信号。
第二简单噪音信号可以是指第一简单噪音信号传播到被测对象的其他位置时的噪音信号。
在一些实施例中,针对每个第一简单噪声信号,可以基于第一简单噪声信号、第一简单噪音信号的输入位置以及简单信号转移关系,确定该第一简单噪声信号在其他位置对应的第二简单噪音信号。可以根据简单信号转移关系以及第一简单噪音信号的输入位置,确定输入位置传播至其他各个位置的传导函数,基于各个传导函数以及第一简单噪音信号,可以确定出其他各个位置的第二简单噪音信号。
步骤1130,合成各个第二简单噪音信号,获得目标噪音信号。
在一些实施例中,针对每一被测对象的位置,可以对该位置的各个第二简单噪音信号进行合成,获得该位置的目标噪音信号。其中,合成方式可以为上述将初始噪音信号处理为第一简单噪音信号的变换方式的逆变换,来获取目标噪音信号。例如,当基于傅里叶变换将初始噪音信号处理为第一简单噪音信号时,可以基于傅里叶逆变换合成各个第二简单噪音信号,获得目标噪音信号。
本说明书的一些实施例通过利用简单噪音信号以及简单信号转移关系,可以比较方便地获取目标噪音信号,使运算得以简化,提高计算效率。
本说明书实施例提供了一种模拟采集脑电信号的生成装置,包括处理器,前述处理器用于执行如本说明书实施例中任一项所述的模拟采集脑电信号的生成方法。
本说明书实施例提供了一种计算机可读存储介质,前述存储介质存储计算机指令,当计算机读取存储介质中的计算机指令后,计算机执行如本说明书实施例中任一项所述的模拟采集脑电信号的生成方法。
本说明书实施例还提供了一种采集脑电信号去噪装置,包括处理器,前述处理器用于执行如本说明书实施例中任一项所述的采集脑电信号去噪方法。
本说明书实施例提供了一种计算机可读存储介质,前述存储介质存储计算机指令,当计算机读取存储介质中的计算机指令后,计算机执行如本说明书实施例中任一项所述的采集脑电信号去噪方法。
应当注意的是,上述有关各个流程的描述仅仅是为了示例和说明,而不限定本说明书的适用范围。对于本领域技术人员来说,在本说明书的指导下可以对各个流程进行各种修正和改变。然而,这些修正和改变仍在本说明书的范围之内。
上文已对基本概念做了描述,显然,对于本领域技术人员来说,上述详细披露仅仅作为示例,而并不构成对本说明书的限定。虽然此处并没有明确说明,本领域技术人员可能会对本说明书进行各种修改、改进和修正。该类修改、改进和修正在本说明书中被建议,所以该类修改、改进、修正仍属于本说明书示范实施例的精神和范围。同时,本说明书使用了特定词语来描述本说明书的实施例。如“一个实施例”、“一实施例”、和/或“一些实施例”意指与本说明书至少一个实施例相关的某一特征、结构或特点。因此,应强调并注意的是,本说明书中在不同位置两次或多次提及的“一实施例”或“一个实施例”或“一个替代性实施例”并不一定是指同一实施例。此外,本说明书的一个或多个实施例中的某些特征、结构或特点可以进行适当的组合。
此外,除非权利要求中明确说明,本说明书所述处理元素和序列的顺序、数字字母的使用、或其他名称的使用,并非用于限定本说明书流程和方法的顺序。尽管上述披露中通过各种示例讨论了一些目前认为有用的发明实施例,但应当理解的是,该类细节仅起到说明的目的,附加的权利要求并不仅限于披露的实施例,相反,权利要求旨在覆盖所有符合本说明书实施例实质和范围的修正和等价组合。例如,虽然以上所描述的系统组件可以通过硬件设备实现,但是也可以只通过软件的解决方案得以实现,如在现有的服务器或移动设备上安装所描述的系统。
同理,应当注意的是,为了简化本说明书披露的表述,从而帮助对一个或多个发明实施例的理解,前文对本说明书实施例的描述中,有时会将多种特征归并至一个实施例、附图或对其的描述中。但是,这种披露方法并不意味着本说明书对象所需要的特征比权利要求中提及的特征多。实际上,实施例的特征要少于上述披露的单个实施例的全部特征。
一些实施例中使用了描述成分、属性数量的数字,应当理解的是,此类用于实施例描述的数字,在一些示例中使用了修饰词“大约”、“近似”或“大体上”来修饰。除非另外说明,“大约”、“近似”或“大体上”表明所述数字允许有±20%的变化。相应地,在一些实施例中,说明书和权利要求中使用的数值参数均为近似值,该近似值根据个别实施例所需特点可以发生改变。在一些实施例中,数值参数应考虑规定的有效数位并采用一般位数保留的方法。尽管本说明书一些实施例中用于确认其范围广度的数值域和参数为近似值,在具体实施例中,此类数值的设定在可行范围内尽可能精确。
针对本说明书引用的每个专利、专利申请、专利申请公开物和其他材料,如文章、书籍、说明书、出版物、文档等,特此将其全部内容并入本说明书作为参考。与本说明书内容不一致或产生冲突的申请历史文件除外,对本说明书权利要求最广范围有限制的文件(当前或之后附加于本说明书中的)也除外。需要说明的是,如果本说明书附属材料中的描述、定义、和/或术语的使用与本说明书所述内容有不一致或冲突的地方,以本说明书的描述、定义和/或术语的使用为准。
最后,应当理解的是,本说明书中所述实施例仅用以说明本说明书实施例的原则。其他的变形也可能属于本说明书的范围。因此,作为示例而非限制,本说明书实施例的替代配置可视为与本说明书的教导一致。相应地,本说明书的实施例不仅限于本说明书明确介绍和描述的实施例。

Claims (10)

1.一种模拟采集脑电信号的生成方法,其特征在于,包括:
基于基础信号和第一参数组,获得无噪音的模拟脑电信号,所述模拟脑电信号用于作为人体中无噪音的脑电信号参与脑电信号的分析与研究;
基于初始噪音信号和第二参数组,获得目标噪音信号,所述目标噪音信号用于作为人体的采集脑电信号中的噪音信号参与所述脑电信号的分析与研究;
获取所述模拟脑电信号的第一权重以及所述目标噪音信号的第二权重;
基于所述模拟脑电信号、所述目标噪音信号、所述第一权重以及所述第二权重,获得目标模拟采集脑电信号,所述目标模拟采集脑电信号用于作为所述采集脑电信号参与所述脑电信号的分析与研究以及减少真实脑电信号的采集。
2.如权利要求1所述的方法,其特征在于,所述第一参数组包括感兴趣区域转移矩阵以及信号源转移矩阵,所述感兴趣区域转移矩阵中元素表征信号在不同感兴趣区域之间进行传播时的转移关系,所述信号源转移矩阵中的元素表征同一感兴趣区域中的信号源与电极点之间进行传播时的转移关系,其中,所述不同感兴趣区域表示人体大脑中的不同脑功能区;
所述基于基础信号和第一参数组,获得无噪音的模拟脑电信号包括:
针对虚拟脑部模型中的每一信号源,基于所述感兴趣区域转移矩阵以及该信号源发出的所述基础信号,确定该信号源的目标基础信号;
针对所述虚拟脑部模型中每一感兴趣区域,基于该感兴趣区域中所述每一信号源的所述目标基础信号以及所述信号源转移矩阵,确定该感兴趣区域对应的初始子模拟脑电信号;
基于各个所述感兴趣区域对应的初始子模拟脑电信号,获得所述模拟脑电信号。
3.如权利要求1所述的方法,其特征在于,所述初始噪音信号包括模拟脑部噪音信号以及电极点噪音信号,所述第二参数组包括感兴趣区域转移矩阵以及电极点转移矩阵,所述感兴趣区域转移矩阵中元素表征信号在不同感兴趣区域之间进行传播时的转移关系,所述电极点转移矩阵中的元素表征信号在不同电极点之间进行传播时的转移关系;
所述基于初始噪音信号和第二参数组,获得目标噪音信号包括:
基于所述感兴趣区域转移矩阵以及所述模拟脑部噪音信号,确定第一子噪音信号;
基于所述电极点转移矩阵以及所述电极点噪音信号,确定第二子噪音信号;
基于所述第一子噪音信号以及所述第二子噪音信号,获得所述目标噪音信号。
4.如权利要求1所述的方法,其特征在于,所述第一参数组基于第一预设参数区间获得,所述第二参数组基于第二预设参数区间获得。
5.一种模拟采集脑电信号的生成系统,其特征在于,所述系统包括:
第一获取模块,用于基于基础信号和第一参数组,获得无噪音的模拟脑电信号,所述模拟脑电信号用于作为人体中无噪音的脑电信号参与脑电信号的分析与研究;
第二获取模块,用于基于初始噪音信号和第二参数组,获得目标噪音信号,所述目标噪音信号用于作为人体的采集脑电信号中的噪音信号参与所述脑电信号的分析与研究;
第三获取模块,用于获取所述模拟脑电信号的第一权重以及所述目标噪音信号的第二权重;
第四获取模块,用于基于所述模拟脑电信号、所述目标噪音信号、所述第一权重以及所述第二权重,获得目标模拟采集脑电信号,所述目标模拟采集脑电信号用于作为所述采集脑电信号参与所述脑电信号的分析与研究以及减少真实脑电信号的采集。
6.如权利要求5所述的系统,其特征在于,所述第一参数组包括感兴趣区域转移矩阵以及信号源转移矩阵,所述感兴趣区域转移矩阵中元素表征信号在不同感兴趣区域之间进行传播时的转移关系,所述信号源转移矩阵中的元素表征同一感兴趣区域中的信号源与电极点之间进行传播时的转移关系,其中,所述不同感兴趣区域表示人体大脑中的不同脑功能区;
所述第一获取模块进一步用于:
针对虚拟脑部模型中的每一信号源,基于所述感兴趣区域转移矩阵以及该信号源发出的所述基础信号,确定该信号源的目标基础信号;
针对所述虚拟脑部模型中每一感兴趣区域,基于该感兴趣区域中所述每一信号源的所述目标基础信号以及所述信号源转移矩阵,确定该感兴趣区域对应的初始子模拟脑电信号;
基于各个所述感兴趣区域对应的初始子模拟脑电信号,获得所述模拟脑电信号。
7.如权利要求5所述的系统,其特征在于,所述初始噪音信号包括模拟脑部噪音信号以及电极点噪音信号,所述第二参数组包括感兴趣区域转移矩阵以及电极点转移矩阵,所述第二参数组包括感兴趣区域转移矩阵以及电极点转移矩阵,所述感兴趣区域转移矩阵中元素表征信号在不同感兴趣区域之间进行传播时的转移关系,所述电极点转移矩阵中的元素表征信号在不同电极点之间进行传播时的转移关系;
所述第二获取模块进一步用于:
基于所述感兴趣区域转移矩阵以及所述模拟脑部噪音信号,确定第一子噪音信号;
基于所述电极点转移矩阵以及所述电极点噪音信号,确定第二子噪音信号;
基于所述第一子噪音信号以及所述第二子噪音信号,获得所述目标噪音信号。
8.如权利要求5所述的系统,其特征在于,所述第一参数组基于第一预设参数区间获得,所述第二参数组基于第二预设参数区间获得。
9.一种模拟采集脑电信号的生成装置,包括处理器,其特征在于,所述处理器用于执行如权利要求1~4任一项所述的模拟采集脑电信号的生成方法。
10.一种计算机可读存储介质,所述存储介质存储计算机指令,其特征在于,当计算机读取存储介质中的计算机指令后,计算机执行如权利要求1~4任一项所述的模拟采集脑电信号的生成方法。
CN202210659471.8A 2022-06-13 2022-06-13 一种模拟采集脑电信号的生成方法和系统 Active CN114742116B (zh)

Priority Applications (4)

Application Number Priority Date Filing Date Title
CN202210910794.XA CN115238745A (zh) 2022-06-13 2022-06-13 一种模拟采集脑电信号的生成方法和系统
CN202210659471.8A CN114742116B (zh) 2022-06-13 2022-06-13 一种模拟采集脑电信号的生成方法和系统
CN202210910733.3A CN115186718A (zh) 2022-06-13 2022-06-13 一种模拟采集脑电信号的生成方法和系统
CN202210666384.5A CN115034266A (zh) 2022-06-13 2022-06-13 一种采集脑电信号去噪方法、系统、装置及介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210659471.8A CN114742116B (zh) 2022-06-13 2022-06-13 一种模拟采集脑电信号的生成方法和系统

Related Child Applications (3)

Application Number Title Priority Date Filing Date
CN202210666384.5A Division CN115034266A (zh) 2022-06-13 2022-06-13 一种采集脑电信号去噪方法、系统、装置及介质
CN202210910794.XA Division CN115238745A (zh) 2022-06-13 2022-06-13 一种模拟采集脑电信号的生成方法和系统
CN202210910733.3A Division CN115186718A (zh) 2022-06-13 2022-06-13 一种模拟采集脑电信号的生成方法和系统

Publications (2)

Publication Number Publication Date
CN114742116A true CN114742116A (zh) 2022-07-12
CN114742116B CN114742116B (zh) 2022-09-02

Family

ID=82287409

Family Applications (4)

Application Number Title Priority Date Filing Date
CN202210910733.3A Pending CN115186718A (zh) 2022-06-13 2022-06-13 一种模拟采集脑电信号的生成方法和系统
CN202210666384.5A Pending CN115034266A (zh) 2022-06-13 2022-06-13 一种采集脑电信号去噪方法、系统、装置及介质
CN202210659471.8A Active CN114742116B (zh) 2022-06-13 2022-06-13 一种模拟采集脑电信号的生成方法和系统
CN202210910794.XA Pending CN115238745A (zh) 2022-06-13 2022-06-13 一种模拟采集脑电信号的生成方法和系统

Family Applications Before (2)

Application Number Title Priority Date Filing Date
CN202210910733.3A Pending CN115186718A (zh) 2022-06-13 2022-06-13 一种模拟采集脑电信号的生成方法和系统
CN202210666384.5A Pending CN115034266A (zh) 2022-06-13 2022-06-13 一种采集脑电信号去噪方法、系统、装置及介质

Family Applications After (1)

Application Number Title Priority Date Filing Date
CN202210910794.XA Pending CN115238745A (zh) 2022-06-13 2022-06-13 一种模拟采集脑电信号的生成方法和系统

Country Status (1)

Country Link
CN (4) CN115186718A (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116636817B (zh) * 2023-07-26 2023-11-03 四川新源生物电子科技有限公司 一种麻醉深度评估方法、系统、装置和存储介质

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060149160A1 (en) * 2004-12-03 2006-07-06 Aspect Medical Systems, Inc. System and method for EEG imaging of cerebral activity using small electrode sets
WO2015196107A1 (en) * 2014-06-20 2015-12-23 Boston Scientific Scimed Inc. Medical devices for mapping cardiac tissue
CN110090017A (zh) * 2019-03-11 2019-08-06 北京工业大学 一种基于lstm的脑电信号源定位方法
CN110859600A (zh) * 2019-12-06 2020-03-06 深圳市德力凯医疗设备股份有限公司 一种生成脑电信号的方法、存储介质及电子设备
CN110879980A (zh) * 2019-11-13 2020-03-13 厦门大学 基于神经网络算法的核磁共振波谱去噪方法
CN111461204A (zh) * 2020-03-30 2020-07-28 华南理工大学 一种用于游戏评估的基于脑电信号的情感识别方法
CN112807000A (zh) * 2021-02-04 2021-05-18 首都师范大学 鲁棒性脑电信号的生成方法及装置
CN114175152A (zh) * 2019-08-01 2022-03-11 杜比实验室特许公司 用于增强劣化音频信号的系统和方法
CN114190953A (zh) * 2021-12-09 2022-03-18 四川新源生物电子科技有限公司 针对脑电采集设备的脑电信号降噪模型的训练方法和系统
CN114521904A (zh) * 2022-01-25 2022-05-24 中山大学 一种基于耦合神经元群的脑电活动模拟方法及系统

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060149160A1 (en) * 2004-12-03 2006-07-06 Aspect Medical Systems, Inc. System and method for EEG imaging of cerebral activity using small electrode sets
WO2015196107A1 (en) * 2014-06-20 2015-12-23 Boston Scientific Scimed Inc. Medical devices for mapping cardiac tissue
CN110090017A (zh) * 2019-03-11 2019-08-06 北京工业大学 一种基于lstm的脑电信号源定位方法
CN114175152A (zh) * 2019-08-01 2022-03-11 杜比实验室特许公司 用于增强劣化音频信号的系统和方法
CN110879980A (zh) * 2019-11-13 2020-03-13 厦门大学 基于神经网络算法的核磁共振波谱去噪方法
CN110859600A (zh) * 2019-12-06 2020-03-06 深圳市德力凯医疗设备股份有限公司 一种生成脑电信号的方法、存储介质及电子设备
CN111461204A (zh) * 2020-03-30 2020-07-28 华南理工大学 一种用于游戏评估的基于脑电信号的情感识别方法
CN112807000A (zh) * 2021-02-04 2021-05-18 首都师范大学 鲁棒性脑电信号的生成方法及装置
CN114190953A (zh) * 2021-12-09 2022-03-18 四川新源生物电子科技有限公司 针对脑电采集设备的脑电信号降噪模型的训练方法和系统
CN114521904A (zh) * 2022-01-25 2022-05-24 中山大学 一种基于耦合神经元群的脑电活动模拟方法及系统

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
ANDREW HAMILTON-WRIGHT 等: "Physiologically Based Simulation of Clinical EMG Signals", 《IEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING》 *
TUOMAS P. MUTANEN 等: "Automatic and robust noise suppression in EEG and MEG: The SOUND algorithm", 《NEUROIMAGE》 *
李肃义 等: "一种综合小波变换的心电信号消噪算法", 《仪器仪表学报》 *
汤璘 等: "融合RFID和脑电信号的机械臂联合控制方法", 《控制工程》 *
陈泽龙 等: "脑电EEG信号的分析分类方法", 《中国医学装备大会暨 2019 医学装备展览会论文汇编》 *
陈灿: "前额柔性干式电极及单通道脑电预处理技术研究", 《中国优秀博硕士学位论文全文数据库(硕士) 信息科技辑》 *

Also Published As

Publication number Publication date
CN115238745A (zh) 2022-10-25
CN115186718A (zh) 2022-10-14
CN114742116B (zh) 2022-09-02
CN115034266A (zh) 2022-09-09

Similar Documents

Publication Publication Date Title
Krasoulis et al. Improved prosthetic hand control with concurrent use of myoelectric and inertial measurements
Lorenz et al. The Automatic Neuroscientist: A framework for optimizing experimental design with closed-loop real-time fMRI
Harrison et al. Multivariate autoregressive modeling of fMRI time series
EP0699413B1 (en) Apparatus and method for analyzing information relating to physical and mental condition
CN105764415B (zh) 将非平稳信号分解成函数分量
CN111714118B (zh) 一种基于集成学习的脑认知模型融合方法
CN111728590A (zh) 基于动态功能连接的个体认知能力预测方法和系统
Li et al. Emotion stimuli-based surface electromyography signal classification employing Markov transition field and deep neural networks
CN112185493A (zh) 人格偏好诊断装置及基于该装置的项目推荐系统
Khaleghi et al. A neuronal population model based on cellular automata to simulate the electrical waves of the brain
CN114742116B (zh) 一种模拟采集脑电信号的生成方法和系统
Gu et al. Nonlinear modeling of cortical responses to mechanical wrist perturbations using the NARMAX method
Nasrolahzadeh et al. Analysis of mean square error surface and its corresponding contour plots of spontaneous speech signals in Alzheimer's disease with adaptive wiener filter
Zhang et al. Objective Extraction of Evoked Event‐Related Oscillation from Time‐Frequency Representation of Event‐Related Potentials
Mouleeshuwarapprabu et al. Nonlinear vector decomposed neural network based EEG signal feature extraction and detection of seizure
Raghu et al. Automated biomedical signal quality assessment of electromyograms: Current challenges and future prospects
Lopes et al. Ensemble deep neural network for automatic classification of EEG independent components
Nunez et al. A tutorial on fitting joint models of M/EEG and behavior to understand cognition
CN113988123A (zh) 一种基于自加权增量rvfl网络的脑电疲劳预测方法
CN111772629B (zh) 一种脑认知技能移植的方法
Lu et al. Time series modeling of surface EMG based hand manipulation identification via expectation maximization algorithm
Abe et al. Paraconsistent neurocomputing and brain signal analysis
Chen et al. Research on AR‐AKF Model Denoising of the EMG Signal
Backenroth et al. Nonnegative decomposition of functional count data
Jaber et al. Elicitation hybrid spatial features from HD-sEMG signals for robust classification of gestures in real-time

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant