IO=imread（‘E:\matlab\program\e6\m1.bmp‘）;
IR=rgb2gray（IO）;
A=im2bw（IR）;%二值化
[lw]=size（A）;
B=zeros（l+2w+2）;
B（2:l+12:w+1）=A;%外圍補零
[mn]=size（B）;
LabelMap=zeros（mn）;%全零矩陣
LabelIndex=0;%標(biāo)簽號
N=0;%已貼標(biāo)簽數(shù)，四連通算法
for?i?=?2?:?m-1
????for?j?=?2?:?n-1???????
??????????if?（LabelMap（ij-1）==0?&&?LabelMap（i-1j）==0?&&?B（ij）==1）%上和左像素值為0
????????????LabelIndex=LabelIndex+1;%標(biāo)簽號+1
????????????LabelMap（ij）=LabelIndex;%在LabelMap矩陣中貼標(biāo)簽
????????????N=N+1;%已貼標(biāo)簽數(shù)+1
??????????end
????????if?（LabelMap（ij-1）~=0?&&?LabelMap（i-1j）==0?&&?B（ij）==1）%已掃描像素標(biāo)簽，上或左不為零
????????????????LabelMap（ij）=LabelMap（ij-1）;%和相鄰的像素點標(biāo)簽相同
????????end
?????????????if（?LabelMap（ij-1）==0?&&?LabelMap（i-1j）~=0?&&?B（ij）==1）
????????????????????LabelMap（ij）=LabelMap（i-1j）;
?????????????end???%已掃描像素標(biāo)簽，上或左同時不為