基于matlab的图像融合源码

%读取并显示原始多光谱和全色图像数据
imgPan=imread（‘spots.bmp‘）;
imgMul=imread（‘dmtm.bmp‘）;
subplot（131）imshow（imgMul）xlabel（‘（a）低分辨率多光谱图像‘）；
subplot（132）imshow（imgPan）xlabel（‘（b）高分辨率全色图像‘）；
%分离TM图像RGB分量
mulR=imgMul（::1）;
mulG=imgMul（::2）;
mulB=imgMul（::3）;
 %对SPOT全色图像进行db13小波三层分解
 [CpanLpan]=wavedec2（imgPan3‘db13‘）;
 imgWH=Lpan（1:）;
 length=imgWH（1）*imgWH（2）;
 %对TM图像三分量进行db13小波三层分解
 [CtmrLtmr]=wavedec2（imgR3‘db13‘）;
Cr=Cpan;Cr