suoyin.m

基于幻方加密算发的全数字全息图制作与加密与解密

% 利用上面的im () 函数, 将RGB 图像转换为索引图像
[XX,map]= im('youdian','bmp') ;
save f0 XX map;
% 加噪并画出含噪声图像
[ggg]=addnoise('youdian','bmp',15);
figure(2)
subplot(221);subimage(ggg,map);title('含噪声图像');
% 强制消噪
[c,s]=wavedec2(g,3,'sym4'); 
%用小波函数sym4 对图像进行3 层小波分解
a3=wrcoef2('a',c,s,'sym4',3);
save f1 a3; % 保存强制消噪后的结果a3 到f1 中
% 默认阈值消噪
% 用ddencmp 获取在消噪过程中的默认阈值
[thr,sorh,keepapp]=ddencmp('den','w v',g) ;
subplot(221);subimage(ggg,map);title('含噪声图像');
g0=wdencmp(gb1,g,'sym4',3,thr,sorh,keepapp);
subplot(222);subimage(g0,map);title('全阈值消噪后的图像');
% 独立阈值消噪
thr_h=[122,124,145];thr_v=[123,148,150];thr_d=[108,134,160];
thr_h=[122,124,145];thr_v=[123,148,150];thr_d=[108,134,160];
thr=[thr_h;thr_d;thr_v];　　　　% 设置水平、垂直、斜线方向阈值
g1=wdencmp('lvd',g,'sym4',3,thr,'h'); % 采取每层不同阈值进行消噪
subplot(223);subimage(g1,map);title('独立阈值消噪后的图像');
save f2 g0; 　
% 保存默认阈值消噪后的结果g0 和独立阈值消噪后的结果g1 到f2