拉普拉斯金字塔变换彩色图像融合matlab代码

以下是一个简单的拉普拉斯金字塔变换彩色图像融合的 Matlab 代码示例：

% 读取两张彩色图像
img1 = imread('image1.jpg');
img2 = imread('image2.jpg');

% 将图像转换为双精度类型
img1 = im2double(img1);
img2 = im2double(img2);

% 构建高斯金字塔
numLevels = 5;
filterSize = 5;
img1GaussianPyramid = buildGaussianPyramid(img1, numLevels, filterSize);
img2GaussianPyramid = buildGaussianPyramid(img2, numLevels, filterSize);

% 构建拉普拉斯金字塔
img1LaplacianPyramid = buildLaplacianPyramid(img1GaussianPyramid);
img2LaplacianPyramid = buildLaplacianPyramid(img2GaussianPyramid);

% 合并两张图像的拉普拉斯金字塔
mergedLaplacianPyramid = mergeLaplacianPyramids(img1LaplacianPyramid, img2LaplacianPyramid);

% 重建图像
mergedImage = reconstructImage(mergedLaplacianPyramid);

% 显示结果
imshow(mergedImage);

% 高斯金字塔构建函数
function [pyramid] = buildGaussianPyramid(img, numLevels, filterSize)
    pyramid = cell(numLevels, 1);
    pyramid{1} = img;
    for i = 2:numLevels
        img = imfilter(img, fspecial('gaussian', filterSize, 0.5), 'symmetric', 'conv');
        img = img(1:2:end, 1:2:end, :);
        pyramid{i} = img;
    end
end

% 拉普拉斯金字塔构建函数
function [pyramid] = buildLaplacianPyramid(gaussianPyramid)
    numLevels = length(gaussianPyramid);
    pyramid = cell(numLevels-1, 1);
    for i = 1:numLevels-1
        img = gaussianPyramid{i} - imresize(gaussianPyramid{i+1}, size(gaussianPyramid{i}), 'bilinear');
        pyramid{i} = img;
    end
end

% 拉普拉斯金字塔合并函数
function [mergedPyramid] = mergeLaplacianPyramids(pyramid1, pyramid2)
    numLevels = length(pyramid1);
    mergedPyramid = cell(numLevels, 1);
    for i = 1:numLevels
        mergedPyramid{i} = pyramid1{i} + pyramid2{i};
    end
end

% 图像重建函数
function [img] = reconstructImage(laplacianPyramid)
    numLevels = length(laplacianPyramid);
    img = laplacianPyramid{numLevels};
    for i = numLevels-1:-1:1
        img = imresize(img, size(laplacianPyramid{i}), 'bilinear') + laplacianPyramid{i};
    end
end

注意：本示例代码仅用于学习和参考，可能不适用于所有情况。在实际应用中，请根据具体需求进行修改和优化。