在Matlab中使用自行实现的霍夫线检测器来可视化线时需要帮助

如何解决在Matlab中使用自行实现的霍夫线检测器来可视化线时需要帮助

我能够使用以下代码在Matlab中重新创建hough（）和houghpeaks（）函数，但是我目前正在努力将它们从参数空间转换为图像空间。

%%% hough() function

corridor = imread('Corridor1.jpg');
grayCorridor = rgb2gray(corridor);
Canny_img = edge(grayCorridor,'Canny');
[rows,cols] = size(Canny_img);
theta_maximum = 90;
rho_maximum = hypot(rows,cols) - 1;
theta_range = -theta_maximum:theta_maximum-1;
rho_range = -rho_maximum:rho_maximum;

Hough = zeros(length(rho_range),length(theta_range));
for j = 1:rows
    for i = 1:cols
        if Canny_img(j,i) ==1
            x = i - 1;
            y = j - 1;
            for T = theta_range
                R = round((x * cosd(T)) + (y * sind(T)));
                R_Index = R + rho_maximum + 1;
                T_Index = T + theta_maximum + 1;
                Hough(R_Index,T_Index) = Hough(R_Index,T_Index) + 1;
            end
        end
    end
end

%%% houghpeaks() function

[X,Y] = find(Hough>0.47*max(Hough(:)))

%Visualising them 

imshow(Hough,[],'XData',theta_range,'YData',rho_range,...
            'InitialMagnification','fit');
xlabel('\theta'),ylabel('\rho');
axis on,axis normal,hold on;

pointsRho = rho_range(X(:));
pointsTheta = theta_range(Y(:));
plot(pointsTheta,pointsRho,'s','color','red');

这是我尝试将它们解码为图像空间。我首先将Rho和Theta配对成一个列表。然后，我创建了一组大小为imageRow x imageColumn x rhotheta对的二进制矩阵。 然后，我将所有x和y像素坐标代入函数x cos（theta）+ y sin（theta）= rho，并且如果left-right == 0，则将标记位置二进制矩阵（x，y）为1，我针对所有theta对对所有x和y值运行此函数，并将它们分别存储在上述二进制矩阵集中。

在第二个for循环中，我遍历了所获得的所有二进制矩阵集，并将它们聚合为一个二进制矩阵。我的问题来了，当我将去毛刺的线彼此靠近时，如何将它们组合成一条线？如果我能够做到这一点，则可以通过降低使用[X，Y] = find（Hough> 0.47 * max（Hough（:)））设置的阈值来可视化更多的霍夫峰。因为现在，如果我将阈值降低甚至0.01，那么我将有更多彼此接近的线。而且，您也许有一种更优雅，更好的dehough方法？

for k = 1:length(pointsRho)
    for i = 1:1280 % Shoulda put 1:size(Canny_img,1) but putting number helps me visualise my logic 
                   % better while I am coding
        for j = 1:960
            r = pointsRho(k);
            t = pointsTheta(k);
            threshold = abs(i*cosd(t)+j*sind(t)-r);
            if threshold ==0
                allLines(j,i,k) = 1;
            end     
        end
    end
end

aggregatedImage = zeros(size(Canny_img));
for k = 1:length(pointsRho)
    for j = 1:1280
        for i = 1:960
            if allLines(i,j,k) ==1
                aggregatedImage(i,j) = 1;
            end
        end
    end
end

figure,imshow(imfuse(aggregatedImage,corridor))

P / S：我想如果我要标记所有满足x cos（theta）+ y cos（ theta）= r仅在一个二进制矩阵中显示（如以下代码所示）。但是，我不确定为什么我做的时候aggregatedImage都全为0。因此，我决定采用上述方法。

aggregatedImage = zeros(size(Canny_img));

for k = 1:length(pointsRho)
    for i = 1:1280
        for j = 1:960
            r = pointsRho(k);
            t = pointsTheta(k);
            threshold = abs(i*cosd(t)+j*sind(t)-r);
            if threshold ==0
                aggregatedImage(j,i) = 1;
            end     
        end
    end
end