实验二
完整代码:https://github.com/Arancew/digtal-image-processing
- 高斯噪声
- 椒盐噪声
1高斯噪声
import cv2
import matplotlib.pyplot as plt
import numpy
import skimage
plt.rc("font", family='Microsoft YaHei')
img = plt.imread('../img/2.1.jpg')
plt.subplot(221), plt.title('原始图像')
plt.imshow(img, 'gray')
noise_img = skimage.util.random_noise(img, mode='gaussian')*255
plt.subplot(222), plt.title('高斯噪声图像')
plt.imshow(noise_img, 'gray')
mean_img = noise_img
for i in range(1, noise_img.shape[0] - 1):
for j in range(1, noise_img.shape[1] - 1):
tmp = 0
for k in range(-1, 2):
for l in range(-1, 2):
tmp += noise_img[i + k][j + l]
mean_img[i][j] = tmp / 9
plt.subplot(223), plt.title('均值滤波后图像')
plt.imshow(mean_img, 'gray')
median_img = noise_img
for i in range(1, noise_img.shape[0] - 1):
for j in range(1, noise_img.shape[1] - 1):
tmpp = []
for k in range(-1, 2):
for l in range(-1, 2):
tmpp.append(noise_img[i + k][j + l])
list.sort(tmpp)
median_img[i][j] = tmpp[4]
plt.subplot(224), plt.title('中值滤波后图像')
plt.imshow(median_img, 'gray')
plt.show()
2.椒盐噪声
import cv2
import matplotlib.pyplot as plt
import numpy
import skimage
plt.rc("font", family='Microsoft YaHei')
img = plt.imread('../img/2.1.jpg')
plt.subplot(221), plt.title('原始图像')
plt.imshow(img, 'gray')
noise_img = skimage.util.random_noise(img, mode='salt')*255
plt.subplot(222), plt.title('椒盐噪声图像')
plt.imshow(noise_img, 'gray')
mean_img = noise_img
for i in range(1, noise_img.shape[0] - 1):
for j in range(1, noise_img.shape[1] - 1):
tmp = 0
for k in range(-1, 2):
for l in range(-1, 2):
tmp += noise_img[i + k][j + l]
mean_img[i][j] = tmp / 9
plt.subplot(223), plt.title('均值滤波后图像')
plt.imshow(mean_img, 'gray')
median_img = noise_img
for i in range(1, noise_img.shape[0] - 1):
for j in range(1, noise_img.shape[1] - 1):
tmpp = []
for k in range(-1, 2):
for l in range(-1, 2):
tmpp.append(noise_img[i + k][j + l])
list.sort(tmpp)
median_img[i][j] = tmpp[4]
plt.subplot(224), plt.title('中值滤波后图像')
plt.imshow(median_img, 'gray')
plt.show()