Извлечение сегментов линий на основе преобразования Хафа
lines = houghlines(BW,theta,rho,peaks)BW связано с определенными ячейками в преобразовании Хафа. theta и rho являются векторами, возвращаемыми функцией hough. peaks - матрица, возвращаемая houghpeaks функция, содержащая координаты строк и столбцов ячеек преобразования Хафа для поиска сегментов линий. Возвращаемое значение lines содержит информацию об извлеченных сегментах линий.
lines = houghlines(___,Name,Value)
Чтение изображения в рабочую область.
I = imread('circuit.tif');Поверните изображение.
rotI = imrotate(I,33,'crop');Создание двоичного изображения.
BW = edge(rotI,'canny');Создайте преобразование Хафа с помощью двоичного изображения.
[H,T,R] = hough(BW); imshow(H,[],'XData',T,'YData',R,... 'InitialMagnification','fit'); xlabel('\theta'), ylabel('\rho'); axis on, axis normal, hold on;
Найдите пики в преобразовании Хафа изображения.
P = houghpeaks(H,5,'threshold',ceil(0.3*max(H(:)))); x = T(P(:,2)); y = R(P(:,1)); plot(x,y,'s','color','white');
Найдите линии и постройте их график.
lines = houghlines(BW,T,R,P,'FillGap',5,'MinLength',7); figure, imshow(rotI), hold on max_len = 0; for k = 1:length(lines) xy = [lines(k).point1; lines(k).point2]; plot(xy(:,1),xy(:,2),'LineWidth',2,'Color','green'); % Plot beginnings and ends of lines plot(xy(1,1),xy(1,2),'x','LineWidth',2,'Color','yellow'); plot(xy(2,1),xy(2,2),'x','LineWidth',2,'Color','red'); % Determine the endpoints of the longest line segment len = norm(lines(k).point1 - lines(k).point2); if ( len > max_len) max_len = len; xy_long = xy; end end
Выделите самый длинный отрезок линии, раскрасив его голубым цветом.
plot(xy_long(:,1),xy_long(:,2),'LineWidth',2,'Color','cyan');
