小波边缘检测，matlab程序

基于方向可调小波的边缘检测，包含matlab源程序，图片，算法介绍，还有一份详细的报告，希望读者都能看懂

%function [allTmpliulei]=test1（IM）
%IM是输入的原图像
tmp=imread（‘vvvvv.jpg‘）;
tmp=tmp（::1）;
IM=double（tmp）;
IM=IM（1:1181:174）;
figure（1）;
imshow（uint8（IM））;
[sizeXsizeY]=size（IM）;

%显示方向可调小波
jiao=-3.1415926/3;
pi=3.1415926;
for f=1:6
    [x1y1]=meshgrid（-5:0.25:5-5:0.25:5）;
    [mn]=size（x1）;
    for i=1:m
        for j=1:n
            x=x1（ij）;
            y=y1（1j）;
            
            z（ij）=-（x/（2*pi）*exp（-（x*x+y*y）/2））*cos（pi/6*（f-1））-（y/（2*pi）*exp（-（x*x+y*y）/2））*sin（pi/6*（f-1））;
        end
    end
    figure（2）;
    subplot（23f）;
    mesh（x1y1z）;
    colormap（[0 0 0]）;
end

l=0;
%x方向变换
for i=1:sizeX
    for j=1:sizeY
        tmp1=0;
        for m=-2:2
            for n=-2:2
                if（i-m）>0&&（i-m）<=sizeX&&（j-n）>0&&（j-n）<=sizeY
                    tmp1=tmp1+IM（i-mj-n）*m*（-1/（2*pi）*exp（-（m*m+n*n）/2））;
                end
            end
        end
        IMtmp1（ij）=tmp1;
    end
end

%y方向变换
for i=1:sizeX
    for j=1:sizeY
        tmp1=0;
        for m=-2:2
            for n=-2:2
                if（i-m）>0&&（i-m）<=sizeX&&（j-n）>0&&（j-n）<=sizeY
                    tmp1=tmp1+IM（i-mj-n）*n*（-1/（2*pi）*exp（-（m*m+n*n）/2））;
                end
            end
        end
        IMtmp2（ij）=tmp1;
    end
end


allTmp=0;
%求6个方向的平均
for pp=2:7
    jiao=pi/6*（pp-2）;
    
    
    IMtmp11=-IMtmp1*cos（jiao）;
    IMtmp22=-IMtmp2*sin（jiao）;
    IMtmp=IMtmp11+IMtmp22;
    
    
    IMtmpAbs=abs（IMtmp）;
    maxv=max（max（IMtmpAbs））;
    minv=min（min（IMtmpAbs））;
    minv1=minv*ones（sizeXsizeY）;
    
    
    IMtmpNom=（IMtmpAbs-minv1）*255/（maxv-minv）;
    liulei（::pp-1）=IMtmpNom;
    allTmp=allTmp+IMtmpNom;
end
%结果输出
allTmp=allTmp/6;
figure（3）;
imshow（uint8（allTmp））;







            

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      70438  2010-04-14 13:00  实验3\p6.bmp

     文件      70438  2010-04-14 13:00  实验3\p8.bmp

     文件       1841  2010-04-14 12:30  实验3\test1.asv

     文件       1852  2010-04-14 12:59  实验3\test1.m

     文件       1872  2010-04-14 18:30  实验3\test8.asv

     文件       1871  2010-04-14 19:00  实验3\test8.m

     文件     157025  2009-06-03 08:22  实验3\vvvvv.JPG

     文件     255210  2010-06-15 12:02  实验3\多媒体实验报告3.docx

     目录          0  2010-06-15 12:02  实验3

----------- ---------  ---------- -----  ----

               560547                    9