写点什么

跟我学 Python 图像处理丨掌握 4 种图像平滑算法

  • 2022 年 5 月 09 日
  • 本文字数:3919 字

    阅读完需:约 13 分钟

本文分享自华为云社区《[Python图像处理] 三.获取图像属性、兴趣ROI区域及通道处理 | 【生长吧!Python】》,作者:eastmount。

 

本篇文章主要讲解 Python 调用 OpenCV 实现图像平滑,包括四个算法:均值滤波、方框滤波、高斯滤波和中值滤波。全文均是基础知识,希望对您有所帮助。


一.图像平滑

1.图像增强

图像增强是对图像进行处理,使其比原始图像更适合于特定的应用,它需要与实际应用相结合。对于图像的某些特征如边缘、轮廓、对比度等,图像增强是进行强调或锐化,以便于显示、观察或进一步分析与处理。图像增强的方法是因应用不同而不同的,研究内容包括:(参考课件和左飞的《数字图像处理》)

2.图像平滑

图像平滑是一种区域增强的算法,平滑算法有邻域平均法、中指滤波、边界保持类滤波等。在图像产生、传输和复制过程中,常常会因为多方面原因而被噪声干扰或出现数据丢失,降低了图像的质量(某一像素,如果它与周围像素点相比有明显的不同,则该点被噪声所感染)。这就需要对图像进行一定的增强处理以减小这些缺陷带来的影响。


简单平滑-邻域平均法

3.邻域平均法

图像简单平滑是指通过邻域简单平均对图像进行平滑处理的方法,用这种方法在一定程度上消除原始图像中的噪声、降低原始图像对比度的作用。它利用卷积运算对图像邻域的像素灰度进行平均,从而达到减小图像中噪声影响、降低图像对比度的目的。


但邻域平均值主要缺点是在降低噪声的同时使图像变得模糊,特别在边缘和细节处,而且邻域越大,在去噪能力增强的同时模糊程度越严重。

首先给出为图像增加噪声的代码。

# -*- coding:utf-8 -*-import cv2import numpy as np
#读取图片img = cv2.imread("test.jpg", cv2.IMREAD_UNCHANGED)rows, cols, chn = img.shape
#加噪声for i in range(5000): x = np.random.randint(0, rows) y = np.random.randint(0, cols) img[x,y,:] = 255
cv2.imshow("noise", img) #等待显示cv2.waitKey(0)cv2.destroyAllWindows()
复制代码

输出结果如下所示:

二.均值滤波

1.原理

均值滤波是指任意一点的像素值,都是周围 N*M 个像素值的均值。例如下图中,红色点的像素值为蓝色背景区域像素值之和除 25。

其中红色区域的像素值均值滤波处理过程为: ((197+25+106+156+159)+ (149+40+107+5+71)+ (163+198+**226**+223+156)+ (222+37+68+193+157)+ (42+72+250+41+75)) / 25

其中 5*5 的矩阵称为核,针对原始图像内的像素点,采用核进行处理,得到结果图像。

提取 1/25 可以将核转换为如下形式:

2.代码

Python 调用 OpenCV 实现均值滤波的核心函数如下:

result = cv2.blur(原始图像,核大小)

其中,核大小是以(宽度,高度)表示的元祖形式。常见的形式包括:核大小(3,3)和(5,5)。

代码如下所示:

#encoding:utf-8import cv2  import numpy as np  import matplotlib.pyplot as plt#读取图片img = cv2.imread('test01.png')source = cv2.cvtColor(img,cv2.COLOR_BGR2RGB)#均值滤波result = cv2.blur(source, (5,5))#显示图形titles = ['Source Image', 'Blur Image']  images = [source, result]  for i in xrange(2):     plt.subplot(1,2,i+1), plt.imshow(images[i], 'gray')     plt.title(titles[i])     plt.xticks([]),plt.yticks([])  plt.show()  
复制代码

输出结果如下图所示:

核设置为(10,10)和(20,20)会让图像变得更加模糊。

如果设置为(1,1)处理结果就是原图,核中每个权重值相同,称为均值。

三.方框滤波

方框滤波和均值滤波核基本一致,区别是需不需要均一化处理。OpenCV 调用 boxFilter()函数实现方框滤波。函数如下:

result = cv2.boxFilter(原始图像, 目标图像深度, 核大小, normalize 属性)

其中,目标图像深度是 int 类型,通常用“-1”表示与原始图像一致;核大小主要包括(3,3)和(5,5),如下所示。

normalize 属性表示是否对目标图像进行归一化处理。当 normalize 为 true 时需要执行均值化处理,当 normalize 为 false 时,不进行均值化处理,实际上为求周围各像素的和,很容易发生溢出,溢出时均为白色,对应像素值为 255。

在图像简单平滑中,算法利用卷积模板逐一处理图像中每个像素,这一过程可以形象地比作对原始图像的像素一一进行过滤整理,在图像处理中把邻域像素逐一处理的算法过程称为滤波器。平滑线性滤波器的工作原理是利用模板对邻域内像素灰度进行加权平均,也称为均值滤波器。

代码如下所示:

#encoding:utf-8import cv2  import numpy as np  import matplotlib.pyplot as plt#读取图片img = cv2.imread('test01.png')source = cv2.cvtColor(img,cv2.COLOR_BGR2RGB)#方框滤波result = cv2.boxFilter(source, -1, (5,5), normalize=1)#显示图形titles = ['Source Image', 'BoxFilter Image']  images = [source, result]  for i in xrange(2):     plt.subplot(1,2,i+1), plt.imshow(images[i], 'gray')     plt.title(titles[i])     plt.xticks([]),plt.yticks([])  plt.show()  
复制代码

代码中使用 5*5 的核,normalize=1 表示进行归一化处理,此时与均值滤波相同,输出结果如下图所示:

下面是图像左上角处理前后的像素结果:

print(source[0:3, 0:3, 0])#[[115 180 106]# [ 83 152  72]# [ 55  58  55]]print(result[0:3, 0:3, 0])#[[92 90 78]# [92 89 77]# [82 80 72]]
复制代码

如果省略参数 normalize,则默认是进行归一化处理。如果 normalize=0 则不进行归一化处理,像素值为周围像素之和,图像更多为白色。

#encoding:utf-8import cv2  import numpy as np  import matplotlib.pyplot as plt#读取图片img = cv2.imread('test01.png')source = cv2.cvtColor(img,cv2.COLOR_BGR2RGB)#方框滤波result = cv2.boxFilter(source, -1, (5,5), normalize=0)#显示图形titles = ['Source Image', 'BoxFilter Image']  images = [source, result]  for i in xrange(2):     plt.subplot(1,2,i+1), plt.imshow(images[i], 'gray')     plt.title(titles[i])     plt.xticks([]),plt.yticks([])  plt.show() 
复制代码

输出结果如下图所示:

上图很多像素为白色,因为图像求和结果几乎都是 255。如果设置的是 2*2 矩阵,只取四个像素结果要好些。result = cv2.boxFilter(source, -1, (2,2), normalize=0)

四.高斯滤波

为了克服简单局部平均法的弊端(图像模糊),目前已提出许多保持边缘、细节的局部平滑算法。它们的出发点都集中在如何选择邻域的大小、形状和方向、参数加平均及邻域各店的权重系数等。

图像高斯平滑也是邻域平均的思想对图像进行平滑的一种方法,在图像高斯平滑中,对图像进行平均时,不同位置的像素被赋予了不同的权重。高斯平滑与简单平滑不同,它在对邻域内像素进行平均时,给予不同位置的像素不同的权值,下图的所示的 3 * 3 和 5 * 5 领域的高斯模板。

高斯滤波让临近的像素具有更高的重要度,对周围像素计算加权平均值,较近的像素具有较大的权重值。如下图所示,中心位置权重最高为 0.4。


Python 中 OpenCV 主要调用 GaussianBlur 函数,如下:

dst = cv2.GaussianBlur(src, ksize, sigmaX)

其中,src 表示原始图像,ksize 表示核大小,sigmaX 表示 X 方向方差。注意,核大小(N, N)必须是奇数,X 方向方差主要控制权重。

代码如下:

#encoding:utf-8import cv2  import numpy as np  import matplotlib.pyplot as plt#读取图片img = cv2.imread('test01.png')source = cv2.cvtColor(img,cv2.COLOR_BGR2RGB)#高斯滤波result = cv2.GaussianBlur(source, (3,3), 0)
#显示图形titles = ['Source Image', 'GaussianBlur Image'] images = [source, result] for i in xrange(2): plt.subplot(1,2,i+1), plt.imshow(images[i], 'gray') plt.title(titles[i]) plt.xticks([]),plt.yticks([]) plt.show()
复制代码

输出结果如下所示:

如果使用 15*15 的核,则图形将更加模糊。

五.中值滤波

1.概念

在使用邻域平均法去噪的同时也使得边界变得模糊。而中值滤波是非线性的图像处理方法,在去噪的同时可以兼顾到边界信息的保留。选一个含有奇数点的窗口 W,将这个窗口在图像上扫描,把窗口中所含的像素点按灰度级的升或降序排列,取位于中间的灰度值来代替该点的灰度值。 例如选择滤波的窗口如下图,是一个一维的窗口,待处理像素的灰度取这个模板中灰度的中值,滤波过程如下:

如下图所示,将临近像素按照大小排列,取排序像素中位于中间位置的值作为中值滤波的像素值。

2.代码

OpenCV 主要调用 medianBlur()函数实现中值滤波。图像平滑里中值滤波的效果最好。

dst = cv2.medianBlur(src, ksize)

其中,src 表示源文件,ksize 表示核大小。核必须是大于 1 的奇数,如 3、5、7 等。

代码如下所示:

#encoding:utf-8import cv2  import numpy as np  import matplotlib.pyplot as plt#读取图片img = cv2.imread('test01.png')#高斯滤波result = cv2.medianBlur(img, 3)
#显示图像cv2.imshow("source img", img)cv2.imshow("medianBlur", result)
#等待显示cv2.waitKey(0)cv2.destroyAllWindows()
复制代码

输出结果如下图所示:

常用的窗口还有方形、十字形、圆形和环形。不同形状的窗口产生不同的滤波效果,方形和圆形窗口适合外轮廓线较长的物体图像,而十字形窗口对有尖顶角状的图像效果好。中值滤波对于消除孤立点和线段的干扰十分有用,尤其是对于二进噪声,但对消除高斯噪声的影响效果不佳。对于一些细节较多的复杂图像,可以多次使用不同的中值滤波。

希望文章对大家有所帮助,如果有错误或不足之处,还请海涵。


本文摘录自 eastmount X 华为云开发者社区联合出品的电子书《从零到一 • Python 图像处理及识别》。

点击免费下载电子书《从零到一 • Python图像处理及识别


点击关注,第一时间了解华为云新鲜技术~​

发布于: 2022 年 05 月 09 日阅读数: 21
用户头像

提供全面深入的云计算技术干货 2020.07.14 加入

华为云开发者社区,提供全面深入的云计算前景分析、丰富的技术干货、程序样例,分享华为云前沿资讯动态,方便开发者快速成长与发展,欢迎提问、互动,多方位了解云计算! 传送门:https://bbs.huaweicloud.com/

评论

发布
暂无评论
跟我学Python图像处理丨掌握4种图像平滑算法_Python_华为云开发者社区_InfoQ写作社区