基于高光谱图像技术的油炸薯片中羧甲基赖氨酸含量检测

王润博，韩文凤，于慧春

（1.河南工业大学漯河工学院，河南 漯河 462000；2.浙江工业职业技术学院鉴湖学院，浙江 绍兴 312000；3.河南科技大学食品与生物工程学院，河南 洛阳 471023）

羧甲基赖氨酸[（Nε-（1-carboxymethyl）-L-lysine，CML]是晚期糖基化终末产物代表之一，当达到一定数量时能促进糖尿病、肾病、动脉粥样硬化等疾病的发展和人体器官的快速衰老[1]。周燕琼[2]发现油炸薯片中CML含量可达22.2 mg/kg，较同类食品高，但实际生活中未将CML含量作为评价油炸薯片的指标[3]。很重要的一方面原因是常见的CML检测方法各有弊端[4-5]，而具有灵敏度高、可利用性好、稳定性较强等特点的液相色谱-质谱法，能对食品中CML进行准确的定量分析，但其操作繁琐、耗时耗力、价格昂贵等。因此，实现快速、高效地检测低剂量CML，对避免食用高含量CML的油炸薯片，保证消费者健康及马铃薯产业的发展具有重要意义。目前，国内外已有对油炸薯片中有害成分[6-7]和吸油量[8-9]的分析，也有对CML的抑制研究[10]，但少有对CML快速检测方法的探讨。

高光谱图像技术具有分辨率高、曲线连续、无损伤、无污染、可对物质特征或成分进行精确探测等特点[11]。近年来，高光谱技术在城市测绘、精确农业、环境监测、军事监视以及计算机视觉等工业领域得到了普遍的认可和关注[12]。并且，由于高光谱光谱反射率与理化参数的相关性，可建立各种理化参数的光谱监测模型，所以在测定食品或农产品中重要营养物质、有害成分含量等方面的研究较多[13-17]，在监测食品或农产品质量方面取得重大进展[18-21]。国内外也有利用高光谱成像快速检测鲜切马铃薯片或块茎的相关研究[22-26]，但鲜见采用高光谱成像技术预测油炸薯片中CML含量的报道。

本文采用高光谱图像技术，以具有不同工艺参数的油炸薯片为研究对象，借助液相色谱-质谱法测定CML含量，构建3种预测模型并进行对比。最终确定最佳快速稳健模型，为寻找高效的CML低剂量检测设备提供依据和手段，为更好地服务于食品企业，保证产品质量。

1 材料与方法

1.1 材料与试剂

1.1.1 材料

试验所用样品为自制油炸薯片，采用如下制作过程：马铃薯粉加水1∶1（质量比），和匀，将其揉成长面团，然后做成厚度为1.0 cm的面块，再熟化、冷却，最后切为厚1.5 mm的小片，油炸。试验采用不同的油炸温度和时间如表1所示，共8个等级。本研究采用粉碎后的油炸薯片作为样品。

表1 油炸温度和时间Table 1 Frying temperature and time

1.1.2 试剂

CML标准品（纯度为98.0%）和CML-D4同位素内标（化学纯度为98.0%，同位素纯度为97.9%）：加拿大TRC公司；乙腈、甲酸（均为质谱纯）：美国Thermo Fisher Scientific公司。

1.2 仪器与设备

高光谱成像采集系统如图1所示，主要由1台高光谱摄像仪（德国 Inno-Spec，IST50—3810型）、4个500 W的光纤卤素灯（德国ESYLUX 90000420108型）、1个传送装置（自制）和1台计算机组成。高光谱摄像仪通过USB 2.0接口连接计算机，软件（SICap-STVR V1.0.x）作为驱动控制摄像仪的平台，电脑实时记录和存贮数据信息。

图1 高光谱图像采集系统Fig.1 Hyperspectral image acquisition system

OasisMCX固相萃取柱（体积3mL，填料质量60mg，填料粒径 30 μm）、Milli－Q Advantage A10 超纯水系统：美国Millipore公司；VacElut SPS 24真空固相萃取装置：美国Agilent公司；N－EVAP－24圆形水浴氮吹仪：美国Organomation公司；Acquity超高效液相色谱仪、UPLC BEH Amide HILIC亲水作用色谱柱（2.1 mm×100 mm，1.7 μm）、XEVO TQ 三重四极杆串联质谱仪：美国Waters公司。

1.3 试验方法

1.3.1 CML含量的测定

液质联用测定油炸薯片中CML含量与高光谱图像采集同步。首先对薯片粉末进行脱脂和水解，然后收集上清液，最后按照文献[1]的方法，采用液相色谱-质谱法测定。

1.3.1.1 液相检测条件

流动相：A相为0.1%甲酸水溶液，B相为乙腈；进行梯度洗脱，梯度洗脱参数如表2所示。柱温：35℃；进样量：5 μL；运行时间 7 min。

表2 梯度洗脱参数Table 2 Gradient elution parameters

1.3.1.2 质谱检测条件

离子源：电喷雾正离子模式（electrospray ionization+，ESI+）；监测模式：多反应监测（multiple reaction monitoring，MRM）；毛细管电压：3.5 kV；锥孔电压：20 V；源温度：150℃；脱溶剂气温度：400℃；脱溶剂气流速：700 L/h；碰撞能量：15 V；MRM 模式：CML 质荷比m/z 205.22～m/z 84.00定量，m/z 205.22～m/z 130.00定性，CML-D4 m/z 209.00～m/z 87.70定性。

1.3.2 高光谱图像的采集

样本高光谱图像采集前，首先采集标准全白图像和全黑图像，然后再将已备好的（40±0.5）g油炸薯片粉末，作为一个待测样本，均匀地平铺在直径为100mm的培养皿中，放置在传送带上，传送装置的速度为1.20mm/s，相机的曝光时间为90ms，物距为350mm。光谱仪的光谱范围为371.05 nm～1 023.82 nm，光谱分辨率为2.8 nm，采样间隔为0.51 nm。采集数据时，摄像头在传送带方向的垂直方向作横向扫描，得到的是所扫区域全部像素点在1 288个波长处的图像数据。同时随着薯片的前进，完成所需区域油炸薯片图像数据的采集。

为获得理想区域的图像，设定高光谱扫描程序中Width和Height的参数都是800。即图像扫描行数为每幅800行，每行扫描的像素点数为800个，得到的高光谱图像分辨率为800像素×800像素。所以，对于每一个油炸薯片样本，采集后最终得到一个大小为800×800×1 288的高光谱图像数据块。最后，为了消除成像仪的暗电流和光源强度分布不均造成的影响，需要对每个油炸薯片样本高光谱图像进行黑白标定[8]。采用标准全白图像和标准全黑图像对采集到的原始高光谱图像数据块按照公式进行标定，得到黑白校正后的高光谱图像数据。黑白校正公式如下。

式中：R为黑白校正后的图像；L为原始反射图像；W为全白反射图像；B为全黑反射图像。

1.3.3 模型建立方法与数据处理方法

本研究采用主成分回归[27]、偏最小二乘回归[28]和BP神经网络[29]3种方法建立模型。平均光谱数据的提取借助软件ENVI5.1；黑白处理、标准散射校正预处理及模型的构建均在MATLAB R2014a中编程实现。

2 结果与分析

2.1 图象的采集与光谱预处理

试验中，每类油炸薯片分别选取50个重复样本，8类样品共400个样本。相机拍摄各类油炸薯片图像如图2所示，在可见光波长658.77 nm下的细节图像如图3所示。

图2 8类样本的可见光图像Fig.2 Visible light images of eight samples

选择标准正态变换（standard normal variate，SNV）作为光谱预处理方法，用来消除油炸薯片固体颗粒大小不一、表面不均匀产生的散射以及光程变化对漫反射光谱造成的影响。400个油炸薯片样本在处理前后的光谱曲线图如图4所示。

图4 400个油炸薯片样本的光谱曲线Fig.4 Spectral curves of 400 samples

图4a是400个样本在1 288个波段的原始平均光谱反射图。由于首尾波段信噪比较低，所以除去371.05 nm～472.86 nm和880.49 nm～1 023.82 nm的数据，得到从第200到1000个波段的信息，即剩余801个波段的光谱值。图4b是所有样本在801个波段的原始光谱值经过黑白校正后得到。图4c是图4b经过标准正态变换后的光谱图。很明显，经过SNV处理过的光谱曲线具有清晰的吸收峰，且每类薯片的光谱各有差别，这有利于预测模型的建立。

2.2 油炸薯片中CML含量的测定结果

每种油炸薯片样品取3个平行进行测定，并确定这3次独立测定结果的绝对值差小于等于其算术平均值的10%，将其平均值作为此类油炸薯片样品的实际CML含量值，结果如表3所示。

表3 8种油炸薯片样品的CML含量测定结果Table 3 Content of CML in eight kinds of fried potato chips

从表3看出，当温度在180℃时，随着时间的延长，CML含量呈现先增多后减少的趋势，160 s达到最大值为12.74 μg/g，而后随着时间的延长CML含量有所下降。当温度220℃、时间240 s时，CML含量比80 s时明显降低，仅为5.17 μg/g。这与韩文凤等[1]的研究结果相一致：反应温度对体系中CML含量的影响较为复杂；在高温加热的反应体系中，其含量随着反应时间的延长呈先增大后减小的趋势。采用此液质联用检测方法，耗时费力，成本高，且有毒。以下采用高光谱图像技术提取数据，借助理化检测方法测定的结果建立不同模型，预测CML含量。

2.3 采用不同方法建立预测模型

随机选择全部样本的70%作为训练集，剩余30%样本作为测试集。本研究选用预测正确率（相对误差不超过真实值的5%为预测正确）、决定系数和均方根误差3个结果指标对3种建模方法的精度进行评价。

2.3.1 主成分回归

前15个主成分的贡献率见表4，不同方法建立预测模型结果对比见图5。

表4 前15个主成分得分贡献率Table 4 Contribution rates of top 15 principal components

由表4可见，前15个主成分得分可以呈现99.999%以上的原始光谱信息。因此，选取前15个主成分作为样本集特征光谱，即将801个波长压缩为15个新变量作为特征变量，明显提升了模型的运行效率。

由图5a可知预测效果不理想，主成分回归预测正确率为30.83%，决定系数为0.85，均方根误差为1.3。

2.3.2 偏最小二乘回归

基于801个波段数的平均光谱反射值建立偏最小二乘回归模型，选择15个主成分，所获得的测试集决定系数为0.84，均方根误差为1.15，预测正确率为55.00%。如图5b所示，与主成分回归对比，预测正确率稍高，决定系数近似。这说明，此模型仍然不能满足预测油炸薯片中CML的含量。

图5 不同方法建立预测模型结果对比Fig.5 Comparison of different prediction models

2.3.3 BP神经网络

运用BP神经网络模型建立光谱值与油炸薯片中CML含量之间的映射关系，首先提取全部波段下图像的平均光谱反射值，然后去除首尾具有噪声的光谱，并对其进行标准正态变量变换光谱预处理。最后对801列光谱数据运用主成分分析融合，选取代表原始信息99.98%以上的前6个主成分作为神经网络的输入变量部分，设计模型结构，确定传递函数和训练函数，设置学习速率、误差、迭代次数等参数，建立具有稳健性和精确度的BP神经网络预测模型。

基于400个样本在第200到1 000个波段下各个高光谱图像的平均光谱反射值所建立的BP神经网络结构，预测模型的结构设为15-10-1，隐含层和输出层的传递函数分别为tansig和logsig，训练函数为traincgf。学习率设定为0.01，误差设定为1×10-6。运行结束，迭代次数共为4 969次，运行时间21 s，误差达到3.12×10-5。预测结果如图5c所示，所获得的测试集决定系数为0.99，均方根误差为0.22，预测正确率为99.67%。

2.3.4 对比结果

以建模集和预测集的决定系数、均方根误差、预测正确率作为3种建模方法，预测CML含量的评价标准，对比结果如表5所示。

表5 不同建模方法的CML含量预测结果对比Table 5 Comparison of CML content obtained in different ways

由表5可知，BP神经网络预测正确率较高，为预测油炸薯片中CML含量最优的模型。

2.3.5 最优预测模型的选择与稳健性分析

对比显示，基于第200到1 000波段下的每个油炸薯片样本图像的平均光谱反射值，预测CML含量建立的BP神经网络模型最优。为验证此模型的稳健性，从油炸薯片样本的高光谱图像数据信息中，随机抽取5组不同的训练集和测试集。保证模型结构及所有参数不变，用这5组数据分别对BP神经网络模型进行训练和测试，结果如表6所示。

表6 不同测试集的测试结果对比Table 6 Comparison of results from different test sets

从表6中可以看出，油炸薯片CML含量预测正确率平均值为96.23%，决定系数平均值为0.99，均方根误差平均值为0.22。这说明，该BP神经网络模型具有可靠的稳健性和精度。

3 结论

以含有不同CML量的油炸薯片为检测对象，用高光谱成像仪采集其在371.05 nm～1 023.82 nm波长下的图像，先提取其第200到1 000波段下图像的平均光谱反射值，然后经过黑白校正和标准正态变量变换光谱预处理，最后对比主成分回归、偏最小二乘回归和BP神经网络3种模型建立方法。结果显示，基于BP神经网络所建立的油炸薯片中CML含量预测模型最为准确。并且经过验证，该预测模型兼备准确性和稳健性。这说明，高光谱图像技术融合BP神经网络能够实现油炸薯片中CML含量的准确预测，可为食品中CML含量快速无损检测提供重要依据。