基于功能磁共振成像的时变约束脑电图或脑磁图溯源方法技术

本发明专利技术提供一种基于功能磁共振成像的时变约束脑电图或脑磁图溯源方法，步骤包括：记录脑电图或脑磁图信号并进行溯源，得到初始的脑电图或脑磁图源强度值，并据以得到脑电图或脑磁图源强度矩阵；记录功能磁共振成像数据并据以得到功能磁共振成像激活图，进而得到功能磁共振成像权重矩阵；根据所述功能磁共振成像权重矩阵和所述脑电图或脑磁图源强度矩阵获得源的协方差矩阵，将该源的协方差矩阵设置为随时间变化；根据所述初始的脑电图或脑磁图源强度值计算交互通讯效应和归一化部分曲线下面积，并据以修正所述源的协方差矩阵；根据所述修正的源的协方差矩阵进行溯源，得到最终的脑电图或脑磁图源强度值，完成溯源。

Time variant constrained electroencephalogram or magnetoencephalogram tracing method based on functional magnetic resonance imaging

The invention provides a variable constrained EEG or MEG tracing method based on functional magnetic resonance imaging, the method comprises the following steps of: recording EEG or MEG signal and traceability, get the initial EEG or MEG source intensity values, and to get the EEG or MEG source intensity matrix recording; functional magnetic resonance imaging according to the data obtained by functional magnetic resonance imaging activation map, and then get the weight matrix of functional magnetic resonance imaging; based on the covariance matrix of the functional magnetic resonance imaging weight matrix and the EEG or MEG source intensity matrix to get the source, will set the source covariance matrix for change over time; according to the calculated area of interactive communication the normalized effect and the initial part of the curve under the EEG or MEG source intensity, and according to the covariance matrix of the source are modified according to the repair; Is the source of the covariance matrix of traceability, get the final EEG or MEG source intensity, complete source.

【技术实现步骤摘要】
基于功能磁共振成像的时变约束脑电图或脑磁图溯源方法
本专利技术涉及脑电图、脑磁图及磁共振成像领域，尤其涉及一种基于功能磁共振成像的时变约束脑电图或脑磁图溯源方法。
技术介绍
脑电图(Electroencephalography，EEG)或脑磁图(Magnetoencephalography，MEG)有很高的时间分辨率，可以达到毫秒量级，并且他们的信号主要来源于突触后电位，是对大脑活动的直接测量。但是由于脑电图电极数或脑磁图传感器数目远远小于大脑内电流源的数目，因此想要唯一地确定脑神经电流源的位置需要加入外加的约束，而这些外加的约束往往是缺乏生理学依据的。采用脑电图或脑磁图信号研究脑神经活动，确定神经电流源的逆问题算法中具有代表性的如最小二次模算法(MinimumNormEstimation，MNE)、深度加权最小范数法(depth-WeightedMinimumNorm，dWMN)等，但这些方法始终不能摆脱空间分辨率低的问题。基于血氧水平依赖(Blood-Oxygen-Level-Dependence，BOLD)的功能磁共振成像(functionalMagneticResonance本文档来自技高网...

【技术保护点】
一种基于功能磁共振成像的时变约束脑电图或脑磁图溯源方法，步骤包括：记录脑电图或脑磁图信号并进行溯源，得到初始的脑电图或脑磁图源强度值，根据所述初始的脑电图或脑磁图源强度值得到脑电图或脑磁图源强度矩阵；记录功能磁共振成像数据并据以得到功能磁共振成像激活图，根据所述功能磁共振成像激活图得到功能磁共振成像权重矩阵；根据所述功能磁共振成像权重矩阵和所述脑电图或脑磁图源强度矩阵获得源的协方差矩阵，将该源的协方差矩阵设置为随时间变化；根据所述初始的脑电图或脑磁图源强度值计算交互通讯效应和归一化部分曲线下面积，并根据该交互通讯效应和归一化部分曲线下面积修正所述源的协方差矩阵；根据所述修正的源的协方差矩阵进行...

【技术特征摘要】
1.一种基于功能磁共振成像的时变约束脑电图或脑磁图溯源方法，步骤包括：记录脑电图或脑磁图信号并进行溯源，得到初始的脑电图或脑磁图源强度值，根据所述初始的脑电图或脑磁图源强度值得到脑电图或脑磁图源强度矩阵；记录功能磁共振成像数据并据以得到功能磁共振成像激活图，根据所述功能磁共振成像激活图得到功能磁共振成像权重矩阵；根据所述功能磁共振成像权重矩阵和所述脑电图或脑磁图源强度矩阵获得源的协方差矩阵，将该源的协方差矩阵设置为随时间变化；根据所述初始的脑电图或脑磁图源强度值计算交互通讯效应和归一化部分曲线下面积，并根据该交互通讯效应和归一化部分曲线下面积修正所述源的协方差矩阵；根据所述修正的源的协方差矩阵进行溯源，得到最终的脑电图或脑磁图源强度值，完成溯源。2.根据权利要求1所述的方法，其特征在于，根据所述脑电图或脑磁图信号、所述修正的源的协方差矩阵进行溯源所采用的线性方程为：其中，是脑电图或脑磁图源强度；R是源的协方差矩阵；L是引导场矩阵；C是噪声协方差矩阵，根据刺激前基线数据来确定；x是记录的脑电图或脑磁图信号。3.根据权利要求2所述的方法，其特征在于，根据所述脑电图或脑磁图信号进行溯源时，使用最小二次模算法、深度加权最小二次模算法得到所述初始的脑电图或脑磁图源强度值。4.根据权利要求1所述的方法，其特征在于，根据所述初始的脑电图或脑磁图源强度值得到所述脑电图或脑磁图源强度矩阵所采用的公式为：其中，是初始的脑电图或脑磁图源强度值；Re是脑电图或脑磁图源强度矩阵。5.根据权利要求1所述的方法，其特征在于，根据所述功能磁共振成像权重矩阵和所述脑电图或脑磁图源强度矩阵获得所述源的协方差矩阵所采用的公式为：R(t)＝RfRe(t)其中，R是源的协方差矩阵；Rf是功能磁共振成像权重矩阵；Re是脑电图或者脑磁图源强度矩阵。6.根据权利要求1所述的方法，其特征在于，所述交互通讯效应的计算公式为：

【专利技术属性】
技术研发人员：高家红，徐静，
申请(专利权)人：北京大学，
类型：发明
国别省市：北京,11