图像分割——边缘检测——均值+Sobel算子（Matlab）

clc;
clear all;
close all;
I=im2double(imread('D:\Gray Files\10-16.tif'));
[M,N]=size(I);

%%
%=============================边缘检测(二)=================================
%对图像先进行均值处理，而后用Sobel算子
%5*5均值算子
n=5;
g_mean=zeros(M,N);
n_l=floor(n/2);
%对原图进行扩展，方便处理边界
I_pad_m=padarray(I,[n_l,n_l],'symmetric');
for i=1:M
    for j=1:N
        %获得图像子块区域
        Block=I_pad_m(i:i+2*n_l,j:j+2*n_l);
        %用拉氏内核对子区域卷积     
        g_mean(i,j)=mean(mean(Block));
    end
end
%Sobel算子
Lx=[-1 -2 -1;
    0 0 0;
    1 2 1];
Ly=[-1 0 1;
     -2 0 2;
     -1 0 1];   
 %检测对角
L_45=[0 1 2;
    -1 0 1;
    -2 -1 0];
L_135=[-2 -1 0;
    -1 0 1;
    0 1 2];
n=3;
%输出图像
g=zeros(M,N);
gx=zeros(M,N);
gy=zeros(M,N);
%门限
th=0.33;
g_th=zeros(M,N);
n_l=floor(n/2);
%对原图进行扩展，方便处理边界
I_pad_s=padarray(g_mean,[n_l,n_l],'symmetric');
for i=1:M
    for j=1:N
        %获得图像子块区域
        Block=I_pad_s(i:i+2*n_l,j:j+2*n_l);
        %用Sobel内核对子区域卷积     
        gx(i,j)=abs(sum(sum(Block.*Lx)));
        gy(i,j)=abs(sum(sum(Block.*Ly)));
        g(i,j)=abs(sum(sum(Block.*Lx)))+abs(sum(sum(Block.*Ly)));
        if g(i,j)>=th
            g_th(i,j)=g(i,j);
        end
    end
end

subplot(221)
imshow(gx)
subplot(222)
imshow(gy)
subplot(223)
imshow(g)
subplot(224)
imshow(g_th)