program_11_15.m

这是基于MATLAB系统上开发的,应用于小波分析方面,很有用的程序模块,书名为"MATLAB在小波分析中的应用"

%载入原始图像
load sinsin;
% 图像编码
nbcol = size(map,1);
cod_X = wcodemat(X,nbcol);
% 画出原始图像
subplot(221)
image(cod_X)
title('原始图像');
colormap(map)
% 使用sym4在第3层对X进行SWT分解
[ca,chd,cvd,cdd] = swt2(X,3,'sym4');
% 输出分解图像
for k = 1:3
    % 对第k层进行图像编码
    cod_ca  = wcodemat(ca(:,:,k),nbcol);
    cod_chd = wcodemat(chd(:,:,k),nbcol);
    cod_cvd = wcodemat(cvd(:,:,k),nbcol);
    cod_cdd = wcodemat(cdd(:,:,k),nbcol);
    decl = [cod_ca,cod_chd;cod_cvd,cod_cdd];
    % 画出第k层图像分解系数
    subplot(2,2,k+1)
    image(decl)
    title(['SWT 分解: 低频系数和高频系数(level. ',num2str(k),')']);
    colormap(map)
end