clc
clear
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%讀取圖像，預處理%%%%%%%%%%%%%%%
he?=?imread（‘C:\Users\lpflying\Desktop\my?doc\程序代碼\BSR\BSDS500\data\images\train\train5.jpg‘）;%讀取圖像
cform?=?makecform（‘srgb2lab‘）;%圖像由RGB轉為lab
lab_he?=?applycform（hecform）;%
lab_he=double（lab_he）;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%設置初值%%%%%%%%%%%%%%%%%%
color=[25000];
thre=0.02;%最終生成分割圖像梯度閾值
m=40;%權值
k=1;%劃分為300個簇
die=20;%kmeans迭代die次
x=size（he）;s=（x（1）*x（2）/k）^0.5;s=ceil（s）;%初始分割網(wǎng)格間距s
r=ceil（x（1）/s）;%網(wǎng)格行數(shù)r
w=ceil（x（2）/s）;%網(wǎng)格列數(shù)w
ct=r*w;
belong=ones（x（1）x（2））;
center=zeros（ct5）;
%初始每個像素點的距離
dist=9999*ones（x（1）x（2））;
%初始中心節(jié)點center
for?i=1:r
????for?j=1:w
????????if?（i????????????x1=（i-1）*s+fix（s/2）;
????????else
????????????x1=（i-1）*s+fix（rem（x（1）s）/2）;
????????end
????????if?（j????????????y1=（j-1）*s+fix（s/2）;
????????else
????????????y1=（j-1）*s+fix（rem（x（2）s）/2）;
????????end
????????z=lab_he（x1y1:）;
????????center（（i-1）*w+j:）=[z（::1）?z（::2）?z（::3）?x1?y1];%初始中心節(jié)點center
????end
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%迭代聚類處理%%%%%%%%%%%%%%%
t1=clock;
move=99999;
for?c=1:die????%進行迭代die次
????if?move<10
????????break;
????end
????move=0;
????c1=zeros（ct）;
????ct_x=zeros（ct）;
????ct_y=zeros（ct）;
????ct_l=zeros（ct）;
????ct_a=zeros（ct）;
????ct_b=zeros（ct）;
????for?i=1:ct
????????%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
????????for?u=center（i4）-s:center（i4）+s
????????????if（u>=1）&&（u<=x（1））
????????????????for?v=center（i5）-s:center（i5）+s
????????????????????if（v>=1）&&（v<=x（2））
????????????????????????dc=（（lab_he（uv1）-center（i1））^2+（lab_he（uv2）-center（i2））^2+（lab_he（uv3）-center（i3））^2）^0.5;
????????????????????????ds=（（u-center（i4））^2+（v-center（i5））^2）^0.5;
????????????????????????d=（（dc）^2+（ds*m/s）^2）^0.5;%計算距離
????????????????????????if?d????????????????????????????dist（uv）=d;
????????????????????????????belong（uv）=i;?
????????????????????????????move=move+1;
????????????????????????end
????????????????????end
????????????????end
????????????end
????????end
????end?
????for?k=1:x（1）
????????for?g=1:x（2）
????????????i=belong（kg）;
????????????c1（i）=c1（i）+1;
????????????ct_x（i）=ct_x（i）+k;
????????????ct_y（i）=ct_y（i）+g;
????????????ct_l（i）=ct_l（i）+lab_he（kg1）;
????????????ct_a（i）=ct_a（i）+lab_he（kg2）;
????????????ct_b（i）=ct_b（i）+lab_he（kg3）;
????????????
????????end
????end
????for?i=1:ct
????????center（i4）=fix（ct_x（i）/c1（i））;
????????center（i5）=fix（ct_y（i）/c1（i））;
????????center（i1）=fix（ct_l（i）/c1（i））;
????????center（i2）=fix（ct_a（i）/c1（i））;
????????center（i3）=fix（ct_b（i）/c1（i））;
????end
end
t2=clock;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%去壞點%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for?i=1:ct
????bw=zeros（x（1）x（2））;
????for?k=1:x（1）
????????for?g=1:x（2）
????????????if?belong（kg）==i
????????????????bw（kg）=1;
?????????

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????6373??2013-10-24?12:34??slic\slic.m
?????文件??????????68??2019-05-05?17:51??slic\說明.txt
?????目錄???????????0??2019-05-05?17:51??slic\