数字图像处理 噪声

import matplotlib.pyplot as plt
import numpy as np

#高斯噪声：均值为0，标准差为64
x1 = np.random.normal(loc=0, scale=64, size=(256, 256))
#瑞利噪声：（2/b）**0.5为1
x2 = np.random.rayleigh(scale=64, size=(256, 256))
#伽马噪声：（b-1）/a为2，放大32倍
x3 = np.random.gamma(shape=2, scale=32, size=(256, 256))
#指数噪声：a=1/32
x4 = np.random.exponential(scale=32, size=(256, 256))
#均匀噪声
x5 = np.random.uniform(low=0, high=1.0, size=(256, 256))
#脉冲噪声
x6 = np.random.random_integers(low=0, high=2.0, size=(256, 256))

for i, x in enumerate([x1, x2, x3, x4, x5, x6]):
    ax = plt.subplot(23 * 10 + i + 1)
    ax.hist(x.reshape(x.size), 64, normed=True)
    ax.set_yticks([])
    ax.set_xticks([])
plt.show()