suanzi.m

用所写算子进行边缘检测

function y=suanzi(image ,th)
[nx,ny]=size(image);
%image=double(image);
h=[1 2 1;0 0 0;-1 -2 -1];
Gx =filter2(h ,image);
Gy = filter2(h' ,image);
F=abs(Gx)+abs(Gy);
th_F=double(F)/255;
for i=1:nx
   for j=1:ny
       if th_F(i,j)<th
         y(i,j)=0;
      else
           y(i,j)=1;
   end
end
end