matlab例程

开发平台：

Matlab

wavelet.m：源码内容

							clear
%装入图像
X = imread('fish2.jpg');
%显示图像
subplot(121);
image(X);
% colormap(map)
title('原始图像');
axis square
disp('压缩前图像X的大小');
% whos('X')
%对图像用小波进行层小波分解
[c,s]=wavedec2(X,2,'bior3.7');
%提取小波分解结构中的一层的低频系数和高频系数
cal=appcoef2(c,s,'bior3.7',1);
%水平方向
ch1=detcoef2('h',c,s,1);
%垂直方向
cv1=detcoef2('v',c,s,1);
%斜线方向
cd1=detcoef2('d',c,s,1);
%各频率成份重构
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);
c1=[a1,h1*10;v1*10,d1*10];
%显示分频信息
subplot(122);
image(c1);
axis square;
title ('分解后低频和高频信息');

						