应用近红外光谱投影模型法分析烟叶的部位特征

近红外分析技术（ NIR） 是近些年发展起来的绿色分析技术 ，具备简单 、高效 、快速等优点 ，适于过程分析和在线分析，是质量控制的理想手段 ，在农业 、食品 、石油等领域特别是烟草行业的应用日益广泛。传统的化学分析等手段需要花费大量的时间和精力 ，难以对整个加工环节进行控制 ，近红外光谱主要是对含氢基团 X—H（X-C，N ，O）振动的倍频和合频吸收 ，其中包含了大多数类型有机化合物的组成信息 ，与烟叶化学成分关联的信息非常丰富；因此利用近红外对烟叶的主要物理指标与化学指标进行检测和控制对提升烟草品质极为有效 ，已经在烟叶工业分级 、复烤在线 、制丝在线 、烟用辅料等多个方面开展了应用 ，对我国中式卷烟品质的稳定与提升起到了重大作用。

 烟叶的质量管理中 ，部位是一个十分重要的要素特征 ，不同的烟叶部位其物理特性 、化学成分以及评吸特性等都有着各自的特点 ，比如从评吸特性上 ，越往顶部香气量越饱满 ，香气浓度越往顶部越高 ，劲头和刺激性也相应增加等。根据烟叶在植株着生位置不同可分为顶叶 、上二棚叶 、腰叶 、下二棚叶和脚叶等 ５ 个部位 ；在实际烟叶收购与质量管理中 ，通常将烟叶分为上部 、中部 、下部三种类型 。应用近红外光谱进行烟叶部位的判别与相似性分析已有相关文献报道，本工作提出应用近红外光谱投影模型及量化最近次近类分值法分析烟叶的部位特征 ，不仅可以实现预测样品的判别分析 ，而且可以给出部位特征程度的量化分值 ，可将预测结果细化为典型上 、上偏中 、中偏上 、典型中 、中偏下 、下偏中 、典型下等 ７ 类或超模型范围样品 ，得到关于预测样品更丰富的特征信息 ，对烟草工业企业实现原料的均值化加工和烟叶配方等提供更细致的分析数据 。

１ 实验部分
1.1 材料
 建模试验样品为2003年 —2012年收集的烤烟样品共5170个，其中上部烟叶1394个，中部烟叶2550个，下部烟叶1226个，分别来自巴西、津巴布韦和湖南、湖北、云南、贵州、四川、重庆、广东、广西、河南、安徽、福建、陕西、黑龙江、山东、吉林、辽宁等 １６ 个地区 。 预测样品为２０１２ 年从各复烤厂收集的烤烟样品共34个

 由表 １ 中可以看出 ，４ ～ ５ 个主成分数之间得到的各类投影均值之间的距离已达到稳定 ，在 PC 为 ４ 和 ５ 情况下得到的二维投影结果情况分别见图 １（a）和（b）（注 ：图 １ 中圆圈表示各类的类内离散度 ，其半径值为类内各样品投影值标准差的 ２ 倍） 。

 遵循相邻主成分数下得到的投影结果没有显著性差异和主成分个数尽量少的原则，推荐主成分数为 ４ 下建立投影分析模型 ，由图 １ 投影模型的分析结果表明 ，在 ２ 倍标准差水平下 ，上 、下部烟叶的近红外光谱特征具有显著差别 ，基本实现完全区分 ；而中部烟叶分别与上部和下部具有一定程度重叠 ，这与部位本身具有连续性特征的实际情况相符合 。
2.2 部位特征模型的预测验证分析
　　应用 １.３ 中的分析方法及 ２. １ 中建立的部位特征模型 ，对于 ２０１２ 在实际复烤生产加工中取得的不同部位 、不同产地的 ３４ 个烟叶样品进行了预测分析 ，各样品的产地 、等级信息 ，以及得到的最近预测类型 、分值 ，次近预测类型 、分值 ，细化的 ７ 类特征描述结果和实际评吸验证评价等情况见表２ 。

　由表 ２ 的分析结果可以看出 ，分析结果中没有出现超模型范围样品 ，也没有出现“上偏下”或“下偏上”等奇异类型品 ，其中 ３０ 个样品的最近预测类型与等级名称描述的部位属性一致 ， 其分析结果整体上合理可靠 。依据最近类预测分值（V １ ）和次近类预测分值（V ２ ）得到的 ７ 类细化特征描述 ，与专家评吸进行的比对验证结果均为准确或合理 。
３ 结论
（１）提出的应用投影模型及量化最近次近类分值法分析烟叶近红外光谱的部位特征 ，不仅可以实现预测样品的判别分析以及部位特征程度的量化分值 ，并且可将预测结果细化为典型上 、上偏中 、中偏上 、典型中 、中偏下 、下偏中 、典型下等 ７ 类或超模型范围样品 ，得到关于预测样品更丰富的属性特征信息 ，可对烟草工业企业实现原料的均值化加工和烟叶配方等提供指导 。
（２）采用的建模数据 ，来源于上海烟草集团公司技术中心 １０ 余年来承担国家烟草局等重大研究课题过程中收集的国内外重要烤烟产地的样品 ，数量大 ，代表性好 ，部位特征信息准确 ，所建模型具有较好的实用性 ，其分析结果对监督和改善目前烟叶收购中部位等级等信息不够准确 ，提升中式卷烟的原料质量水平等可发挥重要作用 。
（３）采用的分析方法 ，对实现其他食品或农产品的产区 、风格 、风味 、口感等类似品质属性特征的分析具有参考价值 。
（４）实现近红外定量或定性分析 ，均是基于大量样本下建立统计模型的分析 ，存在犯错风险 ，尤其是定性所采用的细化属性特征分析 ，其犯错风险更大 ，因此对其分析结果应采用“专家解读” 、“人机交互”的方式进行应用 。同时 ，所采用的最近次近距离法 ，没有考虑特征空间分布的方向性差异 ，有待进一步研究开发更合理的分析方法。