PythonOpenCVHough直线检测算法的原理实现

目录

直线检测原理

OpenCV实现

直线检测原理

核心要点:图像坐标空间、参数空间、极坐标参数空间 -> (极坐标)参数空间表决

给定一个点,我们一般会写成y=ax+b的形式,这是坐标空间的写法;我们也可以写成b=-xa+y的形式,这是参数空间的写法。也就是说,给定一个点,那么经过该点的直线的参数必然满足b=-xa+y这一条件,也就是必然在参数空间中b=-xa+y这条直线上。如果给定两个点,那么这两点确定的唯一的直线的参数,就是参数空间中两条参数直线的交点。

由于上述写法不适合处理水平或垂直的直线,我们可以使用极坐标的形式描述直线,即ρ=xcosθ+ysinθ,其中ρ是从原点到直线的垂直距离,θ是由这条垂直线和水平轴形成的角度(以逆时针方向测量),

如下图所示:

因此,任何垂直线θ=0,水平线θ=90°。那么极坐标参数空间中的曲线交点就是由两个点确定的一条直线,如下图所示。

现在让我们看看Hough变换是如何处理直线的。任何一条线都可以用这两个参数来表示(ρ,θ)。

首先创建一个二维数组,即累加器,用来保存两个参数的值,然后最初将其设置为全0。让行表示ρ,列表示θ。数组的尺寸取决于所需的精度。假设希望角度的精度为1度,则需要180列,枚举0°-179°的所有情况。对于ρ,可能的最大距离是图像的对角线长度。因此,以一个像素的精度计算,行数可以是图像的对角线长度。

枚举所有的点,对于每一个点,将所有经过这一点的直线对应的参数(ρ,θ)在参数空间中找到对应位置,令该位置的累加器加1,即投票。这一过程如下图所示。

枚举完成所有点之后,累加器中值最大的(若干个)参数组合(ρ,θ)就是经过点最多的(若干条)直线,如下图所示,两条直线对应累加器中最亮的两个点。

总的来说,对于多个点,我们可以用(离散)参数空间表决的方法,记录每个点对应的允许的参数组合,求得那些被允许次数最多的参数组合,就是最多点经过的直线。

在图像矫正任务中,我们经过Canny算子检测出了若干边缘点,这些点主要集中在四个边界上,因此我们只需要使用Hough直线检测,求出四条直线,就能确定四个边界。

OpenCV实现

cv.HoughLines()封装了上述步骤,该函数原型为:

cv.HoughLines(image, rho, theta, threshold[, lines[, srn[, stn[, min_theta[, max_theta]]]]]) -> lines

参数:

lines:数组,每一个元素都是一条直线对应的(ρ, θ),ρ以像素为单位,θ以弧度为单位。

image:输入图像,需要是二值图像,所以在应用hough变换之前应用阈值或canny边缘检测。

rho:ρ的精度。

theta:θ的精度。

threshold:阈值,得票数高于该值的线才被认为是线,由于投票数取决于线上的点数,所以它代表了应该被检测到的线的最小点数。

下面是具体代码:

def hough_detect(image_path): # 读取图像并转换为灰度图像 image = cv2.imread(image_path) gray = cv2.cvtColor(image, cv2.COLOR_BGR2GRAY) # 使用Canny算子检测边缘 edges = canny_detect(image_path, False) # 使用Hough检测直线 lines = cv2.HoughLines(edges, 1, np.pi/180, 200) # 绘制直线 for line in lines: rho, theta = line[0] a = np.cos(theta) b = np.sin(theta) x0 = a * rho y0 = b * rho x1 = int(x0 + 1000*(-b)) y1 = int(y0 + 1000*(a)) x2 = int(x0 - 1000*(-b)) y2 = int(y0 - 1000*(a)) cv2.line(image, (x1, y1), (x2, y2), (0, 0, 255), 2) cv2.imshow('line,jpg', image) cv2.waitKey() hough_detect('images/2.webp')

效果:

后面需要调整一下超参数。

到此这篇关于Python OpenCV Hough直线检测算法的原理实现的文章就介绍到这了,更多相关Python OpenCV Hough 内容请搜索易知道(ezd.cc)以前的文章或继续浏览下面的相关文章希望大家以后多多支持易知道(ezd.cc)!

推荐阅读

    电脑十进制算法|十进制的算法教程

    电脑十进制算法|十进制的算法教程,,十进制的算法教程0x10就是十六进制数10,转换为十进制数是16,即10(十六进制) = 16(十进制)。十六进制转换

    伪代码描述算法

    伪代码描述算法,算法,描述,伪代码,自然语言,语言,编程语言,  伪代码是自然语言和类编程语言组成的混合结构。它比自然语言更精确,描述算法很简

    进程调度详细总结|进程调度算法模拟

    进程调度详细总结|进程调度算法模拟,进程,优先级,一、概念: 无论是在批处理系统还是分时系统中,用户进程数一般都多于处理机数、这将导致它

    路由算法区分管理距离和最大跳数

    路由算法区分管理距离和最大跳数,路由,网络,路由器,路由协议,协议,状态,管理距离就是人为指定的一个数字,由这个数字来代表路由协议的优先度,数字

    php hmac_md5 算法函数

    php hmac_md5 算法函数,密钥,下吧,<?php /** * 基于md5的加密算法hmac * * md5已经不是那么安全了,多折腾几下吧 * * @param St

    图解KMP字符串匹配算法+代码实现

    图解KMP字符串匹配算法+代码实现,前缀,字符,  kmp算法跟之前讲的bm算法思想有一定的相似性。之前提到过,bm算法中有个好后缀的概念,而在km

    ORACLE 字段AES算法加密、解密

    ORACLE 字段AES算法加密、解密,密钥,方式,ORACLE 字段AES算法加密、解密(解决中文乱码问题)2014年02月12日 17:13:37 华智互联 阅读数:97971.

    谁能给一个C语言的遗传算法例题

    谁能给一个C语言的遗传算法例题,算法,算子,谁能给一个C语言的遗传算法例题这是一个非常简单的遗传算法源代码,是由Denis Cormier (North Ca