harris 角点检测 matlab代码

以下是一个完整的 Harris 角点检测的 Matlab 代码示例：

% 读入图像
img = imread('image.jpg');

% 转换为灰度图像
gray = rgb2gray(img);

% 计算 Harris 角点响应函数
sigma = 2; % 高斯滤波器的标准差
alpha = 0.04; % Harris 角点响应函数的参数
hsize = ceil(3*sigma)*2 + 1; % 高斯滤波器的尺寸
h = fspecial('gaussian', hsize, sigma);
Ix = conv2(gray, [-1 0 1; -1 0 1; -1 0 1], 'same');
Iy = conv2(gray, [-1 -1 -1; 0 0 0; 1 1 1], 'same');
Ix2 = conv2(Ix.^2, h, 'same');
Iy2 = conv2(Iy.^2, h, 'same');
Ixy = conv2(Ix.*Iy, h, 'same');
R = (Ix2.*Iy2 - Ixy.^2) - alpha*(Ix2 + Iy2).^2;

% 非极大值抑制
corner_mask = imregionalmax(R);
[r, c] = find(corner_mask);
corners = [c, r];

% 可视化结果
imshow(img); hold on;
plot(corners(:,1), corners(:,2), 'r+', 'MarkerSize', 10);

其中，`fspecial` 函数用于生成高斯滤波器，`conv2` 函数用于进行卷积操作，`imregionalmax` 函数用于进行非极大值抑制。最终可得到检测到的角点，并将其可视化在原始图像上。