clear all;
close all;
clc;
im1 = (imread('C:\Users\swr\Desktop\图片6.png'));
im2 = imread('C:\Users\swr\Desktop\图片5.png');
figure(1);
imshow(im1);
title('红外图像');
figure(2);
imshow(im2);
title('可见光图像');
if length(size(im1)) ==3%判断输入图像是否是真彩色图像
im1 = rgb2gray(im1);
end
if length(size(im2)) ==3
im2 = rgb2gray(im2);
end
im1 = double(im1);
im2 = double(im2);
%%
%拉普拉斯滤波器
w = [1 4 6 4 1; 4 16 24 16 4; 6 24 36 24 6; 4 16 24 16 4; 1 4 6 4 1]/256;
G = cell(1,5);
H = cell(1,5);
I = cell(1,5);
G{1} = im1;%第一层为原图像
H{1} = im2;
for k=2:5
G{k} = imfilter(G{k-1},w,'replicate');%滤波
H{k} = imfilter(H{k-1},w,'replicate');%滤波
G{k} = G{k}(1:2:size(G{k},1),1:2:size(G{k},2)); %下采样
H{k} = H{k}(1:2:size(H{k},1),1:2:size(H{k},2)); %下采样
end
for k=5:-1:2
G{k-1} = G{k-1}(1:2*size(G{k},1),1:2*size(G{k},2));%调整图像尺寸
H{k-1} = H{k-1}(1:2*size(H{k},1),1:2*size(H{k},2));%调整图像尺寸
end
for k=1:4
G{k} = G{k} - us(G{k+1},w);%获得拉普拉斯金字塔
H{k} = H{k} - us(H{k+1},w);
a = abs(G{k}) > abs(H{k});%两幅图像的1-4级拉普拉斯金字塔取绝对值大的
I{k} = a.*G{k} + ~a.*H{k};
end
I{5} = H{5};%定义融合后顶层图像尺寸，并为图像边缘赋值，图像中间的值会被比较平均梯度后所取的值替换
gg = g(G{5});
gh = g(H{5});
for i = 2:size(G{5},1)-2
for j = 2:size(G{5},2)-2
if gg(i,j) > gh(i,j)%比较某点邻域平均梯度大小，取大的
I{5}(i,j) = G{5}(i,j);
else
I{5}(i,j) = H{5}(i,j);
end
end
end
for k=4:-1:1
I{k} = I{k}+us(I{k+1},w);%重构图像
end
I{1} = uint8(I{1});%转换为uint8用于显示
figure(3);
imshow(I{1});
title('拉普拉斯金字塔融合后图像');
评论0