這個只需要遍歷一次圖像就能夠完全標記了。我主要參考了WIKI和這位兄弟的博客,這兩個把原理基本上該介紹的都介紹過了,我也不多說什么了。一步法代碼相比兩步法真是清晰又好看matlab在圖像中加標記matlab在圖像中加標記,似乎真的比兩步法要好很多。
代碼如下:
1 clear all;
2 close all;
3 clc;
4
5 img=imread('liantong.bmp');
6 imgn=img>128;
7 s=uint8(1-imgn);
8
9 %{

10 s=[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
11 0 0 1 1 0 0 1 1 0 0 1 1 0 0 1 1 0; %這個矩陣是維基百科中的矩陣
12 0 1 1 1 1 1 1 1 1 0 0 1 1 1 1 0 0;
13 0 0 0 1 1 1 1 0 0 0 1 1 1 1 0 0 0;
14 0 0 1 1 1 1 0 0 0 1 1 1 0 0 1 1 0;
15 0 1 1 1 0 0 1 1 0 0 0 1 1 1 0 0 0;
16 0 0 1 1 0 0 0 0 0 1 1 0 0 0 1 1 0;
17 0 0 0 0 0 0 1 1 1 1 0 0 1 1 1 1 0;
18 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0];

19 %}
20 imshow(mat2gray(s));
21 [m n]=size(s);
22 tmp=zeros(m,n); %標記圖像
23 label=1;
24 queue_head=1; %隊列頭
25 queue_tail=1; %隊列尾
26 neighbour=[-1 -1;-1 0;-1 1;0 -1;0 1;1 -1;1 0;1 1]; %和當前像素坐標相加得到八個鄰域坐標
27

28 for i=2:m-1
29 for j=2:n-1
30
31 if s(i,j)==1 && tmp(i,j) ==0
32 tmp(i,j)=label;
33 q{queue_tail}=[i j]; %用元組模擬隊列,當前坐標入列
34 queue_tail=queue_tail+1;
35
36 while queue_head~=queue_tail

37 pix=q{queue_head};
38 for k=1:8 %8鄰域搜索
39 pix1=pix+neighbour(k,:); if pix1(1)>=2 && pix1(1)<=m-1 && pix1(2) >=2 &&pix1(2)<=n-1
40 if s(pix1(1),pix1(2)) == 1 && tmp(pix1(1),pix1(2)) ==0 %如果當前像素鄰域像素為1并且標記圖像的這個鄰域像素沒有被標記,那么標記
41 tmp(pix1(1),pix1(2))=label;
42 q{queue_tail}=[pix1(1) pix1(2)];
43 queue_tail=queue_tail+1;
44 end end
45 end

46 queue_head=queue_head+1;
47 end
48
49 clear q; %清空隊列,為新的標記做準備
50 label=label+1;
51 queue_head=1;
52 queue_tail=1;
53 end
54
55 end
56 end
57 figure,imshow(mat2gray(tmp))