小波融合matlab代码
时间: 2023-10-27 17:05:21 浏览: 95
以下是使用小波变换实现图像融合的Matlab代码:
```matlab
% 读入需要融合的两幅图像
A = imread('A.jpg');
B = imread('B.jpg');
% 进行小波分解
= dwt2(A,'haar');
= dwt2(B,'haar');
% 对低频分量进行加权平均
w = 0.5;
C = w * AA + (1-w) * AB;
% 对高频分量进行加权平均
HC = w * HA + (1-w) * HB;
VC = w * VA + (1-w) * VB;
DC = w * DA + (1-w) * DB;
% 进行重构
C = idwt2(C,HC,VC,DC,'haar');
% 显示融合后的图像
imshow(uint8(C));
```
这段代码首先读入两幅需要融合的图像A和B,然后利用小波变换对它们进行分解,得到低频分量AA和AB以及高频分量HA、VA、DA和HB、VB、DB。接着对低频分量进行加权平均,得到融合后的低频分量C。然后对高频分量也进行加权平均,分别得到融合后的高频分量HC、VC、DC。最后,利用反小波变换idwt2对融合后的低频分量C和高频分量HC、VC、DC进行重构,得到最终的融合图像C。最后,使用imshow函数显示融合后的图像。
相关问题
小波图像融合matlab代码
以下是一个简单的小波图像融合的MATLAB代码示例:
```matlab
% 读取待融合的图像
image1 = imread('image1.jpg');
image2 = imread('image2.jpg');
% 将图像转换为灰度图像
grayImage1 = rgb2gray(image1);
grayImage2 = rgb2gray(image2);
% 对图像进行小波分解
[cA1, cH1, cV1, cD1] = dwt2(double(grayImage1), 'haar');
[cA2, cH2, cV2, cD2] = dwt2(double(grayImage2), 'haar');
% 将低频分量进行融合
fusedCA = (cA1 + cA2) / 2;
% 将高频分量进行融合
fusedCH = max(abs(cH1), abs(cH2)) .* exp(1i * (angle(cH1) + angle(cH2))) ;
fusedCV = max(abs(cV1), abs(cV2)) .* exp(1i * (angle(cV1) + angle(cV2))) ;
fusedCD = max(abs(cD1), abs(cD2)) .* exp(1i * (angle(cD1) + angle(cD2))) ;
% 逆小波变换得到融合后的图像
fusedImage = idwt2(fusedCA, fusedCH, fusedCV, fusedCD, 'haar');
% 显示融合后的图像
figure;
subplot(1, 3, 1); imshow(image1); title('图像1');
subplot(1, 3, 2); imshow(image2); title('图像2');
subplot(1, 3, 3); imshow(uint8(fusedImage)); title('融合后的图像');
```
请确保在运行代码之前将`image1.jpg`和`image2.jpg`替换为你自己的图像文件名。这段代码使用Haar小波作为小波基函数进行图像融合,并将融合后的图像显示在一个Matlab图形窗口中。
基于小波技术进行三个图像的融合MATLAB代码
以下是一个基于MATLAB实现的基于小波技术进行三个图像融合的示例代码:
```matlab
% 读入三幅图像
img1 = imread('image1.jpg');
img2 = imread('image2.jpg');
img3 = imread('image3.jpg');
% 对三幅图像进行小波变换
[c1,s1] = wavedec2(img1, 2, 'db4');
[c2,s2] = wavedec2(img2, 2, 'db4');
[c3,s3] = wavedec2(img3, 2, 'db4');
% 对小波系数进行融合,这里采用简单的平均值融合规则
c4 = (c1 + c2 + c3) / 3;
% 对融合后的小波系数进行逆小波变换,得到融合后的图像
img4 = waverec2(c4, s1, 'db4');
% 显示融合前后的图像
subplot(2,2,1),imshow(img1),title('Image 1');
subplot(2,2,2),imshow(img2),title('Image 2');
subplot(2,2,3),imshow(img3),title('Image 3');
subplot(2,2,4),imshow(uint8(img4)),title('Fused Image');
```
其中,`imread`函数用于读入图像,`wavedec2`和`waverec2`函数分别用于进行小波变换和逆小波变换,`subplot`和`imshow`函数用于显示图像。这里的融合规则是简单的平均值融合,可以根据实际需求进行修改。
阅读全文
相关推荐












