基于稀疏表示的图像融合代码，matlab

以下是一个基于稀疏表示的图像融合的 Matlab 代码示例：

clear all;
close all;

% 读取两幅图像
I1 = imread('image1.jpg');
I2 = imread('image2.jpg');

% 将图像转换为灰度图
I1 = rgb2gray(I1);
I2 = rgb2gray(I2);

% 将灰度图像转换为 double 类型
I1 = im2double(I1);
I2 = im2double(I2);

% 将图像缩小为原来的 1/4，加快计算速度
I1 = imresize(I1, 0.25);
I2 = imresize(I2, 0.25);

% 设置稀疏表示的参数
lambda = 0.1; % 稀疏度
patchsize = 8; % 图像块大小
overlap = 4; % 图像块重叠区域大小

% 获取图像块
X1 = im2col(I1, [patchsize patchsize], 'sliding');
X2 = im2col(I2, [patchsize patchsize], 'sliding');

% 对每个图像块进行稀疏表示
D1 = dictionary(X1, lambda);
D2 = dictionary(X2, lambda);
alpha1 = sparse_coding(X1, D1, lambda);
alpha2 = sparse_coding(X2, D2, lambda);

% 对每个图像块进行融合
Y = zeros(size(X1));
for i = 1:size(X1, 2)
    y1 = D1 * alpha1(:, i);
    y2 = D2 * alpha2(:, i);
    y = fusion(y1, y2);
    Y(:, i) = y;
end

% 重建图像
I = col2im(Y, [patchsize patchsize], size(I1), 'sliding');
I = imresize(I, 4);

% 显示结果
figure;
subplot(1, 3, 1); imshow(I1); title('Image 1');
subplot(1, 3, 2); imshow(I2); title('Image 2');
subplot(1, 3, 3); imshow(I); title('Fused Image');

% 定义稀疏表示的字典函数
function D = dictionary(X, lambda)
    [U, ~, ~] = svd(X, 'econ');
    D = U(:, 1:round(lambda * size(U, 2)));
end

% 定义稀疏表示的编码函数
function alpha = sparse_coding(X, D, lambda)
    alpha = omp(D' * X, D' * D, lambda);
end

% 定义图像块融合函数
function y = fusion(y1, y2)
    alpha1 = abs(y1) ./ (abs(y1) + abs(y2));
    alpha2 = abs(y2) ./ (abs(y1) + abs(y2));
    y = alpha1 .* y1 + alpha2 .* y2;
end

这个代码使用稀疏表示的方法对两幅图像进行融合。首先将图像缩小为原来的 1/4，然后对每个图像块进行稀疏表示，再对每个图像块进行融合。最后重建图像并显示结果。