血管内脑电图（vECoG）中通过替代性重参考去除Stentrode上心电伪迹同时保留运动特征

《Advanced Science》：Preserving Motor Features by Alternative Re-Referencing to Remove Heart Artifact on the Stentrode

　　

血管内脑电图（vECoG）作为穿透式电极阵列的微创替代方案在神经信号采集方面展现出巨大潜力。本研究以主流vECoG平台Stentrode为对象，考察其信号质量，对比心电（ECG）伪迹残留与低频运动皮层活动保留之间的权衡。典型（重）参考方案，包括单极支架挂载参考、共同平均参考（CAR）与拉普拉斯参考，相较远端单极参考均能显著降低ECG伪迹，但同时静息态β活动也被明显削弱。研究人员采用带限独立成分分析（BL-ICA）——一种允许对各电极噪声进行差异化加权的空间滤波方法——可轻松将ECG伪迹从vECoG记录中分离。该清洗方法重构信号时剔除ECG成分，并评估通道间相关性的降低及静息与执行（go）分布间相对熵的提升。为确认低频运动特征得以保留，研究在源空间与重构信号空间均评估β爆发（beta bursts）特征。结果表明，BL-ICA能有效去除广泛分布的ECG伪迹，同时维持β频段活动中典型的运动相关特征。

研究背景与立项依据

脑机接口（BCI）为严重运动障碍患者恢复独立性带来希望，其家用可靠性取决于所采信号质量。非侵入式脑电图（EEG）因穿越头皮、颅骨等多层组织而受限；侵入式皮层电图（ECoG）与颅内微电极分辨率高但需开颅；以Stentrode为代表的血管内支架电极阵列（vECoG）经颈静脉植入上矢状窦（SSS）毗邻运动皮层，兼顾微创与较高信噪比。然而Stentrode位于SSS内持续受血流搏动影响，vECoG记录常被强心电（ECG）伪迹污染，伪迹能量覆盖0.5–40 Hz，与α、β等皮层节律重叠，传统硬件参考（远端单极参考、支架本地单极参考）与离线重参考（CAR、拉普拉斯）虽压低ECG幅度，却同步削去低频β振荡，导致后续解码特征受损。既有宽带独立成分分析（ICA）去伪迹在vECoG上效果不稳定，且γ频段区分度下降。因此亟需一种在去除ECG同时保留运动皮层β特征的方法。

论文发表于《Advanced Science》，研究对象为COMMAND EFS临床试验（NCT05035823）植入Stentrode的5例有效受试者（P1、P3–P6，P2仅含基线远端参考数据不进入主分析），病因含肌萎缩侧索硬化（ALS）与脑干梗死，电极取阻抗<100 kΩ的9–13通道，采样2 kHz、硬件低通500 Hz。研究人员得出核心结论：远端单极参考虽保留最多β功率却含最强ECG；CAR与拉普拉斯降ECG但β振荡近似被抵消；本地参考抑ECG有限且低频活动受损；提出以远端参考预处理结合带限独立成分分析（BL-ICA）的流程，可在α、β、低γ频段分离出符合生理R-R间期与QRS模板的ECG独立成分，置零对应混合矩阵权重后重构信号，在剔除ECG同时保留β事件相关去同步化（ERD）与β爆发（beta bursts），跨通道零滞后相关下降、静息-执行相对熵上升，意义在于为vECoG提供兼顾伪迹抑制与运动特征保真的离线清洗范式，支撑后续实时BCI特征提取。

主要关键技术方法概括

研究人员基于COMMAND EFS试验6例受试者（5例进入主分析）的Stentrode vECoG数据，硬件端对比远端单极与支架本地单极参考，离线施加CAR与拉普拉斯重参考；预处理经60/120 Hz二阶IIR陷波，按α（8–13 Hz）、β（13–30 Hz）、低γ（30–70 Hz）、高γ（70–200 Hz）五阶巴特沃斯带通；采用FastICA对各频带分别拟合带限独立成分分析（BL-ICA），以BioSPPy校验R-R间期与QRS模板判定ECG成分并置零重构；借助FOOOF分解PSD的振荡与非周期成分估测β振荡功率d_β；以零滞后跨通道相关、静息-执行相对熵、Bycycle周期级β爆发检测量化特征保留情况。

研究结果

3.1 传统参考方案未必保留皮层运动特征

通过对P1–P6基线PSD与ECG峰值比对，远端参考ECG幅度最高但β偏离非周期基线最大（d_β可达79–155 μV²/Hz），CAR与拉普拉斯使ECG降幅67%–88%却将d_β压至近0或个位数；本地参考ECG残存且β峰消失。结论：典型重参考以牺牲低频运动振荡换取伪迹抑制。

3.2 盲源分离实现噪声解耦

对各受试者各频带遍历成分数，取等于良导联数时重构L2误差近零；BL-ICA在α、β、低γ源空间稳定浮现1–2个ECG成分，其R-R分布落生理区间（P5、P6偏快符合ALS窦性心动过速报道），高γ因频带高不易染ECG无须处理。结论：频带限定后ICA源数可控，ECG可作为独立源被隔离。

3.3 源空间保留任务调制运动特征

仅P1、P6有合格远端参考运动数据。其BL-ICA源空间中，α/β成分在go线索后包络持续压抑、低/高γ呈5峰同步化，对应5次尝试动作；最优β成分经Bycycle检出静息多爆发、go期爆发稀疏，ERD与周期一致性保留。结论：任务相关调制不在被删ECG成分上，而在独立神经源中。

3.4 重构降低共享噪声并提升静息/执行可分性

剔除ECG成分重构后，α/β/低γ静息零滞后跨通道相关普遍下降；各通道静息-执行幅度分布相对熵H_R(P_rest‖P_go)上升；观测空间检出的β爆发幅度由异常高值回落至生理区间。结论：BL-ICA去除跨通道共有ECG基底，使真实β爆发显影、两状态分布更易由解码器区分。

讨论与结论翻译

典型参考中CAR假设同构噪声模型在Stentrode圆柱双侧不对称布局下不成立，拉普拉斯受几何限制，本地参考误删广域低频神经活动。远端参考虽噪大，却因ECG源强而利于BL-ICA分离；带限步骤将源数降至通道数内，避免宽带多源混叠。ALS受试者R-R偏短与已知心律失常一致，佐证ECG成分判定可靠。跨通道相关降、相对熵升、β ERD与爆发保真三者一致表明BL-ICA对远端参考vECoG可分离神经与心脏源。局限含远端参考数据量不均、仅P1/P6做源空间演示、离线评估未接解码器；但相对熵与β特征均为解码相关指标。研究者结论：对Stentrode vECoG，采用远端单极参考采集、继以逐带BL-ICA剔除ECG成分，能在大幅抑伪迹的同时保留β频段运动皮层特征；若未来设硬件单极参考，应置于IRTU与支架体之间而非SSS内，以兼顾特征保真与伪迹可分离性。