%?讀取圖像
I?=?imread（‘F:\Camera\111.jpg‘）;
[mnt]?=?size（I）;
%?四鄰域模糊圖像
%?imshow（I）;
I1（:1:n-1:）?=?I（:2:n:）;
I1（:n:）?=?I（:1:）;
I2（:1:）?=?I（:n:）;
I2（:2:n:）?=?I（:1:n-1:）;
I3（1:m-1::）?=?I（2:m::）;
I3（m::）?=?I（1::）;
I4（1::）?=?I（m::）;
I4（2:m::）?=?I（1:m-1::）;
I?=?double（I）;
I1?=?double（I1）;
I2?=?double（I2）;
I3?=?double（I3）;
I4?=?double（I4）;
Ibr?=?（I1+I2+I3+I4+I）?/?5;
%?Ibr=uint8（Ibr）;
%?figureimshow（Ibr）;

%?確定8個顏色空間
R?=?[25500];
G?=?[02550];
B?=?[00255];
RG?=?R?+?G;
RB?=?R?+?B;
GB?=?G?+?B;
White?=?R?+?G?+B;
Black?=?White?*?0;
Color?=?[R;G;B;RG;RB;GB;White;Black];
%?image?=?zeros（40203）;
%?image（11:）?=?[R（1）R（2）R（3）];
%?image（12:）?=?[G（1）G（2）G（3）];
%?image（21:）?=?[B（1）B（2）B（3）];
%?image（22:）?=?[RG（1）RG（2）RG（3）];
%?image（31:）?=?[RB（1）RB（2）RB（3）];
%?image（32:）?=?[GB（1）GB（2）GB（3）];
%?image（41:）?=?[White（1）White（2）White（3）];
%?image（42:）?=?[Black（1）Black（2）Black（3）];
%?image?=?uint8（image）;
%?imshow（image）;

%?講圖像的顏色分布到顏色空間中
for?p=1:m
????for?q=1:n
????????Ipq?=?[Ibr（pq1）Ibr（pq2）Ibr（pq3）];
????????dist1?=?norm（R-Ipq）;
????????dist2?=?norm（G-Ipq）;
????????dist3?=?norm（B-Ipq）;
????????dist4?=?norm（RG-Ipq）;
????????dist5?=?norm（RB-Ipq）;
????????dist6?=?norm（GB-Ipq）;
????????dist7?=?norm（White-Ipq）;
????????dist8?=?norm（Black-Ipq）;
????????Dist?=?[dist1dist2dist3dist4dist5dist6dist7dist8];
????????dist?=?min（Dist）;
????????[m1n1]?=?find（Dist?==?dist）;
????????Ic（pq:）?=?Color（n1:）;
????end
end
%?Ic?=?uint8（Ic）;
%?imshow（Ic）;

%?計算連通區域
%?注：代碼寫的還是有一些問題，如果各個顏色區間是互相隔離的效果最好。
Mask?=?zeros（mn）;

k?=?1;
for?p=1:m
????for?q=1:n
????????if?Mask（pq）?==?0
????????????Mask（pq）?=?k;
????????????%?CColor（k:）?=?Ic（pq:）;
????????????k?=?k?+?1;
????????end
????????%?1
????????if?p-1>0?&&?p-10?&&?q-1????????????if?Ic（p-1q-1:）?==?Ic（pq:）
????????????????if?Mask（p-1q-1）?==?0
????????????????????Mask（p-1q-1）?=?Mask（pq）;
????????????????else
????????????????????Mask（pq）?=?min（Mask（p-1q-1）Mask（pq））;
????????????????end
????????????end
????????end
????????%?2
????????if?p-1>0?&&?p-10?&&?q????????????if?Ic（p-1q:）?==?Ic（pq:）
????????????????if?Mask（p-1q）?==?0
????????????????????Mask（p-1q）?=?Mask（pq）;
????????????????else
????????????????????Mask（pq）?=?min（Mask（p-1q）Mask（pq））;
????????????????end
????????????end
????????end
????????%?3
????????if?p-1>0?&&?p-10?&&?q+1????????????if?Ic（p-1q+1:）?==?Ic（pq:）
????????????????if?Mask（p-1q+1）?==?0
????????????????????Mask（p-1q+1）?=?Mask（pq）;
????????????????else
????????????????????Mask（pq）?=?min（Mask（p-1q+1）Mask（pq））;
????????????????end
????????????end
????????end
????????%?4
????????if?p>0?&&?p0?&&?q-1????????????if?Ic（pq-1:）?==?Ic（pq:）
????????????????if?Mask（pq-1）?==?0
????????????????????Mask（pq-1）?=?Mask（pq）;
????????????????else
?????????????????