尺寸差异对脐橙糖度可见近红外光谱检测模型影响

刘燕德，饶 宇，孙旭东，姜小刚，徐 海，李 雄，王观田，徐 佳

华东交通大学机电与车辆工程学院，江西 南昌 330013

引 言

糖度(SSC)是影响水果口感的一个重要因素，SSC差异在2°Brix以上，水果口感差异明显[1]。近年来，人们对口感品质十分重视，传统外部品质分选难以满足口感品质分选需求，可见/近红外光谱技术作为一种快速无损检测手段，目前已经广泛应用于苹果[2]、梨、桃[3]等水果的内部品质检测。

崔丰娟[4]等研究了苹果动态在线检测的传输速度对模型的精度的影响。曹松涛在他的学位论文中报道了对梨样品相关因素对内部品质普适模型的稳健性影响，样品的尺寸差异、产地差异、品种差异都是主要影响，通过对光谱进行一阶、二阶微分预处理，基本能够消除差异影响，提高了模型的预测精度。Saranwong[5]等对比分析了两种不同近红外漫反射仪器在芒果品质检测上的应用。Greensill[6]等使用自行研发的光学系统，研究了硬皮甜瓜果实不同果肉位置的SSC取样准确度，结果表明，果实赤道部位SSC取样值与光谱建模模型准确度最高，能代表水果的真实糖度值。Mizuki[7]等结合偏最小二乘判别分析，应用苹果分选机顶部和底部的光谱仪，建立了苹果内部褐变发生的判别模型。

本工作主要研究了脐橙的尺寸差异对可见近红外漫透射及两种不同漫反射检测方式模型的影响，对比分析了大果、小果及混合果模型的预测精度，考察了近红外光谱模型及多元散射校正预处理方法对消除尺寸差异的效果。

1 实验部分

1.1 材料

实验所用的脐橙样品购于南昌市某水果市场，将脐橙表面擦净，在25 ℃的室温下静置24个小时，依次采集脐橙样品的横纵径、重量、糖度的理化指标如表1所示。实验共采用480个样品，按照脐橙的横纵径指标划分大小果，其中横径范围在75～90 mm，纵径在80～95 mm为大果共312个，横径范围在60～75 mm，纵径在65～85 mm为小果共168个，大果建模集为234个，预测集为78个，小果建模集为126个，预测集42个，混合果建模集360个，预测集120个。

表1 脐橙的理化指标范围统计结果Table 1 Statistical results of physical andchemical indicators of navel orange

1.2 装置和光谱采集

可见近红外漫反射光谱采集平台如图1所示，实验平台搭建在重量分选装置上的出果口前，考虑实际应用场景，光源发射与接收探头安装在距离果杯15 cm处。采用Ocean Optics公司的QE PRO光谱仪作为实验的光谱采集装置，其波长范围为350～1 150 nm，光源为1个12 V、100 W的卤钨灯，采用稳压电源供电。光纤探头分别为环形发射与接收和多点发射与接收漫反射式光纤探头，均设计为将发射光源与接收光源装置集成在同一探头上，通过卤钨灯提供直射光源，从光源接收端口进入光纤，由探头发射端发射出光照，光信号经过样品内部后，经过吸收、反射等，由探头接收端口接收到漫反射光信号，经过光纤传输并触发光谱仪进行数据保存，每触发一次，光谱采集软件便储存一条光谱。

图1 近红外漫反射在线检测装置1： 光源； 2： 光源接收端口； 3： 光源发射与接收端口；4： QE PRO光谱仪； 5： 电脑； 6： 脐橙样品； 7： 果杯Fig.1 Near-infrared diffuse reflection on-line detecting device

开机预热0.5 h，聚四氟乙烯球作为光谱参比，采集脐橙样品光谱，实验参数为： 积分时间100 ms，速度5个·s-1，光强100 W。光谱采集软件采集光谱，分别标记面光谱，取平均光谱作为该样品可见近红外漫反射光谱。

漫透射光谱采集装置如图2所示，采用海洋光学公司的QE65000光谱仪，波长范围为350～1 150 nm，光源为10个12 V、100 W的卤钨灯，均匀分布在传送带两侧，通过稳压电源供电。光纤探头位于果杯正下方，光信号从样品表面进入，经过透射到达样品下方的光纤探头，探头每接收到光信号一次，光谱仪触发一次，则光谱采集软件保存一条光谱。果杯设计有弹性避光圈与透光孔，避免杂散光与外界干扰光被光纤接收到。

开机预热0.5 h，聚四氟乙烯球作为光谱参比，采集脐橙样品光谱，实验参数为： 积分时间100 ms，速度5个·s-1，光强1 000 W。采用光谱采集软件采集光谱，分别标记面光谱，取平均光谱作为该样品可见近红外漫反射光谱。

图2 近红外漫透射在线检测装置1： 光源； 2： 样品； 3： 遮光圈；4： 光路； 5： 光纤探头； 6： 弹跳爪Fig.2 Near-infrared diffuse transmissionon-line detecting device1： Light source； 2： Sample； 3： Shading ring；4： Optical path； 5： Fiber optic probe； 6： Bouncing claw

1.3 糖度含量判定

脐橙糖度含量由折射式数字糖度计(ATAGO PAL-α，日本ATAGO公司)测定，测量前，用纯净水洗净并擦干糖度计，并进行标定，糖度值为0°Birx，测量时，使用一次性纱布将标记部分约5 mm深果肉滤出少量果汁，滴于糖度计的测试窗口，重复3次，取3次的平均值作为糖度真实值，所得值以质量分数表示。

2.4 数据处理

采用matlab2012a将光谱数据转换为数据矩阵，再通过Unscrambler进行PLS建模。建模前先进行光谱数据预处理，再使用偏最小二乘法(PLS)分别建立漫透射及两种不同漫反射式的糖度模型。偏最小二乘法中的相关系数r和预测集均方根误差RMSEP、建模集均方根误差RMSEC计算公式如式(1)—式(3)所示

(1)

(2)

(3)

2 结果与讨论

2.1 不同检测方式的尺寸差异影响分析

近红外光谱分析是一种应用于水果内部品质无损检测的重要手段，但由于样品的相关因素如样品温度、大小、产地、品种等的影响，模型的稳健性受到较大影响，其中，样品大小是一个重要的影响因素。近红外漫透射和漫反射检测方式不同大小样品光路图如图3所示，其中图3(a)为漫透射检测大、小尺寸的脐橙光路图，在相同光照路径下，光在大果中的光程要比小果中的光程更长，光的能量衰减严重，光纤探头接收到的能量强度更低，而小果由于它的小尺寸，光能快速透过，光纤探头接收到的能量衰减少，强度强，携带的有效信息多。图3(b)为漫反射检测大、小尺寸的脐橙光路图，在相同光照路径下，光在大果中的光程要比小果中的光程要短，可能穿过的果肉组织不多，携带的有效信息少，而在小果中，光程较长，穿过的果肉组织较多，携带的有效信息也多。

图3 漫透射(a)和漫反射(b)上大、小尺寸的脐橙光路示意图

2.2 不同检测方式的光谱对比分析

对样品的尺寸进行测量，每个样品的最大横径测量5次，取平均值作为样品的横径指标，考察尺寸差异对脐橙糖度在线检测精度的影响。由图4可知，无论是漫透射检测还是漫反射检测，小果的光谱能量强度都要比大果强。对同一样品的三种不同检测方式的光谱进行对比分析。两种漫反射光谱的能量要比漫透射光谱强，而环形发射与接收漫反射光谱又要比多点发射与接收漫反射光谱的能量强。漫反射光谱在690 nm左右存在吸收峰，漫透射光谱的峰在700和825 nm左右，从两种检测方式的光谱波形上来看，漫透射光谱要比漫反射光谱多一个吸收峰，且漫透射光谱集中在600～850 nm范围内，漫反射光谱范围则在500～900 nm内，比漫透射光谱范围大一点。对比分析三种不同检测方式，漫透射检测方式中，光源位于样品两侧，接收装置位于样品下方，果杯上的橡胶遮光圈能有效消除样品表面反射光及其他干扰光进入接收器； 在漫反射检测方式中，光源与接收装置集成在同一光纤上，环形发射与接收的漫反射检测方式中，由于它的环形光源的结构特点，光在样品表面形成一个圆形光斑，可以有效的将外界干扰光及样品表面反射光隔离在光斑外面，达到克服杂散光的目的； 多点发射与接收漫反射检测方式中，由于它是点光源，点接收方式，所以它克服样品表面反射光和其他干扰光的能力要比前两种检测方式稍差。

图4 三种检测方式的大小果光谱对比Fig.4 Spectral comparison of fruit with differentsizes using three detection methods

2.3 不同检测方式下尺寸预测模型建立

不同样品尺寸在某一波长范围光下具有一定的衰减规律，实验选用波长范围为900～936 nm建立近红外光谱尺寸预测模型。可见/近红外光谱波长范围320～500 nm是可见光波段，受颜色变化影响较大； 500～900 nm是光谱的有效信息集中波段，用来建立糖度预测模型； 936～1 150 nm波段噪声和无效信息较多，故选用900～936 nm建立尺寸预测模型。不同检测方式下，尺寸模型预测散点图如图5所示。由图5可知，漫透射检测方式下尺寸模型的预测集相关系数为0.60，预测集均方根误差3.95°Brix，截距B=81.27； 多点发射与接收漫反射检测方式下尺寸模型的预测集相关系数为0.97，预测集均方根误差1.46°Brix，截距B=94.3； 环形发射与接收漫反射检测方式下尺寸模型的预测集相关系数为0.96，预测集均方根误差1.73°Brix，截距B=70.73。由不同的检测方式下尺寸预测散点图可知，漫透射检测方式的尺寸模型预测精度较低，两种漫反射检测方式的尺寸模型预测精度较高，可以有效的预测样品的尺寸大小。

图5 不同检测方式下尺寸模型预测结果散点图(a)： 漫透射检测方式尺寸预测散点图；(b)： 多点发射与接收漫反射检测方式尺寸预测散点图；(c)： 环形发射与接收漫反射检测方式尺寸预测散点图Fig.5 Scatter plot of dimensional model predictionresults under different detection modes

2.4 糖度预测模型建立

王动民(本刊34卷9期)利用MSC方法对光谱进行预处理，基于不同波段的光谱数据会得到不同的校正结果，进而影响所得定标模型的可靠性。分别建立三种不同检测方式的不同果型的PLS模型，模型预测结果如表2所示，通过预测集相关系数Rp来评价模型的预测能力。由表2可知，三种不同的近红外光谱检测方式下，小果的PLS模型预测能力要强于大果和混合果，其中漫透射方式的Rp为0.76，多点发射与接收方式的Rp为0.72，环形发射与接收检测方式的Rp为0.72。在相同光照路径下，光在大果中的光程要比小果的长，衰减严重，接收到的光谱能量强度低，因此携带的有效信息比小果光谱的要少很多，所以小果的模型预测精度相对较高。

表2 不同尺寸样品集的不同检测方式PLS建模集预测结果Table 2 PLS modeling prediction results for different size sample sets by using different detection methods

糖度预测模型公式如式(4)所示。

(4)

其中，n为参与建模的光谱变量数，β为糖度相关能量谱强度，b为模型的截距。

2.5 脐橙糖度最佳模型建立和预测

通过以上多种模型的对比分析，漫透射检测方式模型预测效果最好，多元散射校正预处理后，模型预测效果有所优化，可以消除尺寸差异带来的影响。采用多元散射校正预处理漫透射光谱，再利用偏最小二乘(PLS)建立漫透射及漫反射的预测模型如图6所示。由图6可知，经过MSC校正的漫透射光谱预测模型散点距离拟合线较为集中，预测精度较高，建模相关系数为0.88，预测相关系数为0.84，预测集均方根误差RMSEP为0.64，截距b=14.05。回归系数与近红外漫透射大小果光谱曲线对比如图7所示，光谱变量在PLS定量模型中贡献率越大，回归系数越大。所用的光谱变量与回归系数加权求和再加上截距，及PLS模型预测值，实现了对样品的糖度值的预测。

在漫反射检测方式中，先进行尺寸分选，区分出大小果后，再进行大小果的糖度分选方式，也可以消除尺寸差异带来的影响。

图6 漫透射MSC-PLS建模散点图Fig.6 Diffuse transmission MSC-PLS modeling scatter plot

3 结 论

研究了样品尺寸差异对脐橙漫透射及漫反射在线无损检测精度的影响。分析了脐橙的尺寸差异在漫透射和漫反射装置上的不同影响以及脐橙三种不同检测方式的近红外光谱特性。分别建立了三种检测方式下的大果、小果、混合果、经MSC预处理后的混合果糖度PLS模型和尺寸PLS预测模型，其中小果的糖度PLS模型预测效果均最佳，在漫透射糖度预测模型中，经MSC预处理后，混合果预测模型精度有所提高，预测效果要优于小果，模型的预测相关系数为0.84，预测均方根误差为0.64°Brix； 在不同检测方式下的尺寸预测模型中，漫透射检测方式要比两种漫反射检斥方式的尺寸模型预测精度低，在漫反射检测方式下，先采用尺寸模型进行尺寸进行分选，然后再使用糖度模型进行糖度分选，最终实现不同检测方式下不同尺寸脐橙的同时在线分选。论文研究为不同果型的水果糖度同时在线无损检测和分选提供了参考和理论支持。

图7 回归系数与漫透射大小果光谱对比Fig.7 Comparison of regression coefficientand diffuse transmission size