[Python图像处理]四.图像平滑中四种常用的滤波

主要讲解Python调用OpenCV实现图像平滑，包括四个算法：均值滤波、方框滤波、高斯滤波和中值滤波.

给图像增加噪声:

import cv2
import numpy as np

def test10():
    img = cv2.imread("result.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("demo", img)
    if cv2.waitKey(0) == 27:
        cv2.destroyWindow("demo")

test10()

效果如下:

均值滤波: 均值滤波是指任意一点的像素值，都是周围N*M个像素值的均值, result = cv2.blur(图像, 核大小)，其中核大小是以(宽度, 高度)表示的元组形式,常见的形式包括: (3, 3)和(5, 5)

代码如下:

import matplotlib.pyplot as plt
import cv2
def test11():
    img = cv2.imread("result.jpg")
    source = cv2.cvtColor(img, cv2.COLOR_BGR2BGRA)
    # 均值滤波
    result = cv2.blur(source, (5, 5))

    # 显示图像
    title = ["demo1", "demo2"]
    imgs = [source, result]
    for i in range(2):
        plt.subplot(1, 2, i+1)　　　　 plt.imshow(imgs[i], "gray")
        plt.title(title[i])
        plt.xticks([])
        plt.yticks([])
    plt.show()

test11()

效果如下:

 将核值调大会让图像变模糊，例如设置为(20, 20),效果如下:如果设置为(1, 1)处理的结果就是原图，核中每个权重值相同,称为均值。

 方框滤波

方框滤波和均值滤波核基本一致，区别是需不需要均一化处理。OpenCV调用boxFilter()函数实现方框滤波。函数如下：
result = cv2.boxFilter(原始图像, 目标图像深度, 核大小, normalize属性)  其中，目标图像深度是int类型，通常用“-1”表示与原始图像一致；核大小主要包括（3，3）和（5，5）;normalize属性表示是否对目标图像进行归一化处理。当normalize为true时需要执行均值化处理，当normalize为false时，不进行均值化处理，实际上为求周围各像素的和，很容易发生溢出，溢出时均为白色，对应像素值为255。

代码如下:

import cv2
import matplotlib.pyplot as plt
def test12():
    img = cv2.imread("result.jpg")
    source = cv2.cvtColor(img, cv2.COLOR_BGR2BGRA)
    # 方框滤波
    result = cv2.boxFilter(source, -2, (5, 5), normalize=1)
    # 显示图像
    title = ["demo1", "demo2"]
    imgs = [source, result]
    for i in range(2):
        plt.subplot(1, 2, i+1), plt.imshow(imgs[i], "gray")
        plt.title(title[i])
        plt.xticks([])
        plt.yticks([])
    plt.show()

test12()

效果如下:

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

效果如下:

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

 高斯滤波： 高斯滤波让临近的像素具有更高的重要度，对周围像素计算加权平均值，较近的像素具有较大的权重值。dst = cv2.GaussianBlur(src, ksize, sigmaX) 其中，src表示原始图像，ksize表示核大小，sigmaX表示X方向方差。注意，核大小（N, N）必须是奇数，X方向方差主要控制权重

代码如下:

import cv2
import matplotlib.pyplot as plt
def test12():
    img = cv2.imread("result.jpg")
    source = cv2.cvtColor(img, cv2.COLOR_BGR2BGRA)
    # 高斯滤波
    result = cv2.GaussianBlur(source, (3, 3), 0)
    # 显示图像
    title = ["demo1", "demo2"]
    imgs = [source, result]
    for i in range(2):
        plt.subplot(1, 2, i+1), plt.imshow(imgs[i], "gray")
        plt.title(title[i])
        plt.xticks([])
        plt.yticks([])
    plt.show()

test12()

效果如下:

 核数增加，图像会变模糊，例如设置为(17, 17)

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

OpenCV主要调用medianBlur()函数实现中值滤波,图像平滑里中值滤波的效果最好。
dst = cv2.medianBlur(src, ksize)
其中，src表示源文件，ksize表示核大小。核必须是大于1的奇数，如3、5、7等

代码如下：

import cv2
def test13():
    img = cv2.imread("result.jpg")
    # 高斯滤波
    result = cv2.medianBlur(img, 5)
    cv2.imshow("demo1", img)
    cv2.imshow("demo2", result)
    if cv2.waitKey(0) == 27:
        cv2.destroyWindow("demo1")
        cv2.destroyWindow("demo2")
test13()

效果如下:

转自: https://blog.csdn.net/Eastmount/article/details/82216380