%載入圖像
clear?all;
X=imread（‘Lena.bmp‘‘bmp‘）;

%分解
[cs]=wavedec2（X4‘db4‘）;?%用db4作三層小波分解
h=detcoef2（‘h‘cs2）;
v=detcoef2（‘v‘cs2）;??%提取水平細(xì)節(jié)系數(shù)和垂直細(xì)節(jié)系數(shù)
[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;


%求模極大值點(diǎn)
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