使用霍夫变换检测电路图像中的直线方法

I=imread('circuit.tif');%旋转图像并寻找边缘
rotI=imrotate(I,33,'crop');
subplot(221),fig1=imshow(rotI);
BW=edge(rotI,'canny');%使用canny方法检测目标边界
subplot(222),imshow(BW);
%执行Hough变换并显示Hough矩阵
[H,T,R]=hough(BW);
subplot(223),imshow(H,[],'XData',T,'YData',R,'InitialMagnification','fit');
xlabel('\theta'),ylabel('\rho');
axis on,axis normal,hold on;

%在Hough矩阵中寻找前5个大于Hough矩阵中最大值0.3倍的峰值
P=houghpeaks(H,5,'threshold',ceil(0.3*max(H(:))));%hough变换的极值点
x=T(P(:,2));y=R(P(:,1));%由行、列索引转换成实际坐标
plot(x,y,'s','color','white');%在Hough矩阵图像中标出峰值位置
%找出并绘制直线
lines=houghlines(BW,T,R,P,'FillGap',5,'MinLength',7);%合并距离小于5的线段，丢弃所有长度小于7的直线段
subplot(224),imshow(rotI),hold on
max_len=0;
for k = 1:length(lines)%依次标出各条直线段
    xy=[lines(k).point1;lines(k).point2];
    plot(xy(:,1),xy(:,2),'x','LineWidth',2,'Color','green');

%绘制线段端点
    plot(xy(1,1),xy(1,2),'x','LineWidth',2,'Color','yellow');
    plot(xy(2,1),xy(2,2),'x','LineWidth',2,'Color','red');
%确定最长的线段
    len=norm(lines(k).point1-lines(k).point2);
    if(len>max_len)
       max_len=len;