图像分割——边缘检测——Sobel算子——梯度（Matlab）

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

%%
%=============================边缘检测(一)=================================
%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);
ga=zeros(M,N);
%门限
th=0.33;
gth=zeros(M,N);
n_l=floor(n/2);
%对原图进行扩展，方便处理边界
I_pad=padarray(I,[n_l,n_l],'symmetric');
for i=1:M
    for j=1:N
        %获得图像子块区域
        Block=I_pad(i:i+2*n_l,j:j+2*n_l);
        %用Sobel内核对子区域卷积 
        %x方向绝对值
        gx(i,j)=abs(sum(sum(Block.*Lx)));
        %y方向绝对值
        gy(i,j)=abs(sum(sum(Block.*Ly)));
        %梯度强度
        g(i,j)=abs(sum(sum(Block.*Lx)))+abs(sum(sum(Block.*Ly)));
        %角度
        ga(i,j)=atan(gy(i,j)/gx(i,j));
        %带有门限的梯度强度
        if g(i,j)>=th
            gth(i,j)=g(i,j);
        end
    end
end
subplot(221)
imshow(gx)
subplot(222)
imshow(gy)
subplot(223)
imshow(g)
subplot(224)
% imshow(ga)

imshow(gth)