OpenCV-Python官方教程-32-计算摄影学图像去噪

1. 目标：

• 学习使用非局部平均值去噪算法去除图像中的噪音
• 学习函数 cv2.fastNlMeansDenoising()，cv2.fastNlMeansDenoisingColored()等

2. 原理

我们已经学习了很多图像平滑技术，比如高斯平滑，中值平滑等，当噪声比较小时这些技术的效果都是很好的。在这些技术中我们选取像素周围一个小的邻域然后用高斯平均值或者中值平均值取代中心像素。简单来说，像素级别的噪声去除是限制在局部邻域的。

噪声有一个性质。我们认为噪声是平均值为一的随机变量。考虑一个带噪声的像素点，p = p0 + n，其中 p0 为像素的真实值，n 为这个像素的噪声。我们可以从不同图片中选取大量的相同像素（N）然后计算平均值。理想情况下我们会得到 p = p0。因为噪声的平均值为 0。

通过简单的设置我们就可以去除这些噪声。将一个静态摄像头固定在一个位置连续拍摄几秒钟。这样我们就会得到足够多的图像帧，或者同一场景的大量图像。写一段代码求解这些帧的平均值。将最终结果与第一帧图像对比一下。你会发现噪声减小了。不幸的是这种简单的方法对于摄像头和运动场景并不总是适用。大多数情况下我们只有一张带有噪音的图像。

想法很简单，我们需要一组相似的图片，通过取平均值的方法可以去除噪音。考虑图像中一个小的窗口（5x5），有很大可能图像中的其他区域也存在一个相似的窗口。有时这个相似窗口就在邻域周围。如果我们找到这些相似的窗口并取他们的平均值会怎样呢？对于特定的窗口这样做挺好的。
对于特定的窗口这样做挺好的。如下图所示。

上图中的蓝色窗口看起来是相似的。绿色窗口看起来也是相似的。所以我们可以选取包含目标像素的一个小窗口，然后在图像中搜索相似的窗口，最后求取所有窗口的平均值，并用这个值取代目标像素的值。这种方法就是非局部平均值去噪。与我们以前学习的平滑技术相比这种算法要消耗更多的时间，但是结果很好。你可以在更多资源中找到更多的细节和在线演示。

对于彩色图像，要先转换到CIELAB颜色空间，然后对L和AB成分分别去噪。

3.1 OpenCV中的图像去噪

OpenCV 提供了这种技术的四个版本。

• cv2.fastNlMeansDenoising() 使用对象为灰度图。
• cv2.fastNlMeansDenoisingColored() 使用对象为彩色图。
• cv2.fastNlMeansDenoisingMulti() 适用于短时间的图像序列（灰度图像）。
• cv2.fastNlMeansDenoisingColoredMulti() 适用于短时间的图像序列（彩色图像）。

共同的参数有：

• h : 决定过滤器强度。h 值高可以很好的去除噪声，但也会把图像的细节抹去。(取 10 的效果不错)
• hForColorComponents : 与 h 相同，但使用与彩色图像。（与h相同）
• templateWindowSize : 奇数。(推荐值为 7)
• searchWindowSize : 奇数。(推荐值为 21)

3.2 cv2.fastNlMeansDenoisingColored()

和上面提到的一样，它可以被用来去除彩色图像的噪声。（假设是高斯噪声）。代码示例：

import cv2
from matplotlib import pyplot as plt

img = cv2.imread('cell.jpg')
dst = cv2.fastNlMeansDenoisingColored(img,None,10,10,7,21)

plt.subplot(121),plt.imshow(img)
plt.subplot(122),plt.imshow(dst)
plt.show()