OpenCV图像处理入门_opencv alpha和beta-程序员宅基地

技术标签: python  计算机视觉  图像处理  人工智能  OpenCV知识祥解(基础-进阶-项目)  opencv  

欢迎来到本博客
本次博客内容将继续讲解关于OpenCV的相关知识
作者简介️️️目前计算机研究生在读。主要研究方向是人工智能和群智能算法方向。目前熟悉python网页爬虫、机器学习、计算机视觉(OpenCV)、群智能算法。然后正在学习深度学习的相关内容。以后可能会涉及到网络安全相关领域,毕竟这是每一个学习计算机的梦想嘛!
目前更新目前已经更新了关于网络爬虫的相关知识、机器学习的相关知识、目前正在更新计算机视觉-OpenCV。
本文摘要

本文我们将继续讲解OpenCV的相关操作。
在这里插入图片描述

在这里插入图片描述

️1.OpenCV礼帽操作和黑帽操作

图像顶帽操作也叫图像礼帽操作,实际上就是:
结果图像=原图像-图像开运算
那么我们知道图像的开运算就是对图像先腐蚀操作,然后在进行膨胀操作,得到的图像其实就是进行了取噪的一个处理,然后我们所说的礼帽图像操作就是使用原图像-图像的开运算操作。
函数依旧是:

result = cv2.morphologyEx(img, cv2.MORPH_TOPHAT, kernel)

其中img表示原始图像,cv2.MORPH_TOPHAT表示进行礼帽操作,然后kernel表示卷积核,这里我们之前已经讲过。
核心函数是:

k=np.ones((5,5),np.uint8)
r=cv2.morphologyEx(o,cv2.MORPH_TOPHAT,k)

这样我们显示以下我们处理的图像就可以知道:
在这里插入图片描述
处理之后我们得到了相应的噪声结果。

图像黑帽操作就是图像的闭运算-原图像
黑帽结果图像=图像闭运算-原图像
我们之前也介绍过这个闭运算操作就是先对图像进行膨胀操作,然后进行腐蚀操作。
函数是:

result = cv2.morphologyEx(img, cv2.MORPH_BLACKHAT, kernel)

其中cv2.MORPH_BLACKHAT表示的就是黑帽操作的意思,kernel表示卷积核。
然后我们看一下图像处理的核心函数:

k=np.ones((5,5),np.uint8)
r=cv2.morphologyEx(o,cv2.MORPH_BLACKHAT,k)

得到的结果是:
在这里插入图片描述
可以看到我们把图像中的小气泡都取出来了。

️2.Sobel算子理论基础及实际操作

首先我们来了解一个边界的定义,看一个图来理解更加的明白。
在这里插入图片描述
对于如图中的水平梯度,我们看A、B两个地方,右侧像素值减去左侧像素值不为0,那么我们就说在此图中A列和B列是边界,否则不是边界。同样对于垂直梯度我们也是这样定义。
对于Sobel算子,我们先来看x方向他进行了什么操作:
在这里插入图片描述
P5x=(p3-p1)+2*(p6-p4)+(p9-p7)
右侧像素值减去左侧像素值,中间行参数稍大为2。
同样对于y轴方向也做了一个同样的操作,
P5y=(p7-p1)+2*(p8-p2)+(p9-p3)
下一行像素值减去上一行像素值,中间列参数稍大为2
然后我们计算了一个近似梯度值:
G= 根号(2 + 2)
简化版本就是:
G=||+| |
在这里插入图片描述
这中心点的P5的Sobel算子就是按照这个方式进行计算。
对于P5点完成的Sobel算子计算就是这样:
在这里插入图片描述
那么好,我们来看一下在python中对于Sobel是一个什么样的函数:

dst = cv2.Sobel( src , ddepth , dx , dy , [ksize] )

其中src表示原始图像,ddeph表示图像的深度,那么我们一般设定为-1,用来表示和原图像保持一致。dx,dy表示x轴的方向或者y轴的方向,ksize表示核大小。
其中我们在已经知道的256色位图中,白色点像素值255,黑色点像素值0。
这里对于dx,dy方向为什么我们取绝对值进行一下解答,因为加入一块黑块在一个白色背景下,那么我们在水平方向上,左侧那么得到的数值就是一个负数,对于右侧那么我们得到的是就是一个正常的正数,那么对于OpenCV来说,当我们没有加绝对值得时候如果出现了负数,统一处理成0。
在这里插入图片描述
根据图我们就可以知道,左侧的线被OpenCV处理成了0。所以我们这里要加上一个绝对的运算。无论对于水平方向还是垂直方向都是如此。
然后我们要将原始图像处理成256为的色位图。其函数是:

dst = cv2.convertScaleAbs( src [, alpha[, beta]] )

目标图像=调整(原始图像*alpha+beta)
其中alpha和beta可以理解为权重的意思。
当我们想要计算水平方向上的Sobel算子的时候,那么我们就设定函数中的参数dx=1,dy=0.对于垂直方向也是如此。那么这里我们就会想:如果我们两个方向都一起设置这样是不是比计算完x的方向然后计算y方向,然后进行结合这样是不是更加的方便一些呢?也就是说下图中的方式一比方式二更加的简介方便呢?
在这里插入图片描述
事实上不是这样的,而且两个图像有很大的差异,这个我们后续会看到结果。
那么当我们完成了这一步操作之后,我们还要做两个图像的权重和,其函数是:

dst=cv2.addWeighted( src1 , alpha , src2 , beta , gamma )

这里我们说的通俗一点就是,src1表示计算的x轴方向的图像,alpha表示x轴方向上的图像所占的比重,src2表示计算的y轴方向的图像,beta表示y轴方向上的图像所占的比重是多少,最后gamma表示一个修正数值。

dst(I)=saturate(src1(I)*alpha+src2(I)*beta+gamma)

接下来我们来验证一下各个结果:

import cv2
import numpy as np
o = cv2.imread('image\\sobel4.bmp',cv2.IMREAD_GRAYSCALE)
sobelx = cv2.Sobel(o,-1,1,0)
cv2.imshow("original",o)
cv2.imshow("x",sobelx)
cv2.waitKey()
cv2.destroyAllWindows()

这里我们没有加上绝对值计算x轴方向上的Sobel算子,得到的结果是:
在这里插入图片描述
很明显左侧的没有计算出来。

import cv2
import numpy as np
o = cv2.imread('image\\sobel4.bmp',cv2.IMREAD_GRAYSCALE)
sobelx = cv2.Sobel(o,cv2.CV_64F,1,0)
sobelx = cv2.convertScaleAbs(sobelx) # 转回uint8 
cv2.imshow("original",o)
cv2.imshow("x",sobelx)
cv2.waitKey()
cv2.destroyAllWindows()

当我们加上绝对值后的结果是:
在这里插入图片描述
对于垂直方向的结果也是如此。最后我们计算出x,y的方向。进行融合并且和直接将x,y表示成1的结果进行一下对比,我们来看结果:
核心代码:

sobelx = cv2.Sobel(o,cv2.CV_64F,1,0)
sobely = cv2.Sobel(o,cv2.CV_64F,0,1)
sobelx = cv2.convertScaleAbs(sobelx) # 转回uint8 
sobely = cv2.convertScaleAbs(sobely) 
sobelxy = cv2.addWeighted(sobelx,0.5,sobely,0.5,0) 
sobelxy11=cv2.Sobel(o,cv2.CV_64F,1,1)
sobelxy11=cv2.convertScaleAbs(sobelxy11)

在这里插入图片描述
根据结果我们可以得到:使用了dx,dy的情况下无法检测到边。然后我们用lena来测试一下结果。
在这里插入图片描述
在我们相加表示的情况下,边缘检测的线条十分清晰,将分割的地带表达的特别清楚。

️3.Scharr算子简介及相关操作

经过大量的科学家验证,发现3*3的Sobel算子可能并不是特别的精确,于是又提出了相应的Scharr算子,其实就是改变了运算的尺度而已。
在这里插入图片描述
Sobel算子由原来的121变为了3 10 3,这样的一个组合。其余的操作都是和Sobel算子一致我们来简要的说明一下:
dst=Scharr(src, ddpeth, dx, dy)
这里面少了核的大小,但是和Sobel算子依然一直,参数这里就不解释了,唯一不同的就是我们使用Sobel算子可以设置成dx=1,dy=1。那么对于Scharr算子如果这样设置就会报错。这里需要注意一下。满足条件:
dx >= 0 && dy >= 0 && dx+dy == 1

import cv2
import numpy as np
o = cv2.imread('image\\scharr.bmp',cv2.IMREAD_GRAYSCALE)
scharrx = cv2.Scharr(o,cv2.CV_64F,1,0)
scharry = cv2.Scharr(o,cv2.CV_64F,0,1)
scharrx = cv2.convertScaleAbs(scharrx) # 转回uint8 
scharry = cv2.convertScaleAbs(scharry) 
scharrxy = cv2.addWeighted(scharrx,0.5,scharry,0.5,0) 
cv2.imshow("original",o)
cv2.imshow("x",scharrx)
cv2.imshow("y",scharry)
cv2.imshow("xy",scharrxy)
cv2.waitKey()
cv2.destroyAllWindows()

在这里插入图片描述

️4.Sobel算子和Scharr算子的比较

对于上述我们的测试图,其实看不出来Sobel算子和Scharr算子到底有什么区别,看起来是一致的,那么我们这里使用大美女lena进行一次Sobel和Scharr,得到结果进行对比一下可以知道:

sobelx = cv2.Sobel(o,cv2.CV_64F,1,0,ksize=3)
sobely = cv2.Sobel(o,cv2.CV_64F,0,1,ksize=3)
sobelx = cv2.convertScaleAbs(sobelx) # 转回uint8 
sobely = cv2.convertScaleAbs(sobely) 
sobelxy = cv2.addWeighted(sobelx,0.5,sobely,0.5,0) 
scharrx = cv2.Scharr(o,cv2.CV_64F,1,0)
scharry = cv2.Scharr(o,cv2.CV_64F,0,1)
scharrx = cv2.convertScaleAbs(scharrx) # 转回uint8 
scharry = cv2.convertScaleAbs(scharry) 
scharrxy = cv2.addWeighted(scharrx,0.5,scharry,0.5,0)

在这里插入图片描述
很明显,这里我们看一下了lena的肩膀部分,可以看到使用Shcarr算子进行计算的时候,刻画的更加细致。

️5.laplacian算子简介及相关操作

拉普拉斯算子类似于二阶sobel导数。实际上,在OpenCV中通过调用sobel算子来计算拉普拉斯算子。使用的公式为:
∆ = 2/2 + 2/2
使用的卷积核也和Sobel和Shcarr算子的不一致:
在这里插入图片描述

看起来有点像高斯的,哈哈哈。
这里我们进行一下讲解,
P5new=(p2+p4+p6+p8)-4*p5
在这里插入图片描述
当没有处于边界的时候呢,我们可以得到P5的数值是变化不大的,当处于边界的时候,P5的数值是变化非常大的,这是因为什么呢?
在这里插入图片描述
对于拉普拉斯算子的函数估计我们也可以猜到:
dst = cv2.Laplacian( src, ddepth )
自然是没有水平方向和垂直方向的。
实际操作中,计算梯度值可能会出现负数。通常处理的图像是np.uint8类型,如果结果也是该类型,所有负数会自动截断为0,发生信息丢失。所以,通常计算时,使用更高的数据类型cv2.CV_64F,取绝对值后,再转换为np.uint(cv2.CV_8U)类型。
所以我们还是要取一次绝对值的运算:
dst = cv2.convertScaleAbs( src )

import cv2
import numpy as np
o = cv2.imread('image\\laplacian.bmp',cv2.IMREAD_GRAYSCALE)
laplacian = cv2.Laplacian(o,cv2.CV_64F)
laplacian = cv2.convertScaleAbs(laplacian) # 转回uint8 
cv2.imshow("original",o)
cv2.imshow("laplacian",laplacian)
cv2.waitKey()
cv2.destroyAllWindows()

在这里插入图片描述
我们对lena进行一次拉普拉斯算子操作,得到的结果是:
在这里插入图片描述

️6.Canny边缘检测的原理

Canny边缘检测算子是John F. Canny于 1986 年开发出来的一个多级边缘检测算法。更为重要的是 Canny 创立了边缘检测计算理论(Computational theory of edge detection)解释这项技术如何工作。
Canny 的目标是找到一个最优的边缘检测算法,最优边缘检测的含义是:
1.最优检测:算法能够尽可能多地标识出图像中的实际边缘,漏检真实边缘的概率和误检非边缘的概率都尽可能小;
2.最优定位准则:检测到的边缘点的位置距离实际边缘点的位置最近,或者是由于噪声影响引起检测出的边缘偏离物体的真实边缘的程度最小;
3.检测点与边缘点一一对应:算子检测的边缘点与实际边缘点应该是一一对应。
Canny边缘检测的一般步骤:1.去噪 2.梯度运算 3.非极大值抑制 4.滞后阈值。下面我们就分别讲解一下这些操作:
1.去噪
边缘检测容易受到噪声的影响。因此,在进行边缘检测前,通常需要先进行去噪。我们一般使用高斯滤波进行图像去噪处理。也就是让临近的像素具有更高的重要度。对周围像素计算加权平均值,较近的像素具有较大的权重值。
在这里插入图片描述
比较像我们之间讲到的k近邻算法的意思哈。进行去噪处理之后呢,我们对于平滑图像进行一次Sobel算子梯度运算。
2.梯度运算
= 根号下(2 + 2)
对于方向: = −1(/)

在这里插入图片描述
对于每一个点计算之后的方向就是这样。大小和方向在图中都已经进行了标记。梯度的方向一般总是与边界垂直。梯度方向被归为四类:垂直,水平,和两个对角线。
3.非极大值抑制
在获得了梯度和方向后,遍历图像,去除所有不是边界的点。实现方法:逐个遍历像素点,判断当前像素点是否是周围像素点中具有相同方向梯度的最大值。
在这里插入图片描述
点A、点B、点C三点具有相同的方向,梯度方向垂直于边缘。判断点A是否为点A、点B、点C的局部最大值 :如果是,保留该点;否则,它被抑制(归零)。
在这里插入图片描述
这里就表示在垂直方向上进行了抑制,取出最大值。
4.滞后阈值
这里对非极大值抑制之后,我们进行了一次滞后阈值操作,分别取一个最小阈值和一个最大阈值,然后图像像素点的曲线与最大最小阈值相连则保留,如果不相连则抛弃。
在这里插入图片描述
我们来看一下边缘检测的结果图:
在这里插入图片描述

️7.Canny边缘检测的函数及使用

Canny边缘检测的函数是:
edges = cv2.Canny( image, threshold1, threshold2 )
threshold1表示阈值1,threshold2表示阈值2。也就是我们刚刚介绍到的最大最小阈值。

import cv2
import numpy as np
o=cv2.imread("image\\canny.bmp",cv2.IMREAD_GRAYSCALE)
r=cv2.Canny(o,100,200)
cv2.imshow("original",o)
cv2.imshow("result",r)
cv2.waitKey()
cv2.destroyAllWindows()

在这里插入图片描述
然后我们在对lena做一次,结果是:
在这里插入图片描述
调整阈值,可以控制边缘细节。maxVal、minVal变小,有更多的细节信息。
在这里插入图片描述
在这里插入图片描述

支持如果觉得博主的文章还不错或者您用得到的话,可以免费的关注一下博主,如果三连收藏支持就更好啦!这就是给予我最大的支持!

版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。
本文链接:https://blog.csdn.net/m0_37623374/article/details/125267392

智能推荐

android kt框架,GitHub - bitxiao/KtArmor-MVVM: Android快速开发框架, KtArmor 寓意着 为Android 赋予战斗装甲, 方便开发者快速进行An...-程序员宅基地

文章浏览阅读479次。前言学习了Kotlin有一段时间了, 每次写项目/Demo的时候, 总是用到网络请求、MVP、MVVM、常用工具类、通用自定义View, 索性把这些整合到一起, 搭成一个Android的脚手架——KtArmor. 框架是我个人经验的积累, 总结. 如有不妥, 望各位大佬指出.什么是KtArmor ?KtArmor 寓意着 为Android 赋予战斗装甲, 方便开发者快速进行Android 开发。节..._android kt

linux共享动态库中同名对象重复析构-两次析构或多次析构引起的double free解决办法_全局变量跨dll 调用导致析构double free-程序员宅基地

文章浏览阅读2.1k次。原文链接:http://chengxu.org/p/541.htmlLinux 平台下的共享动态库,一般都是通过选项“-fPIC”编译出来。有些应用程序需要链接多个共享库,此时如果在这些共享库中存在相同作用域范围的同名静态成员变量,那么当程序访问完静态成员变量结束析构时,由于内存的 double free 会导致程序 core dump;该问题是由于 Linux 编译器的缺陷造成的,本文就此问题进..._全局变量跨dll 调用导致析构double free

Golang 常见设计模式之装饰模式_goland 装饰模式-程序员宅基地

文章浏览阅读2.4k次,点赞4次,收藏3次。尽管 Go 语言中装饰模式没有 Python 中应用的那么广泛,但是它也有其独到的地方。接下来就一起看下装饰模式在 Go 语言中的应用。_goland 装饰模式

Orthogonal Convolutional Neural Networks_orthognal convolution neural-程序员宅基地

文章浏览阅读1.2k次,点赞2次,收藏2次。文章目录概主要内容符号说明Y=Conv(K,X)Y=Conv(K,X)Y=Conv(K,X)的俩种表示Y=KX~Y=K\tilde{X}Y=KX~Y=KXY=\mathcal{K}XY=KXkernel orthogonal regularizationorthogonal convolutionWang J, Chen Y, Chakraborty R, et al. Orthogonal ..._orthognal convolution neural

Windows7 Server 2008 下安装Oracle 10g提示“程序异常终止,发生未知错误”的解决方法_oracle10程序异常终止一 windows 2008-程序员宅基地

文章浏览阅读7.3k次。我的Oracle 10g版本是10.2.0.1.0,选择高级安装,提示“程序异常终止,发生未知错误”。1.修改Oracle 10G\database\stage\prereq\db\refhost.xml在 后面添加 2.到install目录中找到oraparam.ini文件,把#Windows=4.0,5.0,5.1,5.2修_oracle10程序异常终止一 windows 2008

commonjs与es6 中的模块化机制有什么区别_seajs,es.cjs-程序员宅基地

文章浏览阅读2.1k次。CommonJs导出的是变量的一份拷贝,ES6 Module导出的是变量的绑定(export default 是特殊的)cjs中,对象的传递如同a = b 一样,基本类型值传递,对象则引用传递。而es6 module中,如果不用export default去引用,基本类型也是引用传递的。export default的赋值则与a = b 相同。CommonJs是单个值导出,ES6 Module可..._seajs,es.cjs

随便推点

内存管理之程序内存分布-程序员宅基地

文章浏览阅读377次。在多任务操作系统中的每一个进程都运行在一个属于它自己的内存沙盘中。这个沙盘就是虚拟地址空间(virtual address space)。1 32位虚拟内存布局在32位模式下虚拟地址空间总是一个4GB的内存地址块。这些虚拟地址通过页表(page table)映射到物理内存,页表由操作系统维护并被处理器引用。每一个进程拥有一套属于它自己的页表,但是还有一个隐情。只要虚拟地址被使用,那么...

Bootstrap3 全局样式_bootstrap3 权限选择样式-程序员宅基地

文章浏览阅读423次。Bootstrap 3 的目标是简洁、直观、强悍的前端开发框架,让 Web 开发变得更好、更快、更强壮,我们有必要先了解一下 Bootstrap 底层结构的关键部分。 HTML5文档类型Bootstrap使用了HTML5特定的元素和CSS属性,在使用Bootstrap的时候,所有HTML文档的第一行代码必须是 <!DOCTYPE html>。如:<!DOCTYPE h..._bootstrap3 权限选择样式

linux管道举例理解_linux中的管道举例-程序员宅基地

文章浏览阅读4.7k次,点赞24次,收藏89次。linux管道举例理解一、管道的定义:“|”二、查找2.1统计当前目录下有多少个文件2.2查看当前目录下的前n(3)个文件2.3查看wang.txt文件包含i的字符行2.4查看内存使用情况2.5查询进程三、更改一、管道的定义:“|”一般我们在进行操作的时候,命令很多,但我们只想要其中一部分,那么就可以使用管道了。管道是Linux中很重要的一种通信方式,是把一个前一个结果的输出直接连接到另一个..._linux中的管道举例

hdu 1686 kmp Oulipo-程序员宅基地

文章浏览阅读499次。http://acm.hdu.edu.cn/showproblem.php?pid=1686对于我这种刚开始学的人来说 刚开始我觉的很难 不会做哦 但是看了一下别人的一下子就知道了 还是不是太熟了 需要加强练习哦 呵呵 其实这是一道简单的kmp 的算法 只要你在index_kmp() 将 i 的值一直小于N 就可以了 还是看看我具体

实战篇:Linux 安装 Oracle 11GR2 数据库保姆级教程_liunxoracle11gr2-程序员宅基地

文章浏览阅读2.1w次,点赞204次,收藏603次。没接触Linux的朋友不用害怕,跟着本篇文章一步步操作,安装Oracle如喝水般简单且标准。_liunxoracle11gr2

CSS实现文字折行缩进效果_css 折行的缩进-程序员宅基地

文章浏览阅读1.6k次。移动端开发,经常会实现文字缩进效果。废话不多说,看图例展示。_css 折行的缩进

推荐文章

热门文章

相关标签