小波多尺度边缘检测

Halcom

小波多尺度边缘检测:

1. clc,clear,close all;clear all;
   
2. warning off
   
3. I = imread('lena.jpg');
   
4. % I = imread('lion.jpg');
   
5. if size(I,3)>1
   
6.     I = rgb2gray(I);
   
7. end
   
8. h = [0.125,0.375,0.375,0.125];   % 滤波掩膜
   
9. gx = [0.5,-0.5];   % 水平掩膜
   
10. gy = [0.5;-0.5];   % 垂直掩膜
    
11. 
    
12. imfilter_iter(:,:,:,1) = conv2(h,h,I,'same'); % 图像滤波
    
13. dx(:,:,:,1) = conv2(I,gx,'same');             % X方向卷积
    
14. dy(:,:,:,1) = conv2(I,gy,'same');             % Y方向卷积
    
15. 
    
16. figure('color',[1,1,1])
    
17. subplot(131),imshow(I,[]),title('original image')
    
18. subplot(132),imshow(dx(:,:,:,1),[]),title('X方向卷积')
    
19. subplot(133),imshow(dy(:,:,:,1),[]),title('Y方向卷积')
    
20. 
    
21. imfilter_norm(:,:,:,1) = sqrt(dx(:,:,:,1).^2+dy(:,:,:,1).^2);     % 模
    
22. % 阈值分割
    
23. % threshold = graythresh(imfilter_norm(:,:,:,1));
    
24. % imfilter_bw(:,:,:,1) = im2bw(imfilter_norm(:,:,:,1), threshold);
    
25. imfilter_bw(:,:,:,1) = imfilter_norm(:,:,:,1) > 5;
    
26. 
    
27. S = 3;   % 尺度
    
28. for i=1:S
    
29.     imfilter_iter(:,:,:,i+1) = conv2(h,h,imfilter_iter(:,:,:,i),'same');   % 图像滤波
    
30.     dx(:,:,:,i+1) = conv2(imfilter_iter(:,:,:,i),gx,'same');               % X方向卷积
    
31.     dy(:,:,:,i+1) = conv2(imfilter_iter(:,:,:,i),gy,'same');               % Y方向卷积
    
32.     imfilter_norm(:,:,:,i+1) = sqrt(dx(:,:,:,i+1).^2+dy(:,:,:,i+1).^2);    % 模
    
33.     % 阈值分割
    
34.     %     threshold = graythresh(imfilter_norm(:,:,:,i+1));
    
35.     %     imfilter_bw(:,:,:,i+1) = im2bw(imfilter_norm(:,:,:,i+1), threshold);
    
36.     imfilter_bw(:,:,:,i+1) = imfilter_norm(:,:,:,i+1) > 5;
    
37. end
    
38. 
    
39. figure('color',[1,1,1])
    
40. subplot(321),imshow(I,[]),title('original image')
    
41. subplot(322),imshow(imfilter_norm(:,:,:,1),[]),title('小波多尺度边缘检测--滤波图像模')
    
42. subplot(334),imshow(imfilter_norm(:,:,:,2),[]),title('小波多尺度边缘检测--S=1图像模')
    
43. subplot(335),imshow(imfilter_norm(:,:,:,3),[]),title('小波多尺度边缘检测--S=2图像模')
    
44. subplot(336),imshow(imfilter_norm(:,:,:,4),[]),title('小波多尺度边缘检测--S=3图像模')
    
45. subplot(337),imshow(imfilter_bw(:,:,:,2),[]),title('小波多尺度边缘检测--S=1图像边缘')
    
46. subplot(338),imshow(imfilter_bw(:,:,:,3),[]),title('小波多尺度边缘检测--S=2图像边缘')
    
47. subplot(339),imshow(imfilter_bw(:,:,:,4),[]),title('小波多尺度边缘检测--S=3图像边缘')
    
48. 
    
49. figure('color',[1,1,1])
    
50. subplot(121),imshow(imfilter_norm(:,:,:,4),[]),title('小波多尺度边缘检测--S=3图像模')
    
51. subplot(122),imshow(imfilter_bw(:,:,:,4),[]),title('小波多尺度边缘检测--S=3图像边缘')

复制代码