25 图像去噪声
代码
import cv2 as cv
import cv2 as cv
import numpy as np
def add_salt_pepper_noise(image):
h, w = image.shape[:2]
nums = 10000
rows = np.random.randint(0, h, nums, dtype=np.int)
cols = np.random.randint(0, w, nums, dtype=np.int)
for i in range(nums):
if i % 2 == 1:
image[rows[i], cols[i]] = (255, 255, 255)
else:
image[rows[i], cols[i]] = (0, 0, 0)
return image
def gaussian_noise(image):
noise = np.zeros(image.shape, image.dtype)
m = (15, 15, 15)
s = (30, 30, 30)
cv.randn(noise, m, s)
dst = cv.add(image, noise)
cv.imshow("gaussian noise", dst)
return dst
src = cv.imread("../images/example.png")
cv.imshow("input", src)
h, w = src.shape[:2]
src = gaussian_noise(src)
result1 = cv.blur(src, (5, 5))
cv.imshow("result-1-blur", result1)
result2 = cv.GaussianBlur(src, (5, 5), 0)
cv.imshow("result-2-GaussianBlur", result2)
result3 = cv.medianBlur(src, 5)
cv.imshow("result-3-medianBlur", result3)
result4 = cv.fastNlMeansDenoisingColored(src, None, 15, 15, 10, 30)
cv.imshow("result-4-fastNlMeansDenoisingColored", result4)
cv.waitKey(0)
cv.destroyAllWindows()
实验结果
解释
将图片添加了高斯噪声,并且比较了四种常见的滤波方法的效果
- result1:均值滤波
cv.blur
- result2:高斯滤波
cv.GaussianBlur
- result3:中值滤波
cv.medianBlur
- result4:非局部均值滤波
cv.fastNlMeansDenoisingColored
所有内容均来源于贾志刚老师的知识星球——OpenCV研习社,本文为个人整理学习,已获得贾老师授权,有兴趣、有能力的可以加入贾老师的知识星球进行深入学习。