一种近红外光谱小波融合的食用油种类快速鉴别方法与流程

本发明涉及食用油快速鉴别技术领域，尤其涉及一种近红外光谱小波融合的食用油种类快速鉴别方法。

背景技术：

食用油含有许多人体所需的重要营养成分，在我们日常饮食生活中也是必不可少的，研究表明，部分食用油中富含多种特征成分，包括棕榈酸、油酸、亚油酸、碘值，不同种类的食用油中有着不同的多种特征成分含量组成；另外，由于多种特征成分的含量主要决定了食用油的营养价值，故通常作为评价食用油品质的重要成分指标，同时也是确定其商业价值的重要依据。随着食用油价格不断的提高，许多不法商家为了牟取暴利，将低价值油冒充高价值油投入市场，或将低价值油掺入高价值油中，严重侵害着消费者和合法生产销售企业的利益。因此，有必要研究一种快速鉴别食用油种类的方法，对维护消费者及合法经营者利益、维持食用油市场正常秩序具有重要意义。

技术实现要素：

针对上述技术中存在的不足之处，本发明提供了一种安全可靠、便捷高效的近红外光谱小波融合的食用油种类快速鉴别方法。

本发明解决其技术问题所采用的技术方案是：一种近红外光谱小波融合的食用油种类快速鉴别方法，包括：步骤一、样品选取：选取待鉴别的未知种类的食用油样品；步骤二、一次光谱采集：采集未知种类的食用油样品在1350～1450nm范围的近红外光谱图，得到一次近红外光谱图；步骤三、二次光谱采集：采集未知种类的食用油样品在1700～1800nm范围的近红外光谱图，得到二次近红外光谱图；步骤四、光谱预处理：分别对采集得到的一次近红外光谱图和二次近红外光谱图进行预处理，得到预处理一次近红外光谱图和预处理二次近红外光谱图；步骤五、光谱小波融合：将预处理一次近红外光谱图和预处理二次近红外光谱图进行小波融合，得到小波融合光谱图；步骤六、种类鉴别：根据步骤五中得到的未知种类的食用油样品的小波融合光谱图，采用优化定性模型对所述未知种类的食用油样品进行种类鉴别。

优选的，所述步骤一中选取的未知种类的食用油样品为大豆油、花生油、菜籽油、茶籽油、稻米油、玉米油、葵花籽油、橄榄油中的任意一种。

优选的，所述步骤二中一次光谱采集的条件如下：将食用油样品置于流通池中，设置食用油样品的温度为30±5℃，湿度为40±5％rh，食用油样品的流速为1～1.2cm/s；设置近红外光谱仪的分辨率为4cm^-1，扫描次数为32次，光纤探头的光程为5mm，且光纤探头的入射方向与食用油样品的流向构成一38.38°的入射角；每个食用油样品取三次测量的平均值作为最终采集的一次近红外光谱图。

优选的，所述步骤三中二次光谱采集的条件如下：将食用油样品置于流通池中，设置食用油样品的温度为30±5℃，湿度为40±5％rh，食用油样品的流速为0.4～0.6cm/s；设置近红外光谱仪的分辨率为8cm^-1，扫描次数为64次，光纤探头的光程为10mm，且光纤探头的入射方向与食用油样品的流向构成一38.38°的入射角；每个食用油样品取三次测量的平均值作为最终采集的二次近红外光谱图。

优选的，所述步骤四中对采集得到的一次近红外光谱图分别以1390cm^-1处和1420cm^-1的特征峰强度为基准进行归一化处理，得到预处理一次近红外光谱图。

优选的，所述步骤四中对采集得到的二次近红外光谱图分别以1725cm^-1处和1760cm^-1的特征峰强度为基准进行归一化处理，得到预处理二次近红外光谱图。

优选的，所述步骤五中光谱小波融合的方法如下：将预处理一次近红外光谱图和预处理二次近红外光谱图进行4层小波离散变换，得到13个高低频带，其中包含12个高频带和1个低频带，再根据融合规则分别对12个高频带和1个低频带进行融合，并进行小波逆变换，得到小波融合光谱图。

优选的，所述步骤六中优化定性模型的建立方法如下：收集若干不同已知种类的食用油样品，所述不同已知种类的食用油样品包括：大豆油、花生油、菜籽油、茶籽油、稻米油、玉米油、葵花籽油、橄榄油；采集若干不同已知种类的食用油样品的一次近红外光谱图和二次近红外光谱图，一次近红外光谱图采集条件如步骤二中所述，二次近红外光谱图采集条件如步骤三中所述；分别对采集的若干不同已知种类的食用油样品的一次近红外光谱图和二次近红外光谱图进行预处理，得到若干不同已知种类的食用油样品的预处理一次近红外光谱图和预处理二次近红外光谱图，光谱预处理方法如步骤四中所述；将若干不同已知种类的食用油样品的预处理一次近红外光谱图和预处理二次近红外光谱图进行小波融合，得到若干不同已知种类的食用油样品的小波融合光谱图，光谱小波融合的方法如步骤五中所述；将若干不同已知种类的食用油样品的小波融合光谱图作为定性模型的输入变量，通过多核学习支持向量机分类方法建立若干不同已知种类的食用油样品的定性模型，并采用粒子群优化算法对定性模型中的参数进行优化，得到优化定性模型。

优选的，将若干不同已知种类的食用油样品的小波融合光谱图作为定性模型的输入变量，通过多核学习支持向量机分类方法将小波融合光谱图的数据按特征维数分为10组，这10组均利用高斯核建立单核支持向量机定性模型；并采用粒子群优化算法对定性模型中的惩罚因子矩阵[c]和核函数参数矩阵[g]进行优化，[c]和[g]均是10x8的矩阵，优化时设置种群粒子数为20，每个粒子维数为2，迭代进化次数为100，学习因子初始值设定为c1＝2，c2＝1，得到10组参数([c]，[g])，并对这10组参数([c]，[g])进行加权投票，从而得到优化定性模型。

优选的，经优化，

当[c]等于

[g]等于

且权重为[0.141，0.0976，0.0151，0.043，0.0844，0.1478，0.1490，0.0243，0.1499，0.1478]时，得到优化定性模型。

本发明与现有技术相比，其有益效果是：本发明提供的近红外光谱小波融合的食用油种类快速鉴别方法，通过采用不同波段范围的一次光谱采集和二次光谱采集，并设置两次光谱采集时食用油样品的相同温度、湿度和不同流速，及设置两次光谱采集时近红外光谱仪的不同分辨率、扫描次数、光纤探头光程和相同的入射角，可有效消除外部条件对近红外光谱图的影响；通过小波变换融合一次近红外光谱图和二次近红外光谱图，可有效压缩光谱图中的无用变量并突出有效信息；通过采用多核学习支持向量机分类方法结合粒子群优化算法建立优化定性模型，可显著提升鉴别准确率。

附图说明

图1是近红外光谱小波融合的食用油种类快速鉴别方法的流程示意图。

具体实施方式

下面结合附图对本发明做进一步的详细说明，以令本领域技术人员参照说明书文字能够据以实施。

如图1所示，本发明提供了一种近红外光谱小波融合的食用油种类快速鉴别方法，包括：

步骤一、选取待鉴别的未知种类的食用油样品，所述选取的未知种类的食用油样品为大豆油、花生油、菜籽油、茶籽油、稻米油、玉米油、葵花籽油、橄榄油中的任意一种；

步骤二、将食用油样品置于流通池中，设置食用油样品的温度为30±5℃，湿度为40±5％rh，食用油样品的流速为1～1.2cm/s；设置近红外光谱仪的分辨率为4cm^-1，扫描次数为32次，光纤探头的光程为5mm，且光纤探头的入射方向与食用油样品的流向构成一38.38°的入射角；采集未知种类的食用油样品在1350～1450nm范围的近红外光谱图，每个食用油样品取三次测量的平均值作为最终采集的一次近红外光谱图；

步骤三、将食用油样品置于流通池中，设置食用油样品的温度为30±5℃，湿度为40±5％rh，食用油样品的流速为0.4～0.6cm/s；设置近红外光谱仪的分辨率为8cm^-1，扫描次数为64次，光纤探头的光程为10mm，且光纤探头的入射方向与食用油样品的流向构成一38.38°的入射角；采集未知种类的食用油样品在1700～1800nm范围的近红外光谱图，每个食用油样品取三次测量的平均值作为最终采集的二次近红外光谱图；

步骤四、对采集得到的一次近红外光谱图分别以1390cm^-1处和1420cm^-1的特征峰强度为基准进行归一化处理，得到预处理一次近红外光谱图；对采集得到的二次近红外光谱图分别以1725cm^-1处和1760cm^-1的特征峰强度为基准进行归一化处理，得到预处理二次近红外光谱图；

步骤五、将预处理一次近红外光谱图和预处理二次近红外光谱图进行4层小波离散变换，得到13个高低频带，其中包含12个高频带和1个低频带，再根据融合规则分别对12个高频带和1个低频带进行融合，并进行小波逆变换，得到小波融合光谱图；其中，融合规则如下：

对预处理一次近红外光谱图i和预处理二次近红外光谱图v进行拉普拉斯变换，得到高通滤波图像hi和hv，hi＝i*l，hv＝v*l，其中l是一个11x11的拉普拉斯滤波器；

对经拉普拉斯滤波处理的高通滤波图像hi和hv的系数的绝对值采用高斯滤波处理，构造图像hi和hv的显著性图xi和xv，其中其中g是大小为(2μg+1)(2δg+1)的高斯低通滤波器；参数μg和δg均为5；

当时，取qi＝1否则qi＝0，qv处理方式相同，其中和分别为预处理一次近红外光谱图i和预处理二次近红外光谱图v在像素位置m处的显著性特征值；

对权重图qi和qv分别进行双边滤波获得预处理一次近红外光谱图i和预处理二次近红外光谱图v的低频子带系数ii和iv，融合权重图和其中σu＝10，σv＝10分别是双边滤波中控制空间距离和灰度距离的参数；

步骤六、根据步骤五中得到的未知种类的食用油样品的小波融合光谱图，采用优化定性模型对所述未知种类的食用油样品进行种类鉴别；

其中，优化定性模型的建立方法如下：

收集若干不同已知种类的食用油样品，所述不同已知种类的食用油样品包括：大豆油、花生油、菜籽油、茶籽油、稻米油、玉米油、葵花籽油、橄榄油；

采集若干不同已知种类的食用油样品的一次近红外光谱图和二次近红外光谱图，一次近红外光谱图采集条件如步骤二中所述，二次近红外光谱图采集条件如步骤三中所述；

分别对采集的若干不同已知种类的食用油样品的一次近红外光谱图和二次近红外光谱图进行预处理，得到若干不同已知种类的食用油样品的预处理一次近红外光谱图和预处理二次近红外光谱图，光谱预处理方法如步骤四中所述；

将若干不同已知种类的食用油样品的预处理一次近红外光谱图和预处理二次近红外光谱图进行小波融合，得到若干不同已知种类的食用油样品的小波融合光谱图，光谱小波融合的方法如步骤五中所述；

将若干不同已知种类的食用油样品的小波融合光谱图作为定性模型的输入变量，通过多核学习支持向量机分类方法将小波融合光谱图的数据按特征维数分为10组，这10组均利用高斯核建立单核支持向量机定性模型；并采用粒子群优化算法对定性模型中的惩罚因子矩阵[c]和核函数参数矩阵[g]进行优化，[c]和[g]均是10x8的矩阵，优化时设置种群粒子数为20，每个粒子维数为2，迭代进化次数为100，学习因子初始值设定为c1＝2，c2＝1，得到10组参数([c]，[g])，并对这10组参数([c]，[g])进行加权投票，经优化，

当[c]等于

[g]等于

且权重为[0.141，0.0976，0.0151，0.043，0.0844，0.1478，0.1490，0.0243，0.1499，0.1478]时，得到优化定性模型。

实施例

1、样品选取

收集大豆油、花生油、菜籽油、茶籽油、稻米油、玉米油、葵花籽油、橄榄油8类食用油样品共468份，采用spxy算法按3∶1的比例选取训练集食用油样品351份和预测集食用油样品117份；其中，8类食用油样品的预测集样品和训练集样品的数量分布如下表1。

表1

2、一次光谱采集

将食用油样品置于流通池中，设置食用油样品的温度为30℃，湿度为40％rh，食用油样品的流速为1cm/s；设置近红外光谱仪的分辨率为4cm^-1，扫描次数为32次，光纤探头的光程为5mm，且光纤探头的入射方向与食用油样品的流向构成一38.38°的入射角；采集468份食用油样品在1350～1450nm范围的近红外光谱图，每个食用油样品取三次测量的平均值作为最终采集的一次近红外光谱图。

3、二次光谱采集

将食用油样品置于流通池中，设置食用油样品的温度为30℃，湿度为40％rh，食用油样品的流速为0.5cm/s；设置近红外光谱仪的分辨率为8cm^-1，扫描次数为64次，光纤探头的光程为10mm，且光纤探头的入射方向与食用油样品的流向构成一38.38°的入射角；采集468份食用油样品在1700～1800nm范围的近红外光谱图，每个食用油样品取三次测量的平均值作为最终采集的二次近红外光谱图。

4、光谱预处理

对采集得到的468份食用油样品的一次近红外光谱图分别以1390cm^-1处和1420cm^-1的特征峰强度为基准进行归一化处理，得到预处理一次近红外光谱图；

对采集得到的468份食用油样品的二次近红外光谱图分别以1725cm^-1处和1760cm^-1的特征峰强度为基准进行归一化处理，得到预处理二次近红外光谱图。

5、光谱小波融合

将训练集中的351份食用油样品的预处理一次近红外光谱图和预处理二次近红外光谱图进行4层小波离散变换，得到13个高低频带，其中包含12个高频带和1个低频带，再根据融合规则分别对12个高频带和1个低频带进行融合，并进行小波逆变换，得到小波融合光谱图。

6、优化定性模型建立

将训练集中的351份食用油样品的小波融合光谱图作为定性模型的输入变量，通过多核学习支持向量机分类方法将小波融合光谱图的数据按特征维数分为10组，这10组均利用高斯核建立单核支持向量机定性模型；并采用粒子群优化算法对定性模型中的惩罚因子矩阵[c]和核函数参数矩阵[g]进行优化，[c]和[g]均是10x8的矩阵，优化时设置种群粒子数为20，每个粒子维数为2，迭代进化次数为100，学习因子初始值设定为c1＝2，c2＝1，得到10组参数([c]，[g])，并对这10组参数([c]，[g])进行加权投票，经优化，

当[c]等于

[g]等于

且权重为[0.141，0.0976，0.0151，0.043，0.0844，0.1478，0.1490，0.0243，0.1499，0.1478]时，得到优化定性模型，其对训练集中的351份食用油样品的鉴别准确率为98.58％。

7、种类鉴别

将预测集中的157份食用油样品的小波融合光谱图(融合方法如步骤五中所述)作为定性模型的输入变量，采用上述经优化得到的惩罚因子矩阵[c]和核函数参数矩阵[g]，且权重为[0.141，0.0976，0.0151，0.043，0.0844，0.1478，0.1490，0.0243，0.1499，0.1478]时的优化定性模型对预测集中的157份食用油样品的种类进行预测鉴别，预测鉴别准确率为99.36％。

尽管本发明的实施方案已公开如上，但其并不仅限于说明书和实施方式中所列运用，它完全可以被适用于各种适合本发明的领域，对于熟悉本领域的人员而言，可容易地实现另外的修改，因此在不背离权利要求及等同范围所限定的一般概念下，本发明并不限于特定的细节和这里示出与描述的图例。