步骤 14:SSP 投影(眼电)

import mne

# ====== 方式 A：自动构建眼电 SSP（推荐新手） ======
# MNE 会自动搜索 'eog' 类型通道并构建投影
# 注意：compute_proj_eog 实际返回 (projs, eog_events) 二元组
projs_eog, eog_events = mne.preprocessing.compute_proj_eog(
    raw,
    n_eeg=2,           # 计算几个 EEG 投影（默认 2）
    ch_name=None,      # None = 自动从 eog 类型通道选择
    average=True,      # 把双眼电平均（默认 True）
    reject=None,       # 自动拒绝阈值
)
print(f"构建的 EOG 投影: {projs_eog}")
print(f"检测到的眨眼事件: {len(eog_events)}")

# 加入到 raw.info['projs']
raw.add_proj(projs_eog)

# ====== 方式 B：显式指定双眼电通道 ======
projs_eog = mne.preprocessing.compute_proj_eog(
    raw,
    n_eeg=2,
    ch_name=['HEOG', 'VEOG'],
    average=True,
)
raw.add_proj(projs_eog)

# ====== 验证 ======
print(raw.info['projs'])
for p in raw.info['projs']:
    print(p)

# 验证 1：投影已添加到 raw
print(f"投影数: {len(raw.info['projs'])}")
assert any('eog' in p['desc'].lower() for p in raw.info['projs']), "未添加 EOG 投影"

# 验证 2：可视化投影模式（topomap）
mne.viz.plot_projs_topomap(raw.info['projs'], raw.info)

# 验证 3：施加 vs 未施加的对比
raw_proj = raw.copy().apply_proj()  # 施加投影
# 对比 raw vs raw_proj 的眼电相关时刻差异

# 验证 4：解释方差
# raw.plot_projs_topomap() 不直接显示方差，需结合 ICA 的 plot_properties