图像的二进小波分解matlab源码（保证可运行）

基于二进小波变换的灰度图像分解(Image decomposition based on dyadic wavelet transform) matlab源码 代码说明和具体操作参见压缩包内的说明文档7-4.doc

%非正交的二次样条二进小波
h = [0 0 0.125 0.375 0.375 0.125 0 0].*sqrt（2）;
g = [0 0 0 -1 1 0 0 0].*sqrt（2）;
g1 = [0 0.015625 0.109375 0.34375 -0.34375 -0.109375 -0.015625 0].*sqrt（2）;
l = [1 6 15 84 15 6 1 0]./64;

n = 4;
img=imread（‘1.bmp‘）;
[a d1 d2] = swt1_decomp（imgnhg）;

figure;
imshow（img） title（‘原始矩阵‘）;

row = ceil（sqrt（n））;
col = ceil（n/row）;

% aa = a;
% figure title（‘逼近矩阵‘）;
% for i=1:n
%     aa（::i） = aa（::i）./max（max（aa（::i）））;
%     subplot（rowcoli）;
%     imshow（aa（::i））;
% end
% 
% 
% dd1 = d1;
% figure title（‘水平细节信息矩阵‘）;
% for i=1:n
%     dd1（::i） = dd1（::i）./max（max（dd1（::i）））;
%     subplot（rowcoli）;
%     imshow（dd1（::i））;
% end
% 
% dd2 = d2;
% figure title（‘竖直细节信息矩阵‘）;
% for i=1:n
%     dd2（::i） = dd2（::i）./max（max（dd2（::i）））;
%     subplot（rowcoli）;
%     imshow（dd2（::i））;
% end


aa = a;
for i=1:n
    aa（::i） = aa（::i）./max（max（aa（::i）））;
end
aa = reshape（aa [size（a1） size（a2） 1 size（a3）]）;
figure montage（aa） title（‘逼近矩阵‘）;

dd1 = d1;
for i=1:n
    dd1（::i） = dd1（::i）./max（max（dd1（::i）））;
end
dd1 = reshape（dd1 [size（a1） size（a2） 1 size（a3）]）;
figure montage（dd1） title（‘水平细节信息矩阵‘）;

dd2 = d2;
for i=1:n
    dd2（::i） = dd2（::i）./max（max（dd2（::i）））;
end
dd2 = reshape（dd2 [size（a1） size（a2） 1 size（a3）]）;
figure montage（dd2） title（‘竖直细节信息矩阵‘）;

x = swt1_recon（ad1d2nhg1l）;
% 
xx = x（::[2:n+1]）;
% row = ceil（sqrt（n））;
% col = ceil（n/row）;
% figure title（‘重构矩阵‘）;
for i=1:n
    xx（::i） = xx（::i）./max（[max（max（xx（::i））） 256]）;
%     subplot（rowcoli）;
%     imshow（xx（::i））;
end
xx = reshape（xx [size（xx1） size（xx2） 1 size（xx3）]）;
figure montage（xx） title（‘重构矩阵‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      66616  2005-05-23 20:56  ex7-4\1.BMP

     文件     772096  2006-04-03 10:36  ex7-4\7-4.doc

     文件       1751  2005-05-27 16:20  ex7-4\bat.m

     文件        492  2012-03-28 10:25  ex7-4\psnr.m

     文件       2136  2005-05-27 16:11  ex7-4\swt1_decomp.m

     文件       1891  2005-05-27 16:09  ex7-4\swt1_recon.m

     目录          0  2012-03-28 10:45  ex7-4

----------- ---------  ---------- -----  ----

               844982                    7