Brain MRI Ssegmentation

Brain MRI Ssegmentation自己写的分割脑组织的程序，用FCM方法分割，可以直接运行

clc
clear
img= imread （‘temp.bmp‘）; 
[mn]=size（img）;
subplot（221）imshow（img）;title（‘ 图一 原图像‘）;
hold off;
img=double（img）;
 
for i=1:1000
%归类
     sum1=0;count1=0;
     sum2=0;count2=0;
     sum3=0;count3=0;
     sum4=0;count4=0;
  tic
 for x=1:217
    for y=1:181        
     c1（1）=0;  %背景 
     c2（1）=25; %脑脊液
     c3（1）=83; %脑灰质
     c4（1）=115;%脑白质 
     r1=abs（img（xy）-c1（i））;
     r2=abs（img（xy）-c2（i））;
     r3=abs（img（xy）-c3（i））;
     r4=abs（img（xy）-c4（i））;
     d=[r1r2r3r4];
     min1=min（d）;
     if min1==r1
         sum1=sum1+img（xy）;
         count1=count1+1;
     elseif min1==r2
         sum2=sum2+img（xy）;
         count2=count2+1;
     elseif min1==r3
         sum3=sum3+img（xy）;
         count3=count3+1;
     elseif min1==r4
         sum4=sum4+img（xy）;
         count4=count4+1;
     end
    end
 end
 i=i+1;
 elapsedtime = toc
 % 计算新的聚类中心
 c1（i）=sum1/count1;
 c2（i）=sum2/count2;
 c3（i）=sum3/count3; 
 c4（i）=sum4/count4;
 d1（i）=abs（c1（i）-c1（i-1））;
 d2（i）=abs（c2（i）-c2（i-1））;
 d3（i）=abs（c3（i）-c3（i-1））;
 d4（i）=abs（c4（i）-c4（i-1））;
 if d1（i）<=0.001&&d2（i）<=0.001&&d3（i）<=0.001
         R=c1（i）;
         G=c2（i）;
         B=c3（i）;
         A=c4（i）;
         break;
 end
end  
R
G
B
A
i

%归类显示
 for x=1:217
    for y=1:181        
     r1=abs（img（xy）-R）;
     r2=abs（img（xy）-G）;
     r3=abs（img（xy）-B）;
     r4=abs（img（xy）-A）;
     d=[r1r2r3r4];
     min1=min（d）;
     if min1==r1
         img（xy）=0;       %背景
     elseif min1==r2
         img（xy）=0;       %脑脊液
     elseif min1==r3
          img（xy）=0;      %脑灰质
     elseif min1==r4
          img（xy）=110;    %脑白质
     end
    end
 end
subplot（222）imshow（img）;title（‘ 图二 聚类后脑白质的图像‘） 
img= imread （‘temp.bmp‘）;
 for x=1:217
    for y=1:181        
     r1=abs（img（xy）-R）;
     r2=abs（img（xy）-G）;
     r3=abs（img（xy）-B）;
     r4=abs（img（xy）-A）;
     d=[r1r2r3r4];
     min1=min（d）;
     if min1==r1
         img（xy）=0;         
     elseif min1==r2
         img（xy）=0;      
     elseif min1==r3
          img（xy）=110;      
     elseif min1==r4
          img（xy）=0;      
     end
    end
 end
subplot（223）imshow（img）;title（‘ 图三 聚类后脑灰质的图像‘） 
img= imread （‘temp.bmp‘）;
 for x=1:217
    for y=1:181        
     r1=abs（img（xy）-R）;
     r2=abs（img（xy）-G）;
     r3=abs（img（xy）-B）;
     r4=abs（img（xy）-A）;
     d=[r1r2r3r4];
     min1=min（d）;
     if min1==r1
         img（xy）=0;            
     elseif min1==r2
         img（xy）=110;      
     elseif min1==r3
          img（xy）=0;      
     elseif min1==r4
          img（xy）=0;      
     end
    end
 end
subplot（224）imshow（img）;title（‘ 图四 聚类后脑脊液的图像‘）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       2831  2010-12-02 20:40  Brain MRI Ssegmentation\Kmeans.m

     文件        526  2010-11-02 13:43  Brain MRI Ssegmentation\mask.m

     文件      41006  2010-09-28 14:35  Brain MRI Ssegmentation\MRI.bmp

     文件       1253  2010-12-02 15:19  Brain MRI Ssegmentation\MRI_segmentation.m

     文件        312  2010-09-26 22:04  Brain MRI Ssegmentation\regiongrow.m

     文件      41006  2010-12-02 15:18  Brain MRI Ssegmentation\temp.bmp

     文件       5567  2010-10-27 19:43  Brain MRI Ssegmentation\html\MRI_segmentation.html

     文件       3531  2010-10-27 19:43  Brain MRI Ssegmentation\html\MRI_segmentation.png

     文件      30763  2010-10-27 19:43  Brain MRI Ssegmentation\html\MRI_segmentation_01.png

     文件       1641  2010-10-27 19:43  Brain MRI Ssegmentation\html\MRI_segmentation_02.png

     文件      26531  2010-10-27 19:43  Brain MRI Ssegmentation\html\MRI_segmentation_03.png

     文件       3189  2010-10-27 19:43  Brain MRI Ssegmentation\html\MRI_segmentation_04.png

     文件       3353  2010-10-27 19:43  Brain MRI Ssegmentation\html\MRI_segmentation_05.png

     文件       2331  2010-10-27 19:43  Brain MRI Ssegmentation\html\MRI_segmentation_06.png

     文件       4048  2010-10-27 19:43  Brain MRI Ssegmentation\html\MRI_segmentation_07.png

    ..A.SH.     30208  2011-02-24 14:34  Brain MRI Ssegmentation\html\Thumbs.db

     目录          0  2011-10-27 21:59  Brain MRI Ssegmentation\html

     目录          0  2011-10-27 21:59  Brain MRI Ssegmentation

----------- ---------  ---------- -----  ----

               198096                    18