Hough变换检测直线(Matlab)

Hough变换是直线检测经典算法之一。原理介绍网上很多,不再赘述。参考网上代码,做了一些简单修改,稍微提高了一下速度。

matlab 复制代码
I  = imread('lineSample1.jpg');
HoughLineDetect(I,60);

%% 图像中直线检测
function I_out=HoughLineDetect(Img,thresholdValue)
%input:
%       img:输入图像;
%       thresholdValue:hough阈值;
%output:
%       I_out:检测直线结果,二值图像;
if  ~exist( 'thresholdValue', 'var' )
    thresholdValue = 150;
end

if length(size(Img))>2
    I_gray=rgb2gray(Img);               %如果输入图像是彩色图,需要转灰度图
else
    I_gray=Img;
end
Angle=180;
[xSize,ySize]=size(I_gray);             %图像大小
BW=edge(I_gray,'approxcanny',0.3);      %计算图像边缘

rho_Max=floor(sqrt(xSize^2+ySize^2));   %由图像坐标算出ρ最大值,作为极坐标系最大值
accArray=zeros(rho_Max,Angle);          %初始化极坐标系的数组
Theta=0:pi/Angle:(pi-pi/Angle);         %定义θ数组,范围从0-180度
Angles=1:Angle;

tic;
% hough变换
for n=1:xSize
    for m=1:ySize
        if BW(n,m)==1
            %hough变换方程求ρ值
            rho=abs(m*cos(Theta)+n*sin(Theta));
            rho_Int=floor(rho)+1;
            %在极坐标中标识点,相同点累加(投票)
            inds = sub2ind(size(accArray), rho_Int, Angles);
            accArray(inds)=accArray(inds)+1;
        end
    end
end

% 找出满足条件的直线
find_Ind=find(accArray>=thresholdValue);

% 将直线提取出来,输出图像数组I_out
I_out=zeros(xSize,ySize);
for n=1:xSize
    for m=1:ySize
        if BW(n,m)==1
            rho=abs(m*cos(Theta)+n*sin(Theta));
            rho_Int=floor(rho)+1;
            temp=intersect(find_Ind,rho_Int);
            if(size(temp,1)>=0)
                I_out(n,m)=BW(n,m);
            end
        end
    end
end
toc;

tiledlayout(1,2)
nexttile
imshow(Img);title('输入图像');
nexttile
imshow(I_out);title('Hough变换检测出的直线');

end

有几个函数需要Matalb的版本大于2019b。

简单效果如下图所示(运行速度还是很慢,后期再优化):

相关推荐
有点。1 小时前
C++深度优先搜索(三)
开发语言·c++·深度优先
天空'之城1 小时前
C 语言工业级通用组件手写 25:简易日志系统
c语言·开发语言·日志系统·嵌入式调试·调试打印
程序喵大人1 小时前
【C++进阶】STL算法与函数对象 -【C++进阶】STL算法与函数对象
开发语言·c++·算法
具身新纪元1 小时前
CVPR 2026|新缺陷不断上线,如何让工业视觉检测模型不忘旧账?
人工智能·深度学习·目标检测·机器学习·计算机视觉·目标跟踪·视觉检测
牛马也想出海2 小时前
使用Playwright被检测为机器人的原因及反检测方案
开发语言·网络·人工智能·机器人·php
GrowthDiary0072 小时前
Python 常用函数总结
开发语言·python
热心网友俣先生2 小时前
2026年华数杯C 题 超详细解题思路
c语言·开发语言·人工智能
啄缘之间2 小时前
5.1 uvm‑component 的层级结构管理? new (name,parent) 参数
开发语言·笔记·学习·ic·uvm·sv
旖旎夜光2 小时前
C++(类与对象)(下)
开发语言·c++·学习