内容发布更新时间 : 2024/5/12 22:35:38星期一 下面是文章的全部内容请认真阅读。

1. 经过mriconvert转换之后，得到每个被试的4D数据，b值以及向量 2. fslroi data nodif 0 1

%choose a no weighted imaging,named nodif

从4D图像中选出一张需要的图像,0表示选择从第几个volume开始，1表示选择图像的个数

3. bet nodif nodif_brain -m -f 0.2

%generate binary brain mask（-m）,fractional intensity threshold is 0.2（-f）,result in the file nodif_brain_mask 去除头盖骨

4. eddy_correct data eddy_unwarped 0

ídy_correct <4dinput><4doutput> [] Choose interp from {trilinear,spline} def - trilinear %时间大概是10分钟 5. Dtifit

使用ctrl+alt+t调出终端，输入fsl，选择fdt*下的dtifit，在option中手动加入对应文件，包括data的输入、输出，b值文件，梯度方向（矢量）文件，共4个。运行完成后，就可以得到FA图，MD图等

------------------------------以上为预处理过程，一个人一个人的进行--------------------- --------------------%that processing is going in one subject separately----------------

6. 打开生成的FA文件所在的目录 7. tbss_1_preproc*.nii.gz

将每个人的FA图依次处理为tbss可处理的格式

Ps：可能需要新建两个文件夹（一个为FA，一个为origdata） 8. tbss_2_reg -T

将FA图抽样到1x1x1空间，空间标准化 9. tbss_3_postreg –S

把所有的FA图叠加到一起，生成4D文件mean_FA_skeleton 10. tbss_4_prestats

使骨架上包含所有的点

-------------------------------------以上为TBSS处理过程 ------------------------------------

11. 格式转换（将nii转换到matlab处理）

V=spm_vol(‘input.nii’)

[Y,XYZ]=spm_read_vols(V)； 12. 数据处理

13. 格式转换（将处理结果转换到nii）

V=spm_vol(‘input.nii’) 输入一个文件作为模板

Result=spm_write_vol（V,Yy）

----------------------------------------以上为格式转换过程 ------------------------------------

14. 相关系数的显示

fslmaths -thr 0

将大于0的部分输出

fslmaths -mul -1 -thr 0

将小于0的部分以正数输出 15. 叠加显示

fslview /桌面/fsl_data/p_r –l Red /桌面/fsl_data/n_r –l Green

-------------------------------------以上为显示的处理过程 ------------------------------------

16. 提取mask(ILF,IFOF)

左ILF 图集

/usr/share/data/jhu-dti-whitematter-atlas/JHU/JHU-ICBM-tracts-maxprob-thr25-1mm.nii.gz Intensity:13

Method:fslmaths/usr/share/data/jhu-dti-whitematter-atlas/JHU/JHU-ICBM-tracts-maxprob-thr25-1mm.nii.gz -sub 13 -abs -mul -1 -add 0.5 -thr 0 -bin LILF_mask 右ILF fslmaths

/usr/share/data/jhu-dti-whitematter-atlas/JHU/JHU-ICBM-tracts-maxprob-thr25-1mm.nii.gz -sub 14 -abs -mul -1 -add 0.5 -thr 0 -bin RILF_mask 叠加

fslmathsLILF_mask –add RILF_mask ILF_mask 17. Mask的骨架化

fslmathsILF_mask –mul skeleton_FA ILF_FA 18. Z变换处理（matlab）

------------------------------------以上为ILFmask的处理过程--------------------------------

19. 建立概率分布模型

Bedpost:使用fsl调出fdt，选择bedpostx（时间很长，5h左右） Note :

选择的文件夹中需要包含四个文件：data.nii,nodif_brain_mask, bvecs, bvals (名字必须一致)

20. 配准

使用fsl调出fdt，选择Registration，生成两个转换文件（.mat格式） Note:

需要与T1结构像，对结构像进行格式转换，Bet处理（去除头盖骨，tbss处理用到过）

需要将nodif_brain放在Bedpostx文件夹中 21. 示踪分析

使用fsl调出fdt，选择Probtrackx 1) 选择bedpostx文件夹

2) 选择种子mask（利用roi命令行生成）

-roi

zero outside roi (using voxel coordinates). Inputting -1 for a size will set it to the full image extent for that dimension. Note：

tmin为0，tsize为-1

选择种子点的图集用的是MNI _T1_2mm 3) 选择种子空间为非弥散空间

选择转换矩阵std2diff.mat 4) 点击option

选择距离校正，矩阵输出matrix2，输入tract space mask: data.nii

最后的结果存放在fdt-paths中

--------------------------------------------以上为示踪过程----------------------------------------