图像融合 matlab 源代码 实验报告

图像 融合 matlab 源代码 实验报告

function H=chonggou（h00h01h10h11）
%h00是行低列低，h01是行低列高，h10是行高列低，h11是行高列高
%思路：（（（（（h00列插——低通）+（h01列插——高通））——移位）——行插——低通）+
%（（（h10列插——低通）+（h11列插——高通））——移位）——行插——高通））——移位————————得到重构图像
%======h00偶数行插入0========%
[mn]=size（h00）;
for i=1:m
    H00（2*i-1:）=h00（i:）;
    H00（2*i:）=zeros（1n）;
end
%-----低通----%
for i=1:2*m
    if i==2*m
       F00（i:）=（H00（i:）+H00（1:））/sqrt（2）;
    else
       F00（i:）=（H00（i:）+H00（i+1:））/sqrt（2）;
    end
end
%=====h01偶数行插入0======%
for i=1:m
    H01（2*i-1:）=h01（i:）;
    H01（2*i:）=zeros（1n）;
end
%---高通----%
for i=1:2*m
    if i==2*m
       F01（i:）=（H01（i:）-H01（1:））/sqrt（2）;

    else
       F01（i:）=（H01（i:）-H01（i+1:））/sqrt（2）;
    end
end
%--求和并且移位---%
f00=F00+F01;
 F0（1:）=f00（2*m:）;
for i=1:2*m-1
    F0（i+1:）=f00（i:）;
end
%======h10偶数行插入0========%
for i=1:m
    H10（2*i-1:）=h10（i:）;
    H10（2*i:）=zeros（1n）;
end
%---低通----%
for i=1:2*m
    if i==2*m
       F10（i:）=（H10（i:）+H10（1:））/sqrt（2）;
 
    else
       F10（i:）=（H10（i:）+H10（i+1:））/sqrt（2）;
    end
end
%=====h11偶数行插入0======%
for i=1:m
    H11（2*i-1:）=h11（i:）;
    H11（2*i:）=zeros（1n）;
end
%---高通----%
for i=1:2*m
    if i==2*m
       F11（i:）=（H11（i:）-H11（1:））/sqrt（2）;

    else
       F11（i:）=（H11（i:）-H11（i+1:））/sqrt（2）;
    end
end
%--求和并且移位---%
f11=H10+H11;
F1（1:）=f11（2*m:）;
for i=1:2*m-1
    F1（i+1:）=f11（i:）;
end
%======F0偶数列插入0========%
for i=1:n
    FL0（:2*i-1）=F0（:i）;
    FL0（:2*i）=zeros（2*m1）;
end
%---低通----%
for i=1:2*n
    if i==2*n
       FL（:i）=（FL0（:i）+FL0（:1））/sqrt（2）;   
    else
       FL（:i）=（FL0（:i）+FL0（:i+1））/sqrt（2）;
    end
end
%======F1偶数列插入0========%
for i=1:n
    FHO（:2*i-1）=F1（:i）;
    FH0（:2*i）=zeros（2*m1）;
end
%---高通----%
for i=1:2*n
    if i==2*n
       FH（:i）=（FH0（:i）-FH0（:1））/sqrt（2）;   

    else
       FH（:i）=（FH0（:i）-FH0（:i+1））/sqrt（2）;
    end
end
%--求和并且移位---%
F=FL+FH;
H（:1）=F（:2*n）;
for i=1:2*n-1
    H（:i+1）=F（:i）;
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       2209  2011-11-08 17:12  图像融合\chonggou.m

     文件       1202  2011-11-08 17:11  图像融合\main.m

     文件        959  2011-11-08 17:12  图像融合\nonoise.m

     文件       1410  2011-11-08 17:11  图像融合\pipei.m

     文件        556  2011-11-08 17:11  图像融合\xbbh.m

     文件     193078  2004-05-20 10:38  图像融合\x_1.bmp

     文件     193078  2004-05-20 10:38  图像融合\x_2.bmp

     文件     405504  2011-11-08 17:14  图像融合\实验4.doc

     目录          0  2011-11-08 17:17  图像融合

----------- ---------  ---------- -----  ----

               797996                    9