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）小波融合圖像‘?）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1015??2015-08-14?20:32??DWT\DWT.m

?????文件??????39964??2015-07-04?19:43??DWT\high.jpg

?????文件??????33350??2015-07-04?19:42??DWT\low.jpg

?????目錄??????????0??2015-08-14?20:51??DWT

-----------?---------??----------?-----??----

????????????????74329????????????????????4