5/3小波的分解与重构Matlab

JP2K采用的一种图片压缩方式 --5/3小波，可用于图片压缩，资源内包括分解和重构的代码。

%-----------------------------5/3整数小波3级分解函数---------------------------------
  
function [S11d11]=decompose53（imMatdimn）
% imMat 输入图像的数据变量
%  dim  输入图像的维数
%   n   表示第n次分解
%  S11  表示分解的低频部分
%  d11  表示分解的高频部分

% imMat=imread（‘cameraman.tif‘）;
% dim=size（imMat1）;
% n=1;
imMat=imMat（1:dim/n1:dim/n）;
for m=1:dim/n
    for i=1:2:dim/n-1
        O1（m（（i+1）/2））=imMat（mi）;%得到第一级奇序列
    end
end
for m=1:dim/n
    for i=2:2:dim/n
        E1（mi/2）=imMat（mi）;%得到第一级偶序列
    end
end
E1=[E1zeros（dim/n2）];%再加两列相等行数的零元素（对于行变换）
for m=1:dim/n
    for i=1:dim/（n*2）
        d1（mi）=O1（mi）-E1（mi）/2-E1（mi+1）/2;
    end
end
d_1=[d1zeros（dim/n2）];%再加两列相等行数的零元素（行）
for m=1:dim/n
    for i=1:dim/（2*n）
        S1（mi）=E1（mi）+（d_1（mi）+d_1（mi+2）+2）/4;
    end
end
L1=[S1d1];
for m=1:dim/n
    for i=1:2:dim/n-1
        O11（（（i+1）/2）m）=L1（im）;%得到列变换的第一级奇序列
    end
end
for m=1:dim/n
    for i=2:2:dim/n
        E11（i/2m）=L1（im）;%得到列变换的第一级偶序列
    end
end
E11=[E11;zeros（2dim/n）];%再加两行相等列数的零元素（列）
for m=1:dim/n
    for i=1:dim/（2*n）
         d11（im）=O11（im）-E11（im）/2-E11（i+1m）/2;%各行差值数据；不对细节数据再做变换，只对S1再预测
    end
end
d_11=[d11;zeros（2dim/n）];%再加两行相等列数的零元素（列）
for m=1:dim/n
    for i=1:dim/（2*n）
          S11（im）=E11（im）+（d_11（im）+d_11（i+2m）+2）/4; 
   end
end 

% imshow（S11[]）
% imwrite（uint8（S11）‘cameraman_1.bmp‘）;

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

     文件       1477  2017-05-28 20:14  recompose53.m

     文件       1564  2017-11-23 19:59  decompose53.m

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

                 3041                    2