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,[]);