Реализация детектора углов Harris
Я использую детектор углов Harris для образовательных целей, но я застрял в части ответа Харриса. В основном то, что я делаю, это:
Вычислить градиенты интенсивности изображения в x- и y-направленииРазмытие выхода (1)Вычислить ответ Харриса по выходу (2)Подавить немаксимумы на выходе (3) в 3x3-окрестности и пороговом выходе1 и 2, кажется, работают нормально; тем не менее, я получаю очень маленькие значения в качестве ответа Харриса, и ни одна точка не достигает порога. Вход - стандартная уличная фотография.
[...]
[Ix, Iy] = intensityGradients(img);
g = fspecial('gaussian');
Ix = imfilter(Ix, g);
Iy = imfilter(Iy, g);
H = harrisResponse(Ix, Iy);
[...]
function K = harrisResponse(Ix, Iy)
max = 0;
[sy, sx] = size(Ix);
K = zeros(sy, sx);
for i = 1:sx,
for j = 1:sy,
H = [Ix(j,i) * Ix(j,i), Ix(j,i) * Iy(j,i)
Ix(j,i) * Iy(j,i), Iy(j,i) * Iy(j,i)];
K(j,i) = det(H) / trace(H);
if K(j,i) > max,
max = K(j,i);
end
end
end
max
end
Для примера изображения, max в конечном итоге составляет 6.4163e-018, что кажется слишком низким.