请问如何得到小波边缘检测的模极大值，一定要通过编程实现么，请高手指点，谢谢！

snowyzy

请问如何得到小波边缘检测的模极大值，一定要通过编程实现么，请高手指点，谢谢！

happy

1. clear all;
   
2. load wbarb;
   
3. I = ind2gray(X,map);imshow(I);
   
4. I1 = imadjust(I,stretchlim(I),[0,1]);figure;imshow(I1);
   
5. [N,M] = size(I);
   
6. 
   
7. h = [0.125,0.375,0.375,0.125];
   
8. g = [0.5,-0.5];
   
9. delta = [1,0,0];
   
10. 
    
11. J = 3;
    
12. 
    
13. a(1:N,1:M,1,1:J+1) = 0;
    
14. dx(1:N,1:M,1,1:J+1) = 0;
    
15. dy(1:N,1:M,1,1:J+1) = 0;
    
16. d(1:N,1:M,1,1:J+1) = 0;
    
17. 
    
18. a(:,:,1,1) = conv2(h,h,I,'same');
    
19. dx(:,:,1,1) = conv2(delta,g,I,'same');
    
20. dy(:,:,1,1) = conv2(g,delta,I,'same');
    
21. 
    
22. x = dx(:,:,1,1);
    
23. y = dy(:,:,1,1);
    
24. d(:,:,1,1) = sqrt(x.^2+y.^2);
    
25. I1 = imadjust(d(:,:,1,1),stretchlim(d(:,:,1,1)),[0 1]);figure;imshow(I1);
    
26. 
    
27. lh = length(h);
    
28. lg = length(g);
    
29. 
    
30. for j = 1:J+1
    
31. lhj = 2^j*(lh-1)+1;
    
32. lgj = 2^j*(lg-1)+1;
    
33. hj(1:lhj)=0;
    
34. gj(1:lgj)=0;
    
35. for n = 1:lh
    
36. hj(2^j*(n-1)+1)=h(n);
    
37. end
    
38. 
    
39. for n = 1:lg
    
40. gj(2^j*(n-1)+1)=g(n);
    
41. end
    
42. 
    
43. a(:,:,1,j+1) = conv2(hj,hj,a(:,:,1,j),'same');
    
44. dx(:,:,1,j+1) = conv2(delta,gj,a(:,:,1,j),'same');
    
45. dy(:,:,1,j+1) = conv2(gj,delta,a(:,:,1,j),'same');
    
46. 
    
47. x = dx(:,:,1,j+1);
    
48. y = dy(:,:,1,j+1);
    
49. dj(:,:,1,j+1) = sqrt(x.^2+y.^2);
    
50. 
    
51. I1 = imadjust(dj(:,:,1,j+1),stretchlim(dj(:,:,1,j+1)),[0 1]);figure;imshow(I1);
    
52. end

复制代码

看看是不是这样？