ori_img=imread('C:\Program Files\MATLAB\R2013a\bin\Original_Images\DIP3E_Original_Images_CH03\Fig0342(a)(contact_lens_original).tif');
ori_img=mat2gray(ori_img); % 归一化 转化为double
sobel_v=double([1,2,1;0,0,0;-1,-2,-1]); %Sobel operator 3x3 convolution kernels
% 垂直 Sobel_v kron([1,0,-1]',[1,2,1])
sob_v_img=conv2(ori_img,sobel_v,'same');
% sob_v_img=mat2gray(sob_v_img);
figure,imshow(sob_v_img,[]),title('sobel vetical img');
sobel_h=double([1,0,-1;2,0,-2;1,0,-1]); %Sobel operator 3x3 convolution kernels
sob_h_img=conv2(ori_img,sobel_h,'same');
% sob_h_img=mat2gray(sob_h_img);
figure,imshow(sob_h_img,[]),title('sobel Hor img');
sob_img=abs(sob_h_img)+abs(sob_v_img);
figure,imshow(sob_img,[]),title('sobel img') % 归一化
% 水平 sobel_h kron([1,2,1]',[1,0,-1])
%sobel_h=double([-3,0,3;-10,0,10;-3,0,3]); %更准确
%