小波模极大值算法边缘检测

%载入图像
clear all;
X=imread（‘Lena.bmp‘‘bmp‘）;

%分解
[cs]=wavedec2（X4‘db4‘）; %用db4作三层小波分解
h=detcoef2（‘h‘cs2）;
v=detcoef2（‘v‘cs2）;  %提取水平细节系数和垂直细节系数
[NM]=size（h）;
%求模图
for i=1:N;
    for j=1:M;
        Mo（ij）=sqrt（h（ij）^2+v（ij）^2）;
    end;
end;

%求相角的正切
for i=1:N
    for j=1:M
        Tan（ij）=v（ij）/h（ij）;
    end;
end;


%求模极大值点
edge1=zeros（NM）;
for i=2:N-1;
    for j=2:M-1;
        if（Tan（ij）>=-1-sqrt（2） && Tan（ij）<1-sqrt（2））
            if（Mo（ij）>=Mo（i-1j-1） && Mo（ij）>=Mo（i+1j+1））
                edge1（ij）=Mo（ij）;
            else
                edge1（ij）=0;
            end;
        elseif（Tan（ij）>=1-sqrt（2） && Tan（ij）            if（Mo（ij）>=Mo（ij-1） && Mo（ij）>=Mo（ij+1））
                edge1（ij）=Mo（ij