Python OpenCV形态学处理和图像梯度

形态转换

形态学变换是基于图像形状的一些简单操作。它通常在二值图像上执行。它需要两个输入,一个是我们的原始图像,第二个是决定操作性质的结构元素或内核。两个基本的形态学操作是腐蚀和膨胀。

腐蚀

腐蚀 的基本概念就像土壤侵蚀一样,只侵蚀前景对象的边界(总是尽量保持前景为白色)。

那它有什么作用呢?

内核在图像中滑动(如二维卷积)。只有当内核下的所有像素都为 1 时,原始图像中的像素(1 或 0)才会被视为 1,否则会被侵蚀(变为零)。

所以根据内核的大小,边界附近的所有像素都将被丢弃。因此,前景对象的厚度或大小在图像中减少或只是白色区域减少。它有助于消除小的白色噪音(如我们在"颜色空间"一章中所看到的),分离两个连接的对象等。

python 复制代码
import cv2
import numpy as np
 
# 读取二值图像
image = cv2.imread('binary_image.jpg', 0)
 
# 创建5x5的结构元素(内核)
kernel = np.ones((5, 5), np.uint8)
 
# 执行腐蚀操作
eroded_image = cv2.erode(image, kernel, iterations=1)
 
# 显示原始图像和腐蚀后的图像
cv2.imshow('Original Image', image)
cv2.imshow('Eroded Image', eroded_image)
cv2.waitKey(0)
cv2.destroyAllWindows()
膨胀

膨胀 与腐蚀正好相反。这里,如果内核下至少有一个像素为"1",则像素元素为"1"。所以它会增加图像中的白色区域,或者增加前景对象的大小。

通常情况下,在去除噪音的情况下,腐蚀后会膨胀。因为,腐蚀消除了白噪声,但它也缩小了我们的对象。所以我们扩大它。由于噪音消失了,它们不会再回来,但我们的目标区域会增加到腐蚀之前的状态。它还可用于连接对象的断开部分。

python 复制代码
# 执行膨胀操作
dilated_image = cv2.dilate(image, kernel, iterations=1)
 
# 显示膨胀后的图像
cv2.imshow('Dilated Image', dilated_image)
cv2.waitKey(0)
cv2.destroyAllWindows()
开运算--先腐蚀再膨胀

开运算一般会平滑物体的轮廓、断开较窄的狭颈并消除细的突出物。正如我们上面所解释的,它对消除噪音很有用。在这里,我们使用 cv.morphologyEx()

python 复制代码
# 执行开运算
opened_image = cv2.morphologyEx(image, cv2.MORPH_OPEN, kernel)
 
# 显示开运算后的图像
cv2.imshow('Opened Image', opened_image)
cv2.waitKey(0)
cv2.destroyAllWindows()
闭运算--先膨胀再腐蚀

闭运算同样也会平滑轮廓的一部分。但与开操作相反,它通常会弥合较窄的间断和细长的沟壑,消除小的孔洞,填补轮廓线中的断裂。

python 复制代码
# 执行闭运算
closed_image = cv2.morphologyEx(image, cv2.MORPH_CLOSE, kernel)
 
# 显示闭运算后的图像
cv2.imshow('Closed Image', closed_image)
cv2.waitKey(0)
cv2.destroyAllWindows()
顶帽运算--原始减去开运算
黑帽运算--闭运算减去原始
  • 腐蚀:cv2.erode(src, kernel, iterations) #iterations表示腐蚀几次

  • 膨胀:cv2.dilate(src, kernel, iterations)

  • 开运算:cv2.morphologyEx(src, cv2.MORPH_OPEN, kernel)

  • 闭运算:cv2.morphologyEx(src, cv2.MORPH_CLOSE, kernel)

图像梯度

OpenCv 提供三种类型的梯度滤波器或高通滤波器,Sobel、Scharr 和 Laplacian。

sobel算子

Sobel 算子 是一种边缘检测算子,它结合了高斯平滑与微分运算,从而提高了对噪声的抵抗力。Sobel 算子通过两个方向(水平和垂直)来计算图像的梯度。



计算图像的总梯度(两个方向上)



  • dst = cv2.Sobel(src, ddepth, dx, dy, [ksize])
    • dst:计算结果
    • src:原始图像
    • ddepth :处理结果图像深度
    • dx:x轴方向
    • dy:y轴方向
    • ksize:核大小

ddepth参数说明:

  • 一般,通常直接设置为-1,表示让处理结果与原始图像保持一致
  • 但是此处需要注意,不能直接置为-1

dx dy参数说明:

  • 计算x方向梯度: dx=1,dy=0

  • 计算y方向梯度:dx=0,dy=1

opencv提供了一个计算两幅图像的权重和

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

五个参数分别表示:原图像1,原图像1的系数,原图像2,原图像2的系数,修正值

示例



scharr算子

Scharr 算子 是 Sobel 算子的一个变种,具有更好的旋转不变性和更高的梯度响应。Scharr 算子在某些情况下能够提供更精确的边缘检测,尤其是细节较多的图像。

dst = Scharr(src, ddpeth, dx, dy)

  • src:原始图像
  • ddepth:处理结果图像深度
  • dx:x轴方向
  • dy:y轴方向

ddepth参数:一般在别的图像处理函数中,直接置为-1,表示处理成与原图像一样深度的图像,但是由于此函数处理过程中会出现负数,故与sobel算子一样,将它赋值为cv2.CV_64F

然后再转为unit8正数

python 复制代码
dst = Scharr(src, cv2.CV_64F, dx, dy)

dst = cv2.convertScaleAbs(dst)

sobel和 scharr算子大小是一样的,说明计算起来的工作量是一样的;但是scharr精确度更高,给了临近像素更高的权重


python 复制代码
import cv2
import numpy as np
 
# 读取图像并转换为灰度图
image = cv2.imread('image.jpg', cv2.IMREAD_GRAYSCALE)
 
# Sobel 微分
sobel_x = cv2.Sobel(image, cv2.CV_64F, 1, 0, ksize=3)  # X方向
sobel_y = cv2.Sobel(image, cv2.CV_64F, 0, 1, ksize=3)  # Y方向
sobel_magnitude = cv2.magnitude(sobel_x, sobel_y)  # 计算梯度幅值
 
# Scharr 微分
scharr_x = cv2.Scharr(image, cv2.CV_64F, 1, 0)  # X方向
scharr_y = cv2.Scharr(image, cv2.CV_64F, 0, 1)  # Y方向
scharr_magnitude = cv2.magnitude(scharr_x, scharr_y)  # 计算梯度幅值
 
# 显示结果
cv2.imshow('Original Image', image)
cv2.imshow('Sobel Magnitude', sobel_magnitude)
cv2.imshow('Scharr Magnitude', scharr_magnitude)
cv2.waitKey(0)
cv2.destroyAllWindows()
laplacian算子

拉普拉斯算子类似于二阶sobel导数。实际上,在OpenCV中通过调用sobel算子来计算拉普拉斯算子。

公式:

卷积核:

示例:

函数实现:

  • dst = cv2.Laplacian(src, ddepth)

注意:通常情况下,可以将图像深度ddepth的值设置为-1,让处理结果与原始图像保持一致

但是该函数在计算过程中也会出现负数

python 复制代码
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.destoryAllWindows()
相关推荐
小鸡吃米…16 小时前
Python编程语言面试问题二
开发语言·python·面试
eve杭16 小时前
AI、大数据与智能时代:从理论基石到实战路径
人工智能·python·5g·网络安全·ai
Honmaple17 小时前
中国四级城市联动数据,包含港澳台,内含json , sql , python 脚本
python·sql·json
BoBoZz1917 小时前
Curvatures 曲率的计算、边缘曲率的调整以及曲率、颜色的映射
python·vtk·图形渲染·图形处理
少吃零食多运动17 小时前
【Jupyter notebook修改工作目录】
python·jupyter
Swizard18 小时前
别买树莓派了!3步教你在安卓手机上跑通 CPython + PaddleOCR,打造随身 AI 识别终端
python·ai·移动开发
weixin_4215850119 小时前
PYTHON 迭代器1 - PEP-255
开发语言·python
L.fountain19 小时前
图像自回归生成(Auto-regressive image generation)实战学习(一)
人工智能·深度学习·学习·计算机视觉·图像自回归
phoenix@Capricornus19 小时前
气泡自动计数——数字图像处理设计题
计算机视觉
技术净胜19 小时前
MATLAB进行图像分割从基础阈值到高级分割
opencv·计算机视觉·matlab