资源简介

离散小波变换与离散小波反变换 快速小波变换(Mallat小波分解算法):对一幅图像做2级小波分解(离散小波变换)与合成(离散小波反变换)

代码片段和文件信息

x=imread(‘lena.gif‘);   % 载入图像lena.gif,需要更改路径到指定的文件
figure(1);
imshow(x);
title(‘原始图像‘);
x=double(x);

%图像分解
lpd=[1/2 1/2];hpd=[-1/4 2/3];           % 默认的低通、高通滤波器
[rowcol]=size(x);                      % 读取输入矩阵的大小

%一级分解
for j=1:row                             % 首先对输入矩阵的每一行序列进行一维离散小波分解
    tmp1=x(j:);
    [ca1cd1]=mydwt(tmp1lpdhpd1);
    x(j:)=ca1;                   % 将分解系数序列再存入矩阵x中,得到[L|H]
end
for k=1:col                             % 再对输入矩阵的每一列序列进行一维离散小波分解
    tmp2=x(:k);
    [ca2cd2]=mydwt(tmp2lpdhpd1);
    x(:k)=ca2;                   % 将分解所得系数存入矩阵x中,得到[LLHl;LHHH]
end

%二级分解
for j=1:row/2                             % 首先对输入矩阵左上角四分之一矩阵的每一行序列进行一维离散小波分解
    tmp1=x(j1:col/2);
    [ca1cd1]=mydwt(tmp1lpdhpd1);
    x(j1:col/2)=ca1;               

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

     文件     359424  2010-05-09 21:54  实验二 图像变换.doc

     文件        308  2010-05-08 14:15  downspl.m

     文件       1204  2010-05-08 14:17  myidwt.m

     文件        431  2010-05-08 14:18  upspl.m

     文件        846  2010-05-09 14:54  mydwt.m

     文件       4033  2010-05-09 21:54  2.m

     文件       4804  2010-05-09 21:19  wave2gray.m

     文件       1089  2004-12-16 13:36  wavecut.m

     文件       3528  2004-12-16 13:36  wavework.m

     文件        940  2004-12-16 13:36  wavecopy.m

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

               376607                    10


评论

共有 条评论