图像多层小波变换

clear;
clc;
clf;


%读图像
Data=imread('C:\Users\pc\Desktop\压缩感知PPT\lena.png');




%画原图
subplot(1,3,1)
imshow(Data)
title('原彩图')


Data=rgb2gray(Data);
 
subplot(1,3,2)


imshow(Data)
title('灰度图像')


N=3;
[C,S]=wavedec2(Data,N,'haar');
Y(1:S(2,1),1:S(2,1))=reshape(C(1:S(2,1).^2),S(2,1),S(2,1))


for i=1:N
    sizex=S(i+1,2);
    cH=reshape(C(sizex.^2+1:sizex.^2*2),sizex,sizex);
    cV=reshape(C(sizex.^2*2+1:sizex.^2*3),sizex,sizex);
    cD=reshape(C(sizex.^2*3+1:sizex.^2*4),sizex,sizex);
    Y(sizex+1:sizex*2,1:sizex)=cH;
    Y(1:sizex,sizex+1:sizex*2)=cV;
    Y(sizex+1:sizex*2,sizex+1:sizex*2)=cD;
        
    end 


subplot(1,3,3)
imshow(Y)
title('多层小波变换图像')