例程11-11.m

来自「小波分析Matlab实现」· M 代码 · 共 52 行

% 当前延拓模式是补零
% 装载原始图像
load sinsin;
% X 包含原始图像
% 绘制原始图像
figure(1);
image(X);
colormap(map);
title('原始图像');
% map 包含原始彩图
nbcol = size(map,1);
% 使用db1对X进行单尺度分解
[cA1,cH1,cV1,cD1] = dwt2(X,'db1');
% 绘制各分解系数图像
figure(2);
subplot(2,2,1);
image(cA1);colormap(map);
title('db1低频系数图像');
subplot(2,2,2);
image(cH1);
title('db1水平高频图像');
subplot(2,2,3);
image(cV1);
title('db1垂直高频图像');
subplot(2,2,4);
image(cD1);
title('db1对角高频图像');
% 通过图像量化编码改善图像质量
cod_X = wcodemat(X,nbcol); 
cod_cA1 = wcodemat(cA1,nbcol); 
cod_cH1 = wcodemat(cH1,nbcol); 
cod_cV1 = wcodemat(cV1,nbcol); 
cod_cD1 = wcodemat(cD1,nbcol); 
dec2d = [... 
        cod_cA1,     cod_cH1;    ... 
        cod_cV1,     cod_cD1;    ... 
        ];
% 绘制各分解系数图像
figure(3);
subplot(2,2,1);
image(cod_cA1);colormap(map);
title('编码后低频系数图像');
subplot(2,2,2);
image(cod_cH1);
title('编码后水平高频图像');
subplot(2,2,3);
image(cod_cV1);
title('编码后垂直高频图像');
subplot(2,2,4);
image(cod_cD1);
title('编码后对角高频图像');