空间滤波(五） 综合运用

ori_img=imread('C:\Program Files\MATLAB\R2013a\bin\Original_Images\DIP3E_Original_Images_CH03\Fig0343(a)(skeleton_orig).tif');
stretch_img=intrans(ori_img,'stretch',mean2(im2double(ori_img)),0.9);

w_lap=double([0,1,0;1,-4,1;0,1,0]); %edge_detect   边缘检测  图像锐化，拉普拉斯算子
ori_img=mat2gray(ori_img);
laplace_img=conv2(ori_img,w_lap,'same');
lap_plus_ori_img=mat2gray(ori_img-laplace_img); % 锐化图像
imshow(ori_img),title('ori_img');
figure,imshow(stretch_img),title('stretch img');
figure,imshow(laplace_img,[]),title('laplace img')
figure,imshow(lap_plus_ori_img,[]),title('lap plus ori img')
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=d