clear?all;?
??
s=imread（‘L1.bmp‘）;
?
s=double（s）;
?
?
X1=20;
X2=45;
X3=90;?
?
num1=0;
num2=0;
num3=0;
[MN]=size（s）;
?
R=zeros（MN）;????
?
a=10;
b=20;
counter=0;??????
format?long;
?
?
while（abs（b-a）>0.001）???
????counter=counter+1?
????a
????b=a;?
????a=0;?
????for?i=1:M?
????for?j=1:N?
????????????
????%-------------------------------------------------------
????????????p1=（s（ij）-X1）.^2;?????
????????????p2=（s（ij）-X2）.^2;?
????????????p3=（s（ij）-X3）.^2;???????
????????????min_p=min（[p1p2p3]）;????????
????????????if?（p1==min_p）??
????????????????R（ij）=1;????????
????????????????num1=num1+1;???
????????????elseif?（p2==min_p）?
????????????????R（ij）=2;???????
????????????????num2=num2+1;?
????????????elseif?（p3==min_p）?
????????????????R（ij）=3;????????
????????????????num3=num3+1;?
????????????end????????
????????a=a+sqrt（min_p）;?????????????????????
????????end?
????end?
????
????
???X1=X1*0;
???X2=X2*0;
???X3=X3*0;?????????????????????????????
????a=a/M/N;?
????
????%---------------------------------------------------
????for?i=1:M?
????for?j=1:N?
????????????if?（R（ij）==1）?X1=X1+s（ij）;?
????????????elseif?（R（ij）==2）?X2=X2+s（ij）;?
????????????elseif?（R（ij）==3）?X3=X3+s（ij）;
????????????end?
????????end?
????end?
????
????X1=X1/num1;?num1=0;?????
????X2=X2/num2;?num2=0;?
????X3=X3/num3;?num3=0;?????
end?
?
?
r=zeros（MN）;?
for??i=1:M?
for?j=1:N????????????????
????????if?（R（ij）==1）?????r（ij）=30;???????
????????elseif（R（ij）==2）??r（ij）=100;??????
????????elseif（R（ij）==3）??r（ij）=200;?
????????end?
????end?
end??
r=uint8（r）;?
figure;?
imshow（r）;???????????????????????????
imwrite（r‘r1.bmp‘）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????1542558??2014-01-16?21:36??圖像分析期末設計\2012_圖像分析作業_劉蕾_SC12023170?.pdf

?????文件???????1718??2014-01-17?21:11??圖像分析期末設計\chengxu1.m

?????文件??????66616??2007-03-05?08:43??圖像分析期末設計\L1.bmp

?????文件??????66616??2007-03-05?08:43??圖像分析期末設計\L2.bmp

?????文件??????66616??2007-03-05?08:43??圖像分析期末設計\L3.bmp

?????文件??????66616??2007-03-05?08:43??圖像分析期末設計\L4.bmp

?????文件??????66616??2007-03-05?08:43??圖像分析期末設計\L5.bmp

?????文件??????66616??2007-03-05?08:43??圖像分析期末設計\L6.bmp

?????文件??????66616??2007-03-05?08:43??圖像分析期末設計\L7.bmp

?????文件?????196662??2014-01-17?21:07??圖像分析期末設計\r.bmp

?????文件??????66614??2014-01-17?20:57??圖像分析期末設計\r1.bmp

?????文件??????66614??2014-01-15?20:39??圖像分析期末設計\r2.bmp

?????文件??????66614??2014-01-15?20:40??圖像分析期末設計\r3.bmp

?????文件??????66614??2014-01-15?20:41??圖像分析期末設計\r4.bmp

?????文件??????66614??2014-01-15?20:44??圖像分析期末設計\r5.bmp

?????文件??????66614??2014-01-15?20:44??圖像分析期末設計\r6.bmp

?????文件??????66614??2014-01-15?20:45??圖像分析期末設計\r7.bmp

?????文件???????2689??2014-01-17?21:08??圖像分析期末設計\zongchengxu.m

????..A.SH.???????162??2014-01-16?09:43??圖像分析期末設計\~$鑫SC13023015期末實驗報告.doc

?????文件????2582528??2013-11-19?10:00??圖像分析期末設計\圖像作業\DIA課程實驗2013-2014.ppt

?????文件??????24576??2014-01-15?15:02??圖像分析期末設計\圖像作業\K-means圖像分割.doc

?????文件???????2042??2014-01-15?15:52??圖像分析期末設計\圖像作業\k_means.m

?????文件??????19065??2014-01-15?16:03??圖像分析期末設計\圖像作業\K_means_Image.cpp

?????文件??????65651??2014-01-15?14:45??圖像分析期末設計\圖像作業\K_means算法在遙感分類中的應用.caj

?????文件??????66616??2007-03-05?08:43??圖像分析期末設計\圖像作業\L1.bmp

?????文件??????66616??2007-03-05?08:43??圖像分析期末設計\圖像作業\L2.bmp

?????文件??????66616??2007-03-05?08:43??圖像分析期末設計\圖像作業\L3.bmp

?????文件??????66616??2007-03-05?08:43??圖像分析期末設計\圖像作業\L4.bmp

?????文件??????66616??2007-03-05?08:43??圖像分析期末設計\圖像作業\L5.bmp

?????文件??????66616??2007-03-05?08:43??圖像分析期末設計\圖像作業\L6.bmp

............此處省略46個文件信息