%%%%%%%%%%%?基于聚類法的醫學圖像分割實驗?%%%%%%%%%%%%%%%%
%?10.16.2015?AM-08:22
clcclear?allclose?all;
%?tic
%?Elapsed?time?is?0.520028?seconds.
%%?調入原始圖像
data=imread（‘D:\123\888.jpg‘?）;?
data=double（data）;
%?第1幅圖像
figuresubplot（221）imshow（uint8（data））title（‘原圖像‘）;
%?subplot（231）注解：2*3矩陣中的第1幅圖像-?從左到右，由上至下
%?第2幅圖像
subplot（222）imhist（uint8（data））title（‘原圖像直方圖‘）;
%%?由隨機數產生初始聚類中心
[M?N]=size（data）;
[countx]=imhist（uint8（data））;
mima=find（count~=0）;???%灰度等級?256*1的矩陣
V1=zeros（2）;????%產生2個分組
V2=zeros（2）;
U1=zeros（length（mima）2）;???????%產生一個n*c的矩陣
V1（1）=mima（round（rand（1）*（length（mima）/4））+1）;????%產生隨機數??產生初始聚類中心
V1（2）=mima（round（length（mima）/2+rand（1）*（length（mima）/4））+1）;?
%%?黑白聚類中心的尋找與計算并進行圖像分割
k=0;
while?k==0
?for?j=1:length（mima）
???if?abs（mima（j）-V1（2））????????U1（j1）=0.0;
????????U1（j2）=1.0;
???else
????????U1（j1）=1.0;
????????U1（j2）=0.0;
???end
?end???
???tp1=0.0;
???tp2=0.0;
???tp3=0.0;
???tp4=0.0;
???for?j=1:length（mima）
??????tp1=tp1+count（mima（j））*（mima（j）-1）*U1（j1）;
??????tp2=tp2+count（mima（j））*U1（j1）;
??????tp3=tp3+count（mima（j））*（mima（j）-1）*U1（j2）;
??????tp4=tp4+count（mima（j））*U1（j2）;
????end
????V2（1）=tp1/tp2;
????V2（2）=tp3/tp4;
????if??（abs（V2（1）-V1（1））+abs（V2（2）-V1（2）））>0.0001
???????????????V1（1）=V2（1）;
???????????????V1（2）=V2（2）;
????else
????????k=1;
????end
?end
?subplot（223）
?img1=250.*ones（size（data））;
?for?i=1:M
????for?j=1:N
????????tp=data（ij）;
????????for?l=1:length（mima）
????????????if?tp==mima（l）
????????????????k=l;
????????????end
????????end
????????if?U1（k1）>=U1（k2）
????????????img1（ij）=0;
????????end
????end
?end
%?第3幅圖像
subplot（223）imshow（uint8（img1））title（‘聚類分割圖像結果‘）;
%?第4幅圖像
subplot（224）imhist（uint8（img1））title（‘2D直方圖‘）

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????1933??2015-10-16?09:46??Juleishuiji.m