yasuo.m

小波变换应用MATLAB对于图像处理的压缩程序

% 转载原始图像
load belmont2;
% X包含原始图像
x=X(100:200,100:200);
nbc=size(map,1);
subplot(221),image(x),colormap(map);
axis square;
title('原始图像');
% 使用wdencmp进行图像压缩
% 得到X的小波分解结构
n=5;w='sym2';
[c,l]=wavedec2(x,n,w);
% 小波系数的阈值处理
thr=20;
[xd,cxd,lxd,perf0,perfl2]=wdencmp('gbl',c,l,w,n,thr,'h',l);
% 画出压缩后的图像
subplot(222),image(xd),colormap(map);
title('全局阈值压缩图像');
axis square;
subplot(223),image(x),colormap(map);
axis square;
title('原始图像');
% 定义三个方向的各层不同的阈值
thr_h=[17 18];
thr_d=[19 20];
thr_v=[21 22];
thr=[thr_h;thr_d;thr_v];
% 进行图像压缩
[xd,cxd,lxd,perf0,perfl2]=wdencmp('lvd',x,'sym8',2,thr,'h');
% 画出压缩后的图像
subplot(224),image(xd),colormap(map);
title('独立阈值压缩图像');
axis square;