使用Matlab 进行高光谱数据分析（2）

Preprocess代码 

一. 求导

一阶导数MATLAB代码：

X1st=diff(X,1);%X 为输入光谱矩阵，X1st 为输出一阶导数光谱矩阵

二阶导数MATLAB代码：

X2st=diff(X,2); %X 为输入光谱矩阵，X2st 为输出二阶导数光谱矩阵

 

二. 多元散射校正

me=mean(X); 

[m,~]=size(X);

for i=1:m,

p=polyfit(me,X(i,:),1);

Xmsc(i,:)=(X(i,:)- p(2)*ones(1,n))./(p(1)*ones(1,n));

End

代码中输入光谱矩阵为 X，输出光谱矩阵为 Xmsc。

 

三. 标准正态变量变换 （SNV）

[~,n]=size(X);

rmean=mean(X,2);

dr=X-repmat(rmean,1,n);

Xsnv=dr./repmat(sqrt(sum(dr.^2,2)/(n-1)),1,n);

代码中 X 为光谱输入矩阵，Xsnv 为 SNV 处理后光谱输出矩阵。