%%%%%%%%%%%%%%%%%%%%%
function [ output_args ] = Untitled9( input_args )
%UNTITLED9 Summary of this function goes here
% Detailed explanation goes here
f=imread('lena.bmp');
imshow(f);
[X, map] = gray2ind(f, 256);%将灰度图像转换为索引图像imview(X, map);
nbcol=size(map,1);
[cA1,cH1,cV1,cD1]=dwt2(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);
figure;
imshow(Y,map);
title('重构后的图像')