%-----ImageFusion?based?on?Wavelet------------------------------%
%低頻小波系數絕對值較大融合策略，高頻基于局部顯著性準則
%---2009-03-20@儲彬彬---QQ304150436---E_mail:cbbpaper@126.com---%
%---------------------------------------------------------------%
%function??varargout=myimfus（‘1.bmp‘‘2.bmp‘）
clear??all
clc
x1=imread（‘1.bmp‘）;
x1=rgb2gray（x1）;
x1=double（x1）/255;?
x2=imread（‘2.bmp‘）;
x2=rgb2gray（x2）;
x2=double（x2）/255;?
subplot（221）
imshow（‘1.bmp‘）
title（‘待融合圖像1‘）
subplot（222）
imshow（‘2.bmp‘）
title（‘待融合圖像2‘）
[ca1ch1cv1cd1]=dwt2（x1‘db1‘）;
[ca2ch2cv2cd2]=dwt2（x2‘db1‘）;
[rowcol]=size（ca1）;

fun1?=?inline（‘max（abs（x（:）））‘）
sh1=nlfilter（ch1[3?3]fun1）;
sh2=nlfilter（ch2[3?3]fun1）;
sv1=nlfilter（cv2[3?3]fun1）;
sv2=nlfilter（cv2[3?3]fun1）;
sd1=nlfilter（cd1[3?3]fun1）;
sd2=nlfilter（cd2[3?3]fun1）;

for?i=1:row
????for?j=1:col???????
?%-----低頻融合策略?-----%??
???????%低頻取平均
???????%?ca（ij）=（ca1（ij）+ca2（ij））/2;
????????%低頻取絕對值大????
???????if?abs（ca1（ij））>abs（ca2（ij））
????????????ca（ij）=ca1（ij）;
???????else?
???????????ca（ij）=ca2（ij）;
???????end
??%-----高頻融合策略?-----%????????
???????if?sh1（ij）>=sh2（ij）
????????mhx（ij）=1;
????????mhy（ij）=0;
???????else
????????mhx（ij）=0;
????????mhy（ij）=1;
???????end
????????if?sv1（ij）>=sv2（ij）
????????mvx（ij）=1;
????????mvy（ij）=0;
???????else
????????mvx（ij）=0;
????????mvy（ij）=1;
????????end
????????if?sd1（ij）>=sd2（ij）
????????md（ij）=1;
????????md（ij）=0;
???????else
????????md（ij）=0;
????????md（ij）=1;
????????end
????end
end

fun2?=?inline（‘sum（x（:））‘）
mxx=nlfilter（mhx[3?3]fun2）;?
mvv=nlfilter（mvx[3?3]fun2）;?
mdd=nlfilter（mhx[3?3]fun2）;
for?i=1:row
????for?j=1:col
????????
????????if?mxx（ij）>=5
???????????ch（ij）=ch1（ij）;
????????else
????????????ch（ij）=ch2（ij）;
????????end
????????if?mvv（ij）>=5
???????????cv（ij）=cv1（ij）;
????????else
????????????cv（ij）=cv2（ij）;
????????end??
????????if?mdd（ij）>=5
???????????cd（ij）=cd1（ij）;
????????else
????????????cd（ij）=cd2（ij）;
????????end?????
????end
end
?????
x=idwt2（cachcvcd‘db1‘）;
imwrite（x‘wavefuion2.bmp‘）;
subplot（223）
imshow（x）
title（‘融合后圖像‘）

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件?????196662??2009-02-27?16:41??基于小波變換的圖像融合源代碼\1.bmp

?????文件?????196662??2009-02-27?16:41??基于小波變換的圖像融合源代碼\2.bmp

?????文件??????66614??2011-12-31?15:14??基于小波變換的圖像融合源代碼\3.bmp

?????文件???????2248??2011-12-31?15:24??基于小波變換的圖像融合源代碼\imfusWavelet2.asv

?????文件???????2248??2011-12-31?15:24??基于小波變換的圖像融合源代碼\imfusWavelet2.m

?????文件???????1321??2009-04-03?08:40??基于小波變換的圖像融合源代碼\imgfusionwave1.m

?????文件??????66614??2011-12-31?15:13??基于小波變換的圖像融合源代碼\nDimfus.bmp

?????文件????????942??2011-12-31?15:13??基于小波變換的圖像融合源代碼\nDimFuswav.m

?????文件????????771??2009-04-03?08:36??基于小波變換的圖像融合源代碼\SimpleImgfus.m

?????文件??????66614??2011-12-31?15:25??基于小波變換的圖像融合源代碼\wavefuion2.bmp

?????文件??????66614??2011-12-31?15:02??基于小波變換的圖像融合源代碼\wavefusionV1.bmp

?????目錄??????????0??2012-01-07?18:01??基于小波變換的圖像融合源代碼

-----------?---------??----------?-----??----

???????????????667310????????????????????12