注意:非常抱歉的和大家说我已经不再从事相关学科的研究工作,所以没办法给大家解释此类问题了。欢迎大家移步52脑科学网站等地址寻求帮助,感谢大家一直以来的支持与理解。
近期因为要绘制人脑连接结构的网络连接图,所以接触到了北师大研发的BrainNet Viewer软件,download地址为:http://www.nitrc.org/projects/bnv/
该软件为我们提供了多种图形组合方式,具体的功能可以参见manual,这里我只谈一下我自己的使用感受。
首先是软件启动,比较简单,都是基于matlab开发的,只需要添加路径后在command window中输入'BrainNet'即可。
打开软件后,可以根据我们的选择加载需要的图像,这里我选择全部加载的图像。
1.surface file:这个会自动跳转到系统为大家提供的12个现有surface数据下,我们也可以用同样的方法来绘制自己的surface(感觉意义不大)。
2.Data file(node):这里是输入节点定义信息,我选择使用提供的example中的AAL90节点定义模板,这个模板很棒,直接为我们提供了Lookup table。
3.Data file(edge):这里是输入连接数据信息,也就是edge信息,我选择使用我自己通过PANDA软件生成的确定性determinstric结构连接信息中的FA,FN,Length矩阵进行构建。
4.Volume file:这里输入的是ROI的信息,我们可以将ROI通过他提供的算法映射到surface上,也可以进行ROI绘制。(投射时出现116的问题是因为软件提供的example中的aal.nii是116个分区,也就是AAL116的ROI信息,并非AAL90的信息)。
收入完成后点击OK继续,进行不同的设置,具体设置细节可以参见manual。
因为我对扣带回区域感兴趣,所以我打算研究扣带回结构和其他结构之间的连接关系,这里我就选择提取扣带回(对应于ALL中的Label 31~36)。
这里,我们就需要使用我为大家提供的toolkit来进行ROI连接矩阵的提取,具体代码如下:
%Name:CIN_network_extraction_toolkit
%Function:Extraction CIN(lable 31:36 in AAL 90 atlas set)to show the
%construction
%Author :Changle Zhang e-mail:zhangchanglehit@163.com
clc;clear all;close all;
cd ('E://Anhuishi/00001/Network/Deterministic'); %进入存储目录,后期可以根据需要修改0000i以及确定不确定等
%FA=load ('00001_dti_FACT_35_02_1_0_Matrix_FA_AAL_Contract_90_2MM_90.txt'); %打开FA
%FA=load ('00001_dti_FACT_35_02_1_0_Matrix_FN_AAL_Contract_90_2MM_90.txt');
%%打开FN
FA=load('00001_dti_FACT_35_02_1_0_Matrix_Length_AAL_Contract_90_2MM_90.txt');
%%打开Length
%Upper_Left
for i=1:30 for j=1:30 FA(i,j)=0; end
end
%Upper_Right
for i = 1:30 for j = 37:90 FA(i,j)=0; end
end
%Lower_Left
for i=37:90 for j=1:30 FA(i,j)=0; end
end
%Lower_Right
for i=37:90 for j=37:90 FA(i,j)=0; end
end
%save('FA_AAL.txt','FA','-ascii'); %进行保存
%save('FN_AAL.txt','FA','-ascii'); %进行保存
save('Length_AAL.txt','FA','-ascii'); %进行保存
通过这个toolkit可以很好地实现ROI分割,我计划将自己用到的所有实用toolkit做成GUI的形式提供给BrainNet Viewer的软件制作团队,这样可以方便后续的研究人员。希望大家可以喜欢这个toolkit。
这里附上我的处理结果,明天去和老师汇报~
文章来源: zclhit.blog.csdn.net,作者:zclhit_,版权归原作者所有,如需转载,请联系作者。
原文链接:zclhit.blog.csdn.net/article/details/51088181