diff --git a/ROI_volmue_283_MNI.m b/ROI_volmue_283_MNI.m new file mode 100644 index 0000000..f4fa173 --- /dev/null +++ b/ROI_volmue_283_MNI.m @@ -0,0 +1,40 @@ +%function:统计一个目录下所有分割结果中283个脑区的ROI大小,放到一个txt文档中 +clc;clear all;close all; +pathname_G1 = 'J:\XuanWu\M2_PD分割'; +cd(pathname_G1); +dirs_G1=dir([pathname_G1, '\*.*']); +dircell_G1=struct2cell(dirs_G1)' ; +filenames_G1=dircell_G1(3:length(dircell_G1),1); +%file_path_G1=fullfile(pathname_G1,filenames_G1); %保存各网络参数的全路径 +dirs_G2=dir([pathname_G1, '\*.*']); +dircell_G2=struct2cell(dirs_G2)' ; +filenames_G2=dircell_G2(3:length(dircell_G2),1); +savemat = zeros(length(filenames_G1),283); + for i=1:length(filenames_G1) + cd(char(filenames_G1(i))); + cd output; + clear dirs_G2;clear dircell_G2;clear filenames_G2; + dirs_G2=dir; + dircell_G2=struct2cell(dirs_G2)' ; + filenames_G2=dircell_G2(3:length(dircell_G2),1); + len = length(filenames_G2); + for j = 1:len + + filenamesIS = char(filenames_G2(j)); + jamar = strfind(filenamesIS,'283Labels_M2_MNI.img'); + if length(jamar) ~= 0 + + [hdr,filetype,fileprefix,machine] = load_nii_hdr(filenamesIS); + [img,hdr] = load_nii_img(hdr,filetype,fileprefix,machine); + for k = 1:283 + counter = length(find(img==k)); + savemat(i,k)= counter; + end + end + + end + + + cd(pathname_G1); + end + save wholeBrain_voxel.txt -ascii savemat; diff --git a/ROI_volmue.m b/ROI_volmue_286_MNI.m similarity index 96% rename from ROI_volmue.m rename to ROI_volmue_286_MNI.m index ea5906f..17bb2d8 100644 --- a/ROI_volmue.m +++ b/ROI_volmue_286_MNI.m @@ -1,6 +1,6 @@ %function:统计一个目录下所有分割结果中286个脑区的ROI大小,放到一个txt文档中 clc;clear all;close all; -pathname_G1 = 'G:\韩璎课题组数据\NC\NC_分割结果'; +pathname_G1 = 'G:\韩璎课题组数据\aMCI\aMCI_M2'; cd(pathname_G1); dirs_G1=dir([pathname_G1, '\*.*']); dircell_G1=struct2cell(dirs_G1)' ;