matlab图像融合代码,图像融合+源代码+matlab

【实例简介】

【实例截图】

【核心代码】

小波融合代码：

%%  小波变换图像融合

[imA,map1]

= imread('A.tif');

M1 =

double(imA) / 256;

[imB,map2]

= imread('B.tif');

M2 =

double(imB) / 256;

zt=

4;

wtype =

'haar';

%    M1 - input image A

%    M2 - input image B

%    wtype使用的小波类型

%    Y  -

fused image

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%   小波变换的绝对值大的小波系数，对应着显著的亮度变化，也就是图像中的显著特征。所以，选择绝对值大的小波系数作为我们需要的小波系数。

%%   低频部分系数采用二者求平均的方法

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

[c0,s0] =

wavedec2(M1, zt, wtype);%多尺度二维小波分解

[c1,s1] =

wavedec2(M2, zt, wtype);%多尺度二维小波分解

%%  取绝对值大的小波系数，作为融合后的小波系数

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

KK =

size(c1);

Coef_Fusion

= zeros(1,KK(2));

Temp =

zeros(1,2);

Coef_Fusion(1:s1(1,1))

= (c0(1:s1(1,1)) c1(1:s1(1,1)))/2;

%处理高频系数

MM1 = c0(s1(1,1) 1:KK(2));

MM2 = c1(s1(1,1) 1:KK(2));

mm = (abs(MM1)) > (abs(MM2));

Y  =

(mm.*MM1) ((~mm).*MM2);

Coef_Fusion(s1(1,1) 1:KK(2)) = Y;

%处理高频系数end

%重构

Y = waverec2(Coef_Fusion,s0,wtype);

%显示图像

subplot(2,2,1);imshow(M1);

colormap(gray);

title('

input image A ');

axis

square

subplot(2,2,2);imshow(M2);

colormap(gray);

title('

input image B ');

axis

square

subplot(223);imshow(Y,[]);

colormap(gray);

title('

fused image ');

axis

square;