图像配准技术及MATLAB编程实现_陈显毅_pdf_源码

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）