离散小波变换图像融合MATLAB代码

以下是一个使用MATLAB实现离散小波变换（DWT）图像融合的示例代码：

% 读取两张待融合的图像
image1 = imread('image1.jpg');
image2 = imread('image2.jpg');

% 将图像转换为灰度图像
grayImage1 = rgb2gray(image1);
grayImage2 = rgb2gray(image2);

% 对两张灰度图像执行离散小波变换
[coefficients1, ~] = wavedec2(grayImage1, 2, 'db4');
[coefficients2, ~] = wavedec2(grayImage2, 2, 'db4');

% 将低频系数与高频系数分开
cA1 = appcoef2(coefficients1, [], 'db4', 2);
[cH1, cV1, cD1] = detcoef2('all', coefficients1, 2);

cA2 = appcoef2(coefficients2, [], 'db4', 2);
[cH2, cV2, cD2] = detcoef2('all', coefficients2, 2);

% 对低频系数执行融合
fusedCA = (cA1 + cA2) / 2;

% 对高频系数执行融合
fusedCH = max(cH1, cH2);
fusedCV = max(cV1, cV2);
fusedCD = max(cD1, cD2);

% 重构融合后的图像
fusedCoefficients = [fusedCA; fusedCH, fusedCV, fusedCD];
fusedImage = waverec2(fusedCoefficients, 'db4');

% 显示融合后的图像
figure;
subplot(1, 3, 1); imshow(grayImage1); title('Image 1');
subplot(1, 3, 2); imshow(grayImage2); title('Image 2');
subplot(1, 3, 3); imshow(uint8(fusedImage)); title('Fused Image');

请确保已将图像文件命名为`image1.jpg`和`image2.jpg`，并且这两张图像位于MATLAB当前工作目录下。

这个示例代码使用了`db4`小波作为离散小波变换的基函数。你可以根据需要选择其他小波基函数。代码中还使用了二级分解，你也可以根据需要更改分解的级数。

希望这个示例代码对你有所帮助！如有任何问题，请随时提问。