MATLAB实现PCA影像融合算法

MATLAB实现PCA影像融合算法，其中有三个函数，PCA变换，直方图匹配，PCA变换。输入的影像数据为多光谱影像与全色波段。

function [Data1] = HistogramMatching（data1data2）
    % 进行直方图匹配，第一个参数是待处理影像，第二个参数是参考图像
    data1_tbl = tbl（data1）; % 做表格计算
    data2_tbl = tbl（data2）; 
    
    X1 = data1_tbl（4:）; % 获取累计直方图
    X2 = data2_tbl（4:）; 
       
    % 确定灰度值替代关系
    M = zeros（length（X1）length（X2））;
    for i = 1 :length（X2）
        for j = 1:length（X1）
            M（ij） = abs（X1（j）-X2（i））;
        end
    end
    [~min_location] = min（M）;  % min_location使一个行向量，记录了每一列最小值的位置与下标的关系，
    % 求原图像通过位置信息映射到的新像素的值
    Data1 = [];
    A = size（data1）;
    for i = 1: A（1）
        for j = 1:A（2）
            Data1（ij） = min_location（data1（ij）+1）-1;  % +1是位置，-1是像元
        end
    end
end

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

     文件        797  2019-10-30 09:50  HistogramMatching.m

     文件        608  2019-12-12 21:30  PCA_Fusion.m

     文件        520  2019-12-13 08:46  PCA.m

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

                 1925                    3