给我一段基于拉普拉斯金字塔进行可见光与红外的灰度图像融合的matlab代码

下面是一段基于拉普拉斯金字塔进行可见光与红外的灰度图像融合的Matlab代码：

% 读取可见光图像和红外图像
vis_img = imread('visible.jpg');
ir_img = imread('infrared.jpg');

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

% 构建可见光图像和红外图像的拉普拉斯金字塔
vis_lap_pyr = laplacian_pyramid(vis_gray);
ir_lap_pyr = laplacian_pyramid(ir_gray);

% 构建融合后的拉普拉斯金字塔
fusion_lap_pyr = cell(1, length(vis_lap_pyr));
for i = 1:length(vis_lap_pyr)
    fusion_lap_pyr{i} = (vis_lap_pyr{i} + ir_lap_pyr{i}) / 2;
end

% 重建融合后的图像
fusion_img = reconstruct_image(fusion_lap_pyr);

% 显示融合后的图像
imshow(fusion_img);

% 拉普拉斯金字塔函数
function lap_pyr = laplacian_pyramid(img)
    % 构建高斯金字塔
    gauss_pyr = cell(1, 6);
    gauss_pyr{1} = img;
    for i = 2:6
        gauss_pyr{i} = imresize(gauss_pyr{i-1}, 0.5);
    end

    % 构建拉普拉斯金字塔
    lap_pyr = cell(1, 5);
    for i = 1:5
        lap_pyr{i} = gauss_pyr{i} - imresize(gauss_pyr{i+1}, size(gauss_pyr{i}));
    end
end

% 重建图像函数
function img = reconstruct_image(lap_pyr)
    % 从最小的拉普拉斯金字塔开始重建图像
    img = lap_pyr{5};
    for i = 4:-1:1
        img = imresize(img, size(lap_pyr{i})) + lap_pyr{i};
    end
end

请注意，这段代码仅供参考，并且可能需要根据您的具体需求进行修改。