能否提供一个MATLAB中实现灰度图三层小波变换的详细代码示例

在MATLAB中，你可以使用`wavedec2`函数来进行二维图像的多级小波分解，包括灰度图像的三层小波变换。这是一个简单的代码示例：

% 加载或创建一个灰度图像
img = imread('your_image_gray.jpg'); % 替换为你需要处理的图像文件名

% 确保图像为单通道灰度图像
if size(img, 3) == 3
    img = rgb2gray(img);
end

% 将图像进行三次小波分解
[lvl1, (cA1, cH1, cV1, cD1)] = wavedec2(img, 3, 'db4'); % 使用Daubechies 4级小波

% 分解结果包含四个部分：水平、垂直、对角线和细节系数（cA1-cD1）
% 可以分别查看各级别的系数和重构原始图像
subplot(2, 2, 1), imshow(cA1), title('Approximation Level 1');
subplot(2, 2, 2), imshow(cH1), title(' Horizontal Detail Level 1');
subplot(2, 2, 3), imshow(cV1), title('Vertical Detail Level 1');
subplot(2, 2, 4), imshow(cD1), title('Diagonal Detail Level 1');

% 重建图像
reconstructed_img = waverec2(lvl1, [cA1; cH1; cV1; cD1], 'db4', 'Best');

% 比较原图和重建图
figure, subplot(1, 2, 1), imshow(img, []), title('Original Image');
subplot(1, 2, 2), imshow(reconstructed_img, []), title('Reconstructed Image');