一维信号小波包分解后如何改变分解系数再重构呢?
书上没有看到,然后在网上查,查到的是用于图像的,但是我用其中的write(wp1,'cfs',nd(j),cpjn); 根本没有变化,应该是没有写入,请高手指点,网上查的图像的处理代码如下: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');
页:
[1]