fNIRS脑成像测量方法及分析综述
fNIRS脑成像测量方法及分析综述
功能性近红外光谱成像(functional Near-Infrared Spectroscopy, fNIRS)是一种非侵入性的脑功能成像技术,近年来在认知神经科学、心理学、教育学和医学等领域得到了广泛应用。fNIRS通过近红外光探测脑组织中的血氧变化,提供与脑活动相关的信息。
fNIRS基于近红外光(650-900 nm)在生物组织中的相对透明性,利用近红外光照射头皮,检测反射或透射光的变化,以推测脑血氧浓度的变化,其原理基于近红外光在生物组织中的光学特性和血液中的血红蛋白的吸光特性。大脑活动时,局部血氧饱和度会发生改变,fNIRS通过监测这种变化来间接测量脑功能活动。
fNIRS的基本原理
近红外光的透射和散射:近红外光(650-900 nm)可以穿透生物组织,包括头皮、颅骨和脑组织。fNIRS利用这种特性,将近红外光照射到头皮上,通过散射和透射光的检测,获取脑组织的信息。
血红蛋白的吸光特性:血液中的血红蛋白有两种状态:氧合血红蛋白(HbO)和脱氧血红蛋白(HbR),它们在不同波长的近红外光下具有不同的吸光特性。通过测量这些吸收光谱的变化,fNIRS可以推算出局部脑血氧饱和度的变化。
光路与检测:fNIRS系统通常包括光源(如激光或LED)和探测器(光电二极管)。光源发出的近红外光通过光纤传导到头皮,部分光透过颅骨和脑组织,被组织散射和吸收后,剩余的光被探测器接收。探测器捕捉到的信号反映了血红蛋白的浓度变化。
数据处理:检测到的光信号经过数据处理,可以计算出HbO和HbR的浓度变化。这些浓度变化与脑活动有关,因为当某一脑区活跃时,该区域的血液供应会增加,导致HbO增加和HbR减少。
可测量的数据
氧合血红蛋白(HbO)浓度变化:当脑区域活动增加时,局部血流量和血氧供应增加,导致HbO浓度上升。fNIRS可以实时监测这种变化。
脱氧血红蛋白(HbR)浓度变化:相应地,由于氧气被脑细胞利用,HbR浓度会减少。fNIRS同样可以监测这种变化。
总血红蛋白(HbT)浓度变化:这是HbO和HbR之和,反映了局部血容量的变化。总血红蛋白浓度的变化可以提供关于血流量变化的附加信息。
血氧饱和度(SO2):这是HbO在总血红蛋白中的比例(HbO/(HbO+HbR)),可以反映局部脑组织的氧合状态。
常见fNIRS测量方法
- 持续波法(Continuous Wave, CW)
持续波法是fNIRS中常用的方法,主要依赖于恒定强度的近红外光源。
- 原理:使用恒定强度的光源照射头皮,通过探测器接收反射光或透射光的强度变化,以计算脑组织中血氧和脱氧血红蛋白浓度的相对变化。
- 设备:光源通常是激光或LED,探测器是光电二极管。
- 应用:应用于认知神经科学、心理学和临床医学研究。
- 优点:设备简单、成本低、易于操作。
- 缺点:无法测量绝对血氧浓度,信号容易受到外部噪声干扰。
- 时间分辨法(Time-Resolved, TR)
时间分辨法利用脉冲光源,通过测量光子在组织中传播的时间分布,获取组织的光学特性。
- 原理:发射超短脉冲光,通过探测器测量光子到达的时间分布(时间分辨),计算组织的吸收和散射参数。
- 设备:使用脉冲激光和高时间分辨的光电探测器。
- 应用:高精度空间分辨率的脑功能研究,特别是需要详细的组织光学参数的研究。
- 优点:能够提供组织的光学参数,提高空间分辨率。
- 缺点:设备复杂,成本高,对实验环境要求高。
- 频率调制法(Frequency Domain, FD)
频率调制法通过调制光源的频率,分析光在组织中的相位和振幅变化。
- 原理:使用调制频率的光源,通过测量光在组织中的相位移和振幅变化,计算组织的吸收和散射参数。
- 设备:使用调制激光或LED和锁相放大器。
- 应用:需要同时测量吸收和散射参数的研究,常用于高级脑功能成像。
- 优点:可以同时测量吸收和散射参数,提供更详细的组织信息。
- 缺点:设备较复杂,成本较高,数据分析要求高。
- 高密度fNIRS(High-Density fNIRS, HD-fNIRS)
高密度fNIRS采用高密度的光源和探测器排列,覆盖更大面积的头皮区域,提供高空间分辨率的脑功能图像。
- 原理:通过增加光源和探测器的密度,获取更多的测量点,从而提高空间分辨率。
- 设备:多通道fNIRS系统,配置大量光源和探测器。
- 应用:高级脑功能成像,特别是需要详细空间分布的研究。
- 优点:提供高分辨率的数据,适用于详细的脑功能成像。
- 缺点:数据处理复杂,需要更高的计算能力和存储空间。
- 多模态成像(Multimodal Imaging)
多模态成像结合fNIRS与其他成像技术(如fMRI、EEG)进行联合测量,以提供更全面的脑功能信息。
- 原理:同步使用多种成像技术,通过融合不同技术的数据,获得更全面的脑功能信息。
- 设备:结合fNIRS设备与其他成像设备(如fMRI扫描仪、EEG设备)。
- 应用:跨学科研究,全面分析脑功能和结构。
- 优点:多角度、多层次分析脑功能,弥补单一成像技术的不足。
- 缺点:实验设计复杂,数据融合和分析难度大,设备成本高。
- 移动fNIRS(Mobile fNIRS)
移动fNIRS开发便携式、无线化的fNIRS设备,适用于自然环境下的脑功能研究。
- 原理:使用轻便、无线的fNIRS设备,进行实时动态的脑功能监测。
- 设备:小型化、便携的fNIRS设备,配备无线数据传输模块。
- 应用:日常活动、户外实验、远程医疗等。
- 优点:灵活便捷,适合实时、动态的脑功能监测。
- 缺点:设备稳定性和数据传输的可靠性需要进一步提升。
- 功能连接性分析(Functional Connectivity Analysis)
功能连接性分析通过分析不同脑区之间的协同活动,研究脑功能网络。
- 原理:利用fNIRS数据计算不同脑区之间的相关性或同步性,揭示脑功能网络的结构和功能。
- 设备:多通道fNIRS系统,覆盖广泛的脑区。
- 应用:研究脑功能网络,探索大脑不同区域间的功能连接。
- 优点:揭示复杂的脑功能网络,帮助理解脑功能的整体性。
- 缺点:数据分析复杂,需要高级统计和计算方法。
fNIRS数据采集和处理
数据采集
- 实验设计
- 任务设计:根据研究目标设计合适的任务或刺激,确保实验具有较高的可重复性。例如,认知任务可以包括记忆、注意或语言任务,感知任务可以涉及视觉或听觉刺激。
- 基线测量:在任务开始前进行基线测量,记录受试者在静息状态下的脑血氧水平,作为后续数据分析的参考。
- 设备校准
- 光源校准:确保光源强度和波长稳定,以避免实验过程中光强度的波动影响数据质量。
- 探测器校准:调整探测器的灵敏度,确保其能准确捕捉反射光信号,避免饱和或过低的信号强度。
- 数据采集
- 选择测量位置:根据研究目标选择合适的头皮位置,通常使用国际10-20系统定位,确保覆盖目标脑区。
- 光源和探测器的配置:根据需要配置单通道或多通道系统,确保覆盖目标脑区。光源和探测器之间的距离(通常为2-4厘米)需要根据具体实验要求进行调整。
- 数据记录:使用数据采集软件记录光强度变化,并实时监控数据质量。确保采集环境稳定,避免外部光源干扰和受试者运动伪影。
数据处理
- 预处理
- 去噪:使用滤波器(如带通滤波器)去除心跳、呼吸和运动伪影等噪声。常见的滤波方法包括低通滤波(去除高频噪声)和高通滤波(去除低频漂移)。
- 信号校正:对因探头位置变化或头皮运动引起的信号波动进行校正,可以使用头动校正算法或借助其他运动传感器的数据进行校正。
- 信号分离
- 独立成分分析(ICA):将混合信号分离成独立成分,以去除伪影和提取有用信号。ICA可以有效分离出不同来源的信号,例如心跳伪影和脑活动信号。
- 主成分分析(PCA):通过降维技术减少数据的复杂性,同时保留主要信息。PCA可以帮助识别和去除噪声,提高信号的信噪比。
- 光学密度变化计算
- 原理:根据Beer-Lambert定律,将采集的光强度变化转换为光学密度变化。
- 公式:OD(t) = -log(I(t)/I0)
其中,I(t)为某时刻的光强度,I0为基线光强度。
- 血氧浓度变化计算
- 修正Beer-Lambert定律:利用修正后的Beer-Lambert定律计算氧合血红蛋白(HbO)和脱氧血红蛋白(HbR)的浓度变化。
- 公式:ΔC = (ΔOD/ε) * (1/d) * DPF
其中,ΔC为浓度变化,ΔOD为光学密度变化,ε为吸收系数,d为光程长度,DPF为光程校正因子。
- 统计分析
- 任务块平均:对任务期间的数据进行平均,消除随机噪声的影响。通过分块平均,可以提高数据的稳定性和信噪比。
- 事件相关平均:对多个任务事件的数据进行平均,提高信噪比。这种方法适用于研究快速、瞬时的脑功能变化。
- 统计检验:使用统计方法(如t检验、ANOVA)检验任务与基线之间的差异是否显著。显著性检验可以帮助确定实验结果的可靠性。
- 结果可视化
- 时间序列图:绘制脑血氧浓度随时间变化的图表,直观展示任务相关的脑功能变化。常见的时间序列图包括HbO和HbR浓度随时间变化的曲线图。
- 空间分布图:结合头皮位置绘制脑功能活动的空间分布图,展示不同脑区的活动情况。可以使用头皮顶视图或3D头模图展示数据。
- 三维图像:使用三维重建算法,生成大脑内部的功能活动图像。这种方法可以提供更详细的脑功能分布信息。
fNIRS功能近红外数据处理工具包
- NIRS_SPM
NIRS_SPM基于SPM工具包和Matlab平台,可以识别目前主流fNIRS厂商设备采集的数据格式,它基于一般线性模型GLM对fNIRS进行统计分析和统计(一阶分析和组分析)。此外,NIRS_SPM还具有对数据的简单预处理(比如说滤波等)、通道的定位(如输出每个通道的MNI坐标等)、ROI分析、时间序列分析等功能。
- Homer2/Homer3
Homer是一款基于MATLAB的工具包,是fNIRS数据的分析和显示的一个常用开源平台。其包含有丰富的处理函数可以通过GUI界面灵活选择设置处理流,从而对fNIRS 数据进行预处理。单就运动伪迹去除就提供了多种常用的方法,如样条插值(Spline)、主成分分析(PCA)、小波分析(Wavelet)等。论是Homer3还是Homer2(目前不更新及维护了),但这两个的**功能和函数是不变的。
- FC-NIRS
FC-NIRS 是一款用于近红外光谱数据的功能连通性分析工具,该软件由北京师范大学开发设计。它的功能包括功能连接分析、可视化和网络分析,还具有信号质量控制和批处理功能。该工具包包含Homer,因此可以进行预处理。
- fNIRSOFT
fNIRSOFT 是一个软件包,旨在通过图形用户界面和/或脚本(用于自动化)处理、分析和可视化功能性近红外光谱信号。
该软件是一款付费软件,分为标准版和专业版。标准版具有基本的数据处理功能,如通过 GUI 进行时间可视化、基本降噪、预处理、时间序列分析和空间可视化等。而专业版额外增加了自动信号质量检查,消除饱和和有问题的通道,运动伪影去除算法,用于特征提取的高级信号处理算法,大脑表面图像上的大脑映射和可视化,具有阈值、动画(时间变化)或组/主题/条件平均值的左/右/背视图。
- NeuroDOT_Beta
NeuroDOT_Beta是一个可扩展的Matlab 工具箱,与Matlab2015版兼容,包含 4 个文件夹:Data、Documentation、Functions和 Support_Files。
- Nirstorm
Nirstorm是Brainstorm工具包的一款开源插件,主要用于EEG/MEG的溯源分析,但是借助如Nirstorm可以用于fNIRS数据的分析。该软件具有滤波、降噪、运动矫正、MBLL、窗口平均,使用GLM进行统计分析等多种功能。
- PHOEBE
PHOEBE(Placeing Headgear Optodes Efficiently Before Experiment) 是一款基于MATLAB的开源工具包,具有GUI界面,可实时测量和显示 fNIRS 光极与受试者头皮之间的光耦合。PHOEBE 的主要目标是帮助实验者以高效的方式优化 fNIRS 光极的信号质量。
相关工具链接
- NIRS_SPM:https://www.nitrc.org/projects/nirs_spm/
- Homer2:https://homer-fnirs.org/
- Homer3:https://github.com/BUNPC/Homer3
- nirsLAB:https://nirx.net/nirslab-1
- FC-NIRS:https://www.nitrc.org/projects/fcnirs/
- fNIRSOFT:https://www.biopac.com/product/fnirsoft-professional-edition/?attribute_pa_size=fnir-software-professional-edition
- NeuroDOT_Beta:https://github.com/WUSTL-ORL/NeuroDOT_Beta
- Nirstorm:https://github.com/Nirstorm/nirstorm
- NIRS-KIT:https://www.nitrc.org/frs/?group_id=1423
- Fieldtrip:https://www.fieldtriptoolbox.org/download/
- AtlasViewer(一般与Homer2/Homer3结合使用):https://homer-fnirs.org/
- NIRS-Toolbox:https://bitbucket.org/huppertt/nirs-toolbox/src/master/
- fOLD:https://github.com/nirx/fOLD-public
- PHOEBE:https://bitbucket.org/lpollonini/phoebe/wiki/Home
