PCA 基于PCA的遥感图像融合

基于PCA的遥感图像融合，效果不错。可以作为遥感图像融合方面入门学习使用。

up=imread（‘6.bmp‘）;           %读图像 
low=imread（‘5.bmp‘）;
[low_R]=double（low（::1））; 
[low_G]=double（low（::2））; 
[low_B]=double（low（::3））; 

[MNcolor]=size（up）; 
low_Mx=0; 
for i = 1 : M 
    for j = 1 : N 
         low_S = [low_R（ij）low_G（ij）low_B（ij）]‘;   % 生成由R，G， B组成的三维列向量
         low_Mx = low_Mx + low_S; 
    end 
end 

low_Mx = low_Mx / （M*N）; % 计算三维列向量的平均值
 
low_Cx=0; 

for i = 1 : M 
    for j = 1 : N 
         low_S = [low_R（ij）low_G（ij）low_B（ij）]‘; 
         low_Cx = low_Cx + low_S*low_S‘; 
    end 
end

low_Cx = low_Cx / （M * N）- low_Mx*low_Mx‘;  % 计算协方差矩阵
[low_Alow_latent] = eigs（low_Cx）; %PCA变换的系数矩阵特征值
 
for i = 1 : M 
    for j = 1 : N
       low_X = [low_R（ij）low_G（ij）low_B（ij）]‘; % 生成由R，G， B组成的三维列 
       low_Y = low_A

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件     196662  2010-12-20 21:37  PCA\PCA\5.bmp

     文件     196662  2010-12-20 21:37  PCA\PCA\6.bmp

     文件       1843  2010-12-27 20:39  PCA\PCA\pca.m

    ..A.SH.     13824  2011-03-30 20:00  PCA\PCA\Thumbs.db

     目录          0  2011-03-30 20:00  PCA\PCA

     目录          0  2011-03-30 19:59  PCA

----------- ---------  ---------- -----  ----

               408991                    6