%基于小波变换的多尺度图像边缘检测matlab源代码clearall;loadwbarb;I=ind2gray(X,map);imshow(I);%¼ìË÷ͼת³É»Ò¶ÈͼI1=imadjust(I,stretchlim(I),[0,1]);%µ÷ÕûͼÏñµÄÏñËØÖµ£¬¿ÉÒԸıä¶Ô±È¶ÈºÍÑÕÉ«figure;imshow(I1);[N,M]=size(I);h=[0.125,0.375,0.375,0.125];g=[0.5,-0.5];delta=[1,0,0];J=3;a(1:N,1:M,1,1:J+1)=0;%¸³³õÖµdx(1:N,1:M,1,1:J+1)=0;dy(1:N,1:M,1,1:J+1)=0;d(1:N,1:M,1,1:J+1)=0;a(:,:,1,1)=conv2(h,h,I,'same');%C=conv2(A,B)C=conv2(A,B)%·µ»Ø¾ØÕóAºÍBµÄ¶þά¾í»ýC¡£ÈôAΪma¡ÁnaµÄ¾ØÕó£¬BΪmb¡ÁnbµÄ¾ØÕó£¬ÔòCµÄ´óСΪ(ma+mb+1)¡Á(na+nb+1)dx(:,:,1,1)=conv2(delta,g,I,'same');dy(:,:,1,1)=conv2(g,delta,I,'same');x=dx(:,:,1,1);y=dy(:,:,1,1);d(:,:,1,1)=sqrt(x.^2+y.^2);I1=imadjust(d(:,:,1,1),stretchlim(d(:,:,1,1)),[01]);figure;imshow(I1);lh=length(h);lg=length(g);forj=1:J+1lhj=2^j*(lh-1)+1;lgj=2^j*(lg-1)+1;hj(1:lhj)=0;gj(1:lgj)=0;forn=1:lhhj(2^j*(n-1)+1)=h(n);endforn=1:lggj(2^j*(n-1)+1)=g(n);enda(:,:,1,j+1)=conv2(hj,hj,a(:,:,1,j),'same');dx(:,:,1,j+1)=conv2(delta,gj,a(:,:,1,j),'same');dy(:,:,1,j+1)=conv2(gj,delta,a(:,:,1,j),'same');x=dx(:,:,1,j+1);y=dy(:,:,1,j+1);dj(:,:,1,j+1)=sqrt(x.^2+y.^2);I1=imadjust(dj(:,:,1,j+1),stretchlim(dj(:,:,1,j+1)),[01]);figure;imshow(I1);end