cwq9.m

来自「图象的消噪,图象的压缩,图象的融合,图象的增强」· M 代码 · 共 60 行

load sinsin
[C,S]=wavedec2(X,2,'bior3.7');
%提取分解后的近似和细节系数
cA2=appcoef2(C,S,'bior3.7',2);
[cH2,cV2,cD2]=detcoef2('all',C,S,2);
[cH1,cV1,cD1]=detcoef2('all',C,S,2);
%从系数C重构第二层近似
A2=wrcoef2('a',C,S,'bior3.7',2);
%从系数C重构第一，第二层细节
A1=wrcoef2('a',C,S,'bior3.7',1);
H1=wrcoef2('h',C,S,'bior3.7',1);
V1=wrcoef2('v',C,S,'bior3.7',1);
D1=wrcoef2('d',C,S,'bior3.7',1);
H2=wrcoef2('h',C,S,'bior3.7',2);
V2=wrcoef2('v',C,S,'bior3.7',2);
D2=wrcoef2('d',C,S,'bior3.7',2);

%显示多尺度分解的结果：
colormap(map);

subplot(241);
image(wcodemat(A1,192));
title('近似A1');  

subplot(242);
image(wcodemat(H1,192));
title('水平系节H1'); 

subplot(243);
image(wcodemat(V1,192));
title('垂直细节V1'); 

subplot(244);
image(wcodemat(D1,192));
title('对角细节D1'); 

subplot(245);
image(wcodemat(A2,192));
title('近似A2'); 

subplot(246);
image(wcodemat(H2,192));
title('水平细节H2'); 

subplot(247);
image(wcodemat(V2,192));
title('垂直细节V2'); 

subplot(248);
image(wcodemat(D2,192));
title('对角细节D2'); 
%从多尺度分解后的系数重构原始图象并显示结果

figure(3);
X0=waverec2(C,S,'bior3.7');
image(X0);
colormap(map);
colorbar;