matlab练习程序(二值图像连通区域标记法,一步法)

  这个只需要遍历一次图像就能够完全标记了。我主要参考了WIKI这位兄弟的博客,这两个把原理基本上该介绍的都介绍过了,我也不多说什么了。一步法代码相比两步法真是清晰又好看,似乎真的比两步法要好很多。

  代码如下:

clear all;

close all;

clc;



img=imread('liantong.bmp');

imgn=img>128;

s=uint8(1-imgn);



 %{                     

s=[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;

   0 0 1 1 0 0 1 1 0 0 1 1 0 0 1 1 0;       %这个矩阵是维基百科中的矩阵

   0 1 1 1 1 1 1 1 1 0 0 1 1 1 1 0 0;

   0 0 0 1 1 1 1 0 0 0 1 1 1 1 0 0 0;

   0 0 1 1 1 1 0 0 0 1 1 1 0 0 1 1 0;

   0 1 1 1 0 0 1 1 0 0 0 1 1 1 0 0 0;

   0 0 1 1 0 0 0 0 0 1 1 0 0 0 1 1 0;

   0 0 0 0 0 0 1 1 1 1 0 0 1 1 1 1 0;

   0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0];

%}

imshow(mat2gray(s));

[m n]=size(s);

tmp=zeros(m,n);     %标记图像

label=1;

queue_head=1;       %队列头

queue_tail=1;       %队列尾

neighbour=[-1 -1;-1 0;-1 1;0 -1;0 1;1 -1;1 0;1 1];  %和当前像素坐标相加得到八个邻域坐标



for i=2:m-1

    for j=2:n-1

        

        if s(i,j)==1 && tmp(i,j) ==0           

            tmp(i,j)=label;

            q{queue_tail}=[i j];        %用元组模拟队列,当前坐标入列

            queue_tail=queue_tail+1;

            

            while queue_head~=queue_tail

                pix=q{queue_head};                

                for k=1:8               %8邻域搜索

                    pix1=pix+neighbour(k,:);
            if pix1(1)>=2 && pix1(1)<=m-1 && pix1(2) >=2 &&pix1(2)<=n-1   if s(pix1(1),pix1(2)) == 1 && tmp(pix1(1),pix1(2)) ==0 %如果当前像素邻域像素为1并且标记图像的这个邻域像素没有被标记,那么标记   tmp(pix1(1),pix1(2))=label;   q{queue_tail}=[pix1(1) pix1(2)];   queue_tail=queue_tail+1;   end
            end end queue_head=queue_head+1; end clear q; %清空队列,为新的标记做准备 label=label+1; queue_head=1; queue_tail=1; end end end figure,imshow(mat2gray(tmp))

下面是效果图,就效果而言和上一篇没什么区别的。

matlab练习程序(二值图像连通区域标记法,一步法)原图

matlab练习程序(二值图像连通区域标记法,一步法)效果图

这两篇算是把二值图像连通标记给搞定了。

你可能感兴趣的:(matlab)