%給定一個(gè)原始圖像
load geometry;
nbcol=size(map,1);
colormap(pink(nbcol));
figure(1);
image(wcodemat(X,nbcol));
title('原始圖像');
%================================
%設(shè)置延拓模式為補(bǔ)零方式，對(duì)圖像應(yīng)用小波sym4進(jìn)行3層分解，然后重構(gòu)第3層近似
lev=3;
wname='sym4';
dwtmode('zpd'); 
[c,s]=wavedec2(X,lev,wname);
a=wrcoef2('a',c,s,wname,lev);
figure(2);
image(wcodemat(a,nbcol));
%================================
%改變延拓模式為對(duì)稱延拓且同樣應(yīng)用小波sym4將原始圖像分解到第3層，之后重構(gòu)第3層近似
dwtmode('sym');
[c,s]=wavedec2(X,lev,wname);
aa=wrcoef2('a',c,s,wname,lev);
figure(3);
image(wcodemat(aa,nbcol));
%================================
%設(shè)置平滑填補(bǔ)的延拓模式，應(yīng)用小波sym4對(duì)原始圖像進(jìn)行3層分解，然后重構(gòu)第3層近似
dwtmode('spd');
[c,s]=wavedec2(X,lev,wname);
aaa=wrcoef2('a',c,s,wname,lev);
figure(4);
image(wcodemat(aaa,nbcol));