高光谱图像融合现状
基于Pansharpening的高光谱图像融合方法
基于Pansharpening的高光谱图像融合方法是从多光谱图像与全色图像的融合方法迁移到高光谱图像融合方法的。它主要有以下几种方法:
成分替代
成分替代方法是对多光谱图像进行变换分离其空间部分,利用全色图像代替再进行逆变换。
matlab代码如下(示例):
function I_Fus_PCA = PCA(HS,PAN)
HS = double(HS);
%%%%%%%% Upsampling
ratio1 = size(PAN,1)/size(HS,1);
HSU = interp23tapGeneral(HS,ratio1);
imageLR = double(HSU);
imageHR = double(PAN);
[n,m,d]=size(imageLR);
imageLR = reshape(imageLR, [n*m,d]);
% PCA transform on MS bands
[W,PCAData] = princomp(imageLR);
F = reshape(PCAData, [n,m,d]);
% Equalization
I = F(:,:,1);
imageHR = (imageHR - mean(imageHR(:)))*std2(I)/std(imageHR(:)) + mean2(I);
% Replace 1st band with PAN
F(:,:,1) = imageHR;
% Inverse PCA
I_Fus_PCA = reshape(F,[n*m,d]) * W';
I_Fus_PCA = reshape(I_Fus_PCA,