edickie/ciftify项目：将FreeSurfer recon-all输出转换为HCP风格CIFTI格式的完整指南

概述

edickie/ciftify项目中的ciftify_recon_all工具是一个强大的数据处理管道，它能够将FreeSurfer的recon-all输出转换为符合人类连接组项目(HCP)风格的CIFTI格式文件结构。这一转换过程对于希望使用HCP标准分析流程的研究人员来说至关重要。

核心功能：ciftify_recon_all

基本工作原理

ciftify_recon_all工具基于HCP最小处理流程的核心思想，但进行了重要改进：

1. 它不要求原始数据必须包含HCP标准采集序列
2. 减少了中间文件的生成
3. 自动生成详细的处理日志文件

基本使用方式

运行该工具需要三个关键参数：

1. FreeSurfer的SUBJECTS_DIR：包含所有被试FreeSurfer输出的顶级目录
2. CIFTIFY_WORKDIR：将存放HCP风格输出的工作目录
3. 被试ID：需要处理的被试标识符，必须与FreeSurfer目录中的文件夹名称匹配

示例命令

# 通过环境变量定义路径
export SUBJECTS_DIR=/path/to/freesurfer/outputs
export CIFTIFY_WORKDIR=/path/for/hcp/outputs
ciftify_recon_all Subject001

# 直接通过参数指定路径
ciftify_recon_all --fs-subjects-dir /path/to/freesurfer/outputs \
                 --ciftify-work-dir /path/for/hcp/outputs \
                 Subject001

高级配置选项

表面配准方法选择

工具提供两种表面配准方法：

1. MSMSulc（默认）：

   • 与当前HCP流程保持一致
   • 计算资源密集（可能需要3小时）
   • 需要最新版MSM软件

2. FreeSurfer配准（FS）：

   • 使用FreeSurfer的fsaverage表面配准
   • 计算速度更快
   • 使用--surf-reg FS参数启用

T1w空间32k灰度坐标输出

对于DTI分析等特定应用，可能需要生成T1w空间的32k灰度坐标输出。默认情况下不生成这些文件，需要通过--resample-to-T1w32k参数显式启用。

ciftify_recon_all --resample-to-T1w32k \
                 --fs-subjects-dir /path/to/freesurfer/outputs \
                 --ciftify-work-dir /path/for/hcp/outputs \
                 Subject001

输出结构详解

转换后的输出遵循HCP风格的组织结构，主要包含以下关键目录和文件：

<subject>
├── MNINonLinear
│   ├── fsaverage_LR32k        # 32k表面网格数据（用于fMRI和多模态分析）
│   ├── Native                 # 原生FreeSurfer网格
│   ├── ROIs                   # 皮层下重采样ROI
│   └── xfms                   # T1w空间的FSL变换
└── T1w
    ├── fsaverage_LR32k        # T1w空间的32k表面网格（如启用了--resample-to-T1w32k）
    └── Native                 # 原生FreeSurfer网格

理解"空间"概念

HCP输出数据按"空间"组织，每个空间由两个关键要素定义：

1. 体积空间：数据在三维空间中的位置
2. 表面网格：表面顶点的分布密度

主要空间类型包括：

体积空间	表面网格	典型应用场景
MNINonLinear	32k_fs_LR	fMRI和多模态分析，大多数图谱可用
T1w	32k_fs_LR	DTI分析的种子点，解剖测量
MNINonLinear	164k_fs_LR	高分辨率解剖分析
Native	FreeSurfer	原始FreeSurfer网格

数据可视化与质量控制

cifti_vis_recon_all工具

该工具用于生成质量控制页面，分为两个步骤：

1. snaps步骤：为单个被试生成QC图像

   cifti_vis_recon_all snaps Subject001
   

2. index步骤：创建汇总HTML页面

   cifti_vis_recon_all index
   

QC工作流程建议

1. 浏览组合QC视图，标记异常图像
2. 检查矢状面表面轮廓视图
3. 点击异常图像查看单被试详细视图
4. 记录所有质量问题

常见质量问题示例

• 处理未完成：表面文件缺失
• 极低质量解剖数据：表面皱缩
• 脑掩膜失败：部分脑区异常拉伸
• 枕叶灰质缺失：脑后部看起来分裂

技术细节

神经影像数据类型

HCP/CIFTI格式使用特定文件扩展名区分数据类型：

数据类型	典型内容	CIFTI扩展名	GIFTI扩展名
时间序列	fMRI数据	.dtseries.nii	.func.gii
标量/度量	厚度图,统计图	.dscalar.nii	.shape.gii
标签	ROI,图谱	.dlabel.nii	.label.gii

大脑表面类型

HCP命名约定：<subject>.<hemisphere>.<surface_type>.<surface_mesh>.surf.gii

表面类型	描述	主要用途
white	灰质和白质边界	定义皮层带内侧
pial	灰质和脑外边界	定义皮层带外侧
midthickness	white和pial的中间点	距离和表面积测量
sphere	球面表示的顶点	配准和重采样
inflated	膨胀的midthickness	可视化
very-inflated	更膨胀的midthickness	可视化
flat	展开的皮层	全半球可视化

总结

edickie/ciftify项目提供的工具链为将FreeSurfer输出转换为HCP兼容格式提供了完整解决方案。通过ciftify_recon_all和cifti_vis_recon_all工具，研究人员可以：

1. 标准化数据处理流程
2. 生成多种空间的数据表示
3. 进行全面的质量控制
4. 为后续的HCP风格分析做好准备

这套工具特别适合那些希望利用HCP先进分析方法但数据采集不完全符合HCP标准的研究项目。