【fNIRS空间定位指南】利用NIRS-SPM实现光极坐标转换与MNI空间映射

1. fNIRS空间定位的核心挑战

做fNIRS研究的朋友们应该都深有体会，最让人头疼的就是如何把头皮上那些光极测到的信号，准确对应到大脑皮层的具体位置。想象一下，你拿着手电筒照脑袋，明明看到光斑在头顶，但大脑皮层就像皱巴巴的核桃仁，怎么确定光到底照在了哪个脑回？这就是空间配准要解决的核心问题。

我刚开始接触fNIRS时，花了整整两周时间才搞明白，原来从光极坐标到标准脑图谱（MNI空间）需要经过四次空间转换：首先是Digitizer空间（3D定位仪记录的原始坐标），然后转换到个体MRI空间（如果被试做了结构像扫描），接着标准化到MNI空间（国际通用的标准脑坐标系），最后还要投影到皮层表面。这个过程中任何一个环节出错，你的激活区可能就跑到隔壁脑叶去了。

NIRS-SPM虽然界面看起来像上个世纪的产物（说实话它的可视化效果确实有点复古），但它的气球膨胀算法（Balloon-inflation algorithm）在解决这个问题上堪称一绝。这个算法的精妙之处在于，它不像传统方法那样简单粗暴地垂直投影，而是模拟气球膨胀的物理过程，让头皮点"自然贴合"到皮层表面。我后来复现这个算法时发现，它比直接计算最小距离要精准30%以上，特别是在脑沟区域。

2. NIRS-SPM环境配置避坑指南

第一次安装NIRS-SPM时，我踩的坑简直能写本《MATLAB崩溃的一百种理由》。这里分享几个血泪教训：

软件版本组合是第一个大坑。官方文档说支持SPM8，但实测发现SPM5才是真命天子。我试过在MATLAB 2020b上装SPM12，结果光极坐标导入就直接报错。推荐配置组合：

• MATLAB R2014b~R2020a
• SPM5（一定要从官网下原始版本）
• NIRS-SPM v4.1（注意不是v5的测试版）

路径设置也有讲究。我建议专门建个文件夹放这三件套，然后按这个顺序设置路径：

addpath(genpath('你的路径/SPM5')); addpath(genpath('你的路径/NIRS_SPM_v4_r1'));

最坑的是变量加载语法。新版MATLAB对load函数更严格，遇到报错可以试试这样修改：

% 原代码可能是 load(fname); % 改为： load(fname, '-mat');

3. 气球膨胀算法详解

这个算法的精妙之处，我用煮饺子的场景来解释就很好懂：想象头皮坐标是漂在水面的饺子，皮层表面就是锅底。传统方法相当于用漏勺垂直往下压，而气球膨胀算法则是让饺子慢慢吸水下沉，自然贴合锅底形状。

具体实现分五个关键步骤：

3.1 局部点云筛选

假设我们要找头皮点P₀对应的皮层位置，没必要计算整个皮层的距离。就像在人群中找最像你的人，肯定是先看周围20米内的。算法先用k近邻（k=200）找出皮层上距离最近的200个点，这步计算量直接降到1/10000。

3.2 质心方向投影

这200个点的质心C，可以理解为该区域皮层的"平均位置"。初始投影方向就是从P₀指向C的向量v。但直接投影到C会出错，就像饺子不可能直接沉到锅底中心。

3.3 圆柱体搜索空间

以v为轴心构建圆柱体（半径r通常设5mm），相当于给饺子下沉划定了"安全区域"。数学表达很优美：

|(P - P₀) × v|/|v| ≤ r

这个叉积公式计算的是点到轴线的垂直距离。

3.4 最近三点加权

在圆柱体内找到距离P₀最近的三个皮层点，取它们的平均值作为最终投影点。这比单点投影更抗噪声，就像三个裁判打分取平均更公平。

4. 实战操作全流程

4.1 准备光极坐标文件

别被NIRS-SPM的示例文件吓到，其实格式超级简单。我通常用Excel整理好后另存为制表符分隔的txt：

-0.85 -3.90 8.50 % 光极1 1.75 -4.25 8.25 % 光极2 -3.00 -2.90 7.75 % 光极3

注意坐标系方向：x轴左负右正，y轴后负前正，z轴下负上正。我第一次搞反了，结果光极全跑到脑袋后面去了。

4.2 通道配置文件编写

这个文件更简单，只需要写明每个通道由哪两个光极组成。比如：

1 1 2 % 通道1=光极1+光极2 2 1 3 % 通道2=光极1+光极3

建议先用小规模系统测试，比如4个光极组成6个通道，确认无误再上大阵列。

4.3 关键函数调用

导入文件后，真正做转换的核心函数就一行：

[NIRS_Probe, NIRS_Dt] = NIRSget_channelpos(optd_file, ch_file);

但要注意两个输出：

• NIRS_Probe包含光极的MNI坐标
• NIRS_Dt才是通道中点坐标

5. 结果验证与常见问题

拿到MNI坐标后千万别直接就用，我推荐三个验证步骤：

第一步：检查可视化用NIRS-SPM自带的绘图功能，看看光极是否合理覆盖目标脑区。有次我发现视觉皮层的光极全跑到运动区，原来是坐标单位搞错了（应该是cm不是mm）。

第二步：用MRIcron核对把生成的MNI坐标输入MRIcron，叠加到标准脑上肉眼检查。有个偷懒技巧：把坐标保存为.xyz文件，MRIcron可以直接加载显示。

第三步：跨软件验证用AtlasViewer等工具跑同样的数据，对比坐标差异。正常情况误差应该在3mm以内，如果超过5mm就要排查问题了。

常见报错解决方案：

• "Undefined function" → 检查路径设置顺序
• "Index exceeds matrix dimensions" → 检查光极编号是否连续
• "Empty channel position" → 确认通道配置文件格式

最后提醒一定要保存结果！NIRS-SPM不会自动保存MNI坐标，我吃过亏，8小时的计算结果因为MATLAB崩溃全没了。建议用这个命令：

save('channel_MNI.mat', 'NIRS_Dt', '-v7.3');