10-9.m

来自「matlab图像处理源码」· M 代码 · 共 34 行

load woman;
%打开原始图像
imshow(X,map);
% X包含原始图像信息. 
% map包含打开的色图
nbcol = size(map,1); 
[cA1,cH1,cV1,cD1] = dwt2(X,'db1');
% 对图像X执行单层分解，小波为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); 
%图像编码
figure 
subplot(221)
imshow(cod_cA1,map)
title('近似细节系数')
subplot(222)
imshow(cod_cH1,map)
title('水平细节系数')
subplot(223)
imshow(cod_cV1,map)
title('垂直细节系数')
subplot(224)
imshow(cod_cD1,map)
title('对角细节系数')

Y = idwt2(cA1,cH1,cV1,cD1,'db1',nbcol);
% 对分解的细节系数执行单层重构，小波为db1
figure;
imshow(Y,map);
title('重构后的图像')