% 利用上面的im () 函數(shù), 將RGB 圖像轉(zhuǎn)換為索引圖像
[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'); 
%用小波函數(shù)sym4 對圖像進行3 層小波分解
a3=wrcoef2('a',c,s,'sym4',3);
save f1 a3; % 保存強制消噪后的結(jié)果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];　　　　% 設(shè)置水平、垂直、斜線方向閾值
g1=wdencmp('lvd',g,'sym4',3,thr,'h'); % 采取每層不同閾值進行消噪
subplot(223);subimage(g1,map);title('獨立閾值消噪后的圖像');
save f2 g0; 　
% 保存默認閾值消噪后的結(jié)果g0 和獨立閾值消噪后的結(jié)果g1 到f2