%-----------------基于HSI空間變換的圖像融合方法----------------------------%

clearclose?allclc
PAN?=?imread（‘03_PAN.bmp‘）;
TM?=?imread（‘03_TM.bmp‘）;
PAN0?=?PAN;
TM0?=?TM;

[PC1?PCo?cov]?=?myPCA（TM0）;???%轉換到HSI空間
PC1?=?mat2gray（PC1）;
PAN1?=?rgb2gray（PAN）;

Hist_PC1?=?imhist（PC1）;
[mn]?=?size（PC1）;
Hist_PC1?=?Hist_PC1/（（m*n）+10^-10）;
[PANppTform]?=?hist_op（PAN1‘match‘Hist_PC1）;
%?PANpp?=?mat2gray（PANpp）;

%直接替換亮度法：
PC1_pp?=?double（PANpp）;
FR?=?myIPCA（PC1_ppPCocov）;
FR?=?double2uint8（FR）;


subplot（131）imshow（PAN）;title（‘原始全色圖像‘）
subplot（132）imshow（FR）;title（‘融合結果圖像‘）
subplot（133）imshow（TM）;title（‘原始多光譜圖像‘）

imwrite（FR‘Fusion_Result.bmp‘‘bmp‘）;

B?=?C_PJZB（FRTMPAN）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2015-07-21?21:56??4-PCA變換融合（完）\
?????文件??????786486??2012-07-18?18:00??4-PCA變換融合（完）\01_PAN.bmp
?????文件??????786486??2012-07-18?18:00??4-PCA變換融合（完）\01_TM.bmp
?????文件??????196662??2012-09-06?16:01??4-PCA變換融合（完）\03_PAN.bmp
?????文件??????196662??2012-09-06?16:01??4-PCA變換融合（完）\03_TM.bmp
?????文件?????????753??2013-03-10?13:10??4-PCA變換融合（完）\A_Fusion_PCA.m
?????文件????????8237??2012-09-13?17:07??4-PCA變換融合（完）\C_PJZB.m
?????文件?????????136??2012-09-14?17:02??4-PCA變換融合（完）\double2uint8.m
?????文件??????196662??2013-03-10?13:25??4-PCA變換融合（完）\Fusion_Result.bmp
?????文件????????5547??2012-09-22?10:48??4-PCA變換融合（完）\hist_op.m
?????文件?????????491??2012-07-09?14:10??4-PCA變換融合（完）\mycov.m
?????文件?????????374??2012-07-10?08:32??4-PCA變換融合（完）\myIPCA.m
?????文件?????????400??2012-07-10?08:17??4-PCA變換融合（完）\myPCA.m
?????文件?????????996??2013-03-10?13:19??4-PCA變換融合（完）\ZFTPP.m
?????文件???????38912??2013-03-10?13:14??4-PCA變換融合（完）\單個程序表格.doc