program_11_10.m

这是基于MATLAB系统上开发的,应用于小波分析方面,很有用的程序模块,书名为"MATLAB在小波分析中的应用"

% 装载图像
load belmont2;
% 显示图像
image(X); 
colormap(map);
colorbar;
% 执行图像的单尺度小波分解
[cA1,cH1,cV1,cD1] = dwt2(X,'bior3.7');
% 由系数重构低频和高频部分
A1 = upcoef2('a',cA1,'bior3.7',1);
H1 = upcoef2('h',cH1,'bior3.7',1);
V1 = upcoef2('v',cV1,'bior3.7',1);
D1 = upcoef2('d',cD1,'bior3.7',1);
% 显示低频和高频部分
figure(2);
colormap(map);
nbcol = size(map,1)
subplot(2,2,1); image(wcodemat(A1,nbcol));
title('低频A1')
subplot(2,2,2); image(wcodemat(H1,nbcol));
title('水平高频H1')
subplot(2,2,3); image(wcodemat(V1,nbcol));
title('垂直高频V1')
subplot(2,2,4); image(wcodemat(D1,nbcol));
title('对角高频D1')
% 由小波逆变换恢复原图像信号
Xsyn = idwt2(cA1,cH1,cV1,cD1,'bior3.7');
% 图像的多尺度二维小波分解
[C,S] = wavedec2(X,2,'bior3.7');
% 提取系数的低频和高频部分
cA2 = appcoef2(C,S,'bior3.7',2);
% 重构第二层的低频信号
A2 = wrcoef2('a',C,S,'bior3.7',2);
% 重构第1、2层的高频信号
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);
% 显示多尺度二维分解结果
figure(3)
colormap(map);
subplot(2,4,1);image(wcodemat(A1,nbcol));
title('低频A1')
subplot(2,4,2);image(wcodemat(H1,nbcol));
title('水平高频H1')
subplot(2,4,3);image(wcodemat(V1,nbcol));
title('垂直高频V1')
subplot(2,4,4);image(wcodemat(D1,nbcol));
title('对角高频D1')
subplot(2,4,5);image(wcodemat(A2,nbcol));
title('低频A2')
subplot(2,4,6);image(wcodemat(H2,nbcol));
title('水平高频H2')
subplot(2,4,7);image(wcodemat(V2,nbcol));
title('垂直高频V2')
subplot(2,4,8);image(wcodemat(D2,nbcol));
title('对角高频D2')
% 重构原始图像信号
X0 = waverec2(C,S,'bior3.7');
% 显示重构得到的图像
figure(4)
image(X0); 
colormap(map);