IHS图像变换融合 matlab

%读入已经在ENVI中匹配的EMT影像1238波段
I1=imread（‘G:\10.tif‘‘TIFF‘）;
I2=imread（‘G:\20.tif‘‘TIFF‘）;
I3=imread（‘G:\30.tif‘‘TIFF‘）;
Ipan=imread（‘G:\80.tif‘‘TIFF‘）;
%将RGB图像放大2倍
I1=imresize（I12）;
I2=imresize（I22）;
I3=imresize（I32）;
imgrgb=cat（3I1I2I3）;
%显示RGB图像和pan波段图像
figureimshow（imgrgb）title（‘RGB图像‘）;
figureimshow（Ipan）title（‘pan波段图像‘）;
%对RGB图像进行IHS变换
%A为变换矩阵
A=[1/（3^0.5） 1/（3^0.5） 1/（3^0.5）
    1/（6^0.5） 1/（6^0.5） -2/（6^0.5）
    1/（2^0.5） -1/（2^0.5） 0];


%调整RGB向量
p1=reshape（I11670*670）;
p2=reshape（I21670*670）;
p3=reshape（I31670*670）;
pan=r