count the ROI volume information

ROI_volmue.m

@@ -0,0 +1,40 @@
+%function:统计一个目录下所有分割结果中286个脑区的ROI大小，放到一个txt文档中
+clc;clear all;close all;
+pathname_G1 = 'G:\韩璎课题组数据\NC\NC_分割结果';
+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),286);
+ 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,'286Labels_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:286
+               counter = length(find(img==k));
+               savemat(i,k)= counter;
+            end
+        end
+
+    end
+
+
+    cd(pathname_G1);
+ end
+ save wholeBrain_voxel.txt -ascii savemat;