matlab与数字图像处理--闭操作筛选图像中面积大于一定值的区域

 利用闭操作对图像进行图形元素的筛选，删除规格小于8*8的图形，保留大于8*8的图形

 

 测试用图，图像已经二值化，图像元素为白色，背景为黑色

 

选择结构元素为8*8的矩阵，对测试图像进行闭操作的结果，尺寸小于8*8的图像元素全部为消除。

%手动实现大面积图形查找
img=imread('ces0.bmp');

%关于img：img图像中进行了二值化，非零值（白色1），零值（黑色0）
%其中图像元素为非零值（白色1），图像背景为零值（黑色0）
img=im2bw(img,10/255);
% sc=[0,1,0;1,1,1;0,1,0]; %结构元素
sc=ones(8,8);

%确定结构元素的半径，以便对原图像进行边界的扩展
Hsc=size(sc,1);
Wsc=size(sc,2);
Rh=ceil(size(sc,1)/2); %结构元素的半径
Rw=ceil(size(sc,2)/2); %结构元素的半径
Csc=sum(sc(:)); %结构元素sc中关心的元素个数

Height=size(img,1);
Width=size(img,2);

%构建用于处理的imgin
imgin=zeros(Height+Hsc-1,Width+Wsc-1);
imgin(Rh:Rh+Height-1,Rw:Rw+Width-1)=img;
%对img进行腐蚀

imgout=imgin;
for i=Rh:Rh+Height-1
    for j=Rw:Rw+Width-1
        pitch=imgin(i-(Rh-1):i+Hsc-Rh,j-(Rw-1):j+Wsc-Rw);
       convtmp= pitch.*sc; %用img中的一块与结构元素作卷积
       %所谓腐蚀，指图像中关心的像素（结构元素sc中值为1的部分）有0，则结果赋值为0，没有0，则赋值1
        if sum(convtmp(:))<Csc %只要卷积结果中有零元素，则将输出图像的对应位置赋值0
             imgout(i,j)=0;
        else
            imgout(i,j)=1;
        end
    end
end

%对imgout进行膨胀
imgout1=imgout;
for i=Rh:Rh+Height-1
    for j=Rw:Rw+Width-1
        pitch=imgout(i-(Rh-1):i+Hsc-Rh,j-(Rw-1):j+Wsc-Rw);
       convtmp= pitch.*sc; %用imgout中的一块与结构元素作卷积
       %所谓腐蚀，指图像中关心的像素（结构元素sc中值为1的部分）有1，则结果赋值为1，没有1，则赋值0
        if sum(convtmp(:))>0 %只要卷积结果中有非零元素，则将输出图像的对应位置赋值1
             imgout1(i,j)=1;
        else
            imgout1(i,j)=0;
        end
    end
end
figure
imshow(imgout1,[]);

当选取结构元素为10*10的规格时，被删除掉的元素更多

 

 当选取结构元素为20*20的规格时：