%讀入已經在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