如何使用FreeSurfer 5.3.0处理三维T1加权成像数据,以获得强度不均匀性校正的T1加权图像和白质(WM)掩模
时间: 2024-09-16 14:01:08 浏览: 137
freesuffer
在FreeSurfer 5.3.0中处理三维T1加权成像数据并得到强度不均匀性校正的T1加权图像(通常称为bias-corrected T1)以及白质(WM)掩模,可以按照以下步骤操作:
1. **安装和初始化 FreeSurfer**:
- 首先,确保已经下载并安装了最新版本的FreeSurfer,可以从官网(https://surfer.nmr.mgh.harvard.edu/)获取。
- 运行`recon-all`命令初始化一个新的结构化分析目录,提供你的MRI扫描作为输入。
2. **预处理 MRI 数据**:
- 使用`recon-all`工具的`gunzip_anat`功能解压T1w图像,并对其进行基本的质量检查,如去噪和平滑。
```
gunzip_anat subject_id
```
3. **强度不均匀性校正(Bias Field Correction, BFC)**:
- 通过`apply_bias_field`命令对原始T1w图像进行校正,生成校准后的T1w图像。这一步会去除磁场强度和设备造成的偏移。
```
apply_bias_field subject_id
```
4. **分割和配准**:
- `recon-all`会自动运行`recon-all`, 其中包括了T1脑部分割 (`recon-all -autorecon1`) 和空间标准化 (`recon-all -all`)。这个过程将生成一系列的脑部分割结果,包括WM掩模。
5. **查看和导出 WM 掩模**:
- WM掩模通常保存在`surf/lh.white`和`rh.white`文件中,分别对应左侧和右侧大脑半球。你可以使用`mri_vol2surf`命令将其转换为表面形状文件或导出为nifti格式。
```
mri_vol2surf subject_id lh.white -i > lh.white.surf.gii
mri_vol2surf subject_id rh.white -i > rh.white.surf.gii
```
6. **导出 WM 密度图(if needed)**:
如果需要密度图而非表面掩模,可以使用`mris_anatomical_stats`命令计算平均灰度值。
```
mris_anatomical_stats subject_id lh.white aseg.stats
```
记得替换`subject_id`为你实际的受试者ID。在处理过程中,可能需要根据FreeSurfer文档调整参数,尤其是对于特殊类型的T1加权图像。完成上述步骤后,你将得到校正过的T1图像和WM掩模。
阅读全文