matlab代码：已知有400乘400分辨率的蓝、绿、红、近红外四波段多光谱图象和一张1600乘1600全色图象，参照pansharpening方法和多分辨率分析融合，进行遥感图象融合

以下是 MATLAB 代码实现：

% 读取多光谱图像和全色图像
ms_img_b = imread('ms_blue.tif');
ms_img_g = imread('ms_green.tif');
ms_img_r = imread('ms_red.tif');
ms_img_nir = imread('ms_nir.tif');
pan_img = imread('pan.tif');

% 多分辨率分析融合
ms_img = cat(3, ms_img_b, ms_img_g, ms_img_r, ms_img_nir);
r = 4; % 分辨率等级数
w = 1/2; % 分辨率变换系数
pyr = cell(r, 1);
pyr{1} = ms_img;
for i = 2:r
    pyr{i} = impyramid(pyr{i-1}, 'reduce');
end
for i = r-1:-1:1
    pyr{i} = imresize(pyr{i+1}, size(pyr{i}(1:2:end, 1:2:end, :)), 'bicubic');
    pyr{i} = pyr{i} + imresize(pyr{i-1}, size(pyr{i}(1:2:end, 1:2:end, :)), 'bicubic');
end
fused_img = pyr{1};
for i = 2:r
    fused_img = imresize(fused_img, size(pyr{i}(1:2:end, 1:2:end, :)), 'bicubic');
    fused_img = fused_img + pyr{i};
end

% Pansharpening
fused_img_l = imresize(fused_img, size(pan_img), 'bicubic');
fused_img_h = fused_img - fused_img_l;
fused_img_p = pan_img - fused_img_l;
fused_img_p = imresize(fused_img_p, size(fused_img_h), 'bicubic');
fused_img_ph = fused_img_h + fused_img_p;
fused_img_final = fused_img_l + fused_img_ph;

% 显示结果
figure;
subplot(2,2,1); imshow(ms_img(:,:,1:3)); title('多光谱图像');
subplot(2,2,2); imshow(pan_img); title('全色图像');
subplot(2,2,3); imshow(fused_img); title('多分辨率分析融合结果');
subplot(2,2,4); imshow(fused_img_final); title('Pansharpening后的结果');

代码中用到了 `impyramid` 函数实现图像金字塔，用到了 `imresize` 函数实现图像大小变换。最后使用 `imshow` 函数显示结果。