clc;
clear?all;
close?all;
%?讀入圖像
imgPan?=?imread（‘high.jpg‘）;
imgMul?=?imread（‘low.jpg‘）;

%?顯示源圖像
imshow（imgMul）?xlabel?（‘（a）低分辨率多光譜影像‘）;
figureimshow（imgPan）?xlabel?（‘（b）高分辨率全色圖影像‘）;

mulR?=?imgMul?（::1）;
mulG?=?imgMul?（::2）;
mulB?=?imgMul?（::3）;

%?高分辨率全色圖小波變換
[CpanLpan]?=?wavedec2（imgPan3‘db13‘）;
imgWH?=?Lpan（1:）;
length?=?imgWH（1）*imgWH（2）;

[CtmrLtmr]?=?wavedec2（mulR3‘?db13‘?）;?%?R通道小波變換，得到的低頻替換全色圖的低頻系數
Cr?=?Cpan;?Cr（1:length）?=?Ctmr（1:length）;

[CtmgLtmg]?=?wavedec2（mulG3‘?db13‘?）;
Cg?=?Cpan;?Cg（1:length）?=?Ctmg（1:length）;

[CtmbLtmb]?=?wavedec2（mulB3‘?db13‘?）;
Cb?=?Cpan;?Cb（1:length）?=?Ctmb（1:length）;

%?低頻為多光譜圖的低頻系數，高頻為全色圖的高頻系數，進行小波重構
imgResult（::1）?=?waverec2（CrLpan‘?db13‘?）;
imgResult（::2）?=?waverec2（CgLpan‘?db13‘?）;
imgResult（::3）?=?waverec2（CbLpan‘?db13‘?）;

F=uint8（imgResult）;

figure?imshow（F）?xlabel?（‘?（c）小波融合圖像‘?）;