检测出图片中的直线且标记出图片中的最长线段

clc;
clear all;
I=imread('test.jpg');
figure;
imshow(I);
I=rgb2gray(I);
rot1=imrotate(I,33,'crop');  %旋转图像
bw=edge(rot1,'canny');   %使用边缘检测形成二值图像
[H,T,R]=hough(bw);   %生成点的叠加矩阵和坐标T,R
figure;
imshow(H,[],'XData',T,'YData',R,'InitialMagnification','fit');  %显示统计点的图形
title('变换结果');
xlabel('\theta'),ylabel('\rho');
axis on,axis normal;hold on; %使坐标易于观察
p=houghpeaks(H,5,'threshold',ceil(0.3*max(H(:))));  %提取5个极值点，p中为极坐标的位置
x=T(p(:,2));y=R(p(:,1));
plot(x,y,'s','color','blue');
lines=houghlines(bw,T,R,p,'FillGap',5,'MinLength',7);  %生成这五个点在图像中的直线矩阵
figure;
imshow(rot1);,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(xy(1,1),xy(1,2),'x','LineWidth',2,'Color','red');  %标出起点和终点
    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');  %标记出最大的直线