图像小波包分解后如何改变分解系数再重构呢？

原文地址：图像小波包分解后如何改变分解系数再重构呢？ 作者：杨风暴

目前市场上介绍图像小波包变换Matlab编程的书不多，但也有几本，但是他们不少都是将Matlab中相关的函数的帮助重新叙述而已。因为Matlab帮助中只是说明了分解结果放到一个树的数据结构中，提取查看还行，如何改变分解系数就不知道了。

我们知道，一般小波包分解后就是对小波包系数举行处理，将处理改变后的系数重构，从而得到最后的处理结果，因此，如何得到系数、将改变后的系数再写进树的结构中是关键。

下面介绍图像小波包分解后如何改变分解系数再重构：

p=imread('mwf1.bmp');     %打开某一图像

p1=double(p);             %将图像矩阵变成双精度型

wp1=wpdec2(p1,2,'db2','shannon');    %用db2小波对图像小波包两层分解

nd=leaves(wp1);    %得到分解树wp1的所有也结点

cpj=wpcoef(wp1,nd(j));    %得到第nd(j)个结点对应的小波包分解系数矩阵，例如j=3；

cpjn=cpj.*2;    %改变分解系数矩阵

wp1=write(wp1,'cfs',nd(j),cpjn);    %将改变后的系数矩阵cpjn写入原来的分解树wp1的结构中第nd(j)个结点中。

p2=wprec2(wp1);    %重构改变后的小波分解系数，得到处理的后的图像

figure(1);    %显示处理前后的图像，进行对比

subplot(211);imagesc(p);colormap('gray');

subplot(212);imagesc(uint8(p2));colormap('gray');