怎么给以下程序中加上两种噪声，最后还能恢复出原图像

yanxue0212

程序代码如下：

1. clear all;
   
2. iTimes=input('请输入置换次数iTime：'); %置乱次数
   
3. % 读入水印图像
   
4. M=input('请输入图片像素数');
   
5. file_name=input('请输入图片路径file_name：');
   
6. message=double(imread(file_name));
   
7. %水印图像矩阵的行数与列数
   
8. Mm=size(message,1);
   
9. Nm=size(message,2);
   
10. %对水印图像进行arnold置乱
    
11. if Mm~=Nm
    
12. error('水印矩阵必须为方阵');
    
13. end
    
14. if Mm~=M
    
15. error('必须为图片实际大小,或者修改置乱次数');
    
16. end
    
17. tempImg=message; %图像矩阵赋给tempImg
    
18. for n=1:iTimes %置乱次数
    
19. for u=1:Mm
    
20. for v=1:Nm
    
21. temp=tempImg(u,v);
    
22. ax=mod((u-1)+(v-1),Mm)+1; %新像素行位置
    
23. ay=mod((u-1)+2*(v-1),Nm)+1; %新像素列位置
    
24. outImg(ax,ay)=temp;
    
25. end
    
26. end
    
27. x=1;y=1;
    
28. T=1;m=Mm; %T为要计算的周期
    
29. t=x;x=x+y;y=t+2*y;
    
30. while x~=1&y~=1;
    
31. T=T+1;
    
32. if x>m
    
33. x=mod(x,m);
    
34. end
    
35. if y>m
    
36. y=mod(y,m);
    
37. end
    
38. t=x;x=x+y;y=t+2*y;
    
39. end
    
40. clc
    
41. tempImg=outImg;
    
42. end
    
43. % 显示水印，嵌入水印图像与原始图像
    
44. figure(1)
    
45. subplot(1,3,1);
    
46. imshow(message/255,[]);
    
47. title('原始水印');
    
48. subplot(1,3,2)
    
49. imshow(outImg,[]);
    
50. title('置乱水印');
    
51. % arnold反置乱
    
52. message_arnold=tempImg;
    
53. iTimes1=T-iTimes
    
54. %置乱后水印图像矩阵的行数与列数
    
55. Mo=size(outImg,1);
    
56. No=size(outImg,2);
    
57. for n=1:iTimes1 % 次数
    
58. for u=1:Mo
    
59. for v=1:No
    
60. temp1=tempImg(u,v);
    
61. bx=mod((u-1)+(v-1),Mo)+1;
    
62. by=mod((u-1)+2*(v-1),No)+1;
    
63. outImg1(bx,by)=temp1;
    
64. end
    
65. end
    
66. tempImg=outImg1;
    
67. end
    
68. message=outImg1;
    
69. subplot(1,3,3);
    
70. imshow(message,[]);
    
71. title('反置乱(恢复)水印');
    
72. % 显示反置乱后水印

复制代码

chybeyond

图像加噪声可参考imnoise

梦夕丶

逗 b，自己学吧！这么简单的问题也好意思问，你好意思问我都不好意思回答！