其實寫過一步法二值影象連通區域標記之後我就感覺這個標記和填充基本上是一回事，所以我這裡就用了一步法的那個佇列演算法。也沒什麼好說的，演算法幾乎一樣，就是細節的區別。還有這裡使用了ginput函式。至於堆疊版的實現，看情況吧。

 1 clear all;
 2 close all;
 3 clc;
 4 
 5 img=imread('liantong.bmp');
 6 img=img>128;
 7 img=mat2gray(img);
 8 imshow(img);
 9 
10 [m n]=size(img);
11 [x y]=ginput();
12 x=round(x);
13
 
 y=round(y);
14 
15 tmp=ones(m,n);
16 queue_head=1;       %佇列頭
17 queue_tail=1;       %佇列尾
18 neighbour=[-1 -1;-1 0;-1 1;0 -1;0 1;1 -1;1 0;1 1];  %和當前畫素座標相加得到八個鄰域座標
19 %neighbour=[-1 0;1 0;0 1;0 -1];     %四鄰域用的
20 q{queue_tail}=[y x];
21 queue_tail=queue_tail+1;
22 [ser1 ser2]=size(neighbour);
23 
24 while queue_head~=queue_tail
 
25     pix=q{queue_head};
26     for i=1:ser1
27         pix1=pix+neighbour(i,:);
28         if pix1(1)>=1 && pix1(2)>=1 &&pix1(1)<=m && pix1(2)<=n
29             if img(pix1(1),pix1(2))==1 
30                 img(pix1(1),pix1(2))=0;
31                 q{queue_tail}=[pix1(1
 
) pix1(2)];
32                 queue_tail=queue_tail+1;
33 
34             end      
35         end
36     end
37 
38     queue_head=queue_head+1;
39 end
40 
41 figure(1);
42 imshow(mat2gray(img));

效果圖：

原圖

下面是不同的填充效果：