Matlab实现角点检测：基于Harris算子的焦点定位

 clear; 

ori_im = imread('Couple.bmp');     % 读取图像 

% fx = [5 0 -5;8 0 -8;5 0 -5];          % 高斯函数一阶微分，x方向(用于改进的Harris角点提取算法) 
fx = [-2 -1 0 1 2];                 % x方向梯度算子(用于Harris角点提取算法) 
Ix = filter2(fx,ori_im);              % x方向滤波 
% fy = [5 8 5;0 0 0;-5 -8 -5];          % 高斯函数一阶微分，y方向(用于改进的Harris角点提取算法) 
fy = [-2;-1;0;1;2];                 % y方向梯度算子(用于Harris角点提取算法) 
Iy = filter2(fy,ori_im);              % y方向滤波 
Ix2 = Ix.^2; 
Iy2 = Iy.^2; 
Ixy = Ix.*Iy; 
clear Ix; 
clear Iy; 

h= fspecial('gaussian',[7 7],2);      % 产生7*7的高斯窗函数，sigma=2 

Ix2 = filter2(h,Ix2); 
Iy2 = filter2(h,Iy2); 
Ixy = filter2(h,Ixy); 

height = size(ori_im,1); 
width = size(ori_im,2); 
result = zeros(height,width);         % 纪录角点位置，角点处值为1 

R = zeros(height,width); 

Rmax = 0;                              % 图像中最大的R值 
for i = 1:height