拉普拉斯金字塔图像融合程序.rar_mile47y_图像融合_拉普拉斯

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('拉普拉斯金字塔融合后图像');