凹凸线未到达图像边缘

Question

我试图检测下图中的所有行：

使用以下代码我能够检测到几乎所有行：

%Find the line detected
I = imread('img3.png');
J = histeq(I);
BW = edge(J,'canny');
[H,T,R] = hough(BW);
P = houghpeaks(H,30,'Threshold',0.7*max(H(:)),'nhoodsize',[5 5]);
lines = houghlines(BW,T,R,P );
figure, imshow(I, []), 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','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;
 xy_long = xy;
 end
end
% highlight the longest line segment
plot(xy_long(:,1),xy_long(:,2),'LineWidth',2,'Color','red');

这给出了以下图像

但是，正如您所看到的，线条并没有一直延伸到图像的边缘。我尝试缩小 nhood 大小，但这意味着对某些线条进行双重检测，并且不会将线条延伸到边缘。 还有可能检测到最顶线和最底线吗？我知道这些线由于较短而不会获得那么多选票，我修改了阈值，但最终出现了虚假的对角线，并且仍然没有检测到顶线和底线。

另外，我如何在未知图像上设置 Houghlines 参数。这里非常简单，能够（粗略地）估计我期望看到的行数，然后围绕该数字进行调节。

提前致谢

Answer 1

该问题源于 houghpeaks 属性过于严格，无法找到较小的线段。我使用的方法是：

1. 使用您当前的代码来获取检测到的线所在的角度（我刚刚做了[lines.theta]并发现theta将等于-84)

2. 再次使用您的代码，但仅包含该角度，同时放宽 houghpeaks 属性以允许检测更多线条。

3. 您将检测到更多重叠的线条，我们将通过根据其系数设置相似性阈值来丢弃这些线条。

这是我们收集的实现，

angleList = -84; 
[H,T,R] = hough(BW, 'Theta', angleList);
P = houghpeaks(H,80,'Threshold', 0.1*max(H(:)),'nhoodsize',[1 1]);
lines = houghlines(BW,T,R,P,'FillGap',500,'MinLength',5 );

% just get the eq of the line coef a,b using polyfit (y=a*x+b)
line_coef= @(x) polyfit([lines(x).point1(1) lines(x).point2(1)],...
                        [lines(x).point1(2) lines(x).point2(2)],1);

for n=1:numel(lines)
    p(:,n)=line_coef(n);
end

其中哪条线有相似的线，使得 coef 距离低于阈值

for n=1:numel(lines)-1
     d(:,n) = sum(abs(p(:,n)-p)) <2 ;
end

d 是对称二元相关矩阵，我们找到对角线上方上三角中的所有点，这是线 id#在同一行重复，所以我们丢弃它们。

id=find(sum(triu(d,1)));
lines(id)=[];

%... plot the lines using your code