基于小波变换的图像融合源代码

基于小波变换的图像融合代码，可用于比较实验。

%-----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