DtiStudio读梯度表.bvec的方法步骤

1、DtiStudio是DTI数据处理的一种常用工具，大家在用DtiStudio读入数据时，要经过一个读入梯度表的过程，如下图1所示：GracienlTabs0:0.000,0.000,0ODO-1000.0.0000OOD2.-0.89,0.520,O.OOD3:0.1OS0:0.000,0.000,0ODO-1000.0.0000OOD2.-0.89,0.520,O.OOD3:0.1OS5GE.0S18-0.891-0.35,-0.315斤1JD3-073Gr6776.0.868,-0238,0.367:-0.799,0.370,Q.47S6.0.102,0.967,0ODO9:-0.8GG,

2、-0.129,0.133Hints:Thegradientiadeinedas:zGradient(0A0)meansaBOiinaga.(100.100,100)wllexcludtheimagefromtensorcalcuIatonSelectspre-definedgradenttabls不同品牌MRI机型的原始数据格式是不同，有的厂商如飞利浦的数据的梯度文件.grad可以直接用，用写字板打开，复制粘贴，轻松简单；可其它不少厂商如GE的数据的梯度表不能直接读入，这就不好办了。DTI实验一般扩散方向至少要选择20来个，如果要是手工输入的话那可老费事了，而且不同被试不同实验使用的梯度可能都

3、不同，这样工作量太大了。我花了一些时间上网查资料，总结出了一套读梯度的步骤，分享给大家。一、获取原始梯度表。大家在处理DTI数据之前，一般都要先将原始DICOM文件转化成NIfTI格式文件，它包含四个子文件，后缀分别是.bval.bvec.hdr.img。这个.bvec里面就是梯度表信息，用记事本打开发现是3*N的矩阵，而DtiStudio要求的梯度表却是N*3的矩阵，这样就需要一个矩阵转置的过程。二、用MATLAB对原始梯度表进行转置。大家可能有更好的转置方法，而我就会MATLAB，这里我就展示用MATLAB处理的步骤。首先是读取.bvec文件，在MATLAB命令窗口输入uiopen，选择要

4、处理的.bvec文件，会弹出一个文本，显示一个3*N的矩阵，如下:|1-ld.LdB-4.S34.OLdOd&-4.TQbDliU-OLdM-4.ME-4.EE4.43EOlEM4.E2EOlB-OlKSd.GibSOlM00D.E290.5RI-fl.735-A.23B山570MT-IX.125-0-CLBS2-D.1E2a.7S山M-D-D.D5-fl.075-fl0BP20.74.&43-D.53D.MSQ0Da.Bia-0L31Ga.S?7L4M4.47ED0.4&J丄23-OLfifli0.rD64.192-OlffiO.fM-fl.FEE-fl.flTB-4.EB4-0.366!0

5、L14B-OlTlJ4.202亿圈P3为了转置，首先要将上述梯度表变成矩阵，即套上;如下:3a-LYLa.pB-D.3B4乩口2D-E&3-.735fl.162-D_EtaCl2l2-0.066-D-55&.J3SD.53D.6553-i.MTQ.4.L3Q.舶白-D_462L502T.玄N-0.2&7D-Q4Odd.52B0.5fl5-CLT7W-5.23Bft.370.97-ft?-fl.$15-MMT5U-O.-3EE0l7B0.CI3-0.flEB-O.M-fl.OTflD.SEZOiTI0.41T呵0d.31.9.oaao-Bia-目t.別心(lm阴ol；i嗣yl7ua.2020.爼

6、m-d_M7j|然后复制矩阵赋值给变量a,即输入命令a=复制内容，接着输入b=a，就获得符合要求的N*3转置矩阵b了，见下页图4。但是这样的矩阵直接复制到DtiStudio的梯度表框内是会运行出错的，因为DtiStudio还要求了梯度表输入的格式。三、调整梯度表矩阵格式，使其满足输入要求。DtiStudio的梯度表输入要求是矩阵每行三个变量前要有扩散方向的序号，而且变量之间要有逗号相隔。我就按照这个要求修改矩阵。首先把转置后的梯度表写入Excel，这样就方便编辑了。在MATLAB中，矩阵数据写入Excel的命令如下:xlswrite(F:DTIGradientTable.xls,b,Sheet

7、1,A1:C26)(这里要注意修改文件路径；C26代表有多少个扩散方向，也要注意修改数值。另外，使用这个命令是有必要的，直接复制矩阵到Excel的话，更加不好编辑了，大家可以试试。)然后，就是补入序号，冒号和逗号了。如下图5:-1.0000-0.Sd9O0.1080-0.88400.00300.S6SQ-0.79900.1620-0.S66O0.2120-0.0080-0.55000.43500.69900.52500.6000-0.6530-0.20700.l3O-1.0000-0.Sd9O0.1080-0.88400.00300.S6SQ-0.79900.1620-0.S66O0.2120

8、-0.0080-0.55000.43500.69900.52500.6000-0.6530-0.20700.l3O0.4360-0.46200.SOSO-0.S24O-0.29700.0400O.J0.5fl50-0.3450-0.7360-0.23SO0.37000.9870-0.1290-0.9360-0.S92O-0.5440-0.42200.78000.0300-0.S88O-0.0600-0.0760-0.69900.S22O0.87400.4880-0.53000.34900.3180A10213243546576879B10911LD12111312141315141615171

9、6.817L9LS20L921202221232224232324262b-0.849-0.E660.8180CLR680.162一山S5L00-48B0.3490.31B-0.31500.4J600.47500.23LO一Aid600.1820-0.5310-0.3060ci.sago0.IOS7699a822o,a?4这样就符合要求了，可是复制粘贴过去DtiStudio还是不行啊，会运行出错。我反复试验，原来数值和逗号之间是不能有空格的。这样的话就要把全部空格去掉，空格也挺多的，一个个删太麻烦，这里用WORD的批处理方法就容易了。接下来复制Excel中的梯度表到WORD里（粘贴时选择“只保留文本”），然后调入“查找与替换”选