Торов Собела, LoG и Канни.

f = imread('C: \Users\Шеф\Desktop\Building.tif');

[g_sobel_default, ts] = edge(f, 'sobel');

[g_log_default, tlog] = edge(f, 'log');

[g_canny_default, tc] = edge(f, 'canny');

g_sobel_best = edge(f, 'sobel', 0.05);

g_log_best = edge(f, 'log', 0.003, 2.25);

g_canny_best = edge(f, 'canny', [0.04 0.10], 1.5);

Пример 10.5. Иллюстрация преобразования Хафа.

f = zeros(101, 101);

f(1, 1) = 1; f(101, 1) = 1; f(1, 101) = 1;

f(101, 101) = 1; f(51, 51) = 1;

H = hough(f);

imshow(H, [ ])

[H, theta, rho] = hough(f);

imshow (theta, rho, H, [ ], 'notruesize')

axis on, axis normal

xlabel('\theta'), ylabel('\rho')

Пример 10.6. Использование преобразования Хафа для обнаружения линий и связывания.

f = imread('C: \Users\Шеф\Desktop\Building2.jpg');

[H, theta, rho] = hough(f, 0.5);

imshow(theta, rho, H, [ ], 'notruesize'), axis on, axis normal

xlabel('\theta'), ylabel('\rho')

[r, c] = houghpeaks(H, 5);

hold on

plot(theta(c), rho(r), 'linestyle', 'none',...

'marker', 's', 'color', 'w')

lines = houghlines (f, theta, rho, r, c)

figure, imshow(f), hold on

for k = 1: length(lines)

xy = [lines(k).point1; lines(k).point2];

plot(xy(:, 2), xy(:, 1), 'LineWidth', 4, 'Color3', [.6.6.6]);

end