MATLAB代码

% 载入图像
grayImg = imread('lena256.bmp'); % 替换为你的图像路径


% 选择小波函数和分解级别
waveletFunction = 'db1';
level = 2;

% 执行WPT正向分解
wp= wpdec2(double(grayImg), level, waveletFunction);

% 从小波包分解中重建图像（逆向运算）
reconstructedImg = wprec2(wp);

% 显示原始图像和重建图像
subplot(1, 2, 1);
imshow(grayImg, []);
title('Original Image');

subplot(1, 2, 2);
imshow(uint8(reconstructedImg), []);
title('Reconstructed Image');

输出结果