利用Opencv中Houghline方法实现直线检测


Posted in Python onFebruary 11, 2018

利用Opencv中的Houghline方法进行直线检测—python语言

这是给Python部落翻译的文章,请在这里看原文。

在图像处理中,霍夫变换用来检测任意能够用数学公式表达的形状,即使这个形状被破坏或者有点扭曲。

下面我们将看到利用HoughLine算法来阐述霍夫变化进行直线检测的原理,把此算法应用到特点图像的边缘检测是可取的。边缘检测方法请参考这篇文章?边缘检测。

Houghline算法基础

直线可以表示为y=mx+c,或者以极坐标形式表示为r=xcosθ+ysinθ,其中r是原点到直线的垂直距离,θ是水平轴顺时针方向到垂直线的夹角(这个方向取决于坐标的形式,在OpenCV就是采用这种极坐标形式)。

利用Opencv中Houghline方法实现直线检测

因此任意的曲线都可以用两个参数(r,θ)表示。

HoughLine算法原理:

  • 首先建立一个二维的数组或者累加器(用来保存这两个参数),并初始化为零;
  • 这个二维数组的行代表不同的r,而列代表角度θ;
  • 数组的大小取决于算法的精度。假设所需角度的精度精确到1∘,那么就需要180列(直线的最大角度是180)。
  • 对于r,最大的可能距离是图像的对角长度,因此若需要一个像素的精度,那么需要把行数设为图像对角线的长度。

例子:

假设一幅100x100的图像,在图像中间有一条水平直线。设直线的第一个点的坐标为(x,y),在直线方程中,令参数θ=0,12,⋯,180,观查参数r变化。对每一个参数对(r,θ),在累加器中将(r,θ)对应的单元格中的值递增1,比如现在在累加器中,某个单元(50,90)的值等于1,其它的值也如此。

对于直线上的第二个点,重复上述操作。将得到的参数对(r,θ)的对应值继续递增,然后(50,90)对应的值等于2。实现上我们是对参数对(r,θ)进行投票,对直线上的每一个点重复上述操作,对每一个点,单元格(50,90)对应的值会递增,或者说投票给参数对(50,90),而会或者不会投票给其它参数对。以这种方式,最后单元格(50,90)的值将会是最大的值。然后搜索累加器的最大值,将会找到参数对(50,90)。也就是说,在图像中找到了到原点距离为50,角度为90的一条直线。

利用Opencv中Houghline方法实现直线检测

上述算法的过程被封装成OpenCV函数cv2.HoughLines(),函数返回(r,θ)的一个数组,其中r的单位为像素,θ的单位为弧度。

# Python program to illustrate HoughLine
# method for line detection
import cv2
import numpy as np

# Reading the required image in 
# which operations are to be done. 
# Make sure that the image is in the same 
# directory in which this python program is
img = cv2.imread('xyz.jpg')

# Convert the img to grayscale
gray = cv2.cvtColor(img,cv2.COLOR_BGR2GRAY)

# Apply edge detection method on the image
edges = cv2.Canny(gray,50,150,apertureSize = 3)

# This returns an array of r and theta values
lines = cv2.HoughLines(edges,1,np.pi/180, 200)

# The below for loop runs till r and theta values 
# are in the range of the 2d array
for r,theta in lines[0]:

 # Stores the value of cos(theta) in a
 a = np.cos(theta)

 # Stores the value of sin(theta) in b
 b = np.sin(theta)

 # x0 stores the value rcos(theta)
 x0 = a*r

 # y0 stores the value rsin(theta)
 y0 = b*r

 # x1 stores the rounded off value of (rcos(theta)-1000sin(theta))
 x1 = int(x0 + 1000*(-b))

 # y1 stores the rounded off value of (rsin(theta)+1000cos(theta))
 y1 = int(y0 + 1000*(a))

 # x2 stores the rounded off value of (rcos(theta)+1000sin(theta))
 x2 = int(x0 - 1000*(-b))

 # y2 stores the rounded off value of (rsin(theta)-1000cos(theta))
 y2 = int(y0 - 1000*(a))

 # cv2.line draws a line in img from the point(x1,y1) to (x2,y2).
 # (0,0,255) denotes the colour of the line to be 
 #drawn. In this case, it is red. 
 cv2.line(img,(x1,y1), (x2,y2), (0,0,255),2)

# All the changes made in the input image are finally
# written on a new image houghlines.jpg
cv2.imwrite('houghlines3.jpg', img)

函数(cv2.HoughLines(edges, 1, np.pi/180,200))

  • 第一个参数是输入图像,且必须是二值图像,在进行霍夫变换之前需要采用阈值方法的边缘检测;
  • 第二和第三个参数分别是r,θ对应的精度;
  • 第四个参数是阈值,判定为直线投票数的最小值;
  • 注意,投票数取决于直线上点的个数,因此这个阈值代表了检测到的直线的最短长度。

利用Opencv中Houghline方法实现直线检测

总结

  • 在图像分析应用中,边缘分割点的坐标(即X,Y)是已知的,在直线的参数方程中作为输入,而r,θ是我们要找的未知变量;
  • 对每一个角度θ,求出r的值,也就是将笛卡尔图像空间的点映射到霍夫参数极坐标空间的曲线,这种点到曲线的变换就是直线的霍夫变换;
  • 这种变换是通过量化霍夫参数为有限的间隔或者累加器单元来实现的,随着算法的运行,每一个(X,Y)转换成离散的(r,θ)曲线,曲线上的点对应的累加器(二维数组)的值会递增;
  • 累加器中的峰值就代表了图像中有相应直线的存在。

霍夫变换的应用

  • 用于分离图像中特殊形状的特征;
  • 可以容忍特征边界描述的误差,并且不受噪声的干扰;
  • 广泛用于条形码扫描、验证和识别。

文章作者是Pratima Upadhyay,如果你喜欢GeeksforGeeks,并且愿意分享,可以利用contribute.geeksforgeeks.org写文章,然后发送到contribute@geeksforgeeks.org,在GeeksforGeeks中看到自己的文章,帮助更多的Geeks。

若发现错误欢迎指正,也欢迎评论。

以上就是本文的全部内容,希望对大家的学习有所帮助,也希望大家多多支持三水点靠木。

Python 相关文章推荐
用python删除java文件头上版权信息的方法
Jul 31 Python
python高级特性和高阶函数及使用详解
Oct 17 Python
Python参数解析模块sys、getopt、argparse使用与对比分析
Apr 02 Python
Python判断对象是否为文件对象(file object)的三种方法示例
Apr 26 Python
Python pandas实现excel工作表合并功能详解
Aug 29 Python
python 中值滤波,椒盐去噪,图片增强实例
Dec 18 Python
浅谈python之自动化运维(Paramiko)
Jan 31 Python
Python GUI库PyQt5图形和特效样式QSS介绍
Feb 25 Python
Centos7下源码安装Python3 及shell 脚本自动安装Python3的教程
Mar 07 Python
查看已安装tensorflow版本的方法示例
Apr 19 Python
Python简单实现词云图代码及步骤解析
Jun 04 Python
降低python版本的操作方法
Sep 11 Python
tensorflow输出权重值和偏差的方法
Feb 10 #Python
详解tensorflow实现迁移学习实例
Feb 10 #Python
Python学习之Django的管理界面代码示例
Feb 10 #Python
Tensorflow 自带可视化Tensorboard使用方法(附项目代码)
Feb 10 #Python
tensorflow训练中出现nan问题的解决
Feb 10 #Python
用Eclipse写python程序
Feb 10 #Python
tensorflow建立一个简单的神经网络的方法
Feb 10 #Python
You might like
Zend的MVC机制使用分析(二)
2013/05/02 PHP
dedecms集成财付通支付接口
2014/12/28 PHP
javascript的对话框详解与参数
2007/03/08 Javascript
javascript 星级评分效果(手写)
2012/12/24 Javascript
angularjs实现与服务器交互分享
2014/06/24 Javascript
javascript消除window.close()的提示窗口
2015/05/20 Javascript
js鼠标点击按钮切换图片-图片自动切换-点击左右按钮切换特效代码
2015/09/02 Javascript
JavaScript实现跑马灯抽奖活动实例代码解析与优化(一)
2016/02/16 Javascript
同步文本框内容JS代码实现
2016/08/04 Javascript
AngularJS ngModel实现指令与输入直接的数据通信
2016/09/21 Javascript
Vue.js中轻松解决v-for执行出错的三个方案
2017/06/09 Javascript
Javascript实现的StopWatch功能示例
2017/06/13 Javascript
node.js将MongoDB数据同步到MySQL的步骤
2017/12/10 Javascript
vue学习教程之带你一步步详细解析vue-cli
2017/12/26 Javascript
nodejs实现套接字服务功能详解
2018/06/21 NodeJs
JavaScript实现简单的计算器
2020/01/16 Javascript
Ant Design的Table组件去除
2020/10/24 Javascript
[02:37]2018DOTA2亚洲邀请赛赛前采访 VP.no[o]ne心中最强SOLO是谁
2018/04/04 DOTA
Python中set与frozenset方法和区别详解
2016/05/23 Python
浅谈django中的认证与登录
2016/10/31 Python
解决pyqt中ui编译成窗体.py中文乱码的问题
2016/12/23 Python
python2.7 mayavi 安装图文教程(推荐)
2017/06/22 Python
Django admin.py 在修改/添加表单界面显示额外字段的方法
2019/08/22 Python
20行Python代码实现视频字符化功能
2020/04/13 Python
Python如何实现FTP功能
2020/05/28 Python
Python timeit模块原理及使用方法
2020/10/10 Python
AmazeUI 单选框和多选框的实现示例
2020/08/18 HTML / CSS
总监职责范文
2013/11/09 职场文书
更夫岗位责任制
2014/02/11 职场文书
教师节促销活动方案
2014/02/14 职场文书
颁奖典礼主持词
2014/03/25 职场文书
村党支部对照检查材料思想汇报
2014/09/28 职场文书
2015新学期校长寄语(3篇)
2015/03/25 职场文书
会计工作自我鉴定范文
2019/06/21 职场文书
MySQL系列之六 用户与授权
2021/07/02 MySQL
详解python的异常捕获
2022/03/03 Python