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