写一个matlab 图像融合

图像融合是将两幅或多幅图像合并为一幅图像的过程，常用于图像增强、图像处理等领域。下面给出一个简单的 MATLAB 图像融合程序。

首先，我们需要读入两幅图像，并将它们转换为灰度图像：

% 读入两幅图像
img1 = imread('image1.jpg');
img2 = imread('image2.jpg');

% 将彩色图像转换为灰度图像
gray1 = rgb2gray(img1);
gray2 = rgb2gray(img2);

接下来，我们使用拉普拉斯金字塔对图像进行分解：

% 定义拉普拉斯金字塔层数
levels = 4;

% 构建高斯金字塔
gauss1 = gray1;
gauss2 = gray2;
for i = 1:levels
    gauss1 = impyramid(gauss1, 'reduce');
    gauss2 = impyramid(gauss2, 'reduce');
end

% 构建拉普拉斯金字塔
lap1 = cell(1, levels);
lap2 = cell(1, levels);

lap1{levels} = gauss1;
lap2{levels} = gauss2;
for i = levels-1:-1:1
    lap1{i} = gauss1 - impyramid(gauss1, 'expand');
    lap2{i} = gauss2 - impyramid(gauss2, 'expand');
    gauss1 = impyramid(gauss1, 'reduce');
    gauss2 = impyramid(gauss2, 'reduce');
end

然后，我们对每一层的拉普拉斯金字塔进行融合：

% 定义融合权重
alpha = 0.5;

% 构建融合后的拉普拉斯金字塔
lap_blend = cell(1, levels);
for i = 1:levels
    lap_blend{i} = alpha*lap1{i} + (1-alpha)*lap2{i};
end

最后，我们使用重建图像的方法将融合后的拉普拉斯金字塔合并为一幅图像：

% 重建图像
img_blend = lap_blend{1};
for i = 2:levels
    img_blend = impyramid(img_blend, 'expand');
    img_blend = img_blend + lap_blend{i};
end

完整程序如下：

% 读入两幅图像
img1 = imread('image1.jpg');
img2 = imread('image2.jpg');

% 将彩色图像转换为灰度图像
gray1 = rgb2gray(img1);
gray2 = rgb2gray(img2);

% 定义拉普拉斯金字塔层数
levels = 4;

% 构建高斯金字塔
gauss1 = gray1;
gauss2 = gray2;
for i = 1:levels
    gauss1 = impyramid(gauss1, 'reduce');
    gauss2 = impyramid(gauss2, 'reduce');
end

% 构建拉普拉斯金字塔
lap1 = cell(1, levels);
lap2 = cell(1, levels);

lap1{levels} = gauss1;
lap2{levels} = gauss2;
for i = levels-1:-1:1
    lap1{i} = gauss1 - impyramid(gauss1, 'expand');
    lap2{i} = gauss2 - impyramid(gauss2, 'expand');
    gauss1 = impyramid(gauss1, 'reduce');
    gauss2 = impyramid(gauss2, 'reduce');
end

% 定义融合权重
alpha = 0.5;

% 构建融合后的拉普拉斯金字塔
lap_blend = cell(1, levels);
for i = 1:levels
    lap_blend{i} = alpha*lap1{i} + (1-alpha)*lap2{i};
end

% 重建图像
img_blend = lap_blend{1};
for i = 2:levels
    img_blend = impyramid(img_blend, 'expand');
    img_blend = img_blend + lap_blend{i};
end

% 显示结果
figure;
subplot(2,2,1); imshow(img1); title('Image 1');
subplot(2,2,2); imshow(img2); title('Image 2');
subplot(2,2,[3,4]); imshow(img_blend); title('Blended Image');