CN104605845A - 一种基于diva模型的脑电信号处理方法 - Google Patents
一种基于diva模型的脑电信号处理方法 Download PDFInfo
- Publication number
- CN104605845A CN104605845A CN201510050938.9A CN201510050938A CN104605845A CN 104605845 A CN104605845 A CN 104605845A CN 201510050938 A CN201510050938 A CN 201510050938A CN 104605845 A CN104605845 A CN 104605845A
- Authority
- CN
- China
- Prior art keywords
- brain
- dipole
- data
- fmir
- model
- 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
Links
- 238000003672 processing method Methods 0.000 title claims abstract description 18
- 238000000034 method Methods 0.000 claims abstract description 58
- 238000007781 pre-processing Methods 0.000 claims abstract description 4
- 230000008569 process Effects 0.000 claims description 28
- 239000011159 matrix material Substances 0.000 claims description 25
- 210000004556 brain Anatomy 0.000 claims description 22
- 210000004761 scalp Anatomy 0.000 claims description 16
- 239000000203 mixture Substances 0.000 claims description 14
- 210000003128 head Anatomy 0.000 claims description 13
- 230000004913 activation Effects 0.000 claims description 12
- 238000005457 optimization Methods 0.000 claims description 11
- 238000009826 distribution Methods 0.000 claims description 10
- 230000005684 electric field Effects 0.000 claims description 10
- 230000004927 fusion Effects 0.000 claims description 10
- 238000004088 simulation Methods 0.000 claims description 9
- 230000009466 transformation Effects 0.000 claims description 7
- 238000003062 neural network model Methods 0.000 claims description 5
- PXFBZOLANLWPMH-UHFFFAOYSA-N 16-Epiaffinine Natural products C1C(C2=CC=CC=C2N2)=C2C(=O)CC2C(=CC)CN(C)C1C2CO PXFBZOLANLWPMH-UHFFFAOYSA-N 0.000 claims description 4
- 238000012937 correction Methods 0.000 claims description 4
- 230000005672 electromagnetic field Effects 0.000 claims description 4
- 238000001914 filtration Methods 0.000 claims description 4
- 238000009499 grossing Methods 0.000 claims description 4
- 210000002615 epidermis Anatomy 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 3
- 238000002922 simulated annealing Methods 0.000 claims description 3
- 238000002599 functional magnetic resonance imaging Methods 0.000 claims description 2
- 238000012880 independent component analysis Methods 0.000 abstract description 21
- 238000004458 analytical method Methods 0.000 abstract description 8
- 238000004422 calculation algorithm Methods 0.000 abstract description 8
- 238000012545 processing Methods 0.000 abstract description 8
- 238000011160 research Methods 0.000 abstract description 5
- 210000005036 nerve Anatomy 0.000 abstract description 3
- 230000007547 defect Effects 0.000 abstract description 2
- 238000002474 experimental method Methods 0.000 abstract description 2
- 230000035479 physiological effects, processes and functions Effects 0.000 abstract description 2
- 230000035945 sensitivity Effects 0.000 abstract description 2
- 238000007499 fusion processing Methods 0.000 abstract 1
- 230000006870 function Effects 0.000 description 12
- 238000013507 mapping Methods 0.000 description 11
- 230000000694 effects Effects 0.000 description 7
- 210000000056 organ Anatomy 0.000 description 7
- 230000001537 neural effect Effects 0.000 description 6
- 210000003625 skull Anatomy 0.000 description 3
- 230000002087 whitening effect Effects 0.000 description 3
- 238000007476 Maximum Likelihood Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 210000003710 cerebral cortex Anatomy 0.000 description 2
- 210000001175 cerebrospinal fluid Anatomy 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000004615 ingredient Substances 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000000926 neurological effect Effects 0.000 description 2
- 230000007935 neutral effect Effects 0.000 description 2
- 238000003786 synthesis reaction Methods 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 210000003484 anatomy Anatomy 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 210000001638 cerebellum Anatomy 0.000 description 1
- 230000002490 cerebral effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000004020 conductor Substances 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 210000004704 glottis Anatomy 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 238000002386 leaching Methods 0.000 description 1
- 238000002595 magnetic resonance imaging Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000002156 mixing Methods 0.000 description 1
- 230000036403 neuro physiology Effects 0.000 description 1
- 238000002610 neuroimaging Methods 0.000 description 1
- 210000002569 neuron Anatomy 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 210000005182 tip of the tongue Anatomy 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明公开了一种基于DIVA模型的脑电信号处理方法,该方法通过DIVA模型模拟生成的fMIR数据对脑电信号进行定位分析,并使用独立成分分析方法简化实际计算的复杂度,克服了非侵入式脑电信号分辨率低、干扰大的缺点。本发明使用DIVA模型产生的fMIR数据对EEG数据进行融合处理,克服了脑电信号空间分辨率低、信号干扰大、信噪比很低的问题。通过ICA的预处理,减小了运算的复杂性,在很大程度上克服了等效偶极子定位算法对于噪声的敏感。最后,利用本方法对真实实验数据进行处理,得到的结论符合生理学事实。本方法为汉语神经分析***中的脑电信号处理问题提供了可行的解决方案,为今后汉语语音生成与获取相关研究奠定了基础。
Description
技术领域
本发明公开了一种基于DIVA模型的脑电信号处理方法,涉及脑电信号处理技术领域。
背景技术
脑-机接口(brain computer interface,BCI)是一种基于脑电信号来实现人脑与计算机或其它电子设备进行通讯和控制的***,它是一个不依靠外周神经和肌肉组织等通常的大脑输出通道的通信***。换言之,BCI是在人脑与计算机之间建立的直接的交流和控制通道,通过这种通道,人就可以直接通过脑来表达想法或操纵设备,而不需要语言或肢体动作。
波士顿大学的弗朗克·冈瑟(Frank Guenther)教授带领的科研小组在脑机接口技术基础上研制成功研制成一种神经分析***(Neuralynx System)。神经分析***由脑—计算机接口(BCI)和语音合成***DIVA(Directions IntoVelocities of Articulators)模型两部分组成。该***可以将将人大脑里的思维过程“阅读”出来,然后将其转换为正常语言进行实时表述。
DIVA模型则是一种具有生物学意义的关于语音生成和获取的神经网络。DIVA模型由前馈控制子***、反馈控制子***以及前田(Maeda)模拟声道所组成。模型通过某种规则在以语音作为输入的同时,产生一个发音速度以及器官位置变化的时变序列,应用这个序列,***就可以得到所需要的理想发音。
DIVA模型的一个重要特征就是把模型组件与大脑皮层各相关区域一一对应起来。模型组件与大脑皮层区域的映射关系是在有关神经解剖学和神经生理学研究的基础上根据MNI(Montreal Neurological Institute)标准参照系而确定的。在按照预设参数进行发音任务时,DIVA模型能够产生反应脑区激活状态的fMIR数据。这样,DIVA模型实际上就构成了一个用以解释来自各种研究的相关数据以及对关于语音神经处理过程进行一致性描述的基本框架。
在语音相关研究过程中,普遍采用非侵入式的神经成像术作为脑机之间的接口在人身上进行了实验。虽然这种非侵入式的装置方便佩戴于人体,但是由于颅骨对信号的衰减作用和对神经元发出的电磁波的分散和模糊效应,记录到信号的分辨率并不高。这种信号波仍可被检测到,但信号信噪比较低,对后处理的要求较高。
发明内容
本发明所要解决的技术问题是:针对现有技术的缺陷,提供一种基于DIVA模型的脑电信号处理方法。该方法通过DIVA模型模拟生成的fMIR数据对脑电信号进行定位分析,并使用独立成分分析方法(Independent ComponentCorrelation Algorithm,ICA)简化实际计算的复杂度,克服了非侵入式脑电信号分辨率低、干扰大的缺点。
本发明为解决上述技术问题采用以下技术方案:
一种基于DIVA模型的脑电信号处理方法,利用DIVA神经网络模型产生的功能性核磁共振数据对脑电信号进行融合处理,包括以下步骤:
步骤一、通过非侵入式脑机接口提取发音过程中的脑电信号数据,对采集到的脑电信号数据进行预处理;
步骤二、基于DIVA模型提供的用户接口对声道参数进行设置,使DIVA模型进行模拟发音,获得fMIR数据;
将fMIR数据输入统计绘图工具SPM中进行分析,采用12参数的仿射变换对fMIR数据进行归一化处理,再将处理后的图像与高分辨率的结构图像进行配准,并归一到MNI的空间,然后利用三维高斯函数对经过上述处理的数据做空间平滑处理;
步骤三、对经过步骤二处理的脑电信号数据进行初步的白化或球化处理,去除各观测信号之间的相关性,然后利用ICA方法对脑电信号数据进行处理,得出最优变换矩阵,进而恢复信号源矩阵,从脑电信号数据中分离出有效的成分,每个有效的成分对应一个偶极子成;
步骤四、结合头部模型及其传导系数,对每一个提取的有效的成分求解最优的偶极子配置,使得由上述偶极子产生的头皮电场分布与测量到的头皮电场分布在均方意义下最小,将fMIR数据中的激活点位置作为寻找偶极子位置的种子点,在fMIR数据限定源的可能空间中求解全局最优解;
步骤五、将融合处理后的结果输出。
作为本发明的进一步优选方案,所述步骤三中,利用ICA方法对脑电信号数据进行处理具体如下:
设定记录的头皮电位xi作为观测向量,其中,i=1,2,.....,m,m表示电极的数量,观测向量xi是n个时间上独立的信号源sj的线性混和,其中,j=1,2,...,n,每一个sj都是统计独立并且有其固定的空间信息权重aj,aj是导联场矩阵A的第j列;
在上述设定下,使用ICA方法直接寻找最优的变换矩阵W,并恢复信号源矩阵M;使用估计的信号Y对观测信号X进行重建:
X=W-1Y (1)
X的第j个时间序列可以由下式表示:
其中,是W-1的第n行第m列的元素。
作为本发明的进一步优选方案,所述步骤三中,使用ICA方法对脑电信号数据进行处理,进行下述设定:
301、独立信号源产生的信号是统计独立的;
302、观测到的脑电信号数据是瞬时的线性混和信号;
303、独立信号源的数量不超过电极的数量。
作为本发明的进一步优选方案,所述步骤四的具体过程如下:
设定大脑表皮的电磁场观测信号与大脑内部任何位置的源信号近似地呈线性关系,由下式表示:
X=AS+N (3)
其中,X是头皮电极记录的观测信号,A是导联场矩阵,S是偶极子矢量,N表示每个电极接收信号中的噪声组成的矢量;
导联场矩阵A是偶极子位置、大脑几何形状以及介质传导系数的非线性函数,它的每一列表示位于某一位置的单位强度的源对左端观测信号X的贡献;
优化问题等价于求解下式:
其中,C是信噪比归一化矩阵,λ是正则化参数;
将目标函数与解的空间性质相联系,定义目标函数为:
式中,P项是对解的一种空间约束,即结合fMIR限定源的可能空间位置,由下式定义:
其中,l是偶极子的数量;ri是第i个偶极子的位置向量;si是fMIR激活点的空间信息;c和di均是常数,偶极子的位置与fMIR的激活点位置是相关的,fMIR的激活点位置即是寻找偶极子位置的种子点;
应用最大距离限制拟合偶极子的位置以及偶极矩,这就是fMIR受限的等效偶极子模型,简称FC-ECD模型。FC-ECD的一个基本假设是:所观测到的电场分布可以由一个或者几个偶极子产生,同时一定大小区域的激活可以由单个偶极子表示。
建立目标函数之后,将偶极子的参数看作是高维空间中的一个节点,则采用模拟退火算法求全局最优解。
作为本发明的进一步优选方案,步骤一中所述的预处理过程包括:去除眨眼伪迹、眼球运动伪迹、低通滤波、坏电极复位、平均以及基线校正。
作为本发明的进一步优选方案,步骤二中所述的高斯函数需要满足其半高宽为12mm*12mm*24mm。
作为本发明的进一步优选方案,步骤四中所述的偶极子配置包括位置、方向和大小。
作为本发明的进一步优选方案,步骤五中所述的将融合处理后的结果输出,是通过在MATLAB中使用FIT工具箱进行输出。
作为本发明的进一步优选方案,所述常数di的取值为5mm。
本发明采用以上技术方案与现有技术相比,具有以下技术效果:使用DIVA模型产生的fMIR数据对EEG数据进行融合处理,克服了脑电信号空间分辨率低、信号干扰大、信噪比很低的问题。通过ICA的预处理,减小了运算的复杂性,在很大程度上克服了等效偶极子定位算法对于噪声的敏感。最后,利用本方法对真实实验数据进行处理,得到的结论符合生理学事实。本方法为汉语神经分析***中的脑电信号处理问题提供了可行的解决方案,为今后汉语语音生成与获取相关研究奠定了基础。
附图说明
图1为DIVA神经网络模型的结构示意图;
图2为汉语神经分析***结构示意图;
图3为现有独立成分分析示意图;
图4为理想4层头部模型。
具体实施方式
下面结合附图对本发明的技术方案做进一步的详细说明:
本发明要解决的技术问题是提供一种基于DIVA模型的脑电信号处理方法。包括以下步骤:
步骤1、通过非侵入式脑机接口提取发音过程中的脑电信号(EEG)。对采集到的EEG数据进行预处理,包括去除眨眼伪迹、眼球运动伪迹、低通滤波、坏电极复位、平均以及基线校正;
步骤2、基于DIVA模型提供的用户接口对声道参数进行设置,使模型分别进行模拟发音,获得fMIR数据。然后将模拟产生的数据结果输入统计绘图工具SPM中进行分析,采用12参数的仿射变换对数据进行归一化处理。再将图像与高分辨率的结构图像进行配准,并归一到MNI的空间。然后利用半高宽(FWHM)为12mm*12mm*24mm的三维高斯函数对标准化后的数据做空间平滑处理;
步骤3、对数据进行初步的白化或球化处理,除各观测信号之间的相关性。然后利用ICA方法对EEG数据进行处理,提取有效成分;
在EEG信号的测量中,所观测到的信号实际上是由若干相对独立的源信号叠加而成的。因此,采用ICA方法分解出观测信号的独立成分将有助于我们把握真正有意义的大脑活动。在EEG信号处理中,可以将记录的头皮电位xi,i=1,2,.....,m作为观测向量,m表示电极的数量。假设观测信号矢量xi是n个时间上独立的信号源sj,j=1,2,...,n的线性混和。每一个sj都是统计独立的,并且有其固定的空间信息权重aj。aj是导联场矩阵A的第j列。在上面的表述下,ICA方法直接寻找最优的变换矩阵W,利用ICA的原理恢复信号源矩阵M。因此,从估计的信号Y对观测信号X进行的重建可以通过下式进行:
X=W-1Y (1)
因此,X的第j个时间序列可以由下式表示:
其中,是W-1的第n行第m列的元素。在EEG的数据处理过程中,为了使用ICA方法,一般都作如下假设:
1、信号源产生的信号是统计独立的;
2、观测到的脑电信号是瞬时的线性混和信号;
3、独立信号源的数量不超过电极的数量
如式(2)所示,ICA方法从众多的EEG测量信号中分离出有效的成分,每个有效的成分对应一个偶极子成。
步骤4、对每一个提取有效的成分利用fMIR数据进行融合处理;
人脑的电磁场可以通过半稳态的麦克斯维尔方程解释。根据这一理论,大脑表皮的电磁场观测信号与大脑内部任何位置的源信号近似地呈线性关系,可由下式表示:
X=AS+N (3)
其中,X是头皮电极记录的观测信号,A称为导联场矩阵(lead field matrix),S是偶极子矢量,N表示每个电极接收信号中的噪声组成的矢量。导联场矩阵A是偶极子位置、大脑几何形状以及介质传导系数的非线性函数,它的每一列表示位于某一位置的单位强度的源对左端观测信号X的贡献。
EEG逆向问题是指已知头部模型以及传导系数,寻找最优的偶极子配置(包括位置、方向、大小),使得由这些偶极子产生的头皮电场分布与测量到的头皮电场分布在均方意义下最小。等价于求解下面的优化问题:
其中,C是信噪比归一化矩阵,λ是正则化参数。更一般地,可以把目标函数与解的某种空间性质联系起来,据此可以定义目标函数为:
式中,P项是对解的某种空间约束,也就是结合fMIR限定源的可能空间位置,可以由下式定义:
其中,l是偶极子的数量;ri是第i个偶极子的位置向量;si是fMIR激活点的空间信息;c,di是常数,在本文中,di=5mm。偶极子的位置与fMIR的激活点位置是相关的,fMIR的激活点位置可以看作是寻找偶极子位置的种子点。同时,应用最大距离限制拟合偶极子的位置以及偶极矩,这就是fMIR受限的等效偶极子模型,简称FC-ECD模型。FC-ECD的一个基本假设是:所观测到的电场分布可以由一个或者几个偶极子产生,同时一定大小区域的激活可以由单个偶极子表示。
建立目标函数之后,将偶极子的参数看作是高维空间中的一个节点,采用模拟退火算法求全局最优解。
步骤5、在MATLAB中使用FIT工具箱将融合处理后的结果输出。
下面给出本发明的一个具体实施例:
步骤1、通过非侵入式脑机接口提取汉语元音ɑ发音过程中的脑电信号(EEG)。对采集到的EEG数据进行预处理,包括去除眨眼伪迹、眼球运动伪迹、低通滤波、坏电极复位、平均以及基线校正;
步骤2、基于DIVA模型提供的用户接口将前三个共振峰频率值为805Hz、1265Hz和2770Hz,使模型分别进行汉语元音ɑ模拟发音,获得fMIR数据。然后将模拟产生的数据结果输入统计绘图工具SPM中进行分析,采用12参数的仿射变换对数据进行归一化处理。再将图像与高分辨率的结构图像进行配准,并归一到MNI的空间。然后利用半高宽(FWHM)为12mm*12mm*24mm的三维高斯函数对标准化后的数据做空间平滑处理;
步骤3、对数据进行初步的白化或球化处理,除各观测信号之间的相关性。利用独立成分分析方法(ICA)得出最优变换矩阵,进而恢复信号源矩阵S。从众多的EEG测量信号中分离出有效的成分,每个有效的成分对应一个偶极子成。
步骤4、在已知头部模型以及其传导系数的前提下,对每一个提取的有效的成分求解最优的偶极子配置(包括位置、方向、大小),使得由这些偶极子产生的头皮电场分布与测量到的头皮电场分布在均方意义下最小。将fMIR的激活点位置作为是寻找偶极子位置的种子点,在fMIR限定源的可能空间中求解全局最优解。
步骤5、在MATLAB中使用FIT工具箱将融合处理后的结果输出。
如图1所示,DIVA神经网络模型主要是由前馈控制子***、反馈控制子***和声道所构成的。所述前馈控制子***包括:语音映射集的一部分、小脑、发音器官速率和位置映射集(部分);反馈控制子***包括语音映射集的一部分、听觉误差映射集、听觉状态映射集、体觉误差映射集、体觉位置映射集、发音器官速率和位置映射集;声道采用前田声道模型,把模拟声道分成8个组成部分:舌头的位置、舌头的形状、舌尖、嘴唇高度、唇凸、以及颚、喉咙的高度和声门的张合程度。语音映射集的作用是根据模型中已经存在的映射关系把模型的输入进行映射编码;听觉状态和误差映射集的作用是通过编码来调基频F0和前三个共振峰频率F1~F3的位置来描述当前发音;发音器官速率和位置映射集的作用是来决定DIVA神经网络模型的声道框架中各发音器官的位置;体觉状态和误差映射集的作用是调整发音器官的各个部件的位置和参数来调整当前发音。
如图2所示,汉语神经分析***由两部分组成:脑-计算机接口和语音合成***DIVA模型。BCI中,脑电信号通过使用者佩戴的无线神经电极获取,而检测到的神经信号则被用于驱动语音合成器的连续“运动”,为使用者提供实时语音输出。语音合成器DIVA模型则是一种具有生物学意义的关于语音生成和获取的神经网络。DIVA模型由前馈控制子***、反馈控制子***以及前田(Maeda)模拟声道所组成。模型通过某种规则在以语音作为输入的同时,产生一个发音速度以及器官位置变化的时变序列,应用这个序列,***就可以得到所需要的理想发音。
如图3所示,ICA方法是基于以下假设的:信号源产生的信号是统计独立的;观测到的脑电信号是瞬时的线性混和信号独立信号源的数量不超过电极的数量。
一般情况下,所获得的数据都具有相关性,所以通常都要求对数据进行初步的白化或球化处理,因为白化处理可去除各观测信号之间的相关性,从而简化了后续独立分量的提取过程,而且,通常情况下,数据进行白化处理与不对数据进行白化处理相比,算法的收敛性较好。之后通过解混***将输入信号它们分离开来,输出有效成分。ICA的预处理可以有效地降低问题的复杂度,而且算法简单。ICA算法的研究可分为基于信息论准则的迭代估计方法和基于统计学的代数方法两大类,从原理上来说,它们都是利用了源信号的独立性和非高斯性。基于信息论的方法研究中,各国学者从最大熵、最小互信息、最大似然和负熵最大化等角度提出了一系列估计算法。如FastICA算法,Infomax算法,最大似然估计算法等。基于统计学的方法主要有二阶累积量、四阶累积量等高阶累积量方法。
如图4所示,本专利采用理想的4层头部模型作为真实大脑头部模型的模拟,即将头部分割成4个部分:脑、脑脊液、颅骨及头皮,它提供了大脑体积导体的一个简化模型。很多研究人员采用实际影像数据(CT/MRI)等数据,利用FEM(finite element method)或者BEM(boundary element method)方法,对头部模型进行建模;然后,再通过分割、配准等预处理过程将头部又分为4个部分(脑、脑脊液、颅骨及头皮)。但是,这些过程需要大量的计算和很长的时间。而采用4层模型的一个最大的好处就是在节约计算时间的前提下,能够提供一个对于真实头部模型很好的近似。有研究表明:4层球模型即使与更加真实的有限元模型相比,仍不失为一个有效的模拟工具,它可提供一个误差在10%~20%以内的合理的头表电位估计。如图4所示,本文采用的4层球模型的半径依次是79mm、81mm、85mm及88mm,其传导系数分别为0.461S/m、1.39S/m、0.0058S/m和0.461S/m。可以根据式(3)进行EEG正向问题的计算。正向问题的解不仅能够在后面的仿真实验中提供模拟的头表电势分布,而且能够在求解逆向问题的迭代过程中起到辅助推断的作用。
上面结合附图对本发明的实施方式作了详细说明,但是本发明并不限于上述实施方式,在本领域普通技术人员所具备的知识范围内,还可以在不脱离本发明宗旨的前提下做出各种变化。
Claims (9)
1.一种基于DIVA模型的脑电信号处理方法,其特征在于,利用DIVA神经网络模型产生的功能性核磁共振数据对脑电信号进行融合处理,包括以下步骤:
步骤一、通过非侵入式脑机接口提取发音过程中的脑电信号数据,对采集到的脑电信号数据进行预处理;
步骤二、基于DIVA模型提供的用户接口对声道参数进行设置,使DIVA模型进行模拟发音,获得fMIR数据;
将fMIR数据输入统计绘图工具SPM中进行分析,采用12参数的仿射变换对fMIR数据进行归一化处理,再将处理后的图像与高分辨率的结构图像进行配准,并归一到MNI的空间,然后利用三维高斯函数对经过上述处理的数据做空间平滑处理;
步骤三、对经过步骤二处理的脑电信号数据进行初步的白化或球化处理,去除各观测信号之间的相关性,然后利用ICA方法对脑电信号数据进行处理,得出最优变换矩阵,进而恢复信号源矩阵,从脑电信号数据中分离出有效的成分,每个有效的成分对应一个偶极子成;
步骤四、结合头部模型及其传导系数,对每一个提取的有效的成分求解最优的偶极子配置,使得由上述偶极子产生的头皮电场分布与测量到的头皮电场分布在均方意义下最小,将fMIR数据中的激活点位置作为寻找偶极子位置的种子点,在fMIR数据限定源的可能空间中求解全局最优解;
步骤五、将融合处理后的结果输出。
2.如权利要求1所述的一种基于DIVA模型的脑电信号处理方法,其特征在于,所述步骤三中,利用ICA方法对脑电信号数据进行处理具体如下:
设定记录的头皮电位xi作为观测向量,其中,i=1,2,.....,m,m表示电极的数量,观测向量xi是n个时间上独立的信号源sj的线性混和,其中,j=1,2,...,n,每一个sj都是统计独立并且有其固定的空间信息权重aj,aj是导联场矩阵A的第j列;
在上述设定下,使用ICA方法直接寻找最优的变换矩阵W,并恢复信号源矩阵M;使用估计的信号Y对观测信号X进行重建:
X=W-1Y (1)
X的第j个时间序列可以由下式表示:
其中,是W-1的第n行第m列的元素。
3.如权利要求2所述的一种基于DIVA模型的脑电信号处理方法,其特征在于,所述步骤三中,使用ICA方法对脑电信号数据进行处理,进行下述设定:
301、独立信号源产生的信号是统计独立的;
302、观测到的脑电信号数据是瞬时的线性混和信号;
303、独立信号源的数量不超过电极的数量。
4.如权利要求3所述的一种基于DIVA模型的脑电信号处理方法,其特征在于,所述步骤四的具体过程如下:
设定大脑表皮的电磁场观测信号与大脑内部任何位置的源信号近似地呈线性关系,由下式表示:
X=AS+N (3)
其中,X是头皮电极记录的观测信号,A是导联场矩阵,S是偶极子矢量,N表示每个电极接收信号中的噪声组成的矢量;
导联场矩阵A是偶极子位置、大脑几何形状以及介质传导系数的非线性函数,它的每一列表示位于某一位置的单位强度的源对左端观测信号X的贡献;
优化问题等价于求解下式:
其中,C是信噪比归一化矩阵,λ是正则化参数;
将目标函数与解的空间性质相联系,定义目标函数为:
式中,P项是对解的一种空间约束,即结合fMIR限定源的可能空间位置,由下式定义:
其中,l是偶极子的数量;ri是第i个偶极子的位置向量;si是fMIR激活点的空间信息;c和di均是常数,偶极子的位置与fMIR的激活点位置是相关的,fMIR的激活点位置即是寻找偶极子位置的种子点;
建立目标函数之后,将偶极子的参数看作是高维空间中的一个节点,采用模拟退火算法求全局最优解。
5.如权利要求1或4所述的一种基于DIVA模型的脑电信号处理方法,其特征在于,步骤一中所述的预处理过程包括:去除眨眼伪迹、眼球运动伪迹、低通滤波、坏电极复位、平均以及基线校正。
6.如权利要求1或4所述的一种基于DIVA模型的脑电信号处理方法,其特征在于,步骤二中所述的高斯函数需要满足其半高宽为12mm*12mm*24mm。
7.如权利要求1或4所述的一种基于DIVA模型的脑电信号处理方法,其特征在于,步骤四中所述的偶极子配置包括位置、方向和大小。
8.如权利要求1或4所述的一种基于DIVA模型的脑电信号处理方法,其特征在于,步骤五中所述的将融合处理后的结果输出,是通过在MATLAB中使用FIT工具箱进行输出。
9.如权利要求4所述的一种基于DIVA模型的脑电信号处理方法,其特征在于,所述常数di的取值为5mm。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510050938.9A CN104605845B (zh) | 2015-01-30 | 2015-01-30 | 一种基于diva模型的脑电信号处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510050938.9A CN104605845B (zh) | 2015-01-30 | 2015-01-30 | 一种基于diva模型的脑电信号处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104605845A true CN104605845A (zh) | 2015-05-13 |
CN104605845B CN104605845B (zh) | 2017-01-25 |
Family
ID=53140705
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510050938.9A Expired - Fee Related CN104605845B (zh) | 2015-01-30 | 2015-01-30 | 一种基于diva模型的脑电信号处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104605845B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105054928A (zh) * | 2015-07-17 | 2015-11-18 | 张洪振 | 一种基于bci装置脑电波采集分析的情绪展示设备 |
CN105266804A (zh) * | 2015-11-12 | 2016-01-27 | 杭州电子科技大学 | 一种基于低秩和稀疏矩阵分解的新型脑电信号处理方法 |
CN107961007A (zh) * | 2018-01-05 | 2018-04-27 | 重庆邮电大学 | 一种结合卷积神经网络和长短时记忆网络的脑电识别方法 |
CN108369660A (zh) * | 2015-07-13 | 2018-08-03 | 索邦大学 | 用事件间的时间间隔来表示数值的数据处理装置 |
CN108742602A (zh) * | 2018-03-26 | 2018-11-06 | 校宝在线(杭州)科技股份有限公司 | 一种基于脑电波的结巴检测方法 |
CN111513711A (zh) * | 2020-05-22 | 2020-08-11 | 电子科技大学 | 一种基于参考电极的脑电坏导插值方法 |
CN111616704A (zh) * | 2020-06-24 | 2020-09-04 | 天津大学 | 针对句子听力任务大脑动态功能网络交互模式的研究方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040116791A1 (en) * | 2002-11-27 | 2004-06-17 | Satoru Miyauchi | Method and apparatus for analyzing brain functions |
US20060036153A1 (en) * | 2004-06-14 | 2006-02-16 | Laken Steven J | Questions and control paradigms for detecting deception by measuring brain activity |
CN102880906A (zh) * | 2012-07-10 | 2013-01-16 | 南京邮电大学 | 一种基于diva神经网络模型的汉语元音发音方法 |
CN103310272A (zh) * | 2013-07-02 | 2013-09-18 | 南京邮电大学 | 基于声道动作知识库改进的diva神经网络模型发音方法 |
-
2015
- 2015-01-30 CN CN201510050938.9A patent/CN104605845B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040116791A1 (en) * | 2002-11-27 | 2004-06-17 | Satoru Miyauchi | Method and apparatus for analyzing brain functions |
US20060036153A1 (en) * | 2004-06-14 | 2006-02-16 | Laken Steven J | Questions and control paradigms for detecting deception by measuring brain activity |
CN102880906A (zh) * | 2012-07-10 | 2013-01-16 | 南京邮电大学 | 一种基于diva神经网络模型的汉语元音发音方法 |
CN103310272A (zh) * | 2013-07-02 | 2013-09-18 | 南京邮电大学 | 基于声道动作知识库改进的diva神经网络模型发音方法 |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108369660A (zh) * | 2015-07-13 | 2018-08-03 | 索邦大学 | 用事件间的时间间隔来表示数值的数据处理装置 |
CN105054928A (zh) * | 2015-07-17 | 2015-11-18 | 张洪振 | 一种基于bci装置脑电波采集分析的情绪展示设备 |
CN105266804A (zh) * | 2015-11-12 | 2016-01-27 | 杭州电子科技大学 | 一种基于低秩和稀疏矩阵分解的新型脑电信号处理方法 |
CN105266804B (zh) * | 2015-11-12 | 2017-10-17 | 杭州电子科技大学 | 一种基于低秩和稀疏矩阵分解的脑电信号处理方法 |
CN107961007A (zh) * | 2018-01-05 | 2018-04-27 | 重庆邮电大学 | 一种结合卷积神经网络和长短时记忆网络的脑电识别方法 |
CN108742602A (zh) * | 2018-03-26 | 2018-11-06 | 校宝在线(杭州)科技股份有限公司 | 一种基于脑电波的结巴检测方法 |
CN111513711A (zh) * | 2020-05-22 | 2020-08-11 | 电子科技大学 | 一种基于参考电极的脑电坏导插值方法 |
CN111616704A (zh) * | 2020-06-24 | 2020-09-04 | 天津大学 | 针对句子听力任务大脑动态功能网络交互模式的研究方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104605845B (zh) | 2017-01-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104605845A (zh) | 一种基于diva模型的脑电信号处理方法 | |
CN110090017B (zh) | 一种基于lstm的脑电信号源定位方法 | |
CN105054928A (zh) | 一种基于bci装置脑电波采集分析的情绪展示设备 | |
CN105956624B (zh) | 基于空时频优化特征稀疏表示的运动想象脑电分类方法 | |
CN109965869B (zh) | 基于脑源域空间的mi-eeg识别方法 | |
CN109614885A (zh) | 一种基于lstm的脑电信号快速分类识别方法 | |
Duan et al. | 3D face reconstruction from skull by regression modeling in shape parameter spaces | |
CN114052668B (zh) | 一种基于脑磁图数据的脑功能分析方法 | |
CN109199376B (zh) | 基于oa-wmne脑源成像的运动想象脑电信号的解码方法 | |
EP3912014A1 (en) | Brain-computer interface system, system for brain activity analysis, and method of analysis | |
Wei et al. | Edge sparse basis network: a deep learning framework for EEG source localization | |
Jinliang et al. | EEG emotion recognition based on granger causality and capsnet neural network | |
CN112401905A (zh) | 一种基于源定位和脑网络的自然动作脑电识别方法 | |
Jatoi et al. | BEM based solution of forward problem for brain source estimation | |
Zhang et al. | DCT based multi-head attention-BiGRU model for EEG source location | |
CN108416822B (zh) | 一种基于贝叶斯估计的多层次多尺度层析成像方法 | |
CN102743166B (zh) | 一种事件相关电位源定位方法 | |
CN102855666B (zh) | 基于层次化回归模型的颅面重构方法 | |
Tsai et al. | Mapping single-trial EEG records on the cortical surface through a spatiotemporal modality | |
Ramírez | Source localization | |
CN115813409A (zh) | 一种超低延迟的运动图像脑电图解码方法 | |
KR20090008682A (ko) | 뇌파기반 피질리듬활동 모니터링 시스템 | |
CN114428555A (zh) | 一种基于皮层源信号的脑电运动意图识别方法及*** | |
Hu et al. | A Method of EEG Source Imaging Based on TCN and Attention Network | |
Liu et al. | A semi-supervised classification approach based on restricted Boltzmann machine for fMRI data |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20200512 Address after: 710000 25 / F, block D, Tsinghua Science Park, Keji 2nd Road, Zhangba Street office, hi tech Zone, Xi'an City, Shaanxi Province Patentee after: Xi'an Huinao Intelligent Technology Co.,Ltd. Address before: 210003 Gulou District, Jiangsu, Nanjing new model road, No. 66 Patentee before: NANJING UNIVERSITY OF POSTS AND TELECOMMUNICATIONS |
|
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170125 |