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