CN105203464A - 基于高光谱成像技术检测花生中油酸含量分布的方法 - Google Patents

Abstract

本发明提供基于高光谱成像技术检测花生中油酸含量分布的方法，包括：采集花生样品在特征波长处的光谱图像，将特征波长处经过预处理后的光谱反射值，输入花生油酸含量分布定量模型，得到花生样品油酸含量分布。本发明还提供建立花生中油酸含量分布定量模型的方法，包括采集花生高光谱图像，并利用常规方法测定其油酸含量；高光谱图像经过图像校正与背景删除，提取平均光谱；以预处理后平均光谱为自变量，以油酸含量为因变量，建立全波段油酸含量的回归模型，在此基础上利用回归系数，确定特征波长，建立并验证所述定量模型。本发明快速简便，效率高，不破坏样品，不使用任何化学试剂，测定结果准确，实现了花生油酸含量的可视化。

Description

基于高光谱成像技术检测花生中油酸含量分布的方法

技术领域

本发明涉及一种检测花生中油酸含量的方法，具体地说，涉及基于高光谱成像技术检测花生中油酸含量分布的方法。

背景技术

2013年我国花生产量1697万吨，位居世界第一。我国生产出的花生绝大多数用于加工花生油，花生油中含有大量的单不饱和脂肪酸和多不饱和脂肪酸，其中油酸含量高达35％～72％。花生中油酸含量的高低直接影响花生油的品质和货架期，间接影响企业的效益。传统测定花生中油酸的方法包括：乙酰氯—甲醇甲酯化法和氨水—乙醇提取法，但这些方法存在分析速度慢，操作步骤繁琐，成本高破坏性强，使用试剂污染环境等缺点。因此，急需寻找一种快速，非破坏性的方法为花生油酸含量的测定提供依据。

高光谱成像技术结合了光谱学和成像技术，是一门新兴的快速，无损检测方法。高光谱图像是由一系列连续的波段图像组成的三维图像数据块，其具有某个特定波长下的图像信息，并且针对平面内某个特定像素又具有不同波长下的光谱信息。其原理是利用花生油酸中CH、OH等基团在近红外光谱区的光谱吸收特性，确定光谱与油酸含量之间的定量关系，从而预测花生中油酸含量和分布。

中国专利CN102621077A公布了高光谱反射图像采集系统及基于该系统的玉米种子纯度无损检测方法；中国专利CN1995987公布了基于高光谱图像技术的农畜产品无损检测方法及装置；中国专利CN103636315A公布了一种基于高光谱的种子发芽率在线检测装置及方法。以上发明采用高光谱图像技术检测产品指标，避免了传统方法的局限性。但研究主要集中在种子纯度方面，经检索，到目前为止，国内外还没有用高光谱成像技术检测花生油酸含量分布的报道。

发明内容

为了解决现有技术中存在的问题，本发明的目的是提供基于高光谱成像技术检测花生中油酸含量分布的方法。

为了实现本发明目的，本发明第一方面是提供一种基于高光谱成像技术建立花生中油酸含量分布定量模型的方法，该方法包括以下步骤：

1.1收集具有代表性的花生样品，用高光谱仪扫描获得花生样品中每个像素点在各波长下的图像信息，得到花生样品的原始高光谱三维图像；

优选地，所述高光谱仪扫描的波长范围为900-1700nm，扫描方式为线扫描；

1.2对所述花生样品的原始高光谱三维图像进行校正和背景删除后，提取花生样品图像平均光谱；

优选地，所述校正是指对所述花生样品的原始高光谱三维图像I_raw进行黑白校正，具体方法为对反射率为99％的标准校正板进行采集，得到全白的标定图像I_white，然后关闭镜头采集，得到全黑标定图像I_dark，根据下述公式计算校正后图像I_norm：

$I_{norm}=\left(\frac{I_{raw}-I_{dark}}{I_{white}-I_{dark}}\right)\×100\%$

优选地，所述背景删除具体步骤为：采用主成分分析，确定背景与原料的边界，删除背景，得到花生样品图像；

1.3对所述花生样品图像平均光谱进行中心化、均值方差化结合标准正态变量变换预处理；

1.4采用常规方法检测所述花生样品的油酸含量，得到花生样品的油酸含量；

优选地，所述检测花生样品的油酸含量方法为根据GB5413.27-2010进行，进一步优选为根据GB5413.27-2010中第一法乙酰氯—甲醇甲酯化法进行；

1.5将所述花生样品为校正集，以所述校正集所述预处理后的平均光谱为自变量，以所述花生样品的油酸含量为因变量，通过偏最小二乘法建立所述校正集自变量和因变量的偏最小二乘法回归模型；利用留一法进行验证；即对所述校正集偏最小二乘法回归模型进行验证；

1.6根据所述校正集偏最小二乘法回归模型的回归系数，选择对所述回归模型贡献率绝对值最大的波长为特征波长；并通过偏最小二乘法建立校正集花生中油酸含量分布定量模型；利用留一法进行验证；即对所述定量模型进行验证；

所述留一法是指每次留取一个花生样品作为验证，余下的n-1个花生样品建立验证模型，用所述验证模型对留取的这一个花生样品来进行验证，重复上述的过程，直到对所有花生样品都进行验证。上述步骤1.6所述定量模型表示所述校正集花生样品的油酸含量与所述特征波长处的光谱反射值的定量关系。

优选地，选取所述特征波长分别为：901nm、980nm、1064nm、1147nm、1230nm、1313nm、1397nm、1480nm、1564nm、1648nm；

优选地，所建立的花生中油酸含量分布定量模型如下：

Y_OFA＝10^-3×(–184.75R_901nm+43.84R_980nm–86.76R_1064nm–30.41R_1147nm

+40.89R_1230nm+18.68R_1313nm+133.66R_1397nm+99.63R_1480nm

–92.58R_1564nm–184.26R_1648nm)+791

其中，Y_OFA为花生样品的油酸含量，R_901nm、R_980nm、R_1064nm、R_1147nm、R_1230nm、R_1313nm、R_1397nm、R_1480nm、R_1564nm、R_1648nm分别为花生样品在特征波长901nm、980nm、1064nm、1147nm、1230nm、1313nm、1397nm、1480nm、1564nm、1648nm处经过预处理后的光谱反射值。

本发明对所述花生样品图像平均光谱进行中心化处理可以增加光谱之间的差异，从而提高模型的稳健性和预测能力；均值方差化处理可以给光谱中所有波长变量以相同的权重；进行标准正态变量变换可以消除花生颗粒大小、表面扫射以及光程变化对光谱的影响。

特征波长选取过多或过少都不宜；若特征波长选取过多，则增加计算复杂度；若特征波长选取过少，则会降低检测结果准确度。

本发明进行验证的目的是确保所述定量模型准确性和稳定性。一般地，经验证后若建立的所述定量模型准确、稳定，则可用于检测花生中油酸含量分布；若建立的所述定量模型准确度和稳定性不佳，则需要重新按照上述步骤建立所述回归模型或所述定量模型。

具体地，通过计算所述校正集的相关系数R_cal和留一法的相关系数R_cv以及校正集的标准偏差SEC和留一法的标准偏差SECV来判断所述回归模型和所述定量模型准确度和稳定性。一般地，当相关系数(R_cal或R_cv)≥0.8，标准偏差(SEC或SECV)≤2时，表明所述回归模型或所述定量模型准确度高、稳定性好。

本发明采用下述公式(1)计算相关系数(R_cal或R_cv)；公式(2)计算标准偏差(SEC或SECV)。

$R=\frac{\underset{i=1}{\overset{n}{\Σ}}(x_i-\stackrel{\‾}{x})(y_i-\stackrel{\‾}{y})}{\sqrt{\underset{i=1}{\overset{n}{\Σ}}{(x_i-\stackrel{\‾}{x})}^2\·}\sqrt{\underset{i=1}{\overset{n}{\Σ}}{(y_i-\stackrel{\‾}{y})}^2}}---\left(1\right)$

式(1)中，x_i为第i个样品高光谱方法预测值，是预测值的平均值；y_i为第i个样品常规方法的测定值，是测定值的平均值；n为两个变量的样本值的个数。如果样本为校正集，则R为R_cal；如果为留一法验证，则R为R_cv。

$SEC=\sqrt{\frac{\underset{i=1}{\overset{n}{\Σ}}{(x_i-y_i)}^2}{n-1}}---\left(2\right)$

式(2)中，x_i为校正集第i样品高光谱方法的预测值，y_i为校正集第i样品常规方法的测定值，n为校正集的样品数。如果x_i为留一法验证过程中第i样品高光谱方法的预测值，则公式(2)表示的是SECV。

本发明第二方面是提供上述定量模型在检测花生中油酸含量分布中的应用。

本发明第三方面是提供基于高光谱成像技术检测花生中油酸含量分布的方法，所述方法包括：

1)采集待测花生样品在下列特征波长处的光谱图像：901nm、980nm、1064nm、1147nm、1230nm、1313nm、1397nm、1480nm、1564nm、1648nm；

2)将所述特征波长处经过预处理后的光谱反射值输入花生中油酸含量分布定量模型，得到待测花生样品油酸含量分布；所述花生中油酸含量分布定量模型如下：

Y_OFA＝10^-3×(–184.75R_901nm+43.84R_980nm–86.76R_1064nm–30.41R_1147nm

+40.89R_1230nm+18.68R_1313nm+133.66R_1397nm+99.63R_1480nm

–92.58R_1564nm–184.26R_1648nm)+791

其中，Y_OFA为花生样品的油酸含量，R_901nm、R_980nm、R_1064nm、R_1147nm、R_1230nm、R_1313nm、R_1397nm、R_1480nm、R_1564nm、R_1648nm分别为花生样品在特征波长901nm、980nm、1064nm、1147nm、1230nm、1313nm、1397nm、1480nm、1564nm、1648nm处经过预处理后的光谱反射值。

所述步骤1)采集待测花生样品特征波长处的光谱图像的方法与上述建立花生中油酸含量分布定量模型的方法中获得光谱图像的方法相同。

如无特殊指明，本发明所述预处理是指中心化、均值方差化结合标准正态变量变换预处理。

具体地，所述步骤1)采集待测花生样品特征波长处的光谱图像的方法包括以下步骤：

1.1用高光谱仪扫描获得待测花生样品中每个像素点在各波长下的图像信息，得到待测花生样品的原始高光谱三维图像；

优选地，所述高光谱仪扫描的波长范围为900-1700nm，扫描方式为线扫描；

1.2对所述待测花生样品的原始高光谱三维图像进行校正和背景删除后，提取待测花生样品图像平均光谱；

优选地，所述校正是指对所述花生样品的原始高光谱三维图像I_raw进行黑白校正；具体方法为对反射率为99％的标准校正板进行采集，得到全白的标定图像I_white，然后关闭镜头采集，得到全黑标定图像I_dark，根据下述公式计算校正后图像I_norm：

$I_{norm}=\left(\frac{I_{raw}-I_{dark}}{I_{white}-I_{dark}}\right)\×100\%$

优选地，所述背景删除具体步骤为：采用主成分分析，确定背景与原料的边界，删除背景，得到花生样品图像；

1.3对所述待测花生样品图像平均光谱进行中心化、均值方差化结合标准正态变量变换预处理。

本发明收集了我国主栽地区主栽品种花生，如：白沙1016、海花1号、丰花1号、鲁花11号、鲁花9号，花育19号等，对收集的花生样品无需进行任何的预处理，同时采集高光谱图像和测定油酸含量，并利用偏最小二乘法建立图像中光谱油酸信息与油酸含量的回归模型，在此基础上利用回归系数选择特征波长，应用偏最小二乘法对花生特征波长与油酸含量进行关联研究，确定两者之间的定量关系，即定量模型，测定未知样品的高光谱图像，将图像上的每个像素点下特征波长带入定量模型，计算油酸含量，得到像素点级分辨率的详细油酸空间分布图，实现油酸信息分布的空间可视化。与现有技术相比，本发明具有以下优点和有益效果:

1、本发明通过特征波段实现了快速测定花生油酸含量分布，提高了采集速度，缩短了高光谱数据分析时间，提高了检测效率，并实现了花生油酸信息分布的空间可视化，为真正达到快速无损在线检测提供了理论基础。

2、因样品花生中除含有油酸外，还含蛋白质、碳水化合物、水分等其他物质。这些其他的物质对光谱的影响较大，严重干扰油酸含量检测的准确度。为克服干扰提高精确度，本发明采用科学方法从花生图谱中选择与油酸含量紧密相关的特征波长，填补了高光谱成像技术检测花生油酸含量的空白。

3、花生样品无需进行任何预处理，无破坏性，不使用任何试剂，保护环境，操作快速简单、避免了人为因素的干扰，测定结果更加高效，客观。

4、本发明在大量研究基础上，明确并完善了建立测定花生油酸含量的分析步骤；通过比较不同预处理、不同的建模方法得到模型的优劣程度，确定光谱的最佳预处理为中心化、均值方差化结合标准正态变量变换；最佳建模方法为偏最小二乘法。

5、通过收集近3年来全国花生主栽地区主栽品种，克服了地区的差异、品种的差异和时间的差异，使本发明方法能够涵盖了全国绝大多数品种，适用范围广。

附图说明

图1为实施例1基于高光谱成像技术检测花生中油酸含量分布方法流程图；

图2为实施例1提取花生高光谱图像的平均光谱(未经预处理)；

图3为实施例2待测单个花生品种的花生油酸含量分布图。

图3中“Color”表示用不同的颜色代表花生不同油酸含量；“Amplitude”表示花生中油酸含量的范围；从“18.8至20.2”数字越大表示花生中油酸含量越高。

具体实施方式

以下实施例用于说明本发明，但不用来限制本发明的范围。

以下实施例中所使用的实验方法如无特殊说明，均为常规方法。

以下实施例中所用的材料、试剂等，如无特殊说明，均可从商业途径得到。

如无特殊说明，以下实施例所用高光谱仪ImspectorN17E；感应器TE-cooledInGaAsphotodiodearray；光源10WHalogensidereflector。在20℃下开启高光谱仪，预热10min，设定采集参数，其中曝光时间为5.8s，采集速度为8mm/s，视野范围：200mm，光谱波长范围900-1700nm，分辨率3nm，扫描方式为线扫描。

以下实施例中的高光谱图像的校正、背景删除和光谱的提取、光谱数据的分析和油酸含量空间可视化均由瑞典Umbio公司出售的图像分析软件Evince2.4中完成。

实施例1

本实施例提供一种基于高光谱成像技术建立花生中油酸含量分布定量模型的方法，该方法包括以下步骤：

1.1收集2012、2013和2014年我国主栽省份主栽花生样品96个品种，从每个品种中挑选30粒完整的花生仁，用高光谱仪同时扫描获得花生样品中每个像素点在各波长下的图像信息，重复3遍，取3次扫描的高光谱图像的平均值，得到花生样品的原始高光谱三维图像。每次扫描前，先采集全白的标定图像I_white和全黑标定图像I_dark。

1.2对上述花生样品的原始高光谱三维图像进行校正和背景删除后，提取花生样品图像平均光谱；

所述校正是指对所述花生样品的原始高光谱三维图像I_raw进行黑白校正；具体根据下述公式计算校正后图像I_norm：

$I_{norm}=\left(\frac{I_{raw}-I_{dark}}{I_{white}-I_{dark}}\right)\×100\%$

所述背景删除是采用主成分分析，在主成分1(PC₁)和主成分2(PC₂)下确定背景与花生样品，然后删除背景，确定花生样品为目标区域，将同一品种30粒花生作为整体，从中提取不同品种花生样品平均光谱。

1.3对上述不同品种花生样品图像平均光谱进行中心化、均值方差化结合标准正态变量变换预处理。

1.4采用国家标准GB5413.27-2010中第一法乙酰氯—甲醇甲酯化法测定花生样品的油酸含量，每个品种重复测定三次，取平均值。本步骤利用美国Varian公司GC-450气相色谱仪进行测定。

1.5以上述96个花生品种为校正集，其油酸含量统计见表1，以所述校正集的所述预处理后的花生样品图像平均光谱(具体是指经过预处理后的光谱反射值)为自变量，以所述花生样品的油酸含量为因变量，通过偏最小二乘法建立所述校正集自变量和因变量的偏最小二乘法回归模型(全波段)。然后采用留一法进行验证，采用公式(1)计算校正集的相关系数R_cal和留一法所建立的验证模型的相关系数R_cv，采用公式(2)计算校正集的标准偏差SEC和留一法所建立的验证模型的标准偏差SECV，见表1。

$R=\frac{\underset{i=1}{\overset{n}{\Σ}}(x_i-\stackrel{\‾}{x})(y_i-\stackrel{\‾}{y})}{\sqrt{\underset{i=1}{\overset{n}{\Σ}}{(x_i-\stackrel{\‾}{x})}^2\·}\sqrt{\underset{i=1}{\overset{n}{\Σ}}{(y_i-\stackrel{\‾}{y})}^2}}---\left(1\right)$

式(1)中，x_i为第i个样品高光谱方法预测值，是预测值的平均值；y_i为第i个样品常规方法的测定值，是测定值的平均值；n为两个变量的样本值的个数。如果样本为校正集，则R为R_cal；如果为留一法验证，则R为R_cv。

$SEC=\sqrt{\frac{\underset{i=1}{\overset{n}{\Σ}}{(x_i-y_i)}^2}{n-1}}---\left(2\right)$

式(2)中，x_i为校正集第i样品高光谱方法的预测值，y_i为校正集第i样品常规方法的测定值，n为校正集的样品数。如果x_i为留一法验证过程中第i样品高光谱方法的预测值，则公式表示的是SECV。

表1花生校正集和留一法验证模型的油酸含量统计和模型参数

1.6利用回归系数法(表示波段对油酸含量影响大小的参数)，回归系数绝对值越大表明该波段对油酸含量影响越大，选择对模型贡献率绝对值最大的十点为特征波长，分别为：901nm、980nm、1064nm、1147nm、1230nm、1313nm、1397nm、1480nm、1564nm、1648nm，建立花生中油酸含量分布定量模型，采用公式(1)方法计算校正集的相关系数R_cal，公式(2)计算校正集的误差SEC，见表2，该花生中油酸含量分布定量模型如下：

Y_OFA＝10^-3×(–184.75R_901nm+43.84R_980nm–86.76R_1064nm–30.41R_1147nm

+40.89R_1230nm+18.68R_1313nm+133.66R_1397nm+99.63R_1480nm

–92.58R_1564nm–184.26R_1648nm)+791

其中，Y_OFA为花生样品的油酸含量，R_901nm、R_980nm、R_1064nm、R_1147nm、R_1230nm、R_1313nm、R_1397nm、R_1480nm、R_1564nm、R_1648nm分别为花生样品在特征波长901nm、980nm、1064nm、1147nm、1230nm、1313nm、1397nm、1480nm、1564nm、1648nm处经过预处理后的光谱反射值。

利用留一法对建立的模型进行验证，采用上述公式(1)、(2)分别计算留一法验证模型的相关系数R_cv和标准偏差SECV，见表2。

表2基于特征波长花生油酸含量校正集和留一法验模型参数

选择特征波长能够代表油酸信息，利用本发明方法建立的花生中油酸含量分布定量模型对花生中油酸含量分布进行检测，其检测结果与国标GB5413.27-2010中第一法乙酰氯—甲醇甲酯化法测定花生样品的检测结果呈高度相关，并且能够简化运算分析时间，提高运算速度。

实施例2

本实施例提供一种基于高光谱成像技术检测花生中油酸含量分布的方法，该方法包括以下步骤：

1)采集待测花生样品在下列特征波长处的光谱图像：901nm、980nm、1064nm、1147nm、1230nm、1313nm、1397nm、1480nm、1564nm、1648nm；

具体过程：另取单个花生品种，按与实施例1相同的方法用高光谱仪得到花生样品的原始高光谱三维图像；进而用与实施例1相同的方法提取花生样品图像平均光谱；然后再对该品种花生样品图像平均光谱进行中心化、均值方差化结合标准正态变量变换预处理；最终获得该花生品种样品在上述特征波长处的光谱图像。

2)将上述特征波长处的光谱反射值输入花生中油酸含量分布定量模型，得到待测花生样品油酸含量分布，结果如图3所示；所述花生中油酸含量分布定量模型如下：

Y_OFA＝10^-3×(–184.75R_901nm+43.84R_980nm–86.76R_1064nm–30.41R_1147nm

+40.89R_1230nm+18.68R_1313nm+133.66R_1397nm+99.63R_1480nm

–92.58R_1564nm–184.26R_1648nm)+791

其中，Y_OFA为花生样品的油酸含量，R_901nm、R_980nm、R_1064nm、R_1147nm、R_1230nm、R_1313nm、R_1397nm、R_1480nm、R_1564nm、R_1648nm分别为花生样品在特征波长901nm、980nm、1064nm、1147nm、1230nm、1313nm、1397nm、1480nm、1564nm、1648nm处经过预处理后的光谱反射值。

对比例1

本对比例提供一种基于高光谱成像技术建立花生中油酸含量分布定量模型的方法，与实施例1的区别仅在于选取的偏最小二乘法回归模型的波长不同。本对比例选取十个波长，分别为907nm、947nm、1020nm、1094nm、1217nm、1274nm、1364nm、1423nm、1584nm、1688nm，并基于这十个波长以与实施例1相同方法建立花生中油酸含量分布定量模型，校正集的相关系数R_cal和校正误差SEC，见表3，建立的花生中油酸含量分布定量模型如下：

Y_OFA＝10^-3×(–82.56R_907nm+32.76R_947nm–45.62R_1020nm–42.51R_1094nm

+62.18R_1217nm+32.15R_1274nm+54.68R_1364nm+75.12R_1423nm

–63.88R_1584nm–35.84R_1688nm)+654

利用留一法对建立的模型进行评估和验证，采用上述公式(1)(2)计算留一法验证模型的相关系数R_cv和标准误差SECV，见表3。

表3基于其他波长花生油酸含量校正集和留一法验证模型参数

对比例2

本对比例提供一种基于高光谱成像技术建立花生中油酸含量分布定量模型的方法，与实施例1的区别仅在于选取的偏最小二乘法回归模型的波长不同。本对比例选取十个波长，分别为937nm、994nm、1080nm、1153nm、1254nm、1347nm、1464nm、1484nm、1584nm、1651nm，并基于这十个波长以与实施例1相同方法建立花生中油酸含量分布定量模型，校正集的相关系数R_cal和校正误差SEC，见表4，建立的花生中油酸含量分布定量模型如下：

Y_OFA＝10^-3×(72.16R_937nm+21.54R_994nm–124.23R_1080nm–52.34R_1153nm

+42.64R_1254nm–49.12R_1347nm+35.24R_1464nm+67.12R_1484nm

–32.14R_1584nm–27.32R_1651nm)+712

利用留一法对建立的模型进行评估和验证，采用上述公式(1)(2)计算留一法验证模型的相关系数R_cv和标准误差SECV，见表4。

表4基于其他波长花生油酸含量校正集和留一法验证模型参数

由实施例1-2及对比例1-2的结果来看，特征波长的选取对测定花生中油酸含量，有着重要影响，本发明选择的特征波长建立的模型，相关系数高，误差低，可以用来测定花生中油酸含量。

虽然，上文中已经用一般性说明及具体实施方案对本发明作了详尽的描述，但在本发明基础上，可以对之作一些修改或改进，这对本领域技术人员而言是显而易见的。因此，在不偏离本发明精神的基础上所做的这些修改或改进，均属于本发明要求保护的范围。

Claims (10)

1.一种基于高光谱成像技术建立花生中油酸含量分布定量模型的方法，该方法包括以下步骤：

1.1收集具有代表性的花生样品，用高光谱仪扫描获得花生样品中每个像素点在各波长下的图像信息，得到花生样品的原始高光谱三维图像；

1.2对所述花生样品的原始高光谱三维图像进行校正和背景删除后，提取花生样品图像平均光谱；

1.3对所述花生样品图像平均光谱进行中心化、均值方差化结合标准正态变量变换预处理；

1.4采用常规方法检测所述花生样品的油酸含量，得到花生样品的油酸含量；

1.5将所述花生样品作为校正集，以所述校正集所述预处理后的平均光谱为自变量，以所述花生样品的油酸含量为因变量，通过偏最小二乘法建立所述校正集自变量和因变量的偏最小二乘法回归模型；利用留一法进行验证；

1.6根据所述校正集偏最小二乘法回归模型的回归系数，选择对所述回归模型贡献率绝对值最大的波长为特征波长；并通过偏最小二乘法建立校正集花生中油酸含量分布定量模型；利用留一法进行验证。

2.根据权利要求1所述建立花生中油酸含量分布定量模型的方法，其特征在于，步骤1.6中所述特征波长分别为：901nm、980nm、1064nm、1147nm、1230nm、1313nm、1397nm、1480nm、1564nm、1648nm；所建立的花生中油酸含量分布定量模型如下：

Y_OFA＝10^-3×(–184.75R_901nm+43.84R_980nm–86.76R_1064nm–30.41R_1147nm

+40.89R_1230nm+18.68R_1313nm+133.66R_1397nm+99.63R_1480nm

–92.58R_1564nm–184.26R_1648nm)+791

其中，Y_OFA为花生样品的油酸含量，R_901nm、R_980nm、R_1064nm、R_1147nm、R_1230nm、R_1313nm、R_1397nm、R_1480nm、R_1564nm、R_1648nm分别为花生样品在特征波长901nm、980nm、1064nm、1147nm、1230nm、1313nm、1397nm、1480nm、1564nm、1648nm处经过预处理后的光谱反射值。

3.根据权利要求1或2所述建立花生中油酸含量分布定量模型的方法，其特征在于，所述高光谱仪扫描的波长范围为900-1700nm，扫描方式为线扫描。

4.根据权利要求1或2所述建立花生中油酸含量分布定量模型的方法，其特征在于，步骤1.2所述校正是对指对所述花生样品的原始高光谱三维图像I_raw进行黑白校正，具体方法为对反射率为99％的标准校正板进行采集，得到全白的标定图像I_white，然后关闭镜头采集，得到全黑标定图像I_dark，根据下述公式计算校正后图像I_norm：

$I_{norm}=\left(\frac{I_{raw}-I_{dark}}{I_{white}-I_{dark}}\right)\×100\%.$

5.根据权利要求1或2所述建立花生中油酸含量分布定量模型的方法，其特征在于，步骤1.2所述背景删除具体步骤为：采用主成分分析，确定背景与花生的边界，删除背景，得到花生样品图像。

6.根据权利要求1或2所述建立花生中油酸含量分布定量模型的方法，其特征在于，步骤1.4中所述检测花生样品的油酸含量方法为根据GB5413.27-2010进行；优选为根据GB5413.27-2010中第一法乙酰氯—甲醇甲酯化法进行。

7.根据权利要求1或2所述建立花生中油酸含量分布定量模型的方法，其特征在于，所述留一法是指每次留取一个花生样品作为验证，余下的n-1个花生样品建立验证模型，用所述验证模型对所述留取的这一个花生样品来进行验证，重复上述的过程，直到对所有花生样品都进行验证。

8.权利要求1-7任一项所述方法建立的花生中油酸含量分布定量模型在检测花生中油酸含量分布中的应用。

9.一种基于高光谱成像技术检测花生中油酸含量分布的方法，所述方法包括：

1)采集待测花生样品在下列特征波长处的光谱图像：901nm、980nm、1064nm、1147nm、1230nm、1313nm、1397nm、1480nm、1564nm、1648nm；

2)将所述特征波长处经过预处理后的光谱反射值，输入花生中油酸含量分布定量模型，得到待测花生样品油酸含量分布；所述花生中油酸含量分布定量模型如下：

Y_OFA＝10^-3×(–184.75R_901nm+43.84R_980nm–86.76R_1064nm–30.41R_1147nm

+40.89R_1230nm+18.68R_1313nm+133.66R_1397nm+99.63R_1480nm

–92.58R_1564nm–184.26R_1648nm)+791

其中，Y_OFA为花生样品的油酸含量，R_901nm、R_980nm、R_1064nm、R_1147nm、R_1230nm、R_1313nm、R_1397nm、R_1480nm、R_1564nm、R_1648nm分别为花生样品在特征波长901nm、980nm、1064nm、1147nm、1230nm、1313nm、1397nm、1480nm、1564nm、1648nm处经过预处理后的光谱反射值。

10.根据权利要求9所述基于高光谱成像技术检测花生中油酸含量分布的方法，其特征在于，所述步骤1)采集待测花生样品特征波长处的光谱图像的方法包括以下步骤：

1.1用高光谱仪扫描获得待测花生样品中每个像素点在各波长下的图像信息，得到待测花生样品的原始高光谱三维图像；所述高光谱仪扫描的波长范围为900-1700nm，扫描方式为线扫描；

1.2对所述待测花生样品的原始高光谱三维图像进行校正和背景删除后，提取待测花生样品图像平均光谱；所述校正是指对所述花生样品的原始高光谱三维图像I_raw进行黑白校正；具体方法为对反射率为99％的标准校正板进行采集，得到全白的标定图像I_white，然后关闭镜头采集，得到全黑标定图像I_dark，根据下述公式计算校正后图像I_norm：

$I_{norm}=\left(\frac{I_{raw}-I_{dark}}{I_{white}-I_{dark}}\right)\×100\%$

所述背景删除具体步骤为：采用主成分分析，确定背景与花生的边界，删除背景，得到花生样品图像；

1.3对所述待测花生样品图像平均光谱进行中心化、均值方差化结合标准正态变量变换预处理。