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('對角細節系數');