提取DWI数据的FA和MD

DWI简介

        扩散加权磁共振成像（DWI ）是使用特定的 MRI 序列以所获得的脑成像数据，该成像手段依靠水分子的扩散在 MR 图像中产生对比度。优于组织中的分子扩散不是自由的，而是反映了与许多障碍物（例如大分子，纤维和膜）的相互作用。因此，水分子扩散模式可以揭示有关正常或处于患病状态的组织结构的微观细节。弥散张量成像（DTI）是一种特殊的 DWI，已被广泛用于绘制大脑中的白质束图

FA和MD简介

部分各向异性分数(fractional anisotropy, FA)

        FA测量的是白质的完整性，指水分子各向异性成份占整个弥散张量的比例，它的变化范围从0-1。0代表弥散不受限制，比如脑脊液的 FA 值接近 0；对于非常规则的具有方向性的组织，其FA值大于0，例如大脑白质纤维 FA 值接近 1。

平均弥散率(mean diffusivity, MD)

        MD反映分子整体的弥散水平和弥散阻力的整体情况。MD只表示弥散的大小，而与弥散的方向无关。MD 越大，组织内所含自由水分子则越多。

FA和MD提取流程

1.从HCP官网下载进行了涡流校正之后的DWI数据

2.提取DWI的b0像——使用fslroi函数

#!/bin/bash

DWI_dir=/media/pc02/txt/HCP_DWI

for subj in $(ls $DWI_dir)
do

        dir_temp=$DWI_dir/$subj/${subj:0:6}/T1w/Diffusion_7T
        
        fslroi $dir_temp/data.nii.gz $dir_temp/b0.nii.gz 0 1

done

原始DWI图像

3. 将b0像剥去颅骨——使用bet2函数

#! /bin/bash

DWI_dir=/media/pc02/txt/HCP_DWI

for subj in $(ls $DWI_dir)
do
	
	dir_temp=$DWI_dir/$subj/${subj:0:6}/T1w/Diffusion_7T
	
        bet2 $dir_temp/b0.nii.gz $dir_temp/b0_bet.nii.gz -o -m -f 0.2
        #-o生成叠加图像，-m生成掩模，-f设置分割强度阈值（默认0.5）

done

去颅骨DWI图像

4.创建FA和MD图——使用dtifit函数

#!/bin/bash

DWI_dir=/media/pc02/txt/HCP_DWI

for subj in $(ls $DWI_dir)
do
	
	dir_temp=$DWI_dir/$subj/${subj:0:6}/T1w/Diffusion_7T	
	dtifit -k $dir_temp/data.nii.gz -o $dir_temp/data -m $dir_temp/nodif_brain_mask.nii.gz -r $dir_temp/bvecs -b $dir_temp/bvals --save_tensor
	
done

FA图像

       

MD图像