PCA图像融合源码

%-----------------基于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变换融合（完）1_PAN.bmp
     文件      786486  2012-07-18 18:00  4-PCA变换融合（完）1_TM.bmp
     文件      196662  2012-09-06 16:01  4-PCA变换融合（完）3_PAN.bmp
     文件      196662  2012-09-06 16:01  4-PCA变换融合（完）3_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变换融合（完）FTPP.m
     文件       38912  2013-03-10 13:14  4-PCA变换融合（完）单个程序表格.doc