约束与权衡14 分钟

头皮上的电极,能够还原多少脑内活动?

头皮脑电记录神经电流经过组织和参考方式后形成的电位差,源定位再结合几何与先验推断脑内分布。本文用明确的均匀导体模型计算不可见源分量、增加测点和模式约束的影响,比较干电极接触、2026耳内脑电的任务表现与几何基函数重建。检测到一种信号、准确定位一个峰和恢复完整动态有不同证据要求;合适的先验可以提高特定重建,但不能把缺失信息自动变成已测事实。

从体表读数到脑内图像,中间多了一次推断

头皮EEG首先记录电位差。脑内许多电流源的贡献经过组织叠加,再由电极、放大器与后续处理形成曲线。源定位的任务,则是根据这些曲线及头部模型,推断哪些脑内电流分布可能产生它们。

这两个步骤值得分别理解。一台设备可以稳定检出某种频段变化,但未必有足够空间信息定位其全部来源;一种定位方法也可以在特定验证中表现很好,同时仍依赖源的方向、分布和时间结构等假设。

本篇沿着“神经电流—组织—电极—参考—重建”这条链,解释信息怎样进入记录,又有哪些条件帮助我们把它读出来。

组织把许多电流源混合到每个测点

上一文介绍了神经元的跨膜电流和返回电流。对较远处的测点,一小片有组织的电流常可近似成电流偶极子:它同时具有强度和方向。多个源的电位贡献在通常采用的线性、准静态模型中相加。1 2

因此,一个测点的信号并不专属于它正下方的一块皮层。来源的方向、距离、皮层折叠、组织导电率及其他源,都影响它最后接收到的电位。头部模型常区分脑组织、颅骨、脑脊液与头皮,采用边界元或有限元等方法计算传播后的读数。2

这些方法解决的是正向问题:先给定一组脑内源,计算电极应看到什么。它们需要电极位置与头部几何配准,也需要导电率和边界条件。MRI提供很多结构信息,但不直接告诉我们每一时刻的神经电流。

将许多时刻的关系写成一条式子:

$$ \mathbf y(t)=G\mathbf j(t)+\boldsymbol\varepsilon(t). $$

$\mathbf j(t)$表示候选源的强度,$G$表示从这些源到指定参考下电极读数的映射,$\mathbf y(t)$是实际记录,$\boldsymbol\varepsilon(t)$汇总噪声和模型没有描述的部分。更换源方向、头部模型或电极位置,都会改变$G$。

如果预先固定源的位置与方向,只求每个源的强度,这通常是线性问题;若还要搜索源的位置和数目,问题会更加复杂。测量数大于参数数,也不能单独保证一个非线性定位问题处处具有唯一解。2

参考电极也参与决定你看到了什么

EEG测量需要比较电位。最简单的记录是某一测点减去一个参考;常见的平均参考则是每个测点减去所选测点的平均值:

$$ \mathbf y_{\mathrm{avg}}=R\boldsymbol\phi, \qquad R=I-\frac{\mathbf1\mathbf1^T}{M}. $$

如果所有测点同时增加同样的电位,这个共同部分会被平均参考相消。$M$个平均参考通道之和为零,因此最多有$M-1$项独立电位差信息。这个关系说的是读数的线性结构,不能直接翻译成“大脑只有这么多个活动变量”。

参考还会影响特定源的可见性。假设两耳附近收到相似的某个源成分,相减后它可能被削弱;如果两个位置收到的成分不同,差异便可能保留。相消可以同时作用于干扰和真正的神经信号,是否有利取决于想观察什么。

所以,比较两台设备时,除了电极材料,还需要比较实际测点与参考。将头顶对平均参考的记录改成右侧颞部对左侧颞部,已经同时改变观察位置和比较方式,不能把差异全部归给电极质量。

干电极的关键,是稳定接触与完整测量链

干电极减少了对导电膏或凝胶的依赖,但仍然需要穿过头发建立稳定接触。接触压力、接触面积、头部运动、皮肤状态和电极材料共同影响界面阻抗与噪声。放大器的输入阻抗和自身噪声,又决定这类接口能否提供可用记录。3

单看接触阻抗的大小并不够。用一个只有电阻的简化接口作例子:放大器输入电阻为10亿欧姆,两路接触分别为1万与100万欧姆,单路幅度衰减都很小。然而,如果两路共同接收到1毫伏干扰,由于衰减比例不完全相同,相减后仍留下约0.989微伏的差异。

这只是本库按分压关系计算的示例,不是某款设备的实测规格。它说明为何接触不平衡和不稳定值得关注,也说明更高ADC位数不能自动补回已经被接触与空间相消改变的信息。

一个有实证依据的改进方向是机械结构。2023年的Flower电极研究比较了两种干电极帽:弯曲的接触针在受压时增加接触面积,与传统多针结构相比,参与者报告的佩戴舒适度更好。研究包含20名健康、直发参与者,两种帽在不同日期按随机顺序使用,并比较坐姿和仰卧条件。3

与此同时,平均通道可用比例约在七成上下,部分电极连接在使用中损坏;两种帽的许多群体信号特征接近,也不等于每个通道始终等效。它说明机械舒适、接触稳定和信号可用性需要一起改进,而“干电极”本身不是统一性能等级。该研究直接比较的是两种干电极设计,不能将它改写成全部干、湿设备的对决。

2026年的耳内研究:一种信号成功,不代表每种信号都一样

一项2026年研究比较了19名参与者的同步头皮与耳内脑电。头皮系统有32通道;耳内装置用右耳信号减去左耳参考,另有接地电极,实际只有一个EEG通道。采集板支持8通道,不代表本实验真的记录了8个脑电空间通道。4

较清楚的正向结果来自静息睁眼与闭眼比较。耳内静息记录因缺失片段或噪声排除部分数据后,实际纳入14人,头皮记录纳入19人。两种记录中,闭眼时α功率都增加。

记录方式 静息纳入人数 8—10 Hz闭眼减睁眼(dB) 10—13 Hz闭眼减睁眼(dB)
头皮平均参考,多通道功率汇总 19 2.89 2.79
右耳减左耳 14 0.89 0.96

数值转录自原表1,是各记录方式内部的条件差值;它们不是同一组完整配对下的设备误差,也不是“准确率”。

这支持该配置检出较明显的静息频谱变化,但还不足以证明两台设备等效,或耳内记录能定位α来自哪里。论文没有直接进行设备间差异检验;某种记录显著、另一种不显著,本身也不等于两种记录的差异已经显著。4

听觉事件相关电位ERP部分则有一个实质缺口:正文和表2给出的均值、检验结果不一致,连适配头皮配置的N1是否达到显著标准也不同。本文因此不据这些数字给ERP性能作精确排序。噪声中言语任务的α变化,同样只能按该研究各配置内部的基线比较来理解。

还有一些因素没有被拆开:两系统采样率不同,多通道头皮数据可以做独立成分去伪迹,单通道耳内数据没有同等分解能力;耳内接触阻抗也未量化。因此,研究展示的是这套实际组合的表现,不能单独裁定所有耳内电极、所有干接触或某一种参考方式的能力。

为什么同一组读数可以对应不同脑内源

可以把这个问题算出来。本库在无限、均匀、各向同性导体中放置96个径向点偶极源,并在更外侧安排32个测点。源的位置分布在半径60毫米的球面,测点在半径90毫米的球面;这些球面只规定位置,没有充当导电介质的边界。模型没有颅骨、脑脊液或真实皮层折叠。

单个偶极在测点的电位由其方向和相对位置计算,再对所有源求和。这个简单设置足以展示逆问题的核心,同时不冒充真实个体头部。

均匀导体中的虚拟源、测点与正向映射
均匀导体中的虚拟源、测点与正向映射

模型导电率为自选的0.33 S/m。源系数以nAm表示,读数以微伏表示;这些参数是示例输入,不是人体测量或空间分辨率标准。

计算得到,32测点的平均参考映射秩为31。对于96个自由源系数,便存在65维的变化方向可以不改变这组读数。把其中一项记为$\mathbf z$,有:

$$ G\mathbf z=0, \qquad G(\mathbf j+\mathbf z)=G\mathbf j. $$

也就是说,源分布A和源分布B可以明显不同,电极读数却完全一样。这里的不可见性针对这组测点、参考和源模型,不能扩大成这个分量在任何仪器下都不可观察。

不同源分布及其相同的平均参考读数
不同源分布及其相同的平均参考读数

上两图的源序号只是球面采样点的编号,不代表时间,也不是相邻空间位置的排序。两组源在32个读数上的相对差异约为$4.7\times10^{-16}$,处于数值舍入精度。

把测点扩展为包含原布局的64个,映射秩升为63,前面选出的不可见分量也能在新增测量中留下差异。但96个源系数仍有33维不可见方向。更多空间测量提供了新约束,没有自动变成全部脑内活动的完整录像。

若允许每一时刻的源自由变化,多记录一些时间点也不会自动消去这种空间歧义。已知动力学、跨时刻连续性或其他模态可以提供额外约束,但那需要另行说明加入了什么信息。

先验怎样帮助重建,为什么它不是无价值的“猜测”

如果完全不限制96个源系数,会有很多解。最小范数方法从能解释读数的解中,选择系数平方和较小的一个。它提供了确定的选择规则,但这个平方和并不等于真实大脑的能耗。2

另一种办法是先限制可能的空间形状:假设源由少数平滑模式组合而成。本库用8个低阶球谐模式构造一个源,再加入1000组随机噪声作比较。

当真实源确实属于这8个模式时,受限重建的源误差中位数约为2.68%,96源最小范数约为28.80%。这说明合适的结构知识能够显著帮助推断。

但如果实际源还包含前面那项32测点看不见的分量,传感器收到的数据不变,两种方法的源误差中位数就分别变为约69.76%和74.44%。没有观测依据时,算法无法凭空认出那项分量。

匹配与包含不可见分量两种条件下的源重建误差
匹配与包含不可见分量两种条件下的源重建误差

1000次自选噪声试算的误差分布,不是1000名参与者或设备验证。第一种真值就是按所选模式生成,因而它展示先验匹配时的作用,不是公平评选所有算法。源误差定义为重建与真值差向量的长度,除以真值向量长度。

还有一个容易忽视的现象:96源最小范数在这个例子里可以把含噪测点数据拟合到数值精度,源分布仍然有明显误差。8模式重建留下了一些测点残差,却在匹配条件下更接近真值。因此,“把头皮曲线拟合得更像”与“把脑内活动找得更准”不是同一项评价。

先验的价值来自它与研究对象的对应。几何、稀疏性、平滑性或已知通路都可能提供有用约束;它们需要接受外部数据和适用条件的检验,而不应仅凭有假设就被判无效。2

2026年的几何基函数方法,增加了哪一层信息

Wang等2026年的GBF方法,把每个人MRI提取出的皮层表面分解成几何模式。较低阶模式变化较平缓,较高阶模式包含更细的空间起伏;源活动用这些模式的组合来表示,再通过先验权重抑制不稳定的重建。5

用$A$表示模式矩阵,$\boldsymbol\theta$表示模式系数,便将原问题写成:

$$ \mathbf y=GA\boldsymbol\theta+\boldsymbol\varepsilon. $$

在经过噪声白化的坐标中,可以把求解理解为同时压低两部分:一部分是测量残差,另一部分是先验规定的系数惩罚。

$$ \widehat{\boldsymbol\theta} =\arg\min_{\boldsymbol\theta} \left\{\|\mathbf y-GA\boldsymbol\theta\|^2 +\boldsymbol\theta^T\Lambda\boldsymbol\theta\right\}. $$

这条路线把个体几何引入源的空间形状,区别于只使用几何计算导电路径。本文前面的球谐例子展示相关逻辑,但没有复现GBF的个体皮层网格、具体正则化和完整算法。

作者的模拟基准包含200幅由fMRI元分析空间图生成的源图,经过正向模型再添加噪声。方法在所测条件下表现出重建优势。这里的“真实噪声”主要指按真实EEG估计的空间协方差生成的噪声,不是完整重放全部自然脑电的时间结构。已知模拟源是重要验证,但fMRI统计图构造的源也不是实际逐时刻电流真值。5

其他验证进一步增加了依据,也有不同层级。任务反应与已知解剖分布相符,是有用的外部一致性;不同人群的MEG重建网络与颅内记录的群体网络对应,不等于同一人的每个源时间序列都被直接核验。研究中的皮层下扩展仍属初步,主要结论集中于皮层表面。当前框架逐时刻独立求解,动态先验是后续方向。

刺激定位和术后切除范围也提供临床相关参照,但刺激点、随后活跃的神经网络、切除腔和真正必需的致痫组织并非同一对象。该研究的术后比较选用了24名结局较好的患者,不能直接读作面向所有患者的术前预测能力。其刺激定位的时间窗口在作者稿与公开示例代码之间还有差异,本文未复现全批流程,因而不把展示的传播路径当作已逐段验证的实际神经传导。5

这些边界保留了方法的研究价值:个体几何是可检验的约束,已有若干有利验证;其作用是改善特定条件下的重建,不是让所有可能源在数学上都变得唯一可辨。

验证“真值”时,要问它究竟真在哪里

一种很直接的验证思路,是在已知位置放入电流源,再看头皮反推能否找回来。Mikulan等2020年的Localize-MI数据采用了已知位置的颅内刺激电流,对其在头皮留下的电信号做定位验证。作者明确说明,使用的是刺激伪迹,优点是发生位置清楚;把结果外推到生理信号时则需要额外判断。6

另一类数据记录刺激之后的神经反应。2022年的同步SEEG与高密度头皮EEG研究提供了这类CCEP记录,包含多个脑区对刺激的响应。它们可以帮助研究网络如何被招募,但其全部活跃源不会简单等于刺激电极的位置。7

因此,校准电流、模拟源、颅内神经活动、任务反应和临床结局,各自回答不同的验证问题。它们互相补充;一种参照上的成功,需要经过合理桥梁才能扩展到另一种用途。

怎样理解一台设备或一张重建脑图

可以先确认它实际记录哪些电位差、哪些记录可用,再看目标是检出频谱变化、测量事件波形、定位一个局部峰,还是恢复分布式动态。随后检查空间模型、参考、伪迹处理和先验,最后比较相应任务的外部验证。

目前较稳的认识是,EEG具有明确的电磁测量基础,能够保留有用的时间信息;增加可靠的空间采样、改善接触并加入适当结构约束,可以提高特定任务的表现。但电极数量、ADC位数、重建图的清晰程度和拟合优度,都不能独自充当“还原了多少真实脑活动”的答案。

本篇的整合判断是,源重建值得使用和继续研究,关键是把读数提供的信息与额外约束共同摆在明处。剩余不确定性应落到具体源、场景和参照上,而不是给所有EEG设备统一贴上“读脑成功”或“完全无用”的标签。

计算结果、源分布和测点读数可回查。模型未包含真实分层头部,原临床与硬件数据未重新采集,本篇不提供个人诊断或设备治疗建议。

来源与依据

  1. S1 脑电场与神经电流综述

    The origin of extracellular fields and currents — EEG, ECoG, LFP and spikes

    György Buzsáki; Costas A Anastassiou; Christof Koch · 2012

    · 相关全文已核

    电流/几何/体积导电相关段

    复用T65已取原件,校验同文件。实际读突触/返回电流、几何和叠加,未把单个测点等同局部细胞,不继承所有全物种幅度或1/f解释。

  2. S2 源成像方法与模型综述

    Electrophysiological Source Imaging: a Noninvasive Window to Brain Dynamics

    Bin He; Abbas Sohrabpour; Emery Brown; Zhongming Liu · 2018

    · 相关全文已核

    §2正向/逆向概念,§3 ECD/分布源/最小范数/波束方法

    2018发表,PMC收录日期不当发表年。线性固定源与非线性位置搜索区分;未继承“参数少于传感器便唯一”的过强说法。最小范数是系数平方惩罚,不是脑耗能;2018通道建议不当2026统一监管或设备标准。

  3. S3 原始干电极机械结构与记录比较

    Flower electrodes for comfortable dry electroencephalography

    Indhika Fauzhan Warsito; Milana Komosar; Maria Anne Bernhard; Patrique Fiedler; Jens Haueisen · 2023

    · 相关全文已核

    Methods cap/20 volunteers/randomized days;Results通道可靠性与舒适;Discussion限制

    20健康成人、直发与指定头围,比较两种干电极帽,参考/接地仍是凝胶贴片。不同日期随机顺序,非同时同点干湿比较;通道可靠性按片段/通道,不是整台设备识别准确率。约70%且若干连接损坏,信号群体相似不等每通道等效。摘要equivalent未作等效推断;不继承所有应用或全部发型/人群外推。未取原个体数据或重算混合模型。

  4. S4 原始同步耳内/头皮设备任务比较

    Signal-specific performance of in-ear EEG: strengths and limitations

    Frei, Vanessa ; Mainar, Pablo ; Fritz, Thomas ; Chardon, Jean-Michel ; Giroud, Nathalie · 2026

    · 相关全文已核

    Methods装置/参考/过滤;Results Table1与ERP正文/表2;Discussion

    19招募,静息头皮19/耳内14(正文Results与df13支持,Methods16–19范围遗漏14);ERP17/16、言语16/16,不混分母。耳内右耳−左耳单通道加右耳ground,8通道Cyton不当8实际EEG。Biosemi传统PSD为全通道平均参考汇总,ERP用Cz;适配改为T8−T7,同时变位置/参考而不只是参考。Table1差值dB转录4行与原CI一致;d标记与差值混用未作为标准化效应。ERP正文N1/P2均值与Table2不一致,尤其适配N1正文不显著与表p=.029冲突,作为实质空位,不作精确排序。没有直接系统间检验,不能用显著/不显著推设备显著差异。头皮ICA而耳内单通道不能,接触阻抗未量,眼/心场残留及采样不独立归因。Logitech提供设备与设计支持/资金,作者关联如原文记录;未作原数据复算或全部设备泛化。

  5. S5 原始几何源重建方法,作者全文与期刊附件

    A geometry aware framework enhances noninvasive mapping of whole human brain dynamics

    Song Wang; Kexin Lou; Chen Wei; Zhiyuan Sheng; Jiahao Tang; Kaining Peng; Xinke Shen; Shuhao Mei; Liang Chen; Dongfeng Gu; Quanying Liu · 2026

    作者稿arXiv:2604.25592v1;期刊SI;GitHub ncclab-sustech/GBFs master commit99fc7c6e8afbd2e3e68a700c5f1fa1483b90f49e。当前公开实现参数/窗口未假定等于论文全批运行版。 · 相关全文已核

    arXiv2604.25592v1 Results/Discussion/Methods;期刊SI Notes1/5/6/8、Tables13–14;公开代码

    期刊2026-04-27公开页仅摘要/图/扩展,完整作者稿arXivv1 2026-04-28已取得25805436字节,期刊SI24708658字节,非声称已取期刊付费完整PDF。逐时刻线性几何模式+高斯先验/正则,主要皮层表面;基础模式300/半球实现不同表述,未用固定总数作承重。模拟真值由200fMRI元分析图生成;realistic noise按单人VEP基线协方差抽高斯,非全时间真实噪声。iEEG110患者对HCP80健康人群汇总FC非同步逐人真值。作者报iES35人318session,主文Methods刺激±5ms而当前demo_iES为±50ms且限定已知刺激半球、打印平方距离;未据此宣称重现其Table13均误差25.41mm或证明其错误。原2022来源报36/323、GBF称37初始再35/318,版本/排除差未自行猜解;且原2020Localize-MI是另一电流伪迹数据,不混称本批直接神经真值。24例仅好结局切除距离不当前瞻术前预测,切除腔不等必要EZ;4病例SOZ/重建GC也不确证动态因果。源代码仅核指定版本,不执行整套。

  6. S6 原始已知电流源定位数据说明

    Simultaneous human intracerebral stimulation and HD-EEG, ground-truth for source localization methods

    Ezequiel Mikulan; Simone Russo; Sara Parmigiani; Simone Sarasso; Flavia Maria Zauli; Annalisa Rubino; Pietro Avanzini; Anna Cattani; Alberto Sorrentino; Steve Gibbs; Francesco Cardinale; Ivana Sartori; Lino Nobili; Marcello Massimini; Andrea Pigorini · 2020

    · 相关全文已核

    Background、Electrical stimulation、Preprocessing/Source localization、Usage Notes

    Localize-MI以刺激电流伪迹为已知位置电源,原文明确非生理信号;可以验证导电/定位,不等直接恢复全部神经活动。只用原文说明验证对象,不把它混作GBF引用的2022 CCEP原始数据,不下载11GiB数据或声称复现。

  7. S7 原始颅内/头皮同步神经反应研究的全文选段

    Simultaneous stereo-EEG and high-density scalp EEG recordings to study the effects of intracerebral stimulation parameters

    S Parmigiani; E Mikulan; S Russo; S Sarasso; F M Zauli; A Rubino; A Cattani; M Fecchio; D Giampiccolo; J Lanzone; P D’Orio; M Del Vecchio; P Avanzini; L Nobili; I Sartori; M Massimini; A Pigorini · 2022

    · 相关全文已核

    原PDF页1/3/4/5,经web工具提取;Methods/数据公开/CCEP幅度时程

    原机构/Zenodo PDF本地下载403,但web工具成功读取原12页PDF相关段;落盘为web提取记录,非PDF二进制。36患者379session筛至323公开,185电极位置元数据;N1/N2为神经反应,刺激点不是全活动源真值。没有取得其补充具体伪迹处理或原数据,不推GBF±5ms窗口内到底来自何种信号,也不把数目差自行归因版本。

最后更新 2026-09-11

来源