当前位置：首页 > news >正文

FreeSurfer的脑区表面重建与注意力nii图结合

news 来源：原创 2025/6/14 14:14:54

将FreeSurfer的脑区表面重建与注意力nii图结合，可通过以下步骤实现数据配准、投影和可视化。以下是详细流程及工具使用方法：

一、FreeSurfer预处理原始MRI数据

1. 准备原始MRI数据

确保原始MRI数据为NIfTI格式（.nii/.nii.gz），建议使用T1加权像（T1w MRI），因为FreeSurfer对T1像的皮层分割效果最佳。

2. 运行FreeSurfer的recon-all流程

通过命令行执行完整的预处理和表面重建：

# 假设数据路径为/path/to/mri.nii.gz，subject名为"patient01"
recon-all -i /path/to/mri.nii.gz -s patient01 -all -openmp 8

参数说明：
- -i：输入MRI文件路径。
- -s：指定subject名称，结果会保存在$SUBJECTS_DIR/patient01目录下。
- -all：执行所有处理步骤（颅骨剥离、皮层分割、表面重建等）。
- -openmp：设置并行线程数，加快计算。

3. 生成的关键表面文件

处理完成后，在patient01/surf/目录下会生成：

lh.pial/rh.pial：皮层灰质-脑脊液界面（pia表面）。
lh.white/rh.white：皮层灰质-白质界面（white表面）。
lh.sphere/rh.sphere：标准化球面，用于皮层配准。

二、将注意力nii图配准到FreeSurfer空间

1. 确认注意力nii图与MRI的空间一致性

如果注意力nii图是基于同一MRI数据生成的（如训练时使用的原始MRI），则可能已在同一空间；若不是，需进行空间配准：

# 使用FreeSurfer的mri_coreg工具进行刚体配准
mri_coreg --i /path/to/attention.nii.gz --d /path/to/mri.nii.gz --o registered_attention.nii.gz --lta coreg.lta

若配准精度要求高，可使用FSL的FLIRT/ANTs进行非线性配准。

2. 将注意力图转换到FreeSurfer的Native空间

# 使用mri_convert转换格式并对齐到FreeSurfer空间
mri_convert registered_attention.nii.gz patient01/mri/attention.mgz

三、将注意力值投影到脑表面

1. 使用mri_vol2surf投影体数据到表面

# 投影到左半球pial表面
mri_vol2surf --surf patient01/surf/lh.pial --vol patient01/mri/attention.mgz --o patient01/surf/lh.attention.mgh --method trilinear# 投影到右半球pial表面
mri_vol2surf --surf patient01/surf/rh.pial --vol patient01/mri/attention.mgz --o patient01/surf/rh.attention.mgh --method trilinear

参数说明：
- --method trilinear：使用三线性插值，提高投影精度。
- 输出文件.mgh为FreeSurfer表面数据格式，存储每个顶点的注意力值。

2. （可选）标准化到标准脑模板

若需在标准空间（如MNI）中可视化，可通过球面配准将注意力投影到标准脑表面：

# 左半球投影到标准球面
mri_surf2surf --srcsubject patient01 --srcsurf lh.pial --trgsubject fsaverage --trgsurf lh.pial --sval patient01/surf/lh.attention.mgh --tval fsaverage/surf/lh.attention_std.mgh --hemi lh --regmethod球面# 右半球同理

四、使用Freeview可视化注意力表面

1. 启动Freeview并加载表面和注意力数据

freeview -f patient01/surf/lh.pial:color=white,patient01/surf/rh.pial:color=white \-surf-annot patient01/surf/lh.aparc.a2009s+aseg:name=aparc,patient01/surf/rh.aparc.a2009s+aseg:name=aparc \-surf-Stat patient01/surf/lh.attention.mgh:overlay_color=hot:overlay_threshold=0.1:name=attn_lh,patient01/surf/rh.attention.mgh:overlay_color=hot:overlay_threshold=0.1:name=attn_rh \-viewport 3d

界面操作说明：
- 在Freeview中，通过左侧面板调整表面显示参数：
  - Surface标签页：选择显示pial或white表面。
  - Annotations标签页：加载脑区标注（如aparc.a2009s，显示Broca区、Wernicke区等）。
  - Stats/Overlays标签页：添加注意力数据（.mgh文件），设置颜色映射（如hot色系）和阈值（如0.1，过滤低权重区域）。

2. 自定义可视化参数

颜色映射：在Stats/Overlays中选择Color Table为hot或jet，突出高注意力区域（红色/黄色）。
阈值调整：通过Min/Max滑动条过滤低权重值，仅显示感兴趣的高注意力脑区。
叠加脑区标注：在Annotations中选择aparc+aseg，可查看注意力与解剖脑区的对应关系（如注意力集中在左额下回）。

五、进阶：结合脑图谱分析注意力分布

1. 使用mri_label2surf关联注意力与脑区

# 提取左半球Broca区（BA44）的注意力均值
mri_label2surf --label patient01/label/lh.BA44.label --surf patient01/surf/lh.pial --stat patient01/surf/lh.attention.mgh --out BA44_attention_mean.txt

2. 批量分析各脑区注意力值

# 使用FreeSurfer的mri_segstats工具计算各脑区注意力统计量
mri_segstats --seg patient01/mri/aparc+aseg.mgz --stat patient01/mri/attention.mgz --out stats.txt

注意事项

数据预处理要求：
- 原始MRI需确保无运动伪影，FreeSurfer的recon-all流程建议完整运行（约需数小时）。
- 注意力nii图若为浮点数（如0-1权重），需确认格式兼容性（.mgz支持浮点，.nii.gz需确保数据类型正确）。
配准精度影响：
- 若注意力图与MRI未配准，投影结果可能偏移，建议通过目视检查（如在Freeview中叠加MRI和注意力图）确认配准质量。
可视化优化：
- 高分辨率表面（如512k顶点）可提高注意力分布的显示精度，但会增加计算负担。
- 若注意力值分布不均，可使用对数变换（log1p）或分位数标准化，增强低权重区域的可视性。