opencv图像平滑

学习使用不同的低通滤波器对图像进行模糊
使用自定义的滤波器对图像进行卷积（2D 卷积）

2D 卷积

  对 2D 图像实施低通滤波（LPF），高通滤波（HPF）, LPF去除噪音，模糊图像,HPF找到图像的边缘,锐化图像；
  OpenCV 提供的函数 cv.filter2D() 可以让我们对一幅图像进行卷积操作。下面我们将对一幅图像使用平均滤波器。下面是一个 5x5 的平均滤波器核：

操作如下：将核放在图像的一个像素 A 上，求与核对应的图像上 25（5x5）个像素的和，在取平均数，用这个平均数替代像素 A 的值。重复以上操作直到将图像的每一个像素值都更新一边。代码如下：

import cv2
import numpy as np 
from matplotlib import pyplot as plt

img = cv2.imread('image/lufei.jpeg')
kernel = np.ones((5,5),np.float32)/25
dst = cv2.filter2D(img, -1, kernel)

plt.subplot(121),plt.imshow(img),plt.title('Original')
plt.xticks([]),plt.yticks([]) 
plt.subplot(122),plt.imshow(dst),plt.title('Averaging')
plt.xticks([]),plt.yticks([]) 
plt.show()

结果图：

图像模糊（图像平滑）

  使用低通滤波器可以达到图像模糊的目的。这对与去除噪音很有帮助。其实就是去除图像中的高频成分（比如：噪音，边界）。所以边界也会被模糊一点。（当然，也有一些模糊技术不会模糊掉边界）。 OpenCV 提供了四种模糊技术。

1、均值模糊

  由一个归一化卷积框完成的，用卷积框覆盖区域所有像素的平均值来代替中心元素，使用函数 cv2.blur() 和 cv2.boxFilter() 来完这个任务，需要设定卷积框的宽和高。下面是一个 3x3 的归一化卷积框：

注：如果不想使用归一化卷积框，你应该使用 cv2.boxFilter()，这时要传入参数 normalize=False。

import cv2
import numpy as np 
from matplotlib import pyplot as plt

img = cv2.imread('image/lufei.jpeg')
blur = cv2.blur(img, (5,5))

plt.subplot(121),plt.imshow(img),plt.title('Original')
plt.xticks([]),plt.yticks([]) 
plt.subplot(122),plt.imshow(blur),plt.title('Blurred')
plt.xticks([]),plt.yticks([]) 
plt.show()

结果图：

2 、高斯模糊

  现在把卷积核换成高斯核（简单来说，方框不变，将原来每个方框的值是相等的，现在里面的值是符合高斯分布的，方框中心的值最大，其余方框根据距离中心元素的距离递减，构成一个高斯小山包。原来的求平均数现在变成求加权平均数，全就是方框里的值）。实现的函数是cv2.GaussianBlur()。需要指定高斯核的宽和高（必须是奇数）。以及高斯函数沿 X， Y 方向的标准差。如果只指定了 X 方向的的标准差， Y 方向也会取相同值。如果两个标准差都是 0，那么函数会根据核函数的大小自己计算。高斯滤波可以有效的从图像中去除高斯噪音。
也可以使用函数 cv2.getGaussianKernel() 自己构建一个高斯核。

import cv2
import numpy as np 
from matplotlib import pyplot as plt

img = cv2.imread('image/lufei.jpeg')
#添加噪声点
for i in range(2000): 
    temp_x = np.random.randint(0,img.shape[0])
    temp_y = np.random.randint(0,img.shape[1])
    img[temp_x][temp_y] = 255
blur = cv2.GaussianBlur(img, (5,5),0)
madian = cv2.medianBlur(img, 5)
plt.subplot(121),plt.imshow(img),plt.title('Original')
plt.xticks([]),plt.yticks([]) 
plt.subplot(122),plt.imshow(blur),plt.title('GaussianBlur')
plt.xticks([]),plt.yticks([]) 
plt.show()

结果图：

3 、中值模糊

  顾名思义就是用与卷积框对应像素的中值来替代中心像素的值。这个滤波器经常用来去除椒盐噪声。前面的滤波器都是用计算得到的一个新值来取代中心像素的值，而中值滤波是用中心像素周围（也可以使他本身）的值来取代他,能有效的去除噪声。卷积核的大小也应该是一个奇数。
  给原始图像加上噪声然后再使用中值模糊。
代码：

import cv2
import numpy as np 
from matplotlib import pyplot as plt

img = cv2.imread('image/lufei.jpeg')
for i in range(2000): 
   temp_x = np.random.randint(0,img.shape[0])
   temp_y = np.random.randint(0,img.shape[1])
   img[temp_x][temp_y] = 255
madian = cv2.medianBlur(img, 5)
plt.subplot(121),plt.imshow(img),plt.title('Original')
plt.xticks([]),plt.yticks([]) 
plt.subplot(122),plt.imshow(madian),plt.title('Median')
plt.xticks([]),plt.yticks([]) 
plt.show()

结果图：

4 、双边滤波

  函数 cv2.bilateralFilter() 能在保持边界清晰的情况下有效的去除噪音,但是这种操作与其他滤波器相比会比较慢。我们已经知道高斯滤波器是求中心点邻近区域像素的高斯加权平均值,但高斯滤波器只考虑像素之间的空间关系，而不会考虑像素值之间的关系（像素的相似度），所以这种方法不会考虑一个像素是否位于边界，因此边界也会别模糊掉，而这正不是我们想要。
  双边滤波在同时使用空间高斯权重和灰度值相似性高斯权重。空间高斯函数确保只有邻近区域的像素对中心点有影响，灰度值相似性高斯函数确保只有与中心像素灰度值相近的才会被用来做模糊运算，所以这种方法会确保边界不会被模糊掉，因为边界处的灰度值变化比较大。

import cv2
import numpy as np 
from matplotlib import pyplot as plt

img = cv2.imread('image/1.png')

blur = cv2.bilateralFilter(img, 9, 75, 75)
plt.subplot(121),plt.imshow(img),plt.title('Original')
plt.xticks([]),plt.yticks([]) 
plt.subplot(122),plt.imshow(blur),plt.title('BilateralFilter')
plt.xticks([]),plt.yticks([]) 
plt.show()

结果图：

纹理被模糊掉了，但是边界还在。
注：双边滤波无法去除噪声：

import cv2
import numpy as np 
from matplotlib import pyplot as plt

img = cv2.imread('image/1.png')
for i in range(2000): 
   temp_x = np.random.randint(0,img.shape[0])
   temp_y = np.random.randint(0,img.shape[1])
   img[temp_x][temp_y] = 255
blur = cv2.bilateralFilter(img, 9, 75, 75)
plt.subplot(121),plt.imshow(img),plt.title('Original')
plt.xticks([]),plt.yticks([]) 
plt.subplot(122),plt.imshow(blur),plt.title('BilateralFilter')
plt.xticks([]),plt.yticks([]) 
plt.show()

结果图：

噪声点依然在