边缘检测sobel算法

边缘检测sobel算法：

原理：略（哈哈）

MATLAB实现：

clc;
clear;
imag = imread('1.jpg');  %读取关键帧
imag = rgb2gray(imag);        %转化为灰度图
figure,imshow(imag);
title('原图'); 
gx = [-1 -2 -1;0 0 0;1 2 1];
gy = [-1 0 1;-2 0 2;-1 0 1];
[high,width] = size(imag);   % 获得图像的高度和宽度      
imag1 = double(imag);
uSobel = zeros(high,width);
for i = 2:high - 1   %sobel边缘检测
    for j = 2:width - 1
        Gx = sum(sum(gx.*imag1(i-1:i+1, j-1:j+1))); %先求行和，再求列和
        Gy = sum(sum(gy.*imag1(i-1:i+1, j-1:j+1)));
        uSobel(i,j) = sqrt(Gx^2 + Gy^2); 
    end
end 
uSobel=mat2gray(uSobel);%将matrix值映射到[0,1],归一化
figure,imshow(uSobel);
title('边缘检测后');  %画出边缘检测后的图像
% Matlab自带函数边缘检测
% K为获取得到的关键帧的灰度图
BW3 = edge(imag,'sobel', 0.09);
figure,imshow(BW3,[]);
title('Matlab自带函数边缘检测');