program_13_16.m

这是基于MATLAB系统上开发的,应用于小波分析方面,很有用的程序模块,书名为"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',1);
% 画出压缩后图像
subplot(222),image(xd),colormap(map);
title('全局阈值压缩图像');
axis square;
% 定义三个方向各层不同的阈值
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(223),image(xd),colormap(map);
title('独立阈值压缩后图像');
axis square;