IMAGE FUSION USING LAPLACIAN PYRAMID

%清屏
clear all; 
close all; 
clc;

%读取图像及尺寸
fimg=double(imread(‘1807i_12.bmp‘));
simg=double(imread(‘1807v_12.bmp‘));
[fm,fn]=size(fimg);[sm,sn]=size(simg);

%高斯滤波器
w=1/8*[0,1,0;1,4,1;0,1,0];

fimgn{1}=fimg;
%滤波，下采样
for i=2:5                  
   fimgn{i}=imfilter(fimgn{i-1},w,‘replicate‘);
   fimgn{i}=fimgn{i}(1:2:size(fimgn{i},1)-1,1:2:size(fimgn{i},2)-1); %i-1级近似
%    figure;
%    imshow(uint8(fimgn{i}));
%    imwrite(uint8(fimgn{i}),strcat(‘fimg_‘,num2str(i-1),‘.bmp‘));
end
for i=5:-1:2        %调整图像大小
   fimgn{i-1}=fimgn{i-1}(1:2*size(fimgn{i},1),1:2*size(fimgn{i},2)); 
end

simgn{1}=simg;
%滤波，下采样
for i=2:5                  
   simgn{i}=imfilter(simgn{i-1},w,‘replicate‘);
   simgn{i}=simgn{i}(1:2:size(simgn{i},1)-1,1:2:size(simgn{i},2)-1); %i-1级近似
%    figure;
%    imshow(uint8(simgn{i}));
%    imwrite(uint8(simgn{i}),strcat(‘simg_‘,num2str(i-1),‘.bmp‘));
end
for i=5:-1:2        %调整图像大小
   simgn{i-1}=simgn{i-1}(1:2*size(simgn{i},1),1:2*size(simgn{i},2)); 
end

for i=1:5
    averimgn{i}=(1/2)*(fimgn{i}+simgn{i});
%     figure;
%     imshow(uint8(averimgn{i}));
%     imwrite(uint8(averimgn{i}),strcat(‘averimg_‘,num2str(i-1),‘.bmp‘));
end

 for i=1:4          %获得残差图像，i级预测残差
    averimgn{i}=averimgn{i}-pyexpand(averimgn{i+1},w);    
 end
 
 for i=4:-1:1        %残差图像重构原图像
    averimgn{i}=averimgn{i}+pyexpand(averimgn{i+1},w);
%     figure;
%     imshow(uint8(averimgn{i}));
%     imwrite(uint8(averimgn{i}),strcat(‘averimg_‘,num2str(i-1),‘.bmp‘));
 end
 
 figure;
 imshow(uint8(averimgn{1}));
 imwrite(uint8(averimgn{1}),strcat(‘averimg.jpg‘));