matlab图像处理之图像分割处理

中值滤波

利用中值滤波的方法对图像进行平滑处理，函数如下

fig_med = medfilt2(fig,[m,n]) % 对矩阵[m,n]进行中值滤波

类间方差阈值算法

利用类间方差阈值算法对滤波处理后图像进行分割处理，获取分割图像；
子函数如下所示：

function[threshold]=Ostu(fig)
s=256^2;
N=zeros(256,1);  %-----------各个灰度值出现次数
T=zeros(256,1);

for k=0:255      %-----------各个灰度值出现次数计数
    num=size(find(fig==k),1);
    N(k+1,1)=num;
end

for t=1:256
    w0=0; w1=0; u0=0; u1=0;
    
    for m=1:t
        w0=w0+N(m,1)/s;
        u0=u0+(m-1)*N(m,1)/s;
    end
    for m =t:256
        w1=w1+N(m,1)/s;
        u1=u1+(m-1)*N(m,1)/s;
    end
    u0=u0/w0;
    u1=u1/w1;
    T(t,1)=w0*w1*(u1-u0)^2;
end

threshold=find(T==max(max(T)))-1;
end

图像分割

m=3;
n=3;
fig = imread('dock.png');
fig = rgb2gray(fig);
subplot(1,3,1);
imshow(fig);
title('原图');

fig_med = medfilt2(fig,[m,n]);   %-------------中值滤波器

threshold=Ostu(fig_med);

fig_segment=fig_med;
fig_segment(fig_segment<threshold)=0;
fig_segment(fig_segment>=threshold)=255;

subplot(1,3,2);
imshow(fig_med);
title('中值滤波');

subplot(1,3,3);
imshow(fig_segment);
title('图像分割');

展示效果

笔者做了一个GUI界面，可以输入任意m，n值来改变中值滤波的效果。

原图：

效果如下：