clear?all;
close?all;
%%%%%%%%%%%%%%載入圖像%%%%%%%%%%%%%%
I=imread（‘C:\Users\HP\Desktop\rongchi3.bmp‘）;

%?I=imnoise（I‘gaussian‘00.01）;
X=double（I）;
[IxIy]=size（X）;???????????%圖像長寬
%%%%%%%%%%%%%%二級尺度%%%%%%%%%%%%%%????
m=1.0;
delta=2^m;
%%%%%%%%%高斯函數XY方向的偏導%%%%%%%
N=20;?????????????????????%長度和幅度
A=-1/sqrt（2*pi）;
phi_x=zeros（NN）;
phi_y=zeros（NN）;
for?i=1:N?????????????????%偏導
????for?j=1:N
????????x=i-（N+1）/2;
????????y=j-（N+1）/2;
????????phi_x（ij）=A*（x/delta^2）.*exp（-（x*x+y*y）/（2*delta^2））;
????????phi_y（ij）=A*（y/delta^2）.*exp（-（x*x+y*y）/（2*delta^2））;
????end
end
%%%%%%%%%圖象行、列卷積%%%%%%%%%%%%%%
W_x=conv2（Xphi_x‘same‘）;
W_y=conv2（Xphi_y‘same‘）;
%%%%%%%%%%%%求梯度%%%%%%%%%%%%%%%%%%%
Mod_G=zeros（IxIy）;
for?i2=1:Ix
????for?j2=1:Iy
????????Mod_G（i2j2）=sqrt（W_x（i2j2）*W_x（i2j2）+W_y（i2j2）*W_y（i2j2））;
????end
end
M=fix（Mod_G）;
%%%%%%%%%%%%求幅角%%%%%%%%%%%%%%%%%%%
Arg=zeros（IxIy）;
for?i3=1:Ix
????for?j3=1:Iy
????????angle=atan（W_y（i3j3）/W_x（i3j3））*180/pi;??????%求反正切
????????if?（（W_x（i3j3）>0）&&（W_y（i3j3）>0））????????????%第一象限
????????????angle=angle+0;
????????elseif?（（W_x（i3j3）<0）&&（W_y（i3j3）>0））????????%第二象限
????????????angle=angle+180;
????????elseif?（（W_x（i3j3）<0）&&（W_y（i3j3）<0））????????%第三象限
????????????angle=angle+180;
????????elseif?（（W_x（i3j3）>0）&&（W_y（i3j3）<0））????????%第四象限
????????????angle=angle+360;???????????????????????
????????elseif?（abs（W_y（i3j3）/W_x（i3j3））>100000000）??%大值90度
????????????angle=90;
????????elseif?（abs（W_y（i3j3）/W_x（i3j3））==0）?????????%零值0度
????????????angle=0;
????????end
????????if（angle>180）??????????????????????????????????%幅角限制在0-180????????????????????????
????????????angle=angle-180;
????????end
????????Arg（i3j3）=angle;
????end
end
%%%%%%%%%%%%圖像邊界%%%%%%%%%%%%%%%%%%%%
%沿著梯度方向檢測小波變換系數模的局部極大%
%%%%%%%%值點即可得到圖像的邊緣點%%%%%%%%%
Edge=zeros（IxIy）;
for?i4=2:（Ix-1）
????for?j4=2:（Iy-1）
????????if?（（Arg（i4j4）>=0&&Arg（i4j4）<=22.5））
????????????if?（Mod_G（i4j4）>Mod_G（i4+1j4）&&Mod_G（i4j4）>Mod_G（i4-1j4））%水平
????????????????Edge（i4j4）=Mod_G（i4j4）;
????????????end
????????elseif?（（Arg（i4j4）>=（90-22.5）&&Arg（i4j4）<=（90+22.5）））
????????????if?（Mod_G（i4j4）>Mod_G（i4j4+1）&&Mod_G（i4j4）>Mod_G（i4j4-1））%45度
????????????????Edge（i4