Sobel梯度强度图G
提取梯度强度本质是一种滤波,同高斯模糊滤波一样,只不过所用算子(kernel)不同,sobel梯度图所用的算子如下两个矩阵,分别同图像A做卷积,得到两个大小与A相同的二维浮点矩阵Gx、Gy:
求Gx和Gy两个矩阵的平方和根号,得到二维浮点矩阵:
浮点矩阵G就是图像A的梯度强度图,大小与A相同的。
Matlab实现Sobel边缘提取
%copyright by Institute of Electronic ,Chinese Academy of Sciences
%get the gradient edge of optical image by sobel method
function [ G ] = Sobel( Image )
Image=double(Image);
h=[-1 -2 -1;0 0 0;1 2 1];
Gx=filter2(h,Image); %Gx is the gradient of x axel
Gy=filter2(h',Image); %Gy is the gradient of y axel
G=sqrt(Gx.^2+Gy.^2); %amplitude of gradient
G=255*G./max(G(:));
G=uint8(G); %trans from d