小波图像分解和重构程序每句都带解释.docVIP

小波图像分解程序function coef=mywavedec2(x,N,wname) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % 函数 MYWAVEDEC2() 对输入矩阵 x 进行 dim 层分解，得到相应的分解系数矩阵 y % 输入参数：x —— 输入矩阵 %????????? ????? N —— 分解级数 %????????? ????? wname —— 分解所用的小波函数 % 输出参数：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 ), original : 2007-11-10, modified: 2008-06-04 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % 求出小波函数的滤波器组系数向量 [Lo_D,Hi_D] = wfilters(wname,d); % 画出原始图像 imshow(x);title(Original Image); % 标明图像大小 [r,c]=size(x); xlabel([Size : ,num2str(r),*,num2str(c)]); % 将矩阵x的数据格式转换为适合数值处理的double格式 xd=double(x);?????????????? coef=[]; 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 ??? ?? % 注意细胞矩阵的赋值是用大括号“{}”的，而普通矩阵赋值是用方括号“[]”???? % 细胞矩阵不要求其中的子矩阵的行列数都相同 ??? ?? coef=[outmp;coef];?? % 将细胞矩阵 outmp 存入输出矩阵 coef，coef将由空矩阵变为细胞矩阵 ??? ?? % 注意这里的方括号不能用大括号取代???? % 否则，使用大括号会将初始的coef空矩阵也作为细胞矩阵的子矩阵???? % 而且，在迭代中 coef 将是一个不断嵌套的细胞矩阵，不便于后续处理和读取???????? % 上面这个语句是一种有效的在迭代过程中保存数据的方法???? % 设待存数据为 data，可以是单个数、向量或矩阵???? % 保存数据的矩阵为 mat，初始为空矩阵：mat=[]???? % 则可按以下格式保存迭代过程产生的数据???? % mat=[mat;data];???? % 方括号内的分号“;”表示数据 data 是按“列”排序的方式存入矩阵 mat???? % mat=[mat,data];???? % 方括号内的逗号“,”表示数据 data 是按“行”排序的方式存入矩阵 mat ??? % data 也可以在 mat 前嵌入，即 mat=[data;mat] 或 mat=[data,mat] end % 迭代结束后，矩阵 coef 中保存的是各级分解中的高频系数矩阵 % 故需将迭代后得到的矩阵 cA，即第 dim 级低频矩阵存入矩阵 coef coef=[cA;coef]; % 最后，小波系数矩阵 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} % 画出各级低频、高频系数矩阵 % 首先建立一个名为“Wavelet Decomposition -- Wavelet Type: , Levels: ”的图像窗口 figure(Name,[Wavelet Decomposition -- Wavelet Type: ,wname, , Levels: ,num2str(N)]); % 图像的第1行显示低频系数，置中，左右两个subplot为空 subplot(N+1,3,2); yt=uint8(coef{1}); [yrow,ycol]=size(yt); imshow(yt); title( [Approximation A,num2str(N)]); xlabel([Size : ,num2str(yrow),