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(:)))); % 此处的H(:)是什么意思?
x = T(P(:,2)); y = R(P(:,1)); %这一句是怎么回事儿?一点头绪都没有,抓瞎了。
plot(x,y,'s','color','white');
% Find lines and plot them
lines = houghlines(BW,T,R,P,'FillGap',5,'MinLength',7);
figure, imshow(rotI), hold on
max_len = 0; %这个变量怎么回事儿呢?感觉来的可突然,也不知道是干什么用的。
for k = 1:length(lines) % lines是什么?矩阵吗?什么矩阵?
xy = [lines(k).point1; lines(k).point2];% point1怎么回事儿?lines(k)
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
% highlight the longest line segment
plot(xy_long(:,1),xy_long(:,2),'LineWidth',2,'Color','cyan');
请各位大侠帮帮忙!不胜感激。