本发明涉及及食品质量与安全检测领域，具体涉及在线快速检测鸡肉灰分的方法。
背景技术：
鸡肉作为主要的禽肉产品，因富含蛋白质、氨基酸、脂肪酸、维生素、矿物质且营养丰富比例均衡等特点，其消费量逐年增加，同时鸡肉的食用安全性也越来越受到消费者的重视，其中作为鸡肉重要组成部分的灰分，其即能反映鸡肉的营养价值，也能判断鸡在屠宰过程中受污染的程度。鸡肉经过高温灼烧后的残留物称为总灰分，其总灰分中即包含11种含量大于0.01％的人体常量元素，也包含含量低于0.01％的人体必需的微量元素。正常肉的灰分含量都低于2％，当灰分含量超过2％时说明肉在屠宰过程中受到了污染。一般鸡肉中灰分的检测常使用gb5009.4—2016《食品安全国家标准食品中灰分的测定》方法进行检测，但此方法需用化学试剂对样品进行预处理且步骤繁琐、检测时间长、成本高，满足不了当今肉品行业大规模的在线检测要求。
近年来，高光谱成像系统因能同时检测被测样品的内部品质和外部品质，且还具有高分辨率、快速、无损、对样品无需进行预处理等优点，此技术已成为农畜产品无损检测领域的研究热点，也取得了诸多成果，然而在鸡肉灰分方面的检测鲜有报道。
技术实现要素：
为了解决现有技术的不足，本发明提供了在线快速检测鸡肉灰分的方法。
本发明的技术方案是：在线快速检测鸡肉灰分的方法，采集校正集鸡肉样品的高光谱图像，对获取的光谱图进行预处理再进行目标区域的识别以及光谱图平均光谱数据的提取；将提取的光谱数据代入下式即得y灰分＝0.0331+0.125x900.55nm-0.079x903.845nm+0.0307x908.787nm-0.0566x920.316nm+0.0868x931.844nm+0.0425x943.371nm-0.0226x1042.123nm-0.0218x1162.211nm-0.0247x1213.208nm+0.0245x1257.634nm+0.0124x1506.465nm+0.0798x1695.097nm-0.0767x1700.074nm，其中y灰分为鸡胸肉中灰分的含量，x900.55nm、x903.845nm、x908.787nm、x920.316nm、x931.844nm、x943.371nm、x1042.123nm、x1162.211nm、x1213.208nm、x1257.634nm、x1506.465nm、x1695.097nm、x1700.074nm分别为波长在900.55nm、903.845nm、908.787nm、920.316nm、931.844nm、943.371nm、1042.123nm、1162.211nm、1213.208nm、1257.634nm、1506.465nm、1695.097nm、1700.074nm处的光谱反射率值；上式相关系数为r＝0.943，均方根误差rmse＝0.001。
本发明的进一步改进包括：
对获取的光谱图进行预处理即黑白板校正按照以下公式进行：
其中rc为校正后的图像，rr为原始光谱图像；rd为黑板图像，其反射率为0％，rw为白板图像，其反射率为99.9％。
本发明提供了一种高分辨率、快速、无损、无需对样品进行预处理等优点的高光谱成像技术来检测鸡肉中的灰分，以弥补现有技术在的不足，从而实现鸡肉灰分大规模的在线检测。
本发明从486个全波段内提取13个最优波长时为了剔除大量的冗余信息，提取有用的信息，以降低数据分析的计算量，从而提高偏最小二乘模型的精度，以实现肉品企业大规模的在线生产的需求。与现有技术相比，本发明具有以下有益效果：本发明不需对被测样品进行预处理，仅对样品进行非接触的光谱扫描且没有破坏性；本发明不使用任何化学试剂，即绿色又节约成本；本发明易于操作又节约时间，能实现鸡肉灰分的大规模在线检测。
附图说明
图1是84个校正集样品的平均光谱特征图。
图2是鸡胸肉最优波长的提取图。
图3是鸡胸肉灰分预测值与实测值之间的相关性。
具体实施方式
下面结合实施例对本发明做详细说明。
实施例
本实施例的一种基于近红外高光谱成像技术在线检测鸡肉灰分的方法步骤如下：
1.建立鸡肉灰分含量的全波段的校正模型；
(1.1)将购买的鸡胸肉切割成3*3*1(长*宽*高)的小样品，作为校正集样品，共获取84个样品，将其平均分成7份分别放在一次性带盖的塑料盒里，再放置在4℃的冰箱内进行贮藏0、1、2、3、4、5、6天；
(1.2)在试验之前，提前30min打开高光谱成像系统预热，同时鸡肉样本也提前30min从冰箱内取出待其恢复至室温，首先采集黑白板的图像，再采集校正集样品的高光谱图像，此系统扫描时的速度为6.54mm/s，曝光时间为4.65ms，其成像模式为反射且波长范围为900-1700nm；
(1.3)对采集过光谱图像的校正集样品立即使用gb5009.4—2016《食品安全国家标准食品中灰分的测定》方法检测其灰分含量，其结果如表1；
表184个校正集样品的灰分含量
(1.4)对获取的光谱图首先按照以下公式进行校正：
其中rc为校正后的图像，rr为原始光谱图像；rd为黑板图像，其反射率为0％，rw为白板图像，其反射率为99.9％。
再对校正过的光谱图进行感兴趣区域的识别以及光谱图平均光谱数据的提取，其84个样品的平均光谱图特征如图1。
(1.5)采用偏最小二乘法(plsr)建立步骤(1.3)灰分含量与步骤(1.4)中的光谱数据之间的全波段(486个波长)定量模型，即全波段的偏最小二乘模型，当所建模型的相关系数r越接近于1，均方根误差rmse越小且交叉验证集的相关系数和均方根误差越接近于建模集时说明建模集模型的精度和稳定性就越好，结果如下表2；
表2校正集样品灰分含量的全波段plsr预测模型
从表2中可以得出所建立的建模集plsr模型的相关系数r高达0.943，均方根误差低至0.001，其中交叉验证集的模型相关系和均方根误差数也接近于建模集，表明建模集的模型精度高且较稳定。
2.全波段校正模型的优化
(2.1)上述步骤(1.5)所建的在900-1700nm全波段内偏最小二乘模型中，共有486个波长，而并不是所有的波长都对所建模型具有贡献，其中存在大量的冗余信息，为了剔除冗余信息保留有用信息，通过回归系数法(rc)来提取最优波长，以降低数据的计算量，从而提高计算机的运行速度。结果如图2：
(2.2)使用回归系数法从步骤(1.5)所建的全波段偏最小二乘模型中提取了13个最优波长，分别为900.55、903.845、908.787、920.316、931.844、943.371、1042.123、1162.211、1213.208、1257.634、1506.465、1695.097、1700.074，以提取的13个最优波长作为输入变量来建立优化后的偏最小二乘模型即鸡胸肉灰分含量的校正模型，结果如表3：
表3校正集样品灰分含量的最优波长plsr预测模型
从表3中可得出使用13个最优波长所建立的建模集和交叉验证集的plsr模型相关分别为0.943和0.923，均方根误差均为0.001即两者的差距极小，且与全波段建模集的相关系数和均方根误差的差距也不大，故使用最优波长所建立的plsr模型精度即高又稳定。
(2.3)步骤(2.2)所获得的最优波长的偏最小二乘模型即校正模型公式为：y灰分＝0.0331+0.125x900.55nm-0.079x903.845nm+0.0307x908.787nm-0.0566x920.316nm+0.0868x931.844nm+0.0425x943.371nm-0.0226x1042.123nm-0.0218x1162.211nm-0.0247x1213.208nm+0.0245x1257.634nm+0.0124x1506.465nm+0.0798x1695.097nm-0.0767x1700.074nm，其中y灰分为鸡胸肉中灰分的含量，其中x900.55nm、x903.845nm、x908.787nm、x920.316nm、x931.844nm、x943.371nm、x1042.123nm、x1162.211nm、x1213.208nm、x1257.634nm、x1506.465nm、x1695.097nm、x1700.074nm分别为波长在900.55nm、903.845nm、908.787nm、920.316nm、931.844nm、943.371nm、1042.123nm、1162.211nm、1213.208nm、1257.634nm、1506.465nm、1695.097nm、1700.074nm处的光谱反射率值。利用此模型来检测待测鸡胸肉灰分的含量。
3.测试
(3.1)采集28个待测鸡肉样品的高光谱图像，对获取的光谱图像进行预处理及感兴趣区域的识别和光谱数据的提取；
(3.2)将获取的每个待测样品的900.55nm、903.845nm、908.787nm、920.316nm、931.844nm、943.371nm、1042.123nm、1162.211nm、1213.208nm、1257.634nm、1506.465nm、1695.097nm、1700.074nm处下的光谱反射率值输入到步骤(2.3)的校正模型内，便可获得28和待测鸡胸肉灰分的含量。将得到的鸡胸肉的灰分含量与使用gb5009.4—2016《食品安全国家标准食品中灰分的测定》方法所测得的灰分含量值进行关联，其相关系数高达0.921，均方根误差为0.0012，其真实值与预测值之间的相关很好。结果如图3。
以上显示和描述了本发明的基本原理和主要特征和本发明的优点。本行业的技术人员应该了解，本发明不受上述实施例的限制，上述实施例和说明书中描述的只是说明本发明的原理，在不脱离本发明精神和范围的前提下，本发明还会有各种变化和改进，这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。