program_11_05.m

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

% 当前延拓模式是补零
% 装载原始图像
load sinsin; 
% X 包含装载的图像
subplot(3,2,1);
image(X);
colormap(map);
title('原始图像');
% 使用db4对X进行尺度为2的分解
[c,s] = wavedec2(X,2,'db4');
% 由系数重构第1层低频和细节部分
% 也可以使用wrcoef2完成，或者使用下面的步骤等价实现
% Step 1: 从分解结构[c,s]中提取系数
% Step 2: 使用upcoef2进行重构
siz = s(size(s,1),:); 
ca1 = appcoef2(c,s,'db4',1); 
a1 = upcoef2('a',ca1,'db4',1,siz);
% 绘制低频重构图
subplot(3,2,3);
image(a1);
colormap(map);
title('低频重构图');
chd1 = detcoef2('h',c,s,1); 
hd1 = upcoef2('h',chd1,'db4',1,siz); 
% 绘制高频水平部分重构图
subplot(3,2,4);
image(hd1);
colormap(map);
title('高频水平部分重构图');
cvd1 = detcoef2('v',c,s,1); 
vd1 = upcoef2('v',cvd1,'db4',1,siz);
% 绘制高频垂直部分重构图
subplot(3,2,5);
image(vd1);
colormap(map);
title('高频垂直部分重构图');
cdd1 = detcoef2('d',c,s,1); 
dd1 = upcoef2('d',cdd1,'db4',1,siz);
% 绘制高频对角部分重构图
subplot(3,2,6);
image(dd1);
colormap(map);
title('高频对角部分重构图');