原始图
matlab软件中的参考代码
clc;
clearvars;
close all;
I = imread('C:\Users\xiaosa\Desktop\IMG_0638.JPG');
I = I(:,:,1);
rotI = I;
thresh = graythresh(rotI);
bwTemp = im2bw(rotI, thresh);
[H,T,R] = hough(bwTemp);
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','red');
lines = houghlines(bwTemp,T,R,P,'FillGap',5,'MinLength',7);
% lines = houghlines(bwTemp,T,R,P);
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');