光谱数据重采样

matlab代码，用于高光谱、多光谱数据重采样，内容清晰

%%输入文件
clear all
load ‘D:\工具\GER光谱数据_重采样\all.txt‘;
[rowcol] = size（all）;
band=all（:1）;
all（:1）=[]; %%%%%%%%%%%%%%%%将第一列波长从矩阵中去掉
col=col-1; %%%%%%%%%%%%%%%%矩阵的采样点数目，即光谱曲线数目
% newband=[400:1:750];
load ‘D:\工具\GER光谱数据_重采样\newband.txt‘;

%%Piecewise Cubic Hermite Interpolating Polynomial （PCHIP）内插
for point=1:col
    %point;
    spec=all（:point）;
Interp_spec = pchip（bandspecnewband）; %%%%PCHIP既能保持光谱的形状特征，又能够不会过拟合，避免某些直线的部分拟合成曲线。pchip（xyf）  x表示内插前的x变量，y表示内插前的y变量，f是内插后的x变量。
Interp_Spec（:point）=Interp_spec;      
end


%%输出
plot（bandall（:1）‘o‘newbandInterp_Spec（:1）‘-.‘）
legend（‘data‘‘pchip‘4）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件       84924  2017-05-15 15:26  GER光谱数据_重采样\all.txt
     文件         748  2017-05-15 15:26  GER光谱数据_重采样\Interpchip.m
     文件        3006  2017-05-15 15:26  GER光谱数据_重采样\newband.txt