%讀取并顯示原始多光譜和全色圖像數據
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