为什么散射校正要放在首位
食品与饲料行业中粉末、谷物和粕类的漫反射 NIR 光谱,除化学吸收信号外,还叠加了随样品而变化的物理效应。Barnes 等将粒度和光散射共同造成的乘性效应描述为解释此类光谱的主要障碍,并指出粒度占光谱方差的大部分,而化学组成的贡献相对较小 [1]。Rinnan 等将常用的 NIR 预处理方法分为散射校正(MSC、SNV、归一化)和光谱导数两类,并在实际数据上比较了它们的效果 [3]。本指南介绍第一类;导数方法见Savitzky-Golay 指南。
多元散射校正(MSC)
MSC 由 Geladi、MacDougall 和 Martens 于 1985 年针对肉类的 NIR 反射光谱提出,是一种将化学吸收与物理散射分离的多波长方法,其中散射被建模为加性分量与乘性分量之和 [2]。每条光谱 xi 在全部波长上对参考光谱 xref 进行回归:xi ≈ ai + bi·xref;校正后的光谱为 (xi − ai)/bi。实践中参考光谱取定标集的平均光谱,因此 Dhanoa 等称 MSC 为依赖于样品集的变换 [4]。在该肉类研究中,MSC 与校正 log(1/R) 或逆 Kubelka-Munk 等线性化变换结合,得到了最佳的脂肪预测结果和更简单的定标模型 [2]。随后,Isaksson 和 Næs 在五种食品上用主成分回归对蛋白质、脂肪、水分和碳水化合物检验了 MSC,发现它在每一个回归中都降低了预测误差,相对于未校正的吸光度光谱降低 7 % 至 68 % [5]。
这种对样品集的依赖带来一个实际后果:参考光谱是定标模型的一部分,必须原样应用于每条新光谱;若在新的样品总体上重新定标后更换参考光谱,则每条校正后的光谱都会随之改变 [4]。
标准正态变量(SNV)与去趋势
Barnes、Dhanoa 和 Lister 于 1989 年提出 SNV,作为对每条光谱单独施加的变换:减去该光谱自身在全部波长上的均值,再除以其自身的标准差,从而在不参照任何其他样品的情况下消除散射和粒度的乘性效应 [1]。在同一篇论文中,他们还提出了去趋势:对经 SNV 校正的光谱拟合一个关于波长的二次多项式并将其减去,以消除他们所述在粉状或紧密填充样品中常见的基线偏移和弯曲;他们将这一结果作为导数光谱的替代方案,认为导数光谱的形状更难解释 [1]。由于 SNV 只使用待校正的光谱本身,它与样品集无关:无需选择、存储或传递任何参考光谱。
MSC 与 SNV 的关系与差异
Dhanoa 等证明,正如理论所预言,依赖样品集的 MSC 变换与不依赖样品集的 SNV 变换之间存在线性关系;这一关系只需要样品集平均光谱的均值和标准差,以及单条光谱与该平均光谱之间的相关系数,而 MSC 对样品集的依赖正是通过这三个量引入的 [4]。他们的结论是,二者是降低粒度效应的两条可互换的替代途径 [4]。对于均匀的定标集,二者通常给出相似的光谱和相似的模型,因此常被视为可以互换。
然而二者并不完全相同。Fearn 等在光谱空间中考察了两种变换的几何性质,报告了两种不同的行为:SNV 可能在处理后光谱的得分图中诱发弯曲结构,而 MSC 则倾向于产生离群点;两种效应均以 NIR 成像相机数据加以说明并作了几何解释 [6]。因此,SNV 之后 PCA 得分图中出现的弯曲条带,或 MSC 之后出现的少数极端样品,可能是预处理造成的假象而非化学特征,在解释或据此剔除样品之前应先加以核查 [6]。二者还有一个共同的陷阱:它们都假定化学信号只占光谱的一小部分,且物理效应在各波长上相同。当被测组分主导光谱变异时,校正可能会去除部分被测组分信号;Martens、Nielsen 和 Engelsen 在小麦面筋与淀粉粉末混合物的透射光谱上用标准 MSC 恰好发现了这种情况,它去除了过多的被测组分信息 [7]。
扩展 MSC
扩展乘性信号校正(EMSC)由 Martens 和 Stark 于 1991 年提出,旨在将先验知识引入预处理:它消除不受控光程和光散射的影响,其配套方法光谱干扰扣除则去除已知干扰物;两者均以苯–二甲苯混合物中的甲苯、并施加了显著的模拟散射效应进行了演示 [8]。Martens、Nielsen 和 Engelsen 后来扩展了该模型,可同时校正加性、乘性和波长依赖的散射效应;经 EMSC 处理后,他们的面筋–淀粉粉末光谱表现得与透明溶液光谱非常相似,并遵循 Beer 定律 [7]。当普通 MSC 或 SNV 明显扭曲被测组分谱带,或已知散射效应随波长变化时,应考虑采用 EMSC。
实践中如何选择
不存在唯一正确的选择;文献支持的是一个简短的流程,而非某种默认设置:
- 由模型而非肉眼决定。将 MSC、SNV、SNV 加去趋势以及 EMSC 作为候选步骤,分别与拟采用的导数设置组合,并按交叉验证误差进行比较,正如 Rinnan 等对常用方法所做的那样 [3]。
- 将参考光谱与模型一并保存。MSC 参考光谱在预测时必须原样应用;若定标样品总体可能发生变化,SNV 可完全避免对样品集的依赖 [4]。
- 校正后检查得分图。SNV 之后的弯曲结构和 MSC 之后新出现的离群点是变换本身已知的行为,未必反映样品 [6];若被测组分主导光谱,应考虑 EMSC [7]。
caliX 中的散射校正
在 caliX 软件套件中,MSC、SNV、去趋势和导数滤波器在模型开发的预处理步骤中配置:在确定 PLS 模型之前,通过交叉验证比较各候选组合,MSC 参考光谱随模型一并存储。同一工作流程适用于来自 MasterLine17EC 谷物和饲料原料近线分析仪的光谱,也适用于来自 ProLine17EC 在线分析仪的光谱,后者的粒度和填充状态随工艺而变化。哪种校正适合特定物料,需通过使用您自己样品的可行性研究来确定。
参考文献
- [1] Barnes, R.J., Dhanoa, M.S., Lister, S.J. (1989). Standard Normal Variate Transformation and De-trending of Near-Infrared Diffuse Reflectance Spectra. Applied Spectroscopy, 43(5), 772–777. doi:10.1366/0003702894202201
- [2] Geladi, P., MacDougall, D., Martens, H. (1985). Linearization and Scatter-Correction for Near-Infrared Reflectance Spectra of Meat. Applied Spectroscopy, 39(3), 491–500. doi:10.1366/0003702854248656
- [3] Rinnan, Å., van den Berg, F., Engelsen, S.B. (2009). Review of the most common pre-processing techniques for near-infrared spectra. TrAC Trends in Analytical Chemistry, 28(10), 1201–1222. doi:10.1016/j.trac.2009.07.007
- [4] Dhanoa, M.S., Lister, S.J., Sanderson, R., Barnes, R.J. (1994). The Link between Multiplicative Scatter Correction (MSC) and Standard Normal Variate (SNV) Transformations of NIR Spectra. Journal of Near Infrared Spectroscopy, 2(1), 43–47. doi:10.1255/jnirs.30
- [5] Isaksson, T., Næs, T. (1988). The Effect of Multiplicative Scatter Correction (MSC) and Linearity Improvement in NIR Spectroscopy. Applied Spectroscopy, 42(7), 1273–1284. doi:10.1366/0003702884429869
- [6] Fearn, T., Riccioli, C., Garrido-Varo, A., Guerrero-Ginel, J.E. (2009). On the geometry of SNV and MSC. Chemometrics and Intelligent Laboratory Systems, 96(1), 22–26. doi:10.1016/j.chemolab.2008.11.006
- [7] Martens, H., Nielsen, J.P., Engelsen, S.B. (2003). Light Scattering and Light Absorbance Separated by Extended Multiplicative Signal Correction. Application to Near-Infrared Transmission Analysis of Powder Mixtures. Analytical Chemistry, 75(3), 394–404. doi:10.1021/ac020194w
- [8] Martens, H., Stark, E. (1991). Extended multiplicative signal correction and spectral interference subtraction: new preprocessing methods for near infrared spectroscopy. Journal of Pharmaceutical and Biomedical Analysis, 9(8), 625–635. doi:10.1016/0731-7085(91)80188-F