function?h=EdgeDetect（filename）
f=imread（filename）;
figureimshow（f）
[widthheight]=size（f）;
g=zeros（widthheight）;
h=zeros（widthheight）;
df=im2double（f）;

wx=[-1?0;1?0];
wy=[-1?1;0?0];
for?i=2:width-1
????for?j=2:height-1
????????gw=[df（ij）?df（ij+1）;df（i+1j）?df（i+1j+1）];
????????g（ij）=sqrt（（sum（sum（wx.*gw）））^2+（sum（sum（wy.*gw）））^2）;
????end
end

T=0.25*max（g（:））;
h=g>=T;
figureimshow（h）