matlab练习程序（条件膨胀）

顾名思义，条件膨胀是带有限定条件的膨胀。限定条件通常有两个。

1.膨胀始终在原始集合内。如下图：

原始集合图像：

集合内部的子集图像：

如果对后一个图进行规则1的条件膨胀，那么最后的结果就会得到第一幅图。所以这里第一幅图即是原始集合图像又是最终条件膨胀图像。

2.膨胀之后各自粒子互不连接。如下图：

原图：

按规则2条件膨胀需首先标记连通区域：

然后再进行互补侵犯的膨胀：

根据具体应用环境，可以使用单一的条件或两个条件同时使用。

分水岭分割的最后筑坝那一步就是两个条件同时使用到了。

单独规则1的matlab代码：

clear all;
close all;
clc;

mask=imread('p1.jpg');
mask=mask>128;

img=imread('p2.jpg');
img=img>128;            

[m n]=size(img);
imgn=zeros(m,n);
preimg=imgn;
while sum(sum(preimg-img))~=0
    preimg=img;
    
    for i=2:m-1
        for j=2:n-1
            if mask(i,j)==1     %在模板限定的范围之内膨胀，即为条件膨胀
                imgn(i,j)=max(max(img(i-1:i+1,j-1:j+1)));
            end
        end
    end
    img=imgn;

end

figure;
imshow(img,[]);

单独规则2的matlab代码：

clear all;
close all;
clc;

img=imread('l1.jpg');
img=img>128;            %二值化图像

img=liantong(img);      %连通标记，实现见先前的博文
[m n]=size(img);
imshow(img,[])

imgn=zeros(m,n);
preimg=imgn;
while sum(sum(preimg-img))~=0   %判断是否完全膨胀完
    preimg=img;
    
    for i=2:m-1
        for j=2:n-1
            if img(i,j)==0      %只对背景膨胀。
                flag=0;         %3*3区域是否含有两部分的标记变量，0为有一个前景一个背景，1为有多个前景一个背景 
                tmp=img(i-1:i+1,j-1:j+1);   %3*3区域
                ma=max(tmp(:));
                mi=min(tmp(:));
                for k=1:9           %判断3*3区域中是否含有两部分区域，如果有，则当前像素不再膨胀
                   if tmp(k)~=ma && tmp(k)~=mi;
                       flag=1;
                       break;
                   end
                end
                
                if flag==0          %有一个前景一个背景情况正常膨胀，否则不膨胀。不膨胀也就是筑坝了
                    imgn(i,j)=ma;
                end
                
            end
        end
    end
    
    img=imgn;  
end


for i=2:m-1         %将区块边缘的黑像素去掉
   for j=2:n-1 
       if img(i,j)==0
            imgn(i,j)=max(max(img(i-1:i+1,j-1:j+1)));
       end
   end
end

figure;
imshow(imgn,[]);