wtest.m

用dwt函数实现了一层小波变换

info=imfinfo('E:\照片\图片\背景\lena.bmp');
[source,map]=imread('E:\照片\图片\背景\lena.bmp');
source
subplot(121);
imshow(source);
title('原始图像');
smodify=im2double(source);
%%这里用double或im2double都可以，dct时是因为变换矩阵是0到1的数所以那时需要用im2double
[ca1,ch1,cv1,cd1]=dwt2(smodify,'db1');
%[m,n]=size(ca1)   用dwt变换是
ca1
mcol=size(map,1);
cod_source=wcodemat(smodify,mcol);       %根据值的整体范围，将它的值变成对应的灰度值
cod_a=wcodemat(ca1,mcol);
cod_h=wcodemat(ch1,mcol);
cod_v=wcodemat(cv1,mcol);
cod_d=wcodemat(cd1,mcol);
subplot(122);
imshow(cod_source,map);
title('编码后的原始图像');
figure
subplot(221)
imshow(cod_a,map);
title('近似小波系数');
subplot(222)
imshow(cod_h,map);
title('水平细节系数');
subplot(223)
imshow(cod_v,map);
title('垂直细节系数');
subplot(224)
imshow(cod_d,map);
title('对角细节系数');