sobel算子，matlab实现

1. A=imread('D:\image\Body.jpg');    
2. I = rgb2gray(A);  
3.   
4. subplot(2,2,1);   
5. imshow(I);    
6. title('原图');   
7.   
8.   
9.  hx=[-1 -2 -1;0 0 0 ;1 2 1];%生产sobel垂直梯度模板  
10.  hy=hx';                             %生产sobel水平梯度模板  
11.    
12.  gradx=filter2(hx,I,'same');  
13.  gradx=abs(gradx); %计算图像的sobel垂直梯度  
14.  subplot(2,2,2);  
15.  imshow(gradx,[]);  
16.  title('图像的sobel垂直梯度');  
17.    
18.    
19.  grady=filter2(hy,I,'same');  
20.  grady=abs(grady); %计算图像的sobel水平梯度  
21.  subplot(2,2,3);  
22.  imshow(grady,[]);  
23.   title('图像的sobel水平梯度');  
24.     
25.   grad=gradx+grady;  %得到图像的sobel梯度  
26.   subplot(2,2,4);  
27.   imshow(grad,[]);  
28.   title('图像的sobel梯度'); 
    

结果为：

可以看出：

Sobel算子中的垂直模板得到的梯度图，由于梯度方向与边缘走向垂直，所以该梯度图对水平边缘有较强的响应，从而水平细节信息非常清晰；

sobel算子中的水平模板得到的梯度图，它对垂直边缘有较强的响应，垂直细节非常清晰。

Sobel算子水平和垂直方向叠加的梯度图，水平和垂直细节都非常清晰。

转自 http://blog.youkuaiyun.com/zhufanqie/article/details/8709910