CIN 计算智能和神经科学 1687 - 5273 1687 - 5265 Hindawi出版公司 656092年 10.1155 / 2009/656092 656092年 评论文章 脑电图/梅格源成像:方法,挑战,和开放的问题 Wendel的 卡特里娜飓风 1 Vaisanen Outi 1 Malmivuo 1 Gencer Nevzat G。 2 Vanrumste 巴特 3、4 4 Durka 彼得亚雷 5 Magjarević 6 Supek 塞尔玛 7 Pascu Mihail卢西恩 8 Fontenelle Hugues 9 严重的德佩拉尔塔梅内德斯 罗兰多 10、11 11 Astolfi 劳拉 1 生物医学工程系 坦佩雷理工大学 33101年坦佩雷 芬兰 tut.fi 2 电气电子工程系 中东技术大学 06531年土耳其 土耳其 metu.edu.tr 3 谈到工程系(计) 鲁汶Katholieke项目 Kasteelpark阿伦伯格10 b - 3001鲁汶 比利时 kuleuven.be 4 MOBILAB、生物科学和技术部门 Katholieke Hogeschool Kempen Kleinhoefstraat 4 b - 2440 Geel 比利时 khk.be 5 生物物理系 华沙大学ul。Hoża 69 00 - 681年华沙 波兰 uw.edu.pl 6 学院电气工程和计算 萨格勒布大学 10000年萨格勒布 克罗地亚 unizg.hr 7 物理系 理学院 萨格勒布大学 10000年萨格勒布 克罗地亚 unizg.hr 8 国家激光研究所 等离子体物理和辐射,激光 077125年布加勒斯特 罗马尼亚 inflpr.ro 9 医学物理学部门 医学院的 佩特雷大学 26504年佩特雷 希腊 upatras.gr 10 电神经影像组 临床神经科学部门 日内瓦大学医院 1211年日内瓦 瑞士 hug-ge.ch 11 神经动力学实验室 精神病学和临床心理生物学 巴塞罗那大学 08035年巴塞罗那 加泰罗尼亚 西班牙 ub.es 2009年 19 07年 2009年 2009年 25 11 2008年 31日 03 2009年 29日 04 2009年 2009年 版权©2009 这是一个开放的文章在知识共享归属许可下发布的,它允许无限制的使用,分布和繁殖在任何媒介,提供最初的工作是正确的引用。

我们提出的四个关键领域research-preprocessing,容积导体,逆问题的解决提出的问题,影响脑电图和梅格源成像的性能。我们确定在每一个关键领域著名的途径和方法,开放问题需要进一步调查在社区内,挑战与某些相关技术和算法需要澄清他们的影响。提供明确的答案我们的目标是识别多源定位的重要追求开放的问题。

1。介绍

脑电图(EEG)、脑磁图描记术(MEG)代表两种无创性脑功能成像方法,测量电势颅外记录的差异和极弱磁场产生的神经细胞的电活动,分别。这些录音提供直接、实时监测自发和诱发的大脑活动,让时空定位潜在的神经发电机。脑电图和梅格具有以下特点: ( 1 ) 他们有一个毫秒的时间分辨率; ( 2 ) 潜在的差异和磁场源的优点和非线性函数的线性函数的源支持(例如,偶极子位置); ( 3 ) 它们是由相同的神经生理学的事件,也就是说,电流从同步激活神经组织通常被称为初级或外加电流源密度,因此都可以同样用于神经发电机的本地化。

三维重建的神经活动通常被误解了 断层摄影术,这是定义( 1)“ 任何技术,使对象的详细的预先确定的平面部分而模糊的图像其他飞机。”物理治理电磁场的传播取决于容积导体的构成,这意味着源活动在预定的飞机也会影响传感器的读数躺在飞机。所以,实际的程序应用于断层和逆解逆转:在断层扫描我们重建三维图像结合分别获得片,然而,逆解计算整个三维分布,可后来切片。在层析技术(如CT、PET、核磁共振等)与适定的数学问题,大脑的无创性评估活动本质上是一个不适定问题的无数的解决方案。

在后续部分中我们对一些重要问题的理解的话,选择,和评价源成像方法;因此,我们强调的一般方法,而不是特定的解决方案。这些部分反思我们的小组讨论NeuroMath举行研讨会。作为一个团体我们承认我们有不同的意见关于方法的选择,我们面临的各种挑战作为单独的研究中心,和我们不同的关键问题,由于我们的不同的利益。因此,我们试图写一篇文章,新手,将源图像的不同部分社区,主要需要关注几个开放的问题。

1.1。理论

头内的来源之间的关系和外部测量 d 被描述为 d = l j , 在哪里 l 是代表领导领域的线性算子(也称为增益模型或直接模型),然后呢 j 代表了来源。两个数学性质( 1)反映了属性的物理大小。首先,同质性属性指出一个放大的图像源 k * j 是一个放大测量 k * d ,其次,附加属性两个来源的总和 j = j 1 + j 2 产生一个测量等于每个测量单独的总和 d = d 1 + d 2 。在这两个属性遵循叠加原理,也就是说, l ( k 1 * j 1 + k 2 * j 2 ) = k 1 * d 1 + k 2 * d 2 ,在那里 d 1 = l j 1 d 2 = l j 2

这个问题的不适定的本质来自于一个事实,就是两个不同的来源 j 1 j 2 可能会产生相同的测量 d ,也就是说, d = l j 1 = l j 2 ,这是非常相当于说,存在一个“沉默”的来源 h 这样 l h = 0 。为了看到等价,注意,如果 d = l j 1 = l j 2 因此,沉默的来源 h = j 1 - - - - - - j 2 满足 l h = l j 1 - - - - - - l j 2 = d - - - - - - d = 0 。如果我们假设,在另一个方向 l h = 0 和源的存在 j 1 这样 l j 1 = d ,我们可以建立一个新的来源 j 2 = j 1 + h 收益率相同的数据,即 l j 2 = l ( j 1 + h ) = l j 1 + l h = l j 1 = d

那就是说,我们可以建立一个脑电图/梅格头皮分布的主要属性(地图)。虽然类似的头皮地图不能排除的可能性在底下的来源分布不同,不同地图必然是由于不同的源分布。重要的是,我们不需要诉诸任何逆方法得出。

基于线性和在缺乏先验信息为借口,否则,我们可以代表的解决方案( 1)和一个线性算子 G “估计” j 如下: j 美国东部时间 = G * d 替换 d 中定义它的值( 1)产生一个线性算子的基本方程估计和原始源分布 j 美国东部时间 = G * l * j = R * j , 在哪里 R = G * l 是解析操作符( 2, 3]。在实践中传感器和几何形状是由离散的测量,因此可以假定 l , G , R 有限维矩阵近似连续(积分)操作符。

2。预处理

在本节中,我们讨论一些有关问题相关数据的准备识别一些有用的预处理和避免的事情。一般来说,预处理的哲学是准备解决的信号。通常,这些步骤分解复杂的信号和减少噪声传感器以及其他不良的来源。

(图的脑电图和梅格的逆问题 1绿色箭头)从时间序列(图开始 1在头皮上的传感器记录)。因此,头皮的本地化基于分布振幅在单一时间瞬间可能提高了信号处理的应用技术来测量时间序列(图 1蓝色箭头)。特别是输入噪声可以减少选择性和敏感提取相关活动的脑电图/梅格数据。这可以通过本地化组件提取信号的盲源分离(例如,ICA [ 4])。其他方法依赖于来自时频表示的信息,对应相关的现象我们要本地化(例如,睡眠纺锤波 5])。类似但更敏感和选择性预处理提出了在 6使用多道匹配追踪算法。总的来说,大多数预处理算法有望受益逆解的质量和准确性。

源成像的关键部分。预处理准备记录信号求解逆问题。反问题试图定位记录测量数据的来源,而提出的问题假设源定义为了计算电位分布地图。

2.1。时代

我们应该权衡利弊的时代扮演的角色与事件相关电位的录音。到目前为止还没有标准的试验,抖动,平均振幅,或单个试验的合理性分析。例如,信噪比(信噪比)与试验的数量增加,也就是说,时代的数量;然而,习惯可以影响一些研究的结果。我们建议一个文档概述这些类别将有利于未来的研究比较和正规化。

neuroelectric信号被埋在自发的脑电图与信噪比低至5分贝。为了降低噪音水平,找到一个 模板诱发电位(EP)信号,一个总体均值(EA)获得使用大量的重复测量( 7]。这种方法( 8将背景EEG视为附加噪声和EP为不相关的信号。的大小和延迟EP波形显示大型个人间的差异和变化取决于心理生理因素对于一个给定的个人( 9]。因此,一个方法论的EP研究目标是开发技术来提取的真正的EP波形 单扫描

临床评估,模板EP信号或可能的振幅和/或延迟变化单一使用扫描( 9]。观察这些变化,可以确定具体的功能从模板参考/ EP基于各种评估方法。因为有相对严格约束相关的记录时间或与主题,使用EA(作为参考EP信号)通常是不切实际的。这导致了另一种信噪比改善的发展基于加法模型方法。有些算法加权平均的方法,子空间平均方法,参数滤波,自适应滤波,维纳滤波。在所有这些方法中,假设EP(即。整个实验、信号)是静止的。然而,这种假设也在质疑一些报道描述的事件相关电位叠加相位调制有节奏的活动可能与大脑的不同认知过程( 10]。

2.2。事情,以避免

相反大多数预处理算法的好处,有一些算法之前,我们应该避免源定位算法的应用。特别是,以下选择威胁逆的完整性解决方案。

( 1 ) 基线校正。不同值的单个电极通过“任意”基线转移或扩展因素改变了表面地图,从而估计来源。虽然线性逆解,而稳定的(连续性的数据),基线校正的应用程序两个条件(相比,将来源)的基础上可以产生人工差异引起的修正而不是真正的来源。

( 2 ) 人工地图生产的总平均数据或分割算法。统计平均水平(例如,意味着)屈服值,通常没有出现在原始数据。,这并不让人感到意外,如果平均地图不存在任何主题的平均水平。此外,这种效应可以被放大的延迟差异科目。

( 3 ) 传感器的使用非常高的密度也可能危害源分析由于不同类型的噪声在不同的传感器。此外,没有明显的信息添加在大约128个电极由于噪音水平。最后,一些传感器可能会比其他人更工件测量由于其位置附近活跃的肌肉。

3所示。容积导体

头部模型作为一个容积导体在源定位是一个关键元素。容积导体的配置直接影响到解决方案和逆解的解决方案。三个几乎同样重要的区域是头几何,组织导率,和电极位置。

3.1。几何和分割

的开创性研究,德里斯科尔 11)使用三个同心球体,而当代研究实现现实的模型。我们发现在社区内不同的模型是必要的,但如何每种类型的几何模型有助于源定位的目标?球面模型回答一般问题的理论提供几厘米的脑电图定位精度,而现实的模型试图查明确切位置但实际上提高偶极子定位几厘米( 12- - - - - - 14]。另一方面,大多数MEG-based数值模拟的球面几何学是充分的。只有本地化的深源的头骨底部附近额颞叶和额叶区域需要一个realistically-shaped-head-volume MEG-based导体模型模拟( 15]。

几何成像形态直接相关,计算机断层扫描(CT)或磁共振成像(MRI)和分割的质量。自然,我们将问题形式segment-CT, MRI或融合CT-MRI图片( 16]。我们鼓励建模者理解的意义边界定义为一个特定的模式,不以任何方式暗示医学界提供任何病人不安全的和不必要的辐射。我们必须记住,我们选择影响分割的形式由于其灵敏度相应困难或软组织。此外,有多少组织,组织,和蛀牙的模型应该包括?我们预见到我们正在接近一个高原定位精度的改善作为细分解决方案增加的包含太多的小组织区域。研究大道之一可能受益于head-model几何学的发展是人体测量和craniometric数据的集成 17]。这条路能证明个人模型代表一个族群和否定研究声称,歪曲一个确定分组人口统计学意义。此外,它可以引导我们建立的形状和大小的统计意义的个人和组内模型的几何特性。

3.2。电导率值

许多生物组织的电特性不均匀,各向异性,色散和非线性。头组织如头骨,头皮,肌肉,脑脊液,灰色和白色物质有不同的导率 σ ,介电系数 ε 和磁渗透率 μ (在大多数情况下,它被认为是等于水的渗透,这是接近自由空间的渗透性 μ o )。头骨以及头皮有多层结构,每一层具有不同的电特性。这一事实导致要么多层建模的几何组织( 18)或者将非齐次性质的组织,也就是说,张量分配 σ = σ ( x , y , z ) ε = ε ( x , y , z ) 。的值和分布是一个更严重的非均质性问题在病理过程的患者群体可能会显著影响导率影响的大脑区域。可能存在一个等效混合各向同性模型代表多个各向异性层?明显这样近似如何影响源定位在健康个体相比,患者头部病态?

任何模型的电导率值影响的主要领域提出问题和逆问题的解决方案。因此,它是至关重要的,我们必须分配尽可能准确的电导率值报告从以前的文献研究和推断和插入。作为一个社区建立了电性质对于大多数头部组织电导率的范围 σ 和介电常数 ε ;然而,我们必须确定个体的实际导电性分布的头。由于这些范围,许多历史研究值平均分配给他们的组织 15, 18- - - - - - 26]。使用平均值可能不是适合个性化模型,因为这些模型可能会导致不准确的解决方案由于位置的函数 27年龄的]或[ 28, 29日]。然而,研究患者( 30., 31日)表明,用近似导率比率的头部(即一个精确的几何描述。基于主题的MRI)可能产生合理,皮层和深部脑电图来源证实的结果。不过,未来的模型可能受益于使用特定导率。我们推测,年龄的电导率值的应用程序应用于tissues-most尤其是trilayer头骨组织主要受益青春的模型,其骨化中心在前两年迅速变化和高原电导率值在18到20岁的颅顶的骨化过程完成时( 32]。

为了巩固我们的动机凸显颅骨电导率的意义,我们必须简要地深入研究它的历史。的开创性工作 11]介绍了一个标准的电导率比1:80:1 brain-to-skull-to-scalp的,这是一个历史价值仍然被一些研究者在四十年后使用。在[ 33]在尸检报告测量尸体产生brain-to-skull电导率值的比例15:1。三年后( 28]给出了电导率值现场组织低至4:1的比例。随后,万德尔,Malmivuo 29日)相关 后期活组织测量作为一种整合和评价过去的数据由于缺乏在活组织测量。他们以前的工作使用的缩放比率0.33到0.4从生活适应电导率的变化 后期组织根据电导率的死亡记录组织样本( 34, 35]。摘要之前他们提出了一个开放的问题社区更测量现场组织samples-most尤其是过头骨样品正常的体温和水分,至今仍然作为一个开放的问题。因此,它是我们歧视的条件相关的组织电导率和介电常数的值,并将被收购。值通过体内或体外测量应优先于生活 后期测量。在的情况下 后期测量,采集的时间和温度应指定的自组织特性细胞死后迅速变化。

3.3。收购电导率值

在过去的二十年里,许多方法已经提出了人体形象的导电性。在传统应用当前电阻抗断层扫描(ACEIT) low-frequency-sinusoidal电流通过电极应用于身体表面( 36]。在感应电流电阻抗断层扫描(ICEIT)时变磁场不同spatial-field模式应用于体内诱导电流。在这两种情况下,表面电极用于使电压测量。

最近,提出了两种新方法,利用磁测量电导率分布。在磁感应成像(MII),发射线圈是由正弦电流提供时变磁场( 37, 38]。当附近的身体带来这些线圈时,涡流诱导体内。这些洋流的分布是身体的电导率分布的函数。这些电流产生了二次磁场,测量接收机中的电动势感应线圈。在磁共振电阻抗断层成像(EIT)先生,低频电流应用从身体表面,以及由此产生的磁场测量系统[先生使用一个 39, 40]。因为体内磁场测量,可以测量高分辨率图像。注意,所有方法仍在调查阶段,和没有人可以提供高分辨率的要求电导率源定位所需信息。

3.4。电极蒙太奇

脑电图是传统测量使用标准的10 - 20电极系统包括只有21个测量电极。它被广泛承认的空间分辨率现代大脑研究[10 - 20系统是不够的 41- - - - - - 44]。提高脑电图的空间分辨率的第一步是提高脑电图电极的数量,市场反应与商用系统包括256个电极。

在过去的二十年中几项研究调查增加脑电图电极的数量的好处。对精度的影响提出解决方案和逆解决方案评估。在几篇文章,增加电极的数量至少有128已经被证明可以提高结果的准确性( 45- - - - - - 50]。

不同的因素影响到适当数量的电极。这些措施包括,例如,头骨的广泛讨论价值的相对电导率,这对逆解的准确性有很大的影响。此外,尤其是密集的脑电图系统的空间分辨率(128 - 512)电极测量噪声非常敏感。因此,对于不同的脑电图测量在不同的环境中进行的,适当的电极数量可能相差很大( 48]。使用活性电极将减少噪音。

4所示。提出问题

1969年,德里斯科尔的研究( 11在脑电图电极的灵敏度又迎来了源定位的新时代。他们的工作分析解决了麦克斯韦方程映射的字段,它只可能至少有椭圆对称。当代模型由一个组合的复杂几何和/或电气参数,从而迫使数字解决方案,比如边界元法(BEM),有限元方法(FEM),有限差分法(FDM)(表 1)。在本节中,我们的目标是识别的一些并发症,这些数值方法的利弊。我们希望读者通过以下解释收益的理解差异,采用一个或多个特定于他/她的需求,适当的方法,指的是引用的特定信息。

比较四种方法求解泊松方程的真实头模型:边界元法(BEM),有限元方法(FEM),各向同性有限差分法(iFDM),和各向异性有限差分方法(aFDM)。

有限元法 iFDM aFDM
计算点的位置 表面 体积 体积 体积
自由选择的计算点 是的 是的 没有 没有
系统矩阵 完整的 稀疏的 稀疏的 稀疏的
解决者 直接/迭代 迭代 迭代 迭代
数量的隔间
需要镶嵌 是的 是的 没有 没有
处理各向异性 没有 是的 没有 是的

大多数模特都是无法获得直接的解决方案,所以他们依靠迭代解决如逐次超松弛(SOR)、共轭梯度(CG)、预处理共轭梯度法(PCG),和解决代数多重网格(AMG)。虽然这些方法已经开发了正则线性系统,也可以应用在我们的半定的情况。在一致的右边,semiconvergence可以保证对法师和CG (P),而AMG理论结果更加复杂( 51]。摘要给出了每种方法基于[ 52)第一三种方法( 53, 54最后一个方法)。

第一本和有限元法的区别或FDM计算域的解决方案。本解决方案中计算均匀各向同性隔间之间的界限在有限元方法和有限差分法计算提出问题的解决方案在整个体积。随后,有限元方法和FDM导致比本更大数量的计算点。另一方面,可以确定任意点的电位与有限元法和有限差分法的插值计算点的附近,而对于本有必要重新应用巴纳德公式( 55)和数值积分。

另一个重要方面是计算效率。本,一个完整的矩阵( - - - - - - C ),在 V = C V + V 0 , 需要倒。当头皮电位需要知道另一个偶极子, V 0 在( 4)需要重新计算并乘以已经可用 ( - - - - - - C ) - - - - - - 1 。因此一旦倒矩阵,只需要一个矩阵乘法来获取头皮电位。这个有限的计算负载是一个有吸引力的特性解决逆问题时,在大量的评估需要执行。另外,加速BEM方法增加速度明显只通过计算 (例如,the number of electrodes) rows of the corresponding inverse, whereas, the normal inversion process requires a lot more time due to the dimensionality of the matrix as n × n (例如, n 等于节点的数量) 56, 57]。投影的方法( 58)参数表示的基础上的表面也允许大幅减少计算负荷。

有限元法和FDM,直接反演的大型稀疏矩阵是不可能由于矩阵的维数。通常至少有500 000个计算点被认为是因此导致系统矩阵与500 000 500 000方程的未知数,而不能以直接的方式解决当前可用的电脑。然而,矩阵中有限元方法和有限差分法可以为一个给定源配置或倒右边项,利用迭代解决如逐次超松弛法(SSOR)、共轭梯度(CG)方法( 59),或代数多重网格(AMG)方法 60, 61年]。迭代解决的缺点是,对于每个源配置解决者必须重新应用。有限元方法和有限差分法将计算效率低下,当迭代解算器需要用于每一个偶极子。为了克服这种低效率使用互易定理( 62年]。

介绍了大量进行隔间时,需要大量的边界采样本。这导致大量完整的系统矩阵,从而降低计算效率。在有限元方法和有限差分法建模,真实头模型的异构特性会降低刚度矩阵稀疏,严重制约。此外,整合各向异性导率将降低刚度矩阵的稀疏。这可能会导致一个不稳定的系统或如果使用迭代方法收敛非常缓慢。获得一个快速收敛或一个稳定的系统,应该使用预处理。预处理变换的系统方程 一个 x = b 预处理系统 - - - - - - 1 一个 x = - - - - - - 1 b 有相同的解决方案及系统。 是一个预处理矩阵或预调节器,其目的是为了减少条件数(比例最大的特征值最小的特征值)的刚度矩阵的最优值1。基本预处理的形式可以使用雅可比,高斯-赛德尔,逐次超松弛(SOR),对称连续超松弛(SSOR)。这些很容易实现 63年]。更先进的方法使用不完全LU分解和多项式预处理[ 63年, 64年]。

对FDM BEM和有限元法相比,计算点谎言固定在各向同性的立方体中心方法和多维数据集上各向异性的方法。在有限元方法和边界元计算点,四面体的顶点和三角形,分别可以更自由地选择。因此,有限元法能更好地代表了比FDM不规则之间的接口不同的隔间,相同数量的节点。然而,分割医学图像用于获得现实的容积导体模型构造出立方体素。直接生成结构化网格用于FDM从这些分割图像。有限元法和边界元,额外的镶嵌算法( 65年)需要使用获取四面体元素和表面三角形,分别虽然立方和直角棱镜元素可能在像FDM有限元法。

最后,众所周知,一些组织在人类头上的导率各向异性头骨和白质组织等。各向异性可以引入有限元法( 66年]在FDM [ 67年),但不是本。

5。反问题

虽然前面部分集中在不同步骤前逆的应用程序,也就是说,头部几何近似,电导率,几何轮廓,电导率值的准确性,等等,本节讨论一些开放的问题,包括记录方式的选择、源模型,可能后期处理来提高逆解估计的鲁棒性。

5.1。记录方式:梅格和脑电图

超导量子干涉器件的引入(鱿鱼)成为可能测量非常低磁场引起大脑的电活动,称为脑磁图描记术,梅格。

在生物磁效应研究的开始,有很多希望生物磁效应信号包括信息独立的生物信号。如Plonsey所述 68年),这一事实根据亥姆霍兹定理标量和矢量势字段可以选择独立被认为是作为证据的独立电和磁测量。另一方面,考虑到生物的起源电流得出散度和旋度的主要电流不能任意分配。进一步的实验中描述( 68年)指出,有关贡献电场和磁场辅助源文档。因此,虽然我们不能声称的生物或生物磁效应字段就足以定义其他 69年),我们不应指望重要差异的信息记录。

结论,电和磁测量提供可比的信息已被证实的理论和仿真。使用灵敏度一半体积的新概念,Malmivuo et al。 70年]表明,脑电图和梅格记录电活动以一种非常相似的方式,也就是说,脑电图之间的差异和梅格敏感性一半量的大小和形式的敏感性分布非常小。进一步评估皮质来源的空间分辨率脑电图的球形模型导致更好的结果( 71年]。使用模拟[ 72年]证实也略有优势的脑电图在许多源位置和方向与最好的结果结合脑电图/梅格测量。最近[ 73年)模式识别技术应用于解码手运动方向同时脑电图/梅格测量,得出结论的推理运动方向对这两种技术也同样有效。

因此,它可能是有益的考虑也影响成本的记录方式。梅格仪器成本约20倍的脑电图仪器相同数量的渠道。因此,为提高逆解的准确性可能首先有利于提高脑电图技术的各个方面,也就是说,许多渠道,电极位置的准确性,头部几何模型,和组织电阻率的准确性,等等,因为改善所有这些成本远低于梅格仪器。

总之我们可以证实读者,除了成本差异,这两种技巧提供类似大脑信息来源对源定位精度,时空分辨率,和解码或预测能力。我们想强调,虽然类似的信息发现,EEG和梅格测量敏感性是正交的。EEG主要检测电力来源与足够遥远的径向头皮表面电极和切向分量时,导致彼此附近(图 2);然而,梅格主要感官磁电流所产生的电在径向方向(图源 3)。

脑电图的灵敏度分布。(左)脑电图设置测量的切向分量neuroelectrical活动,每个双领先位置相对较近。(右)脑电图设置测量neuroelectric活动的径向分量,测量电极的位置远离参比电极。箭头在这两个数字代表macrocolumns蜂窝架构不是偶极源。

梅格的灵敏度分布。(左)一个梅格设置测量的切向分量neuroelectrical活动,使用一个轴向梯度仪。(右)一个梅格设置测量的切向分量neuroelectric活动,使用平面梯度仪。箭头在这两个数字代表macrocolumns蜂窝架构不是偶极源。

5.2。源模型

有巨大的文献回顾的方法用于解决所谓的bioelectromagnetic-inverse问题处理的估计电活动(即。,来源)的分布在给定外部电场和磁场的测量,例如,( 44, 74年, 75年]。然而,应用逆解决方案之前我们必须决定类型的来源及其可能的分布(即。在头部、位置)。

逆解估计不同源建模假设。通过比较的数量未知的源模型的数据量,我们可以区分两种主要类型的问题(即。的解决方案)。首先,过多的问题(例如,偶极解决方案),更多的数据比未知数可以不同最小化算法及其效率逃离局部最小值,拟合优度指标以及使用生理和/或数学约束通常要求在解决方案评估/选择过程。这些模型需要假设关于大脑来源的数量和位置建模为点电流偶极子给一个独特的解决方案确定提供全球最低。这样的方法需要一个模型顺序搜索除了源参数优化( 76年]。数值模拟研究已经证明,一个准确的估计的时间动态偶极模型是非常依赖于解决能力和准确定位所有活跃的大脑区域( 77年]。虽然有一系列的生理和解剖原因,动物研究以及已经融合证据从人类血液动力学和/或代谢功能磁共振成像和宠物研究表明,感觉和认知过程可以被视为一个分布式网络焦点活动;扩展激活神经组织的可能性在某些情况下不能被忽视。深入的电极记录使用后者假设的主要证据,证明在广泛的脑区活动。然而,即使这样的录音,总结贡献的主要来源贡献和体积电流是可以预期的,应该考虑和逆模型,而不是采取这种措施作为扩展大脑激活的强有力的证据。

另一方面,我们有欠定的问题(如分布式逆解)未知数比数据与linear-minimum-norm相关方法,提出( 75年, 78年, 79年)的情况下当焦点源的假设是不合理的。这种方法是具有挑战性的,因为它可能需要进一步的权重和正则化来补偿深度偏差,选择实施数学标准或生理的。

为了帮助读者作出正确的选择,我们描述四个主要来源模型通过限制源类型和/或他们的位置和他们的主要假设。

< inline-formula > < mml:数学xmlns: mml = " http://www.w3.org/1998/Math/MathML " id = " M58 " > < mml:莫弹性=“false”> (< / mml:莫> < mml: mn > 1 < / mml: mn > < mml:莫弹性= "假" >)< / mml:莫> < / mml:数学> < / inline-formula >等效电流偶极子模型

它假设测量是由于单一集中源。它主要是强大和有效空间有限的来源(例如,一些局灶性癫痫)或观察到的来自远方测量表面。这可能是更有用的总结测量领域比源本身,这是一个特殊情况下下列源模型。

< inline-formula > < mml:数学xmlns: mml = " http://www.w3.org/1998/Math/MathML " id = " M59 " > < mml:莫弹性=“false”> (< / mml:莫> < mml: mn > 2 < / mml: mn > < mml:莫弹性= "假" >)< / mml:莫> < / mml:数学> < / inline-formula >偶极模型用于多因素决定的问题

这些模型考虑到测量领域是由于少量的来源未知位置和方向。他们非常适合煤产生的数据过滤和averaged-evoked反应( 80年, 81年]。

< inline-formula > < mml:数学xmlns: mml = " http://www.w3.org/1998/Math/MathML " id = " M60 " > < mml:莫弹性=“false”> (< / mml:莫> < mml: mn > 3 < / mml: mn > < mml:莫弹性= "假" >)< / mml:莫> < / mml:数学> < / inline-formula >皮层模型

极端的假设下深源不会导致头部的外部领域,它假设主要来源位于皮质只在壁炉架方向约束。它可能是非常适合测量与分析一些主要的皮质的激活( 82年]。

先前的模型可以看作是数据驱动的,他们只能用非常具体和严格的实验条件下,将不会接受脑电图作为一般模型和梅格来源。此外,缺乏实验证据支持的偶极子。事实上一个偶极子将意味着一个无限期地增加潜在当我们接近它的位置。希望这从未被报告过,因为在这个位置对应于一个未定义的潜力。然而,一个更完整的源模型必须包含特定的情况下,之前的源模型,同时将这些元素的讨论到目前为止,如下。

< inline-formula > < mml:数学xmlns: mml = " http://www.w3.org/1998/Math/MathML " id = "路" > < mml:莫弹性=“false”> (< / mml:莫> < mml: mn > 4 < / mml: mn > < mml:莫弹性= "假" >)< / mml:莫> < / mml:数学> < / inline-formula >内的电位分布

附近的电磁测量在/头皮是由于大脑内的电位分布。这些(颅内)潜力代表的主要来源在生成,至少,不仅整个灰质和大脑皮层。这个源模型是兼容所有先前的几何约束,包括偶极子作为一种特定的情况下。源模型重要的是,这意味着重大理论和数值简化,解决的问题焦点和扩展来源,由于潜在的总是一个连续函数定义的所有点。

定义适当的预处理和源模型后对我们的数据,我们面临的问题选择逆过程。以下问题可能在这个阶段有关。

5.3。偶极子定位错误

overdetermined-dipolar模型的评价似乎主张采用一个直截了当的解决办法通过比较目标与措施和估计源偶极子定位错误。不幸的是,这些措施不能直接外推到我要说的分布式解决方案。这或许就是为什么分布式解决方案的评估仍然是作为一个开放的问题。显然,这可能会影响逆选择的解决方案。虽然我们不想告诉读者他/她应该做什么/使用,我们想避免讨论一些事情。

有人建议,零偶极子定位可能选择逆解的方法。这可能是出于真正的应用程序的数据是由单一焦来源(如癫痫焦点定位)的长期经验积累以及超定的(偶极)模型。这可能是一种滥用语言,这让人们相信“如果我们正确定位每一个源,然后叠加的原则我们应该正确定位任何来源的组合”。这个声明有两个明显的错误。

在这种情况下,正确定位只意味着最大的模量电流源与目标站点密度一致。这忽略了振幅将几乎总是线性方法和多个来源,错估由于不可避免的非对角元素分辨率矩阵( 3)。

定义清楚,偶极子定位误差(DLE)估计的模量电流源密度,这意味着半夜不是数据的线性函数 d ,因此叠加原理不成立。因此,线性系统理论,描述系统的响应(δ)输入冲动,不能引起。

鉴于之前的理论缺陷,也就不足为奇了半夜无法预测的性能存在多个来源的逆解。事实上它可以证明,单一来源的正确定位是一个微不足道的财产简单而强大的方法(见工作严重·德·佩拉尔塔等人对于这个问题),我们坚持,只适用于如果集中单一源假说成立。

5.4。逆解和空间过滤器

声音逆问题解的物理卷的方法是空间滤波器的估计,这“过滤器”活动,起源于一个特别的位置,在试图抑制从所有其他活动。这些方法现在已经重新出现的名字“beamformers”非常感激,除此之外,因为每个解决方案的解决方案可以独立计算点。继续原来的描述这些方法( 83年, 84年),很明显,减少串扰(即。,distance to the ideal resolution matrix) between sources does not necessarily imply an optimal resolution kernel. Nevertheless, current applications suggest that the solution provided by these methods is not affected by the crosstalk.

有很好的理由选择Backus-Gilbert(即。,beamformer) method such as its adaptive properties to deal with specific noise structures [ 85年]。然而,我们不能强调足够的唯一方法评估线性反演过程提供的估计是看解决内核( 2, 3, 86年]。我们为每个点单独建立一个独立的估计并不意味着这个估计不污染的影响同时活跃的来源。

为了达成这个问题的逆过程的选择积极的一面,我们提到,有一个良好的理论方法选择,更重要的是,建立一个逆解。足够的注意,infinitely-many-linear-inverse方法可以用方程描述 G = C * l * ( l * C * l ) + 。源 j 美国东部时间 用这种方法估计将属于空间横跨的列 C * l 无噪声和噪声情况。另一方面,很明显的 3)的唯一方法改变分辨率矩阵的行(即。,the resolution kernels) is by right transformations of the lead field. These two procedures together yield meaningful source estimates, when C 根据声音先验信息,当选择一个适当的右手领先领域的转换是由( 87年]。

5.5。健壮的方法分析EEG /梅格来源

估计脑电图/梅格的问题来源可以解释如下。提供精确的测量数据(噪声)但当地信息。为了更了解整个系统(即。,the brain), we need to ascend to qualitatively higher levels corresponding to the surface maps and the 3D distribution of sources. By doing this we obtain a more complete global descriptor but probably also with a higher incertitude (if compared with the sensor data).

它也是如此的fMRI信号( 88年),一般来说我们不能依靠提供的振幅逆解比较神经活动在两个不同的位置。出于同样的原因,鬼和失去来源出现在每一个重建与真正的混合来源。因此,区分真正的来源从工件几乎是不可能的,除非我们知道真正的分布。因此,我们可以说,获得的源分布从一个地图可能是最不精确的图片,我们的大脑活动。

我们能做什么来增加这些功能图像的可靠性?至于部分答案,我们建议以下几点。

选择您的逆解记住前面讨论的关于空间过滤器和零偶极子定位误差和小心考虑源分布估计从单个地图(生产,例如,通过分割算法)。

使用减少underdetermination源模型的反问题。优先身体声音转换减少标量场的问题,估计提高分辨率的内核。

计算大小或数据的基于时间的大脑而不是当地瞬时振幅的来源和使用措施,独立于颅内的比例因子信号相关系数。

评估对比实验条件或prestimulus与poststimulus条件以减少系统的幽灵和失去了源效应( 31日, 89年]。

计算之间的相关性大小源自大脑活动的时间进程和行为测量反应时间( 90年]。

6。结论

有许多关键领域,严重影响源定位的准确度和精密度。在本文中,我们讨论了四个关键领域的脑电图/梅格源成像,即预处理,容积导体,和逆问题。尽管这些广泛的组件,我们着重引导社会分配注意这些关键的未决问题。首先,电导率同样影响解决方案从而保证和逆解实际电导率测量现场组织需要填补空虚的这些关键参数。这些未来的研究应该准确记录他们的测量setups-most特别是湿度和温度。其次,未来的研究应该把如何建模病理改变一个正常的,健康的头模型。最后,关键是选择源模型和逆过程基于良好的理论和实验基础。

最终,我们应该做出明智的决策优化模型元素,获得最精密度和准确度源成像和抑制那些贡献最小收益源定位。这种优化后,如何这些未来的模型表示他们的生理,也就是说,人类的头?作为我们继续前进作为一个社区,我们应该记得要突出未来研究反映新的导率的缺点,病态,源模型,等等,这些模型,以防止任何进一步的误解,而集体建设的贡献模型。

确认

这项工作已经部分支持的欧洲项目BACS fp6 -是- 027140。本文只反映了作者的观点,和资助机构不承担任何使用这可能是由本文所包含的信息。这项工作也支持的成本行动BM0601“NeuroMath”。无论是办公室成本的代理也没有人负责使用,这可能是出版物中包含的信息。成本办公室不负责外部网站本刊物。

ven D。 泰伯的百科全书的医学词典 2005年 20 美国费城,宾夕法尼亚州 f·a·戴维斯公司 戴尔 a . M。 塞里诺 m . I。 改善大脑皮层活动的本地化结合脑电图和梅格MRI大脑皮层表面重建:线性方法 认知神经科学杂志》 1993年 5 2 162年 176年 严重的德佩拉尔塔梅内德斯 R。 冈萨雷斯Andino s . L。 Lutkenhoner B。 数据的价值比较分布式线性逆解 大脑地形 1996年 9 2 117年 124年 10.1007 / BF01200711 忠实的 a·j·R。 迪亚斯 答:我。 维埃拉 j . P。 分析gelastic发作癫痫的脑电图动态活动在独立的组件使用分解 临床神经生理学 2006年 117年 7 1595年 1601年 10.1016 / j.clinph.2006.03.020 安德利果汁 P。 Kloesch G。 格鲁伯 G。 低分辨率的大脑额叶和顶叶电磁层析成像显示同时活跃在人类睡眠纺锤体来源皮层 神经科学 2001年 103年 3 581年 592年 10.1016 / s0306 - 4522 (01) 00028 - 8 Durka P。 Matysiak 一个。 蒙特斯 e . M。 索萨 p V。 Blinowska k·J。 多道匹配追踪和EEG逆解决方案 神经科学杂志》上的方法 2005年 148年 1 49 59 10.1016 / j.jneumeth.2005.04.001 戴维拉 c, E。 Srebro R。 子空间平均稳态视觉诱发电位 IEEE生物医学工程 2000年 47 6 720年 728年 O。 中村 M。 Nagamine T。 Shibasaki H。 参数化建模的躯体感觉诱发电位使用离散余弦变换 IEEE生物医学工程 2001年 48 11 1347年 1351年 10.1109/10.959331 皮克顿 t·W。 手册的脑电图与临床神经生理学:人类与事件相关电位 1988年 阿姆斯特丹,荷兰 爱思唯尔 詹森 b . H。 阿加瓦尔 G。 对冲基金 一个。 加利 N . N。 相位同步的脑电图和听觉EP的一代 临床神经生理学 2003年 114年 1 79年 85年 10.1016 / s1388 - 2457 (02) 00327 - 9 年代。 德里斯科尔 d . A。 脑电图电极sensitivity-an互惠的应用 IEEE生物医学工程 1969年 16 1 15 22 罗斯 b . J。 Balish M。 Gorbach 一个。 佐藤 年代。 情况如何three-sphere模型预测位置的偶极子实际形状的头? 脑电图与临床神经生理学 1993年 87年 4 175年 184年 10.1016 / 0013 - 4694 (93)90017 - p Crouzeix 一个。 Yvert B。 贝特朗 O。 Pernier J。 偶极子的评价重建精度梅格的球形和真实头模型 临床神经生理学 1999年 110年 12 2176年 2188年 10.1016 / s1388 - 2457 (99) 00174 - 1 Cuffin b . N。 脑电图定位精度使用实际形状的头部模型的改进 IEEE生物医学工程 1996年 43 3 299年 303年 10.1109/10.486287 Hamalainen m . S。 Sarvas J。 现实的电导率的几何模型解释的人头neuromagnetic数据 IEEE生物医学工程 1989年 36 2 165年 171年 10.1109/10.16463 Wendel的 K。 紫檀 n G。 Hannula M。 Kauppinen P。 Malmivuo J。 CSF对脑电图的影响灵敏度分布的多层模型 IEEE生物医学工程 2008年 55 4 1454年 1456年 10.1109 / TBME.2007.912427 Wendel的 K。 Osadebey M。 Malmivuo J。 将人体测量和craniometric数据纳入realistically-shaped容积导体的头模型 学报》第11届世界大会在医学物理学和生物医学工程 2009年9月 德国慕尼黑 Haueisen J。 拉蒙 C。 Eiselt M。 布劳尔 H。 诺瓦克 H。 组织电阻率对neuromagnetic领域和电动电位的影响研究的有限元模型 IEEE生物医学工程 1997年 44 8 727年 735年 10.1109/10.605429 施万 h·P。 Nastuk w . L。 生物阻抗的测定 物理技术在生物研究 1963年 纽约,纽约,美国 学术出版社 盖伯瑞尔 C。 盖伯瑞尔 年代。 Corthout E。 生物tissues-I的介电性能:文献调查 物理学在医学和生物学 1996年 41 11 2231年 2249年 10.1088 / 0031 - 9155/41/11/001 戈德斯 l。 贝克 l E。 生物作为数据的一种概略的比电阻的生物医学工程师和生理学家 医学和生物工程 1967年 5 3 271年 293年 10.1007 / BF02474537 Radich b . M。 巴克利 k . M。 脑电偶极子定位范围和映射算法模型具有参数不确定性 IEEE生物医学工程 1995年 42 3 233年 241年 10.1109/10.364509 van den Broek s P。 Reinders F。 Donderwinkel M。 彼得斯 m·J。 体积传导影响脑电图和梅格 脑电图与临床神经生理学 1998年 106年 6 522年 534年 10.1016 / s0013 - 4694 (97) 00147 - 8 Huiskamp G。 Vroeijenstijn M。 冯·迪 R。 Wieneke G。 范Huffelen a . C。 需要正确的现实几何的EEG逆问题 IEEE生物医学工程 1999年 46 11 1281年 1287年 10.1109/10.797987 Vanrumste B。 范·霍利 G。 Van de机器人瓦力 R。 D 'Have M。 Lemahieu 我。 P。 偶极子定位误差源analyssis脑电图由于容积导体模型错误 医学和生物工程和计算机 2000年 38 5 528年 534年 Gencer n G。 Acar c, E。 脑电图的敏感性和梅格组织电导率测量 物理学在医学和生物学 2004年 49 5 701年 717年 10.1088 / 0031 - 9155/49/5/004 Ollikainen j . O。 Vauhkonen M。 Karjalainen p。 Kaipio j . P。 当地的头骨尺度对脑电图的影响源估计 医学工程学和物理学 1999年 21 3 143年 154年 10.1016 / s1350 - 4533 (99) 00038 - 7 Hoekema R。 Wieneke g . H。 Leijten f·S·S。 电导率的测量头骨,癫痫手术期间暂时移除 大脑地形 2003年 16 1 29日 38 10.1023 /:1025606415858 Wendel的 K。 Malmivuo J。 联系生活和事后剖析颅骨电导率测量 学报的IEEE年度国际会议在医学和生物学社会工程 2006年 1 4285年 4288年 米歇尔 c . M。 Lantz G。 卢卡雷利 l 严重的德佩拉尔塔梅内德斯 R。 兰迪斯 T。 Seeck M。 128 -通道脑电图在癫痫源成像:临床产量和定位精度 临床神经生理学杂志 2004年 21 2 71年 83年 10.1097 / 00004691-200403000-00001 冈萨雷斯Andino s . L。 严重的德佩拉尔塔梅内德斯 R。 Khateb 一个。 兰迪斯 T。 a·J。 电生理学相关的情感盲视 科学杂志 2009年 44 2 581年 589年 10.1016 / j.neuroimage.2008.09.002 摩尔 k . L。 Dalley 答:F。 面向临床解剖学 2005年 5日 美国费城,宾夕法尼亚州 Lippincott Williams &威尔金斯 Oostendorp t F。 Delbeke J。 斯蒂门 d F。 人类头骨的导电率:体内和体外测量的结果 IEEE生物医学工程 2000年 47 11 1487年 1492年 Kosterich j . D。 福斯特 k·R。 波拉克 s R。 电介质的介电常数和电导率流体饱和的骨头 IEEE生物医学工程 1983年 30. 2 81年 86年 Kosterich j . D。 福斯特 k·R。 波拉克 s R。 介电性能的流体—骨头。浸液的电导率的变化的影响 IEEE生物医学工程 1984年 31日 4 369年 374年 Goncalves 我美国。 de Munck j . C。 Verbunt j . p。 Bijma F。 Heethaar r·M。 达席尔瓦 f . L。 体内使用EIT-based大脑和头骨的电阻率测量方法和现实的模型 IEEE生物医学工程 2003年 50 6 754年 767年 10.1109 / TBME.2003.812164 格里菲思 H。 斯图尔特 w·R。 咳嗽 W。 磁感应断层。生物组织的测量系统 纽约科学院上 1999年 873年 335年 345年 10.1111 / j.1749-6632.1999.tb09481.x Gencer n G。 Tek m . N。 通过非接触式测量电导率成像 IEEE医学成像 1999年 18 7 617年 627年 N。 基于电流密度电阻抗断层成像 1992年 加拿大多伦多 多伦多大学 考虑 y Z。 Onart 年代。 狮心王 w·r·B。 独特性和重建在磁性resonance-electrical阻抗断层扫描(MR-EIT) 生理测量 2003年 24 2 591年 604年 10.1088 / 0967 - 3334/24/2/368 Gevins 一个。 H。 史密斯 m E。 J。 R。 映射的认知脑功能与现代高分辨率脑电图 神经科学的趋势 1995年 18 10 429年 436年 10.1016 / 0166 - 2236 (95)94489 - r Gevins 一个。 J。 H。 McEvoy l·K。 史密斯 m E。 由模糊变清晰 临床神经生理学杂志 1999年 16 3 204年 213年 10.1097 / 00004691-199905000-00002 Babiloni F。 Cincotti F。 性格外向 F。 罗西尼 p . M。 Babiloni C。 脑电图数据的空间增强表面拉普拉斯算子的估计:磁共振成像头模型的使用 临床神经生理学 2001年 112年 5 724年 727年 10.1016 / s1388 - 2457 (01) 00494 - 1 米歇尔 c . M。 穆雷 M . M。 Lantz G。 冈萨雷斯 年代。 卢卡雷利 l 严重的德佩拉尔塔梅内德斯 R。 脑电图源成像 临床神经生理学 2004年 115年 10 2195年 2222年 10.1016 / j.clinph.2004.06.001 Gevins 一个。 Brickett P。 B。 J。 发音 B。 除了地形映射:对functional-anatomical成像与124 -通道脑电图和三维核磁共振成像 大脑地形 1990年 3 1 53 64年 Srinivasan R。 Nunez p . L。 塔克 d . M。 Silberstein r B。 Cadusch p . J。 空间采样和滤波花键拉普拉斯算子的脑电图评估皮质潜力 大脑地形 1996年 8 4 355年 366年 Babiloni F。 Babiloni C。 性格外向 F。 高分辨率的脑电图:一种新的依赖于模式的空间去模糊方法使用realistically-shaped MR-constructed主题的头模型 脑电图与临床神经生理学 1997年 102年 2 69年 80年 10.1016 / s0921 - 884 x (96) 96508 - x Ryynanen o . r . M。 Hyttinen j·a·K。 Malmivuo J。 测量噪声的影响和电极在皮层电位分布的空间分辨率和密度不同的电阻率值的头骨 IEEE生物医学工程 2006年 53 9 1851年 1858年 10.1109 / TBME.2006.873744 Srinivasan R。 塔克 d . M。 Murias M。 估计的空间尼奎斯特人的脑电图 行为研究方法、仪器和电脑 1998年 30. 1 8 19 Y。 B。 计算机模拟的研究从头皮电位皮质成像 IEEE生物医学工程 1998年 45 6 724年 735年 10.1109/10.678607 伯曼 一个。 Plemmons R。 非负矩阵在数学科学 1994年 美国费城,宾夕法尼亚州 暹罗 新闻 w·H。 Teukolsky 美国一个。 Vetterling w·T。 弗兰纳里 b P。 C中的数值食谱 1995年 英国剑桥 剑桥大学出版社 Ruge j·W。 Stuben K。 麦考密克 美国F。 代数多重网格(AMG) 多重网格方法 1987年 3 美国费城,宾夕法尼亚州 暹罗 73年 130年 前沿应用数学 布里格斯 w . L。 亨森 诉E。 麦考密克 美国F。 多栅的教程 2000年 美国费城,宾夕法尼亚州 暹罗 巴纳德 a·c·L。 i M。 林恩 m . S。 电磁理论的应用electrocardiography-I:积分方程的推导 生物物理学杂志 1967年 7 443年 462年 Akalin-Acar Z。 Gencer n G。 一个先进的边界元法(BEM)实现远期电磁源成像问题 物理学在医学和生物学 2004年 49 21 5011年 5028年 10.1088 / 0031 - 9155/49/21/012 Ataseven Y。 Akalin-Acar Z。 Acar c, E。 Gencer n G。 EMSI加速BEM方法的并行实现人类的大脑 医学和生物工程和计算机 2008年 46 7 671年 679年 10.1007 / s11517 - 008 - 0316 - 0 冈萨雷斯Andino s . L。 严重的德佩拉尔塔梅内德斯 R。 比斯开湾 R。 吉梅内斯 j . C。 帕斯卡 r D。 Lemagne J。 射影方法磁直接问题 生物磁学的发展 1989年 纽约,纽约,美国 达塔 b . N。 数值线性代数和应用程序 1995年 美国加州太平洋格罗夫 布鲁克斯/科尔 布里格斯 w . L。 多栅的教程 1987年 美国费城,宾夕法尼亚州 暹罗 Hoekema R。 文纳 K。 Struijk J·J。 Holsheimer J。 多栅的解决方案建模的势场电神经刺激 计算机和生物医学研究 1998年 31日 5 348年 362年 10.1006 / cbmr.1998.1486 Laarne P。 Hyttinen J。 Dodel 年代。 Malmivuo J。 Eskola H。 两个偶极逆算法精度互惠申请提出的计算 计算机和生物医学研究 2000年 33 3 172年 185年 10.1006 / cbmr.1999.1538 萨德 Y。 稀疏线性系统的迭代方法 2003年 2日 美国费城,宾夕法尼亚州 暹罗 尼尔森 l。 尤夫 M。 科尔 z . J。 效率计算方法准确地解决脑电图提出问题一头细离散模型 临床神经生理学 2005年 116年 10 2302年 2314年 10.1016 / j.clinph.2005.07.010 汤普森 j·F。 索尼 b K。 Weatherrill n P。 手册的网格生成 1998年 美国佛罗里达州波卡拉顿 CRC的新闻 Ottosen N。 彼得森 H。 介绍了有限元法 1992年 恩格尔伍德悬崖,新泽西,美国 新世纪 Saleheen h . I。 T。 新的有限差分公式一般不均匀各向异性生物问题 IEEE生物医学工程 1997年 44 9 800年 809年 10.1109/10.623049 Plonsey R。 的本质来源的生物和生物磁效应字段 生物物理期刊 1982年 39 3 309年 312年 Hamalainen M。 Ilmoniemi r . J。 Sarvas J。 Hario H。 相互依赖的信息传达magnetoencephalogram和脑电图 逆的理论和应用问题 1988年 纽约,纽约,美国 约翰威利& Sons Malmivuo J。 Suihko V。 Eskola H。 脑电图和梅格的灵敏度分布测量 IEEE生物医学工程 1997年 44 3 196年 208年 10.1109/10.554766 Malmivuo J。 Suihko 诉E。 头骨电阻率对脑电图的空间分辨率的影响,梅格 IEEE生物医学工程 2004年 51 7 1276年 1280年 10.1109 / TBME.2004.827255 答:K。 戴尔 a . M。 Belliveau j·W。 蒙特卡罗模拟研究脑电图和MEG定位精度 人类大脑图谱 2002年 16 1 47 62年 10.1002 / hbm.10024 Waldert 年代。 Preissl H。 Demandt E。 手的运动方向解码来自梅格和脑电图 神经科学杂志》上 2008年 28 4 1000年 1008年 10.1523 / jneurosci.5171 - 07.2008 Hamalainen M。 哈里 R。 Ilmoniemi r . J。 Knuutila J。 Lounasmaa o . V。 Magnetoencephalography-theory、仪器和应用无创性研究人类大脑的工作 现代物理学的评论 1993年 65年 2 413年 497年 10.1103 / RevModPhys.65.413 严重的德佩拉尔塔梅内德斯 R。 Gonzalez-Andino s . L。 的关键分析线性逆neuroelectromagnetic逆问题的解决方案 IEEE生物医学工程 1998年 45 4 440年 448年 Supek 年代。 皇家艺术 c·J。 多个偶极子neuromagnetic源定位:模型的模拟研究秩序和源分辨率的极限 IEEE生物医学工程 1993年 40 6 529年 540年 10.1109/10.237672 Supek 年代。 皇家艺术 c·J。 时空建模neuromagnetic我:多源位置和时间进程估计精度 人类大脑图谱 1997年 5 3 139年 153年 10.1002 / (SICI) 1097 - 0193 (1997) 5:3 < 139:: AID-HBM1 > 3.0.CO; 2 - 4 Hamalainen m . S。 Ilmoniemi r . J。 解释大脑的磁场:估计最低标准 医学和生物工程和计算机 1994年 32 1 35 42 霍克 O。 保持简单:古典最小准则估计还是可以使用的,脑电图和梅格数据分析 科学杂志 2004年 21 4 1612年 1621年 10.1016 / j.neuroimage.2003.12.018 Scherg M。 冯Cramon D。 唤起人类听觉皮层偶极子源的潜力 脑电图与临床神经生理学 1986年 65年 5 344年 360年 Scherg M。 皮克顿 t·W。 分离和鉴定与事件相关电位的脑电源分析组件 脑电图与临床神经生理学。补充 1991年 42 24 37 Babiloni F。 Babiloni C。 Locche l Cincotti F。 罗西尼 p . M。 性格外向 F。 高分辨率脑电图:源估计Laplacian-transformed体感觉诱发电位使用现实主题头部模型由磁共振图像 医学和生物工程和计算机 2000年 38 5 512年 519年 女人气的男人 J。 高分辨率的频率-波数 IEEE学报》 1969年 57 8 1408年 1418年 巴克斯 g . E。 吉尔伯特 j·F。 地球总值数据的分辨能力 地球物理《皇家天文学会杂志》上 1968年 16 169年 205年 Sekihara K。 美国年代。 Poeppel D。 马兰士 一个。 性能的梅格adaptive-beamformer技术的存在相关的神经活动:对信号强度和时间进程的影响估计 IEEE生物医学工程 2002年 49 12我 1534年 1546年 10.1109 / TBME.2002.805485 Lutkenhoner B。 严重的德佩拉尔塔梅内德斯 R。 resolution-field概念 脑电图与临床神经生理学 1997年 102年 4 326年 334年 10.1016 / s0013 - 4694 (96) 96590 - 6 严重的德佩拉尔塔梅内德斯 R。 穆雷 M . M。 米歇尔 c . M。 Martuzzi R。 冈萨雷斯Andino s . L。 电神经成像基于生物物理限制 科学杂志 2004年 21 2 527年 539年 10.1016 / j.neuroimage.2003.09.051 Logothetis n K。 万德尔 b。 解释的信号 年度回顾的生理 2004年 66年 735年 769年 10.1146 / annurev.physiol.66.082602.092845 冈萨雷斯Andino s . L。 严重的德佩拉尔塔梅内德斯 R。 Khateb 一个。 a·J。 Thut G。 兰迪斯 T。 一窥你的视力 人类大脑图谱 2007年 28 7 614年 624年 10.1002 / hbm.20302 冈萨雷斯Andino s . L。 米歇尔 c . M。 Thut G。 兰迪斯 T。 严重的德佩拉尔塔梅内德斯 R。 响应速度的预测先行高频振荡(伽玛乐队)在人类的大脑 人类大脑图谱 2005年 24 1 50 58 10.1002 / hbm.20056