work 基于IHS融合算法的多光谱与全色图像融合代码

基于IHS融合算法的多光谱与全色图像融合代码，使用于融合ETM+遥感图像

imgPan1=imread（‘pan6.bmp‘）;
imgMul=imread（‘multi6.bmp‘）;
imgPan=rgb2gray（imgPan1）;
subplot（141）imshow（imgPan）xlabel（‘（a）高分辨率灰度图像‘）;
subplot（142）imshow（imgMul）xlabel（‘（b）低分辨率彩色图像‘）;
[rowscolscolor]=size（imgMul）;
for i=1:rows
m=（i-1）*cols+1;
mulRGB（1m:m+cols-1）=imgMul（i:1）;
mulRGB（2m:m+cols-1）=imgMul（i:2）;
mulRGB（3m:m+cols-1）=imgMul（i:3）;
end
IHS=[1/sqrt（3）1/sqrt（3）1/sqrt（3）;1/sqrt（6）1/sqrt（6）-2/sqrt（6）;1/sqrt（2）-1/sqrt（2）0];
imgl=IHS*double（mulRGB）
imgl_m

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1017  2010-04-17 12:22  work.m

----------- ---------  ---------- -----  ----

                 1017                    1