program_11_15.m

来自「小波分析理论与matlab2007配套程序」· M 代码 · 共 27 行

%载入原始图像
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