Demo entry 6768083

userjzb

   

Submitted by anonymous on Nov 01, 2018 at 08:04
Language: Matlab. Code size: 2.5 kB.

clear all;
close all;
clc;
%matlab练习程序(二值图像连通区域标记法,一步法)
img=imread('rice.png');
figure,imshow(img);
%imgn=img>128;
%用ostu方法获取二值化阈值,进行二值化并进行显示
level=graythresh(img);
s=im2bw(img,level);
figure,imshow(s);
%s=uint8(1-imgn);

 %{                     
s=[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;       %这个矩阵是维基百科中的矩阵
1 1 1 1 1 1 1 1 0 0 1 1 1 1 0 0;
0 0 1 1 1 1 0 0 0 1 1 1 1 0 0 0;
0 1 1 1 1 0 0 0 1 1 1 0 0 1 1 0;
1 1 1 0 0 1 1 0 0 0 1 1 1 0 0 0;
0 1 1 0 0 0 0 0 1 1 0 0 0 1 1 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];
%}
%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
    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)>=1 && pix1(1)<=m && pix1(2) >=1 &&pix1(2)<=n
                    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
%在连通域中心标记数字。
qq=zeros(2,label);  
for i=1:label
    % 计算目标区域中心,用于显示编号的位置
    [p,q]=find(tmp==i); %返回满足条件的行和列
    temp=[p,q];       
    [x,y]=size(temp); 
    qq(1,i)=sum(p)/x;%中心点
    qq(2,i)=sum(q)/x;
end
for i=1:label
    figure(2);
    text(qq(2,i),qq(1,i),int2str(i),'color','red')
end
figure,imshow(tmp)

%提取显示特定连通域
extract_tmp=tmp;
extract_tmp(extract_tmp ~=3 & extract_tmp ~=2) = 0;%%这边进行区域的选择,例如只保留2、3.
figure,imshow(extract_tmp);
%缺点,时间问题,图像边界没处理过。

This snippet took 0.01 seconds to highlight.

Back to the Entry List or Home.

Delete this entry (admin only).