function R=fharris(image,k)
W=fspecial(*gaussian*,[7,7],2);
dx=[-1,1;-1,1];
dy=[-1,-1;1,1];
image=filter2(W,image);
IX=filter2(dx,image);
IY=filter2(dy,image);
IX=filter2(W,IX);
IY=filter2(W,IY);
[m,n]=size(image);
for i=1:m
for j=1:n
M=[IX(i,j).^2,IX(i,j).*IY(i,j);IX(i,j).*IY(i,j),IY(i,j).^2];
R(i,j)=dec(M)-k*(trance(M).^2);
end
end
W=fspecial(*gaussian*,[7,7],2);
dx=[-1,1;-1,1];
dy=[-1,-1;1,1];
image=filter2(W,image);
IX=filter2(dx,image);
IY=filter2(dy,image);
IX=filter2(W,IX);
IY=filter2(W,IY);
[m,n]=size(image);
for i=1:m
for j=1:n
M=[IX(i,j).^2,IX(i,j).*IY(i,j);IX(i,j).*IY(i,j),IY(i,j).^2];
R(i,j)=dec(M)-k*(trance(M).^2);
end
end
