请用matlab写出利用拉普拉斯金字塔分解图像实现可见光和红外光图像融合的程序

以下是利用拉普拉斯金字塔分解图像实现可见光和红外光图像融合的 MATLAB 程序：

% 读取可见光图像和红外光图像
vis_img = imread('Visible_Image.png');
ir_img = imread('IR_Image.png');

% 将图像转换为灰度图像
vis_img_gray = rgb2gray(vis_img);
ir_img_gray = rgb2gray(ir_img);

% 对灰度图像进行拉普拉斯金字塔分解
n_levels = 5;
vis_lap_pyramid = laplacian_pyramid(vis_img_gray, n_levels);
ir_lap_pyramid = laplacian_pyramid(ir_img_gray, n_levels);

% 创建拉普拉斯金字塔混合图像
mix_lap_pyramid = cell(1, n_levels);
for i = 1:n_levels
  mix_lap_pyramid{i} = (vis_lap_pyramid{i} + ir_lap_pyramid{i}) / 2;
end

% 重建混合图像
mix_img = reconstruct_laplacian_pyramid(mix_lap_pyramid);

% 显示混合图像
imshow(mix_img);

% 拉普拉斯金字塔分解函数
function lap_pyramid = laplacian_pyramid(img, n_levels)
  lap_pyramid = cell(1, n_levels);
  for i = 1:n_levels
    if i == 1
      % 第一层为原始图像
      lap_pyramid{i} = img;
    else
      % 上采样前一层并减去下采样后一层
      prev_level = lap_pyramid{i-1};
      downsampled = imresize(prev_level, 0.5);
      upsampled = imresize(downsampled, size(prev_level));
      lap_pyramid{i} = prev_level - upsampled;
    end
  end
end

% 拉普拉斯金字塔重建函数
function img = reconstruct_laplacian_pyramid(lap_pyramid)
  n_levels = length(lap_pyramid);
  img = lap_pyramid{n_levels};
  for i = n_levels-1:-1:1
    % 上采样后一层并加上当前层
    upsampled = imresize(img, 2);
    img = lap_pyramid{i} + upsampled;
  end
end

请注意，此程序仅仅是实现了基于拉普拉斯金字塔的图像融合方法，实际应用中可能需要根据具体情况进行调整和优化。