一种抛光粮食的快速无损鉴别方法
【技术领域】
[0001 ]本发明涉及农产品质量安全领域，具体地说，涉及一种基于振动光谱的矿物油抛 光粮食快速无损检测方法。
【背景技术】
[0002] 农产品是迄今为止人类赖以生存的食物重要来源。随着生产力的发展和社会的进 步，对农产品质量安全的关注日益提高。农产品质量安全包含两个方面的含义，其一是质 量，即营养成分的含量是否达标、风味是否良好;其二是安全，即是否存在危害或潜在危害 人类健康的物理、化学、生物及其他各方面可能的因素。质量良好的农产品可以给人类提供 优质的能量、营养来源，对保障人类生存、繁衍具有重要意义。然而近年来，不法商贩对劣质 农产品以次充好、造假制假等行为屡见不鲜。以大米为例，已有报道称不法商贩用矿物油抛 光陈旧大米，使其"油光可鉴"，以冒充优质大米出售给消费者。食用上述矿物油抛光大米对 人存在这极大的安全隐患，甚至直接危害消费者身体健康。并且，大米是我国大宗农产品， 拥有巨大的消费人群和消费市场。因此，针对矿物油抛光大米的快速无损检测成为亟待解 决的问题。
[0003] 矿物油是石油化工产品，以脂肪烃为主。由于其中几乎不含有芳香烃，因此采用传 统的荧光法很难对其进行检测。采用传统的感官评价方法则具有较强的主观因素，容易形 成误判。现代大型仪器方法虽然可以得到精准的检测结果，但是操作繁琐、对检测人员的专 业性要求高、检测过程需要使用化学试剂、检测周期长、检测效率低、难以适应大宗农产品 的快速检测需求。
[0004] 振动光谱是一种基于物质分子的振动、转动能级跃迀对电磁波产生吸收的光谱。 根据量子力学原理，当电磁波的能量和物质分子的振动、转动能级跃迀能量差相等，且满足 振动光谱发生条件时，会产生振动光谱吸收。常见的振动光谱有:近红外光谱、中红外光谱、 拉曼光谱、太赫兹波谱等。振动光谱技术具有采样速度快、分析效率高、绿色无污染等特点， 目前已在生命科学、医学及生理学、药物学、农学、化学等诸多领域具有广泛应用。
[0005] 本发明提出一种基于振动光谱的矿物油抛光粮食快速无损检测方法。该方法具有 快速、准确、高效、绿色无污染等特点，可以对矿物油抛光粮食进行快速无损鉴别，并且对保 障农产品质量安全、提高工农业生产效率、促进农产品市场健康发展、快速无损检测仪器的 研发等方面具有积极作用。
【发明内容】
[0006] 本发明针对矿物油抛光粮食的快速无损鉴别难题，提供一种基于振动光谱的矿物 油抛光粮食快速无损检测方法，其目的是对矿物油抛光粮食进行快速无损检测与鉴别。
[0007] 本发明提供的一种抛光粮食的快速无损检测方法，所述方法的流程可参考图1所 7Jn 〇
[0008] 具体而言，所述方法包括以下步骤：
[0009] (1)用已知的抛光粮食和未抛光粮食为样品，构建样本总组数为η的训练集;所述 训练集中，50,抛光样品与未抛光样品的组数之比为1~3.5:1;
[0010] (2)采集所述训练集中每组样品在8900(^1~WOOOcnf1波长范围内的近红外吸收 光谱，每组样品的近红外吸收光谱以lcnf1~3CHT1为间隔包括m个光谱数据；
[0011] (3)以所述训练集中全部样品的近红外吸收光谱数据构建nXm阶矩阵，对所述训 练集中的抛光样品和未抛光样品分别赋予不同的数值，依据所述数值用回归算法提取所述 矩阵的若干组特征向量并获得与所述特征向量一一对应的特征值，以特征值最高的两组特 征向量作为第1特征向量和第2特征向量；
[0012] (4)依据训练集中每组样品的近红外吸收光谱数据，求出分别与所述第1特征向量 和第2特征向量对应的权重值，作为各训练集样品的得分值1和得分值2,并以其分别为横、 纵坐标，绘制训练集样品的二维得分散点图（一组样品对应图中的一个点）；在该图上标记 一条抛光样品与未抛光样品的分界线，所述分界线两侧分别为抛光区域和未抛光区域；
[0013] (5)用与步骤(2)相同的方法采集待测样品的近红外吸收光谱数据，据此数据求出 分别与所述第1特征向量和第2特征向量对应的权重值，作为待测样品的得分值1和得分值 2；
[0014] (6)以待测样品的得分值1和得分值2分别为横、纵坐标，在所述二维得分散点图中 绘出对应的待测点，根据所述待测点所在的区域，鉴别待测样品是否为抛光粮食。
[0015] 本发明所述抛光粮食是指经矿物油抛光处理后的粮食，所述未抛光粮食是指未经 过任何矿物油抛光处理的粮食。具体而言，本发明所述未抛光样品中的矿物油含量为〇;所 述抛光样品中的矿物油含量为〇.4~0.6%。
[0016] 所述矿物油指石油裂解成分中的脂肪烃类物质，优选液体石蜡。
[0017] 所述粮食为优选为谷物，进一步优选为大米。
[0018] 本发明所述方法中，样品的单位为组，即训练集中包括η组样品。每组样品的取样 方法为:取相同品种、产地以及采收时间的粮食不少于3kg，采用堆锥四分缩分法获得一组 样品，每组样品的质量为200 ±5g。
[0019] 所述步骤（1)中，抛光样品和未抛光样品的组数可以相同，也可以不同；具体而言， 所述抛光样品数量与未抛光样品组数之比为1~3.5:1。为了确保检测结果的准确性，所述 训练集中，优选所述抛光样品和未抛光样品的组数均不小于40。本发明所述方法中，为了确 保检测结果的客观性，选取抛光样品和未抛光样品时，应确保其代表性;所述代表性具体是 指样品的品种、产地以及采收时间全面覆盖客观情况。
[0020] 本发明用于采集所述近红外吸收光谱数据的光谱仪分辨率为lcnf1~64CHT1，优选 为4cm- 1 ~16cm-1。
[0021] 所述步骤(2)中，每组样品近红外吸收光谱所含的光谱数据之间的间隔相等，具体 为lcnf1~3CHT 1，优选为2CHT1。所述步骤(2)还可以包括对所述近红外吸收光谱数据进行处 理，以使检测结果更加准确。所述处理包括平滑、微分、基线校正、数据标准正态化中的一种 或几种;优选为数据标准正态化。
[0022]具体而言，所述数据标准正态化处理，是针对每组样品的近红外吸收光谱数据，即 每一条光谱数据自身所做的标准正态化变换，其计算公式如公式(I)所示。
[0023]
(T)
[0024] 所述公式（I)中，X中是处理前的光谱数据，Xsnv是处理后的光谱数据;每组样品的 近红外吸收光谱数据中，即每一条光谱数据中，包含m个数据，经过处理后，得到m个标准正 态化处理后的数据，即SNV处理后的一条光谱数据。
[0025] 其中，ΧΜ是该条光谱数据处理前的平均值，其计算公式如公式（II)所示；S是该条 光谱数据处理前的标准差，其计算公式(III)所示。
[0026]
[0027]
[0028] 本发明步骤(3)提取所述特征向量的回归方法选自多元线性回归、主成分回归、偏 最小二乘回归、人工神经网络回归、支持向量机回归;优选为偏最小二乘回归。
[0029] 在偏最小二乘回归算法下，从训练集数据中求取若干组特征向量并获得与所述特 征向量一一对应的特征值，以特征值最高的两组特征向量作为第1特征向量和第2特征向 量，即依据训练集中每组样品的近红外吸收光谱数据，求出分别与所述第1特征向量 和第2特征向量对应的权重值，作为各训练集样品的得分值1和得分值2，即Ti、T 2。其中，h、 T2U2和Xsnv之间的关系为：
[0030] Xsnv = TiPi+T2P2+E (IV)
[0031] 公式（IV)中，E是光谱残差。
[0032] 针对待测样品，在相同实验条件下采集待测样品的近红外光谱数据，经过相同的 数据预处理，即SNV处理后所得待测样品预处理后的近红外光谱数据X P_SNV，求取XP_SNV在Pi、 P 2的权重值，即待测样品的得分值1(TPJ)和待测样品的得分值2(Tp_2):
[0033] Tp_i=Xp_snvPi
[0034] Tp-2=Xp-snvP2
[0035] 以上均为线性代数领域的常规数据处理方法，本发明不做特殊限定。
[0036] 针对训练集样品，分别以!^、!^为横、纵坐标，绘制二维得分散点图。针对每个待测 样品，分别以TPJ和Τ Ρ_2为横、纵坐标，在二维得分散点图上绘制对应的待测点，根据所述待 测点所在的区域，鉴别待测样品是否为抛光粮食。
[0037] 作为本发明的一种优选方案，所述步骤(3)中，对抛光样品赋值0.1~5,优选为1， 对未抛光样品赋值-5~-0.1，优选为-1;步骤(4)所述分界线为穿过坐标轴原点且与横坐标 正方向呈135°~170°、优选为150°的直线，所述抛光区域位于所述分界线的右上方。
[0038] 本发明提出的一种基于振动光谱的矿物油抛光粮食快速无损检测方法，扩展了振 动光谱在农产品质量安全快速无损检测鉴别方面的用途。该方法具有快速、准确、工作效率 高、绿色无污染等特点，不仅可以为矿物油抛光农产品的快速无损检测鉴别提供有力保障 和技术支持，而且对保障农产品质量安全、提高工农业生产效率、促进农产品市场健康发 展、快速无损检测仪器的研发等方面具有积极作用。
【附图说明】
[0039] 图1为本发明所述一种基于振动光谱的矿物油抛光粮食快速无损检测方法操作流 程图；
[0040] 图2为本发明所述矿物油抛光大米和非抛光大米的近红外吸收光谱图；
[0041] 图3为本发明所述一种基于振动光谱的矿物油抛光粮食快速无损检测方法实施例 1的散点图；
[0042] 图4为本发明所述一种基于振动光谱的矿物油抛光粮食快速无损检测方法对比例 1的散点图。
【具体实施方式】
[0043] 为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例 中的附图，对本发明实施例中的技术方案进行清楚、完整地描述。显然，所描述的实施例是 本发明的一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人 员在没有做出创造性劳动的前提下所获得的所有其他实施例，都属于本发明保护的范围。 [0044]以下各实施例中，所用仪器为:傅里叶变换中红外光谱仪，型号：Spectrum 400,美 国Perkin Elmer公司；附件:积分球漫反射附件；
[0045] 参数:光谱分辨率:8cnf1，光谱范围：12800cnf1~4000cnf1，累加次数:64次，以积分 球空光路作参比。
[0046]样品:未抛光大米，矿物油抛光大米。
[0047]数据处理:采用数据标准正态化预处理方法对每条光谱进行数据预处理。所谓数 据标准正态化，是针对每组样品的近红外吸收光谱数据，即每一条光谱数据自身所做的标