《基于功能磁共振成像的时变约束脑电图或脑磁图溯源方法.pdf》由会员分享，可在线阅读，更多相关《基于功能磁共振成像的时变约束脑电图或脑磁图溯源方法.pdf（11页珍藏版）》请在专利查询网上搜索。
1、(19)中华人民共和国国家知识产权局 (12)发明专利申请 (10)申请公布号 (43)申请公布日 (21)申请号 201710576924.X (22)申请日 2017.07.14 (71)申请人 北京大学 地址 100871 北京市海淀区颐和园路5号北 京大学 (72)发明人 高家红 徐静 (74)专利代理机构 北京君尚知识产权代理事务 所(普通合伙) 11200 代理人 刘召民 (51)Int.Cl. A61B 5/055(2006.01) A61B 5/0488(2006.01) (54)发明名称 基于功能磁共振成像的时变约束脑电图或 脑磁图溯源方法 (57)摘要 本发明提供一种基于功。
2、能磁共振成像的时 变约束脑电图或脑磁图溯源方法， 步骤包括： 记 录脑电图或脑磁图信号并进行溯源， 得到初始的 脑电图或脑磁图源强度值， 并据以得到脑电图或 脑磁图源强度矩阵； 记录功能磁共振成像数据并 据以得到功能磁共振成像激活图， 进而得到功能 磁共振成像权重矩阵； 根据所述功能磁共振成像 权重矩阵和所述脑电图或脑磁图源强度矩阵获 得源的协方差矩阵， 将该源的协方差矩阵设置为 随时间变化； 根据所述初始的脑电图或脑磁图源 强度值计算交互通讯效应和归一化部分曲线下 面积， 并据以修正所述源的协方差矩阵； 根据所 述修正的源的协方差矩阵进行溯源， 得到最终的 脑电图或脑磁图源强度值， 完成溯源。
3、。 权利要求书2页 说明书5页 附图3页 CN 107550493 A 2018.01.09 CN 107550493 A 1.一种基于功能磁共振成像的时变约束脑电图或脑磁图溯源方法， 步骤包括： 记录脑电图或脑磁图信号并进行溯源， 得到初始的脑电图或脑磁图源强度值， 根据所 述初始的脑电图或脑磁图源强度值得到脑电图或脑磁图源强度矩阵； 记录功能磁共振成像数据并据以得到功能磁共振成像激活图， 根据所述功能磁共振成 像激活图得到功能磁共振成像权重矩阵； 根据所述功能磁共振成像权重矩阵和所述脑电图或脑磁图源强度矩阵获得源的协方 差矩阵， 将该源的协方差矩阵设置为随时间变化； 根据所述初始的脑电图或。
4、脑磁图源强度值计算交互通讯效应和归一化部分曲线下面 积， 并根据该交互通讯效应和归一化部分曲线下面积修正所述源的协方差矩阵； 根据所述修正的源的协方差矩阵进行溯源， 得到最终的脑电图或脑磁图源强度值， 完 成溯源。 2.根据权利要求1所述的方法， 其特征在于， 根据所述脑电图或脑磁图信号、 所述修正 的源的协方差矩阵进行溯源所采用的线性方程为： 其中，是脑电图或脑磁图源强度； R是源的协方差矩阵； L是引导场矩阵； C是噪声协方 差矩阵， 根据刺激前基线数据来确定； x是记录的脑电图或脑磁图信号。 3.根据权利要求2所述的方法， 其特征在于， 根据所述脑电图或脑磁图信号进行溯源 时， 使用最小。
5、二次模算法、 深度加权最小二次模算法得到所述初始的脑电图或脑磁图源强 度值。 4.根据权利要求1所述的方法， 其特征在于， 根据所述初始的脑电图或脑磁图源强度值 得到所述脑电图或脑磁图源强度矩阵所采用的公式为： 其中，是初始的脑电图或脑磁图源强度值； Re是脑电图或脑磁图源强度矩阵。 5.根据权利要求1所述的方法， 其特征在于， 根据所述功能磁共振成像权重矩阵和所述 脑电图或脑磁图源强度矩阵获得所述源的协方差矩阵所采用的公式为： R(t)RfRe(t) 其中， R是源的协方差矩阵； Rf是功能磁共振成像权重矩阵； Re是脑电图或者脑磁图源强 度矩阵。 6.根据权利要求1所述的方法， 其特征在于。
6、， 所述交互通讯效应的计算公式为： 其中， 是交互通讯效应； G是线性逆算子； L是引导场矩阵； 所述归一化部分曲线下面积的计算公式为： 权 利 要 求 书 1/2 页 2 CN 107550493 A 2 其中， 是归一化部分曲线下面积； 是第i个位置得到的初始的脑电图或脑磁图源强度 值； q是通过所述初始的脑电图或脑磁图源强度值的T检验来确定的向量。 7.根据权利要求6所述的方法， 其特征在于， 所述向量q通过对任务状态和基线状态下 的初始的脑电图或脑磁图源强度值做T检验， 选用统计检验参数p值和多重比较校正方法； 所述统计检验参数p0.05， 所述多重比较校正方法采用邦费罗尼校正。 8.。
7、根据权利要求6所述的方法， 其特征在于， 根据所述计算交互通讯效应和所述归一化 部分曲线下面积修正所述源的协方差矩阵的方法是， 根据所述计算交互通讯效应和所述归 一化部分曲线下面积修正所述功能磁共振成像权重矩阵， 进而修正所述源的协方差矩阵， 采用的公式为： R(t)R, Re(t) 其中， R是源的协方差矩阵； R, 是修正的功能磁共振成像权重矩阵， 由交互通讯效应 和归一化部分曲线下面积 两个参数来确定。 9.根据权利要求8所述的方法， 其特征在于， 针对所述修正的功能磁共振成像权重矩阵 R, ， 根据所述交互通讯效应和所述归一化部分曲线下面积的计算公式计算在功能磁共振 成像激活区域内的平。
8、均值 ave和 ave， 在功能磁共振成像激活区域外的延展的激活区域应满 足条件： i ave并且 i ave。 10.根据权利要求9所述的方法， 其特征在于， 根据所述功能磁共振成像激活图得到所 述功能磁共振成像权重矩阵的方法是， 所述功能磁共振成像权重矩阵的对角元素在功能磁 共振成像激活区域和延展的激活区域内则置为1， 其余位置则置为0.1， 非对角元素设置为 0。 权 利 要 求 书 2/2 页 3 CN 107550493 A 3 基于功能磁共振成像的时变约束脑电图或脑磁图溯源方法 技术领域 0001 本发明涉及脑电图、 脑磁图及磁共振成像领域， 尤其涉及一种基于功能磁共振成 像的时变。
9、约束脑电图或脑磁图溯源方法。 背景技术 0002 脑电图(Electroencephalography， EEG)或脑磁图(Magnetoencephalography， MEG)有很高的时间分辨率， 可以达到毫秒量级， 并且他们的信号主要来源于突触后电位， 是 对大脑活动的直接测量。 但是由于脑电图电极数或脑磁图传感器数目远远小于大脑内电流 源的数目， 因此想要唯一地确定脑神经电流源的位置需要加入外加的约束， 而这些外加的 约束往往是缺乏生理学依据的。 采用脑电图或脑磁图信号研究脑神经活动， 确定神经电流 源的逆问题算法中具有代表性的如最小二次模算法(Minimum Norm Estimat。
10、ion， MNE)、 深 度加权最小范数法(depth-Weighted Minimum Norm， dWMN)等， 但这些方法始终不能摆脱空 间分辨率低的问题。 0003 基于血氧水平依赖(Blood-Oxygen-Level-Dependence， BOLD)的功能磁共振成像 (functional Magnetic Resonance Imaging， fMRI)方法自1992年Ogawa等小组在人体上成 功后得到了空前的发展。 BOLD功能磁共振成像现在已经成为认知神经科学中重要的研究手 段。 由于血氧动力学效应， 功能磁共振成像信号的时间分辨率低于脑电图或脑磁图信号， 但 是功能磁共。
11、振成像的空间分辨率要优于脑电图或者脑磁图， 可以达到毫米的量级。 0004 采用功能磁共振成像约束脑电图或脑磁图溯源方法确定大脑内电流源位置， 不但 有效克服了脑电图或者脑磁图数据空间分辨率不足的问题,而且在分析中同时引入了神经 活动电磁信息和代谢信息， 提高了分析结果的可信度。 0005 1998年， liu等第一次将功能磁共振成像的激活结果用于约束脑电图或脑磁图的 源定位， 并提出了功能磁共振成像约束最小二次模算法(fMRI-weighted minimum norm estimation， fMNE)。 但是由于脑电图或脑磁图信号反映的是大脑神经活动的电磁特征， 功 能磁共振成像信号反映。
12、的是大脑神经活动的血氧代谢信息， 脑电图或脑磁图电流源和功能 磁共振成像激活区域间关系并不可知。 用功能磁共振成像约束脑电图或者脑磁图溯源的这 种方法依赖于两种假设： 1、 功能磁共振成像激活区位置存在可被脑电或者脑磁记录到的神 经活动； 2、 神经活动可以引起对应的位置产生可探测到的血氧信号。 但是并不是所有情况 都满足这两种假设， 举例说明： 如果一个区域内有血氧代谢信号， 但是不能记录到对应的电 磁信号， 这种情况可能是因为神经活动产生的是封闭场， 故而违反了第一个假设， 会造成功 能磁共振成像额外源的存在； 如果能在头皮记录到某个神经活动产生的电磁信号但是却测 量不到血氧代谢信号， 这。
13、种情况可能因为这些神经元放电时间太短而造成很少的血氧消 耗， 这种情况违反了第二个假设， 可能会造成功能磁共振成像丢失源的存在。 0006 因此， 我们需要一种更好的方法来解决两者位置不匹配所产生的溯源结果不准确 的问题。 说 明 书 1/5 页 4 CN 107550493 A 4 发明内容 0007 为了解决上述不足， 本发明提出一种基于功能磁共振成像的时变约束脑电图或脑 磁图溯源方法， 将源的协方差矩阵认为是随时间变化， 以此为先验条件来约束脑电图或脑 磁图溯源， 本方法适用于所有的功能磁共振成像约束脑电图或脑磁图溯源分析。 0008 为了解决上述技术问题， 本发明采用以下技术方案： 0。
14、009 一种基于功能磁共振成像的时变约束脑电图或脑磁图溯源方法， 步骤包括： 0010 记录脑电图或脑磁图信号并进行溯源， 得到初始的脑电图或脑磁图源强度值， 根 据所述初始的脑电图或脑磁图源强度值得到脑电图或脑磁图源强度矩阵； 0011 记录功能磁共振成像数据并据以得到功能磁共振成像激活图， 根据所述功能磁共 振成像激活图得到功能磁共振成像权重矩阵； 0012 根据所述功能磁共振成像权重矩阵和所述脑电图或脑磁图源强度矩阵获得源的 协方差矩阵， 将该源的协方差矩阵设置为随时间变化； 0013 根据所述初始的脑电图或脑磁图源强度值计算交互通讯效应(cross talk effect)和归一化部分。
15、曲线下面积(normalized partial area under curve matrix)， 并 根据该交互通讯效应和归一化部分曲线下面积修正所述源的协方差矩阵； 0014 根据所述修正的源的协方差矩阵进行溯源， 得到最终的脑电图或脑磁图源强度 值， 完成溯源。 0015 进一步地， 根据所述脑电图或脑磁图信号、 所述修正的源的协方差矩阵进行溯源 所采用的线性方程为： 0016 0017其中，是脑电图或脑磁图源强度； R是源的协方差矩阵； L是引导场矩阵； C是噪声 协方差矩阵， 根据刺激前基线数据来确定； x是记录的脑电图或脑磁图信号。 其中R和C是未 知量， x和L是已知量。 00。
16、18 进一步地， 根据所述脑电图或脑磁图信号进行溯源时， 使用最小二次模算法、 深度 加权最小二次模算法得到所述初始的脑电图或脑磁图源强度值。 0019 进一步地， 根据所述初始的脑电图或脑磁图源强度值得到所述脑电图或脑磁图源 强度矩阵所采用的公式为： 0020 0021其中，是初始的脑电图或脑磁图源强度值； Re是脑电图或脑磁图源强度矩阵。 0022 进一步地， 根据所述功能磁共振成像权重矩阵和所述脑电图或脑磁图源强度矩阵 获得所述源的协方差矩阵所采用的公式为： 0023 R(t)RfRe(t) (3) 0024 其中， R是源的协方差矩阵； Rf是功能磁共振成像权重矩阵； Re是脑电图或者。
17、脑磁图 源强度矩阵。 说 明 书 2/5 页 5 CN 107550493 A 5 0025 进一步地， 所述交互通讯效应的计算公式为： 0026 0027 其中， 是交互通讯效应； G是线性逆算子(inverse operator)； L是引导场矩阵。 0028 所述归一化部分曲线下面积的计算公式为： 0029 0030其中， 是归一化部分曲线下面积； 是第i个位置得到的初始的脑电图或脑磁图源 强度值； q是通过所述初始的脑电图或脑磁图源强度值的T检验来确定的向量。 0031 进一步地， 所述向量q通过对任务状态和基线状态下的初始的脑电图或脑磁图源 强度值做T检验， 选用统计检验参数p值和多。
18、重比较校正方法。 0032 进一步地， 所述统计检验参数p0.05， 所述多重比较校正方法采用邦费罗尼校正 (Bonferroni corrected)。 0033 进一步地， 根据所述计算交互通讯效应和所述归一化部分曲线下面积修正所述源 的协方差矩阵的方法是， 根据所述计算交互通讯效应和所述归一化部分曲线下面积修正所 述功能磁共振成像权重矩阵， 进而修正所述源的协方差矩阵， 采用的公式为： 0034 R(t)R, Re(t) (6) 0035 其中， R是源的协方差矩阵； R, 是修正的功能磁共振成像权重矩阵， 由交互通讯效 应 和归一化部分曲线下面积 两个参数来确定。 0036 进一步地，。
19、 针对所述修正的功能磁共振成像权重矩阵R, ， 首先根据所述交互通讯 效应和所述归一化部分曲线下面积的计算公式计算在功能磁共振成像激活区域内的平均 值 ave和 ave， 然后在功能磁共振成像激活区域外满足条件： i ave并且 i ave， 满足此条 件的区域定义为延展的激活区域。 0037 进一步地， 根据所述功能磁共振成像激活图得到所述功能磁共振成像权重矩阵的 方法是， 所述功能磁共振成像权重矩阵的对角元素在实际功能磁共振成像激活区域和延展 的激活区域内则置为1， 其余位置则置为0.1， 非对角元素设置为0。 0038 本方法通过功能磁共振成像约束脑电图或脑磁图信号溯源， 由功能磁共振成。
20、像权 重矩阵和脑电图或脑磁图源强度矩阵两部分共同决定源的协方差矩阵， 能够解决功能磁共 振成像额外源的影响， 并且在一定程度上降低了功能磁共振成像丢失源的影响； 通过交互 通讯效应和归一化部分曲线下面积修正功能磁共振成像权重矩阵， 进而修正源的协方差矩 阵， 能够更好地解决功能磁共振成像丢失源的影响。 在上述的约束条件下， 通过消除两方面 的影响， 达到溯源结果的准确定位。 本方法适用于不论神经电活动位置和功能磁共振成像 激活区域位置的匹配或者不匹配等所有情况。 附图说明 0039 图1是本发明的溯源方法的流程图。 0040 图2是单侧视觉刺激下得到的事件相关电位曲线图。 说 明 书 3/5 。
21、页 6 CN 107550493 A 6 0041 图3是单侧视觉刺激下得到的功能磁共振激活图。 0042 图4是MNE、 fMNE和FITC对图2中脑电信号溯源的溯源结果示意图。 具体实施方式 0043 为使本发明的上述特征和优点能更明显易懂， 下文特举实施例， 并配合所附图作 详细说明如下。 0044 应指出的是， 本发明公开的仅是具有代表性的一种较佳实施例。 显然， 本发明并不 局限于所描述的具体结构、 功能、 器件和方法， 也可以具有其他实施方式， 或者是其他实施 方式的组合。 本发明中所描述的元素数目也可以设置为多个， 除非明确限制为单数。 此外， 为避免其他例与本发明发生混淆， 对。
22、于本领域中众所周知的一些技术特征和细节未进行描 述。 0045 本实施例提出了一种基于功能磁共振成像的动态约束脑电图或脑磁图溯源方法， 亦称之为基于功能磁共振成像的时变约束法(fMRI informed time-variant constrains， FITC)。 本例以单侧棋盘格视觉刺激下得到的脑电图或脑磁图信号和功能磁共振成像数据 为例， 详细说明溯源的具体方法， 如图1所示。 值得强调的是， 该方法并不局限于视觉刺激， 其他类型刺激任务可以按照相似的步骤处理。 0046 脑电图或脑磁图信号和大脑内电流源的关系可以用线性方程(1)来进行溯源， 利 用最小二次模算法或其他溯源算法得到初始的。
23、脑电图或脑磁图电流源强度值根据功能 磁共振成像数据利用现有技术得到功能磁共振成像激活图， 再进而根据功能磁共振成像激 活图得到功能磁共振成像权重矩阵Rf。 由于功能磁共振成像激活区域位置与脑电图或脑磁 图电流源位置可能会出现不匹配的情况， 为了解决这样由于这种不匹配情况产生的溯源问 题， 我们加入了一些约束条件。 0047 首先解决功能磁共振成像额外源的影响， 因为电流源的活动是随着时间发生变化 的， 它在某一时间的活动位置可能只是功能磁共振成像激活区域的一个子集， 因此在这一 时间点其他功能磁共振成像激活区域则表现为功能磁共振成像额外的源。 为了解决这一问 题， 本方法的源的协方差矩阵R由功。
24、能磁共振成像权重矩阵Rf和脑电图或脑磁图源强度矩 阵Re两部分决定， 并且该源的协方差矩阵R是随时间变化的。 从公式(3)就可以看出， 本方法 主要强调了那些在功能磁共振成像激活区内并且神经电流源强度很高的位置， 所以可以解 决功能磁共振成像额外源的影响， 并且在一定程度上降低了功能磁共振成像丢失源的影 响。 0048 为了更好地解决功能磁共振成像丢失源的影响， 功能磁共振成像权重矩阵Rf通过 公式(4)和公式(5)中的参数 和 来决定。 首先计算在功能磁共振成像激活区域内的两个参 数的平均值 ave和 ave， 然后在功能磁共振成像激活区域外满足： i大于 ave并且 i小于 ave， 满足。
25、这两个条件的位置定义为延展的激活区域。 修正的功能磁共振成像权重矩阵R, 的对 角线元素在功能磁共振成像激活区域和延展的激活区域内设置为1， 其余设置为0.1， 非对 角元素设置为0。 0049 在计算归一化部分曲线下面积 时， 通过对任务状态和基线状态下的初始的脑电 图或脑磁图源强度值做T检验， p0.05， 邦费罗尼校正。 此处需要说明的是， 选用统计检验 参数p值和多重比较校正方法不一定非要选择0.05和邦费罗尼校正， 只要是计算向量q时采 说 明 书 4/5 页 7 CN 107550493 A 7 用统计检验的方式， 都在本发明的保护范围内。 0050 将经过修正后的源的协方差矩阵R。
26、, 再重新带入公式(1)中进行溯源， 得到最终的 脑电图或脑磁图源强度值完成溯源。 0051 本发明所述的脑电图或脑磁图信号具体是指脑电图事件相关电位(event related potential， ERP)或脑磁图事件相关场(event related field， ERF)， 图2是受试者 在2赫兹棋盘格的刺激下得到的脑电事件相关电位曲线图， 该视觉刺激只出现在视野的右 下侧。 从图中可以清楚看到与视觉刺激相关的N75、 P100和N135等成分， 这些成分是脑电波 的峰， 大概发生在69毫秒、 107毫秒和128毫秒。 0052 图3是受试者在相同视觉刺激下得到的功能磁共振激活图， 从。
27、图中也可以看到包 括初级视觉皮层V1以及V2、 V3区域的激活。 这些区域是人类视觉系统最常见的区域， 现有技 术中的相关研究指出在观看棋盘格刺激时会导致这三个区域的激活， 本图结果与相关研究 结论一致， 由此可见本方法是可行的。 0053 图4是对图2中脑电信号用三种方法在三个时间点下的溯源得到的源强度值的示 意图， 第一行是用最小二次模算法(MNE)得到的溯源结果， 虽然可见视觉区的激活， 但是激 活区域很宽泛， 定位准确性差； 第二行是用功能磁共振成像约束最小二次模算法(fMNE)溯 源得到的溯源结果， 相对于最小二次模算法来说， 功能磁共振成像约束最小二次模算法得 到的激活区更加集中，。
28、 但是在128ms时出现了错误的激活区。 第三行是本发明阐述的方法 (FITC)得到的溯源结果， 相对于最小二次模算法和功能磁共振成像约束最小二次模算法， 利用本发明的方法可以得到更好的定位准确性。 0054 综上所述， 利用本发明的方法处理功能磁共振成像约束脑电图溯源结果能够得到 更好的溯源结果。 此外， 本方法适用于不论神经电活动位置和功能磁共振成像激活区域位 置的匹配或者不匹配等所有情况， 根据同步记录脑电图或者脑磁图和功能磁共振成像数据 的广泛应用范围， 本发明还可运用到认知神经科学中(比如语言、 记忆等)以及精准医疗或 者与大脑活动相关的个体差异研究中。 说 明 书 5/5 页 8 CN 107550493 A 8 图1 说 明 书 附 图 1/3 页 9 CN 107550493 A 9 图2 图3 说 明 书 附 图 2/3 页 10 CN 107550493 A 10 图4 说 明 书 附 图 3/3 页 11 CN 107550493 A 11 。