matlab dwt2 小波分解,小波分解矩阵 Matlab.doc

小波分解矩阵 Matlab.doc

%----------------------------------------------------------%

小波图像分解 Matlab 程序 - V2.0版

/chenyusiyuan/archive/2008/06/05/2513865.aspx

?

小波图像重构 Matlab 程序 - V2.0版

/chenyusiyuan/archive/2008/06/05/2514119.aspx

%----------------------------------------------------------%

?

?

%----------------------------------------------------------%?

???? 小波分解矩阵 Matlab 程序 - V3.0版

%----------------------------------------------------------%

?

function [coef,scf]=mywavedec2(x,N,wname)

%----------------------------------------------------------%

% 函数 MYWAVEDEC2() 对输入矩阵 x 进行 dim 层分解，得到相应的分解系数矩阵 y

% 输入参数：x —— 输入矩阵

%????????? N —— 分解级数

%??????? ??wname —— 分解所用的小波函数

% 输出参数：scf —— 存储各级分解系数矩阵的大小以及原始(图像)矩阵的大小

%?????????? coef —— 分解系数矩阵，其结构如下：

% coef = {cA_N;cV_N;cH_N;cD_N;cV_N-1;cH_N-1;cD_N-1;……;cV_1;cH_1;cD_1}

%

% Copyright by Zou Yuhua ( chenyusiyuan )

% Version: 3.0, Date: 2008-07-08

%----------------------------------------------------------%

?

?

% 求出小波函数的滤波器组系数向量

[Lo_D,Hi_D] = wfilters(wname,'d');

lf=length(Lo_D);

% 画出原始图像

imshow(x);title('Original Image');

% 标明图像大小

[r,c]=size(x);

xlabel(['Size : ',num2str(r),'*',num2str(c)]);

?

?

% 将矩阵x的数据格式转换为适合数值处理的double格式

xd=double(x);??????????????

[rx,cx]=size(x);

[o2sa,f1sa,rsx]=sizcoef([rx,cx],lf,N);? a=[o2sa,f1sa,rsx]

?

coef=[];

scf=[rx,cx];

for i=1:N

??? [cA,cV,cH,cD]=mydwt2(xd,Lo_D,Hi_D);???? % 第 i 级小波分解???

??? xd=cA;????????????? % 将第 i 级分解得到的低频系数矩阵作为第 i+1 级分解的源矩阵

??? outmp={cV;cH;cD};????? % 将第 i 级分解得到的高频系数矩阵cV,cH,cD存入细胞矩阵 outmp

???

??? scf=[size(cV);scf]; % 将各级分解系数矩阵的大小存入矩阵 scf

??? coef=[outmp;coef];?????? % 将细胞矩阵 outmp 存入输出矩阵 coef，coef将由空矩阵变为细胞矩阵???

end

% 迭代结束后，矩阵 coef 中保存的是各级分解中的高频系数矩阵

% 故需将迭代后得到的矩阵 cA，即第 dim 级低频矩阵存入矩阵 coef

coef=[cA;coef];

scf=[size(cA);scf];

?

?

plotcoef(N,wname,coef);

plotcoef2one(N,wname,coef,rsx);

?

function plotcoef2one(N,wname,coef,rsx)

% 画出小波分解系数的塔式结构图

rsx=rsx(end:-1:1,:);

tmpcoef=[];

tA=wkeep(coef{1},rsx(1,:),'c');

tmpcoef=tA;

tA=uint8(tA); tA(end,:)=255; tA(:,end)=255;

for j=1:N

??? tV=wkeep(coef{(j-1)*3+2},rs