一种基于DIVA模型的脑电信号处理方法技术

本发明专利技术公开了一种基于DIVA模型的脑电信号处理方法，该方法通过DIVA模型模拟生成的fMIR数据对脑电信号进行定位分析，并使用独立成分分析方法简化实际计算的复杂度，克服了非侵入式脑电信号分辨率低、干扰大的缺点。本发明专利技术使用DIVA模型产生的fMIR数据对EEG数据进行融合处理，克服了脑电信号空间分辨率低、信号干扰大、信噪比很低的问题。通过ICA的预处理，减小了运算的复杂性，在很大程度上克服了等效偶极子定位算法对于噪声的敏感。最后，利用本方法对真实实验数据进行处理，得到的结论符合生理学事实。本方法为汉语神经分析系统中的脑电信号处理问题提供了可行的解决方案，为今后汉语语音生成与获取相关研究奠定了基础。

【技术实现步骤摘要】
一种基于DIVA模型的脑电信号处理方法
本专利技术公开了一种基于DIVA模型的脑电信号处理方法，涉及脑电信号处理

技术介绍
脑-机接口(braincomputerinterface,BCI)是一种基于脑电信号来实现人脑与计算机或其它电子设备进行通讯和控制的系统，它是一个不依靠外周神经和肌肉组织等通常的大脑输出通道的通信系统。换言之，BCI是在人脑与计算机之间建立的直接的交流和控制通道，通过这种通道，人就可以直接通过脑来表达想法或操纵设备，而不需要语言或肢体动作。波士顿大学的弗朗克·冈瑟(FrankGuenther)教授带领的科研小组在脑机接口技术基础上研制成功研制成一种神经分析系统(NeuralynxSystem)。神经分析系统由脑—计算机接口(BCI)和语音合成系统DIVA(DirectionsIntoVelocitiesofArticulators)模型两部分组成。该系统可以将将人大脑里的思维过程“阅读”出来，然后将其转换为正常语言进行实时表述。DIVA模型则是一种具有生物学意义的关于语音生成和获取的神经网络。DIVA模型由前馈控制子系统、反馈控制子系统以及前田(Maeda)模拟声道所组成。模型通过某种规则在以语音作为输入的同时，产生一个发音速度以及器官位置变化的时变序列，应用这个序列，系统就可以得到所需要的理想发音。DIVA模型的一个重要特征就是把模型组件与大脑皮层各相关区域一一对应起来。模型组件与大脑皮层区域的映射关系是在有关神经解剖学和神经生理学研究的基础上根据MNI(MontrealNeurologicalInstitute)标准参照系而确定的。在按照预设参数进行发音任务时，DIVA模型能够产生反应脑区激活状态的fMIR数据。这样，DIVA模型实际上就构成了一个用以解释来自各种研究的相关数据以及对关于语音神经处理过程进行一致性描述的基本框架。在语音相关研究过程中，普遍采用非侵入式的神经成像术作为脑机之间的接口在人身上进行了实验。虽然这种非侵入式的装置方便佩戴于人体，但是由于颅骨对信号的衰减作用和对神经元发出的电磁波的分散和模糊效应，记录到信号的分辨率并不高。这种信号波仍可被检测到，但信号信噪比较低，对后处理的要求较高。
技术实现思路
本专利技术所要解决的技术问题是：针对现有技术的缺陷，提供一种基于DIVA模型的脑电信号处理方法。该方法通过DIVA模型模拟生成的fMIR数据对脑电信号进行定位分析，并使用独立成分分析方法(IndependentComponentCorrelationAlgorithm,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限定源的可能空间位置，由下式定义：>λ---(6)]]>其中，l是偶极子的数量；ri是第i个偶极子的位置向量；si是fMIR激活点的空间信息；c和di均是常数，偶极子的位置与fMIR的激活点位置是相关的，fMIR的激活点位置即是寻找偶极子位置的种子点；应用最大距离限制拟合偶极子的位置以及偶极矩，这就是fMIR受限的等效偶极子模型，简称FC-ECD模型。FC-ECD的一个基本假设是：所观测到的电场分布可以由一个或者几个偶极子产生，同时一定大小区域的激活可以由单个偶极子表示。建立目标函数之后，将偶极子的参数看作是高维空间中的一个节点，则采用模拟退火算法求全局最优解。作为本专利技术的进一步优选方案，步骤一中所述的预处理过程包括：去除眨眼伪迹、眼球运动伪迹、低通滤波、坏电极复位、平均以及基线校正。作为本专利技术的进一步优选方案，步骤二中所述的高斯函数需要满足其半高宽为12mm*1本文档来自技高网...

【技术保护点】
一种基于DIVA模型的脑电信号处理方法，其特征在于，利用DIVA神经网络模型产生的功能性核磁共振数据对脑电信号进行融合处理，包括以下步骤：步骤一、通过非侵入式脑机接口提取发音过程中的脑电信号数据，对采集到的脑电信号数据进行预处理；步骤二、基于DIVA模型提供的用户接口对声道参数进行设置，使DIVA模型进行模拟发音，获得fMIR数据；将fMIR数据输入统计绘图工具SPM中进行分析，采用12参数的仿射变换对fMIR数据进行归一化处理，再将处理后的图像与高分辨率的结构图像进行配准，并归一到MNI的空间，然后利用三维高斯函数对经过上述处理的数据做空间平滑处理；步骤三、对经过步骤二处理的脑电信号数据进行初步的白化或球化处理，去除各观测信号之间的相关性，然后利用ICA方法对脑电信号数据进行处理，得出最优变换矩阵，进而恢复信号源矩阵，从脑电信号数据中分离出有效的成分，每个有效的成分对应一个偶极子成；步骤四、结合头部模型及其传导系数，对每一个提取的有效的成分求解最优的偶极子配置，使得由上述偶极子产生的头皮电场分布与测量到的头皮电场分布在均方意义下最小，将fMIR数据中的激活点位置作为寻找偶极子位置的种子点，在fMIR数据限定源的可能空间中求解全局最优解；步骤五、将融合处理后的结果输出。...

【技术特征摘要】
1.一种基于DIVA模型的脑电信号处理方法，其特征在于，利用DIVA模型产生的功能性核磁共振数据对脑电信号进行融合处理，包括以下步骤：步骤一、通过非侵入式脑机接口提取发音过程中的脑电信号数据，对采集到的脑电信号数据进行预处理；步骤二、基于DIVA模型提供的用户接口对声道参数进行设置，使DIVA模型进行模拟发音，获得fMRI数据；将fMRI数据输入统计绘图工具SPM中进行分析，采用12参数的仿射变换对fMRI数据进行归一化处理，再将处理后的图像与高分辨率的结构图像进行配准，并归一到MNI的空间，然后利用三维高斯函数对经过上述归一化处理的数据做空间平滑处理；步骤三、对经过步骤一处理的脑电信号数据进行初步的白化或球化处理，去除各观测信号之间的相关性，然后利用ICA方法对脑电信号数据进行处理，得出最优的变换矩阵，进而恢复信号源矩阵，从脑电信号数据中分离出有效的成分，每个有效的成分对应一个偶极子；步骤四、结合头部模型及其传导系数，对每一个提取的有效的成分求解最优的偶极子配置，使得由上述偶极子产生的头皮电场分布与测量到的头皮电场分布在均方意义下最小，将fMRI数据中的激活点位置作为寻找偶极子位置的种子点，在fMRI数据限定源的可能空间中求解全局最优解；步骤五、将融合处理后的结果输出。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的第m行第n列的元素。3.如权利要求2所述的一种基于DIVA模型的脑电信号处理方法，其特征在于，所述步骤三中，利用ICA方法对脑电信号数据进行处理，进行下述设定：301...

【专利技术属性】
技术研发人员：张少白，陈彦霖，
申请(专利权)人：南京邮电大学，
类型：发明
国别省市：江苏;32