基于PDS算法的不同光照下模型传递研究

基于PDS算法的不同光照下模型传递研究作者：张文君唐红来源：《湖北农业科学》2017年第05期摘要：针对在室外光照对样品使用近红外光谱检测带来误差的问题，提出基于模型传递来减少检测误差的方法。

以圆黄梨为样品，分析样品在室内、室外阴影下的近红外光谱，建立室内光谱的PLS模型。

采用分段直接校正（Piecewise Direct Standardization，PDS）算法，减小室内外光谱差距，使得室内PLS模型能预测室外光谱。

结果表明，在室内建立的模型能预测经PDS算法传递后的室外光谱，预测决定系数（R2）和标准差（Root Mean Square Error of Prediction，RMSEP）分别为0.64和0.565 27，能有效地解决室外光照对光谱检测影响的问题。

关键词：近红外光谱；分段直接校正；光照影响；模型传递；糖度中图分类号：O657.33 文献标识码：A 文章编号：0439-8114（2017）05-0969-04DOI：10.14088/ki.issn0439-8114.2017.05.046Research on Model Transfer Based on the Piecewise Direct Standardizationat Different IlluminationZHANG Wen-jun， TANG Hong（College of Metrology Measurement Engineering， China Jiliang University， Hangzhou 310018，China）Abstract： For an error problem caused by near infrared spectrometer was used to detect sample at outdoor， a method was proposed based on model transfer to reduce the measurement error. Using the Wonhuwang pears samples， which analysis near infrared spectrum of samples of indoor and outdoor shade by near infrared spectroscopy（NIRS） and using an indoor spectrum by establishing a partial least squares（PLS） model. This experimental was using the Piecewise Direct Standardization （PDS） algorithm to reduce the gap between the indoor and outdoor shade spectrum. The model could effectively predict the outdoor spectrum. The experimental results show that the coefficient of determination R-squared of the prediction set was 0.64 and the root mean of square error of prediction（RMSEP） of the prediction was 0.565 27， it could effectively solve the measurement error using near infrared spectrometer at outdoor.Key words： NIRS； PDS； illumination； model transfer； sugar中国是水果产量大国，但水果出口率仍然较低，而导致出口率低下的主要因素是对水果处理技术的不成熟，水果质量的管控没有得到很好的保障。

目前，水果检测技术采用较多的是近红外光谱技术，该技术具有快速、绿色、无损等特点，但由于影响近红外光谱仪检测准确性的因素较多，更多的应用在室内恒定的条件下，在室外应用的结果不理想。

在室外检测水果，太阳光的照射会改变入射光谱，使得水果漫反射光谱发生改变，从而影响检测结果。

为了解决室外的光照对光谱仪检测的影响，国内外学者做了很多研究。

Saranwong等[1]使用便携式近红外光谱仪在果园对芒果进行了光谱检测及分析，发现黑色袋子能消除部分光照的影响，但该方法适用范围有限，且对实际应用带来不便。

吴方龙等[2]在室内利用灯光来模拟室外光照对梨进行了光谱检测，用预处理方法消除了光照影响，但该试验不能很好地反应实际应用情况。

近年来，近红外光谱仪在使用中经常会遇到测量环境发生改变时，模型无法共用的问题。

很多学者提出采用模型传递来解决这个问题。

模型传递可以减少由于测试环境的变化需要重建模型的麻烦。

如，林振兴等[3]将PDS算法用于检测不同温度下的喷漆燃料光谱，有效地解决了温度对喷漆燃料光谱检测带来的影响。

本研究以圆黄梨为样本，对其糖度进行了室内、外不同光照条件下的近红外光谱检测，采用PDS算法对室外光谱进行模型传递，使得室内模型能很好地预测室外光谱，以期为近红外光谱仪在室外的应用提供参考。

1 模型传递原理与PDS算法随着化学计量学及光谱学的快速发展，近红外检测技术随之进入一个新的时代，它已成为一种快速和高效的分析技术，成功应用在食品[4]、药物[5]、工业[6]等行业。

这种分析技术为多元校正方法，即对样品的自变量和因变量之间的建立函数关系（建立数学模型），该过程为校正。

预测过程是通过自变量数据和已建模型来计算因变量数据[7，8]。

该分析技术在光谱分析中也被称为模型传递（Model Transfer），基于光谱模型传递的校正过程有DS（Direct Standardization，直接校正）算法、PDS（Piecewise Direct Standardization，分段直接校正）算法和Shenk’s算法[9]。

以下介绍PDS算法和传递模型稳健性的评价。

1.1 模型传递模型传递是为了解决样品所建模型在不同仪器上不能共用的问题，即在某一仪器（称源机，Master）上建立一个多元校正模型，而该模型在另外一台与源机相同型号的仪器（称目标机，Slaves）上无法使用，或者结果产生较大的偏差，通过模型传递可以消除仪器之间的偏差。

同理，在不同光照环境下，对光谱分析也可以采用模型传递来消除光照环境对样品检测的影响。

首先，利用近红外光谱仪测量样品的光谱。

即室内光谱采用主仪器光谱（Master），室外所测光谱采用从仪器光谱（Slaves）进行模型传递。

然后，采用PDS算法来消除或减小主、从仪器间光谱的差距。

1.2 PDS算法PDS算法是一种多元全光谱标准化方法[10]。

它通过转移矩阵F将室内光谱矩阵Rm（主仪器光谱）转换成室外光谱矩阵Rs（从仪器光谱），减小Rm与Rs之间的差异，使二者匹配。

其中矩阵F为室内、外光谱下的光谱响应差异。

PDS算法的具体步骤：1）对应于Rms，i，在从仪器标样光谱矩阵Rss选取窗口为（j+k+1）大小的光谱段Rss，j+k+1，（from（i-j）th to （i+k）th wavelengths），组成矩阵Xs，i=[Rss，i-j，Rss，i-j+1，…，Rss，j+k-1，Rss，j+k]，i代表波长点，j、k分别代表i波长点左右不同的窗口宽度。

2）将Rms，i与Xs、i进行关联，Rms，i=Xsbi，转换系数可由PCR或PLS方法求出。

3）循环i，求出所有的bi，i=1，2，…，m。

4）对未知样品光谱rs，un经固定窗口分段，由转换系数bi循环得到rm，un相一致的光谱r ■■。

对于rs，un两端，窗口大小的波长范围（i～j和m-k～m）不能转换，一般舍去，也可通过外推法获得。

在本研究中，数据点两端通过样条插值获得。

1.3 KS（Kennard-Stone）算法在模型传递效果中，为了克服测量环境不同而引起的误差，必须使转换集在主仪器和从仪器所测信号中包含足够多的信息；同时，在主仪器上选择转换集的效果最好。

KS算法是由Kennard-Stone提出，选择标准样品的一种选取转换集的有效和广泛应用的方法，它是根据光谱主成分之间欧式距离d来决定的，从而在光谱特征中均匀地选取标准样品。

dij=■式中，dij为样本i和样本j之间的欧式距离，p为样品光谱的波长个数。

1.4 模型传递的评价为了比较建模效果质量的好坏，可以通过决定系数R2（相关系数R）、校正均方根误差（RMSEC）和预测均方根误差（RMSEP）等来判断。

R2越大越接近1，均方根误差越小，越接近0建模效果越好。

R2=1-■RMSEC=■RMSEP=■其中，n为校正集样品的数目，f为主成分数目，Ip为预测集样品的数目，yi为样品测量的实际值，y为样品的平均值，■为预测样品集的预测值。

2 材料与方法2.1 材料与处理试验样品均采自浙江省杭州市果园，选择大小、颜色相近的圆黄梨150个，样品按2∶1分为校正集和预测集。

将样品擦拭干净并按序标记，在每个样品赤道处均匀取4个点作为测量位置，将样品保存在室内恒定温度下，直至样品表面温度与室内温度相同。

表1是校正样本和预测样本的糖度值。

从表1可知，样品的糖度校正集15.11 °Brix和预测集15.69 °Brix近似相等，样品预测值在正常的校正样品范围内。

2.2 试验仪器MCS600型德国蔡司光谱仪，并携带90 °垂直积分球进行测量，光谱范围300～1 700 nm，分辨率为3 nm，光谱采集软件为仪器自带软件Aspect plus 2.0；理化值测定使用ATAGOPR-101α型数字式糖度计，测定范围为0～45 °Brix，误差为±0.1 °Brix；光照度计使用希玛AR823数字照度计，测量范围是1～200 klx，分辨率为1 klx；数据处理及建模软件为MATLAB2012b。

2.3 试验设计将样品进行预处理后对样品进行不同光照条件下的光谱检测。

为了减小温度、湿度等因素对检测的影响，先在室内保存24 h以上再进行检测。

在室外检测时，为了避免太阳长时间照射改变样品表面温度，快速检测样品的近红外光谱。

在室内恒定条件下对样品采用漫反射光谱检测，对每个样品采集点进行采集，取均值作为该样品的最终光谱；在室外光照条件下进行上述试验，作为该样品光照下的光谱。