自编程序实现matlab对图像矩阵贴标签求连通域

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;%标签号
N=0;%已贴标签数，四连通算法
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;%标签号+1
            LabelMap（ij）=LabelIndex;%在LabelMap矩阵中贴标签
            N=N+1;%已贴标签数+1
          end
        if （LabelMap（ij-1）~=0 && LabelMap（i-1j）==0 && B（ij）==1）%已扫描像素标签，上或左不为零
                LabelMap（ij）=LabelMap（ij-1）;%和相邻的像素点标签相同
        end
             if（ LabelMap（ij-1）==0 && LabelMap（i-1j）~=0 && B（ij）==1）
                    LabelMap（ij）=LabelMap（i-1j）;
             end   %已扫描像素标签，上或左同时不为