function?[?Result?]?=?Image_Fusion（?K_PanK_Mul?）
%?此函數用于圖像融合
%?K_Pan和K_Mul分別為待融合的全色圖像與多光譜圖像
%?Result為融合之后的結果
?
A1=double（K_Pan）;
A2=double（K_Mul）;
[RowColThi]=size（A2）;
Result=zeros（RowColThi）;

%?選擇要執行的操作
str={‘Brovey變換‘‘乘積變換‘‘PCA變換‘‘HSI變換‘};
[selok]=listdlg（‘Liststring‘str‘PromptString‘‘請選擇融合算法?‘?‘SelectionMode‘‘single‘）;??

%?brovey變換
if?sel==1
???B=zeros（RowCol）;
???for?i=1:Thi
???????B=A2（::i）+B;
???end
???for?i=1:Thi
???????C（::i）=A1.*A2（::i）./（B+eps）;
???end
???C（B==0）=0;
???
???for?i=1:Thi
???????d=max（max（C（::i）））;
???????Result（::i）=C（::i）/d;
???end

%?乘積變換
elseif?sel==2
???for?i=1:Thi
???????C（::i）=A1.*A2（::i）;
???end
???for?i=1:Thi
???????d=max（max（C（::i）））;
???????Result（::i）=C（::i）/d;
??