OpenCV中Canny边缘检测的实现
作者:Allen Chou 发布时间:2022-10-17 10:10:19
边缘检测一般是识别目标图像中亮度变化明显的像素点. 因为显著变化的像素点通常反映了图像变化比较重要的地方.
1. Canny 边缘检测理论
Canny 是一种常用的边缘检测算法. 其是在 1986 年 John F.Canny 提出的.
Canny 是一种 multi-stage 算法,分别如下:
具体流程:
高斯滤波:平滑图像,消除噪声
梯度和方向计算:利用Sobel算子计算每个像素点的梯度和方向
非极大值抑制:消除边缘检测带来的杂散相应
双阈值:检测真正和潜在的边缘
抑制弱边缘:通过抑制弱边缘来完成边缘检测
1.1、高斯滤波
高斯滤波最重要的还是卷积核核,通常使用高斯平滑滤波器卷积降噪,这里以size=3的高斯内核为例:这里做了归一化处理(元素和为 1)
高斯去噪其实就是一个低通滤波器,滤除高频噪声。
1.2、Sobel算子计算梯度和方向
计算方法:
这里 是指水平方向的掩码模板,
是指垂直方向的掩码模板。根据上面的模板可以计算出图像梯度幅值和方向。
1.3、非极大值抑制(定位准确的边缘同时可缩小边缘线宽)
非极大值抑制是进行边缘检测的重要步骤,通俗的来说,就是获取局部的最大值,将非极大值所对应的灰度值设置为背景像素点。像素邻近区域满足梯度值的局部最优值判断为该像素的边缘,对非极大值相关信息进行抑制。利用这个准则可以剔除大部分的非边缘点。
简单的说呢?就是保留梯度大的像素点点,对于那些在边缘旁边的杂散点,梯度相对较小,利用非极大值抑制就可以很好的去除杂散点。
1.4、双阈值检测
这里的双阈值并不是说介于阈值之间的像素保留,外面的的去除。这里的阈值检测有所不同。
分析:
梯度大于maxVal的像素点保留,视为边缘
梯度小于minVal的像素点弃用,不视为边缘
梯度介于两者之间的,分情况判断:
四周有包含大于maxVal的像素点,视为边缘
四周没有大于maxVal的像素点,不视为边缘
根据上面的分析,我们可以得出来:A, D点位边界; B, C点不是边界。
注意:
具体这两个值怎么设置,我们就要分析两个值变化对图像的影响。
maxVal: 带来最明显的差异,增大maxVal无疑会导致原来的边界点可能会直接消失。但这种消失时是成片消失。
minVal: 增大minVal,会导致有些待定像素点被弃用,也就是靠近边界像素点的介于双阈值之间的被弃用。导致的现象就是边界出现破损,这种非成片消失。只是边界信息不完整。
下面以 video = cv2.Canny(img, 80, 250) 为例:分别增大minVal和maxVal。
增大minVal: (边界出现缺损)
增大maxVal: (边界出现成片消失,边界信息完整)
总结:
在实际应用中,观察梯度图像,如果边界信息缺损,那么适当的减小minVal;如果有不想要的区域出现,那么适当的增加MaxVal。
2. OpenCV 之 Canny 边缘检测
OpenCV 提供了 cv2.canny() 函数.
edge = cv2.Canny(image, threshold1, threshold2[, edges[, apertureSize[, L2gradient ]]])
参数 Image - 输入图片,必须为单通道的灰度图
参数 threshold1 和 threshold2 - 分别对应于阈值 minVal 和 maxVal
参数 apertureSize - 用于计算图片提取的 Sobel kernel 尺寸. 默认为 3.
参数 L2gradient - 指定计算梯度的等式. 当参数为 True 时,采用 1.2 中的梯度计算公式,其精度更高;否则采用的梯度计算公式为:
. 该参数默认为 False.
具体代码:
import cv2
import numpy as np
import matplotlib.pyplot as plt
img = cv2.imread('test.jpg', 0)
edges = cv2.Canny(img, 100, 200)
plt.subplot(121), plt.imshow(img, cmap='gray')
plt.title('Original Image'), plt.xticks([]), plt.yticks([])
plt.subplot(122),plt.imshow(edges,cmap = 'gray')
plt.title('Edge Image'), plt.xticks([]), plt.yticks([])
plt.show()
带 minVal 和 maxVal 滑动条的使用:
import numpy as np
import cv2
def nothing(x):
pass
img=cv2.imread('test.jpg',0)
cv2.namedWindow('res')
cv2.createTrackbar('min','res',0,25,nothing)
cv2.createTrackbar('max','res',0,25,nothing)
while(1):
if cv2.waitKey(1) & 0xFF == 27:
break
maxVal=cv2.getTrackbarPos('max','res')
minVal=cv2.getTrackbarPos('min','res')
canny=cv2.Canny(img,10*minVal,10*maxVal)
cv2.imshow('res',canny)
cv2.destroyAllWindows()
来源:https://blog.csdn.net/Vermont_/article/details/108471205


猜你喜欢
- 一个else语句可以使用if语句结合起来。如果在if语句中的条件表达式解析为0或false值,那么else语句包含代码执行。el
- 引言书接上回 Gradio机器学习模型快速部署工具【quickstart】翻译,讲到多输入输出,其实很简单,就是把多个组件包装到列表,inp
- 在使用Pytorch进行神经网络训练时,有时会遇到训练学习率不下降的问题。出现这种问题的可能原因有很多,包括学习率过小,数据没有进行Norm
- 需求:根据医保中心的文档和提供的dll动态库调用相关接口下载医保中心的账单。文档:对调用dll动态库的描述,调用哪个dll文件,同时了解清楚
- 每次写博客都是源于纳闷,python解析pcap这么常用的例子网上竟然没有,全是一堆命令行执行的python,能用吗?玩呢?pip安装sca
- 1.MySQL 5.5命令行里面set global log_slow_queries = on; &nb
- 对于一些数据量较大的系统,数据库面临的问题除了查询效率低下,还有就是数据入库时间长。特别像报表系统,每天花费在数据导入上的时间可能会长达几个
- 我们最近的项目中需要使用谷歌机器人验证,这个最主要的就是要有vpn,还需要有公司申请的google账号(自己申请的没用)用于商用的,利用这个
- 数据增强的必要性深度学习在最近十年得以风靡得益于计算机算力的提高以及数据资源获取的难度下降。一个好的深度模型往往需要大量具有label的数据
- 本篇文章主要介绍了python OpenCV学习笔记之绘制直方图的方法,小编觉得挺不错的,现在分享给大家,也给大家做个参考。一起跟随小编过来
- 目录原理先行实践篇完整代码原理先行作为一个资深的小说爱好者,国内很多小说网站如出一辙,什么 🖊*阁啊等等,大都是 get 请求返回 html
- 本文实例为大家分享了js鼠标经过表格某行时此行变色的具体代码,供大家参考,具体内容如下表格表头为蓝色,表主体为白色,当鼠标放到表的非表头行时
- 问题背景 基于PyQt5开发了一个可以用于目标跟踪的软件,在开发过程中遇到一个问题,就是如何在PyQt5的组件QLable中自主选定目标框
- 一、MySQL的下载 1.登陆MySQL的官网下载适用于64位系统的ZIP压缩包(https://dev.mysql.com/do
- 本文通过实例解析了 SQL Server 数据库扩展存储过程,实现远程备份与恢复的方法和步骤实例说明: 环境:win2k+sqlserver
- 本文要实现的功能是:根据下拉列表的选项将数据库中对应的内容显示在页面,选定要排除的选项后,提交剩余的选项到数据库。为了方便前后台交互,利用了
- 从过往MySQL数据库生产环境的维护工作中,总结的一些小经验和知识,未必有多深奥,但是对我们消除隐患,确保MySQL数据库生产环境四个9的作
- 一、统计语句1、--统计当前【>当天00点以后的数据】SELECT * FROM 表 WHERE CONVERT(Nvarchar,
- 如何用ASP来识别操作系统是vista的?我在网上找了个函数,但是不能判断是vista系统,希望大家帮忙. 这个是我在网上找的函数: Fun
- 二、XML的定义 XML是一个精简的SGML,它将SGML的丰富功能与HTML的易用性结合到Web的应用中。XML保留了SGML的可扩展功能