diff --git a/img2txt_wholebrain&hippo.m b/img2txt_wholebrain&hippo.m new file mode 100644 index 0000000..883b749 --- /dev/null +++ b/img2txt_wholebrain&hippo.m @@ -0,0 +1,54 @@ +%Name:Img2txt +%Function:将img文件读取并转换为一列扫描格式的txt文件(用于unity输入) +%Author:Changle Zhang, zhangchanglehit@163.com +clc;clear all;close all; +[hdr,filetype,fileprefix,machine] = load_nii_hdr('p01_2808584_286L.hdr'); +[img,hdr] = load_nii_img(hdr,filetype,fileprefix,machine); +%pic=double(img(:,:,60)); +%imshow(pic,[0 255]); +% savemat=zeros(1,216*256); +% counter=1; +%抛出XXXXX区域之外的结果 +% for i=1:181 +% for j=1:217 +% for k=1:181 +% if(img(i,j,k)==0||img(i,j,k)==169||img(i,j,k)==170||(img(i,j,k)>=254&&img(i,j,k)<=260)) +% img(i,j,k)=0; +% else +% img(i,j,k)=256; +% end +% end +% end +% end +%figure; +%imshow(img,[0 255]); +%选择海马区的结果 +for i=1:181 + for j=1:217 + for k=1:181 + if(img(i,j,k)~=75&&img(i,j,k)~=76&&img(i,j,k)~=284&&img(i,j,k)~=282) + img(i,j,k)=0; + else + img(i,j,k)=256; + end + end + end +end + + +matsize=size(img); +numl=matsize(1)*matsize(2)*matsize(3);%length of txt; +savemat=zeros(numl,4); +line=1; +for x=1:matsize(1) + for y=1:matsize(2) + for z=1:matsize(3) + savemat(line,1)=x; + savemat(line,2)=y; + savemat(line,3)=z; + savemat(line,4)=img(x,y,z); + line=line+1; + end + end +end +save hippo_img.txt -ascii savemat; \ No newline at end of file