光谱预处理算法(matlab版本)

系列文章目录

近红外光谱分析技术属于交叉领域，需要化学、计算机科学、生物科学等多领域的合作。为此，在（北京邮电大学杨辉华老师团队）指导下，近期准备开源传统的PLS，SVM，ANN，RF等经典算和SG，MSC，一阶导，二阶导等预处理以及GA等波长选择算法以及CNN、AE等最新深度学习算法，以帮助其他专业的更容易建立具有良好预测能力和鲁棒性的近红外光谱模型。

---

---

前言

NIRS是介于可见光和中红外光之间的电磁波，其波长范围为（1100∼2526 nm。 由于近红外光谱区与有机分子中含氢基团（OH、NH、CH、SH）振动的合频和 各级倍频的吸收区一致，通过扫描样品的近红外光谱，可以得到样品中有机分子含氢 基团的特征信息，常被作为获取样本信息的一种有效的载体。 基于NIRS的检测方法具有方便、高效、准确、成本低、可现场检测、不 破坏样品等优势，被广泛应用于各类检测领域。但 近红外光谱存在谱带宽、重叠较严重、吸收信号弱、信息解析复杂等问题，与常用的 化学分析方法不同，仅能作为一种间接测量方法，无法直接分析出被测样本的含量或 类别，它依赖于化学计量学方法，在样品待测属性值与近红外光谱数据之间建立一个 关联模型(或称校正模型，Calibration Model) ，再通过模型对未知样品的近红外光谱 进行预测来得到各性质成分的预测值。现有近红外建模方法主要为经典建模 （预处理+波长筛选进行特征降维和突出，再通过pls、svm算法进行建模）以及深度学习方法（端到端的建模，对预处理、波长选择等依赖性很低）

---

本篇主要讲述基于matlab语言版本的光谱预处理方法,

一、预处理算法

%% 标准化Autoscales
function [ax,mx,stdx] = auto(x)
% Autoscales matrix to zero mean and unit variance
%
% [ax,mx,stdx] = auto(x)
%
% input:
% x 	data to autoscale
%
% output:
% ax	autoscaled data
% mx	means of data
% stdx	stantard deviations of data
%
% By Cleiton A. Nunes
% UFLA,MG,Brazil



[m,n] = size(x);
mx    = mean(x);
stdx  = std(x);
ax    = (x-mx(ones(m,1),:))./stdx(ones(m,1),:);

%# 数据中心化(Mean centering)
function [mcx,mx] = center(x)
% Mean center scales matrix to zero mean
%
% [mcx,mx] = center(x)
%
% input:
% x 	data to mean center
%
% output:
% ax	mean center data
% mx	means of data
%
% By Cleiton A. Nunes
% UFLA,MG,Brazil

[m,n] = size(x);
mx    = mean(x);
mcx   = (x-mx(ones(m,1),:));

%#归一化Normalize
function [nx] = normaliz(x)
% Normalize matrix rows dividing by its norm
% 
% [nx] = normaliz(x)
%
% input:
% x (samples x variables)   data to normalize
%
% output:
% nx (samples x variables)  normalized data
%
% By Cleiton A. Nunes
% UFLA,MG,Brazil

[m,n]=size(x);
nx=x;
nm=zeros(m,1);
for i = 1:m
nm(i)=norm(nx(i,:));
nx(i,:)=nx(i,:)/nm(i);
end
% z=[1