写一段两张大小不同的图像融合的MATLAB 程序

以下为一段将两张大小不同的图像进行融合的MATLAB程序：

% 读入两张不同大小的图像
img1 = imread('image1.jpg');
img2 = imread('image2.jpg');

% 将两张图像进行大小调整
img2 = imresize(img2, size(img1));

% 将两张图像进行融合
alpha = 0.5; % 设置融合系数
blended_img = alpha*double(img1) + (1-alpha)*double(img2);
blended_img = uint8(blended_img);

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

在程序中，首先通过`imread`函数分别读入两张不同大小的图像，并通过`imresize`函数将第二张图像的大小调整为与第一张图像相同。然后，通过设置融合系数`alpha`，将两张图像进行融合。最后，通过`imshow`函数将融合后的图像显示出来。

需要注意的是，在进行融合时需要将图像转换为`double`类型，完成后再转回`uint8`类型。同时，在进行图像融合时需要保证两张图像的大小相同，否则会出现错误。

相关问题

写一段图像配准、图像拼接、图像融合的MATLAB程序

以下是一个简单的图像配准、图像拼接、图像融合的MATLAB程序示例：

% 读取待拼接的两张图片
image1 = imread('image1.jpg');
image2 = imread('image2.jpg');

% 提取图片1和图片2的SIFT特征点及其描述子
[f1, d1] = vl_sift(single(rgb2gray(image1)));
[f2, d2] = vl_sift(single(rgb2gray(image2)));

% 对SIFT特征点进行匹配
[matches, scores] = vl_ubcmatch(d1, d2);

% 提取匹配点对应的坐标
x1 = f1(1,matches(1,:));
y1 = f1(2,matches(1,:));
x2 = f2(1,matches(2,:));
y2 = f2(2,matches(2,:));

% 计算图像1到图像2的变换矩阵
H = estimateGeometricTransform([x1', y1'], [x2', y2'], 'affine');

% 使用变换矩阵将图像1与图像2进行拼接
transformedImage = imwarp(image1, H);

% 对拼接后的图像进行融合
blendedImage = imfuse(transformedImage, image2, 'blend', 'Scaling', 'joint');

% 显示结果
figure;
imshow(blendedImage);
title('Blended Image');

能给我一段基于拉普拉斯金字塔的图像融合matlab代码吗

当然可以，以下是基于拉普拉斯金字塔的图像融合MATLAB代码：

% 读取两张待融合的图像
img1 = imread('img1.jpg');
img2 = imread('img2.jpg');

% 将两张图像转换成灰度图
gray1 = rgb2gray(img1);
gray2 = rgb2gray(img2);

% 将灰度图像转换成双精度浮点型数据类型
gray1 = im2double(gray1);
gray2 = im2double(gray2);

% 对两张灰度图像分别进行拉普拉斯金字塔分解
levels = 3; % 金字塔层数
pyr1 = cell(levels,1);
pyr2 = cell(levels,1);
pyr1{1} = gray1;
pyr2{1} = gray2;
for i = 2:levels
    pyr1{i} = impyramid(pyr1{i-1}, 'reduce');
    pyr2{i} = impyramid(pyr2{i-1}, 'reduce');
end

% 对每一层金字塔进行权值计算，得到权值金字塔
wPyr = cell(levels,1);
for i = 1:levels
    [M,N] = size(pyr1{i});
    w = ones(M,N);
    w(:,1:round(N/2)) = 0;
    wPyr{i} = w;
end

% 对两张图像的拉普拉斯金字塔进行图像融合
blendedPyr = cell(levels,1);
for i = 1:levels
    blendedPyr{i} = pyr1{i}.*wPyr{i} + pyr2{i}.*(1-wPyr{i});
end

% 对融合后的图像金字塔进行重建
blended = blendedPyr{levels};
for i = levels-1:-1:1
    blended = impyramid(blended, 'expand');
    blended = blended + blendedPyr{i};
end

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

这段代码实现了基于拉普拉斯金字塔的图像融合，其中`img1.jpg`和`img2.jpg`是待融合的两张图像。代码中首先将两张图像转换成灰度图，然后对灰度图进行拉普拉斯金字塔分解，得到两个金字塔。接着根据金字塔的层数计算每一层的权值金字塔，然后对两张图像的拉普拉斯金字塔进行融合。最后将融合后的图像金字塔进行重建，得到融合后的图像并显示出来。