1. from pylab import *
2. from matplotlib import *
3. import numpy
4.  
5. def rgb2gray(rgb):
6.         r,g,b=rgb[:, :, 0], rgb[:,:,1], rgb[:,:,2]
7.         gray=0.2989*r+0.5870*g+0.1140*b
8.         return gray
9.  
10. image = imread("oko.jpg")
11. gimage=rgb2gray(image)
12. imshow(gimage, cmap="gray")
13. title("Obraz podstawowy")
14. show()
15.  
16. szum=rand(gimage.shape[0],gimage.shape[1])
17. szum=szum*100
18. imshow(szum, cmap="gray")
19. title("Szum o rozkladzie rownomiernym")
20. show()
21.  
22. gimageszum = gimage + szum
23. imshow(gimageszum, cmap="gray")
24. title("Obraz z szumem o rozkladzie rownomiernym")
25. show()
26.  
27. szumGauss=randn(gimage.shape[0],gimage.shape[1])
28. szumGauss=szumGauss * 100
29. imshow(szumGauss, cmap="gray")
30. title("Szum Gaussa")
31. show()
32.  
33. gimageszumGauss= gimage + szumGauss
34. imshow(gimageszumGauss, cmap="gray")
35. title("Obraz z szumem Gaussa")
36. show()
37.  
38. sol_pieprz = gimage
39. for i in range(8000):
40. sol_pieprz[randint(0,sol_pieprz.shape[0]),randint(0,sol_pieprz.shape[1])] =1
41. sol_pieprz[randint(0,sol_pieprz.shape[0]),randint(0,sol_pieprz.shape[1])] =0
42. imshow(sol_pieprz, cmap="gray")
43. title("Obraz z szumem pieprz i sol")
44. show()