检测苹果可溶性固体含量的方法

摘要

本发明提供一种检测苹果可溶性固体含量的方法，其特征在于，包括以下步骤：准备苹果样品步骤；调节高光谱成像仪步骤；黑白版校正步骤；测量苹果样品的光谱反射率步骤；测量其他样品的光谱反射率步骤；用国标方法测定苹果样品可溶性固体含量步骤；建立数学回归模型步骤；以及苹果可溶性固体含量的测算程序化步骤。根据本发明所提供的检测苹果可溶性固体含量的方法，只需采用高光谱成像仪测得苹果的光谱反射率信息，带入线性回归模型中，即可计算得出该苹果的可溶性固体含量，因此该方法无需破坏苹果就可以快速准确地测算出苹果的可溶性固体含量。

说明书

技术领域

本发明涉及可溶性固体含量的测定，具体涉及一种基于光谱分析 的无损检测苹果可溶性固体含量的方法。

背景技术

中国是苹果生产大国，苹果质量等级的准确判定一直是中国非常 重视的问题，这也是中国苹果产业化生产的关键。苹果的可溶性固体 含量是衡量苹果品质的重要因素之一，而且苹果中的可溶性固体含量 直接影响到苹果制品的美味。

目前，国内外检测苹果可溶性固体含量的方法，大多需要将苹果 损坏。国家标准的测量苹果可溶性固体含量的方法是：首先将苹果切 碎，然后将苹果碎块放入高速组织匀浆机中捣成浆，最后提取浆液， 采用糖度计测定苹果的可溶性固体含量。这种方法虽然检测结果比较 准确，但是操作非常复杂耗时，而且只能进行抽样检测，从而判断某 类苹果的可溶性固体含量的大致范围，根本无法实现测量每个苹果样 品的可溶性固体含量。

另外，目前还没有适用于市场上检测苹果可溶性固体含量的技 术，经验人士往往通过观察外观来判断苹果的可溶性固体含量，这样 的判断结果存在很大的误差，也成为苹果产业发展的瓶颈。

因此，发展一种快速无损检测苹果可溶性固体含量的方法是当今 市场的迫切需求。

发明内容

本发明是针对上述问题进行的，目的在于提供一种不需要破坏苹 果，即可快速准确地测定苹果可溶性固体含量的方法。

本发明为实现上述目的，采用了以下的技术方案：

本发明提供一种检测苹果可溶性固体含量的方法，其特征在于， 包括以下步骤：准备苹果样品步骤，从一批苹果中抽取若干个苹果样 品，将该苹果样品清洗干净，在0℃下保存至少12小时；调节高光 谱成像仪步骤，将一个苹果样品放置在高光谱成像仪的水平位移台 上，调节高光谱相机的高度和光源的强度以及矩形聚光透镜照射在苹 果样品上的角度，使苹果样品表面处在矩形聚光透镜的焦平面处，调 节矩形聚光透镜的景深，使高光谱相机能够拍摄到苹果样品整个曲面 的清晰图像；黑白板校正步骤，移开苹果样品，保持光源的强度不变， 将白板放置在焦平面处，记录白板在光波长为400～1000nm范围内的 光谱反射率r₁，然后将矩形聚光透镜用黑板封闭，记录黑板在光波长 为400～1000nm范围内的光谱反射率r₂；测量苹果样品的光谱反射率 步骤，移开黑板和白板，将苹果样品放置在焦平面处，用黑板封闭矩 形聚光透镜，记录样品黑板光谱反射率r₀，然后移开黑板，通过移动 水平位移平台，扫描苹果样品的检测面，得到苹果样品在光波长为 400～1000nm范围内的光谱反射率r_s，按如下公式校正苹果样品的光 谱反射率r_s：式中R即为苹果样品校正后的光谱反射率； 测量其他样品的光谱反射率步骤，重复测量苹果样品的光谱反射率步 骤，测量其他的苹果样品在光波长为400～1000nm范围内的光谱反射 率，计算校正后的光谱反射率，计算所有的苹果样品在光波长 400～1000nm范围内各个波长所对应的光谱反射率平均值，作光波长 -光谱反射率平均值曲线；用国标方法测定苹果样品可溶性固体含量 步骤，采用国标ISO2173规定的水果和蔬菜制品可溶性固体含量的 测定方法测定所有的苹果样品的可溶性固体含量；建立数学回归模型 步骤，根据光波长-光谱反射率平均值曲线在400～1000nm范围内选 取至少37个特征波长，将该37个特征波长下所有的苹果样品校正后 的光谱反射率与采用国标方法测得的相应的苹果样品可溶性固体含 量进行关联，运用偏最小二乘回归法建立数学回归模型，以选取37 个特征波长下的光谱反射率为例，得到数学回归模型如下：

Y＝11.9607537294788-0.143642520620794·X₁-0.105863600481352 ·X₂+0.109246444282623·X₃-0.254882848332586·X₄-2.78691407831185 E-02·X₅-7.86541002189555E-02·X₆-7.24996914024644E-02·X₇-0.1443 04439911315·X₈-5.07006523779758E-02·X₉+5.95223777462032E-02·X ₁₀-5.62302521068629E-02·X₁₁-2.35601365841085E-03·X₁₂-5.211636820 57966E-02·X₁₃+1.83828562670345E-02·X₁₄-0.140737132450369·X₁₅-0. 049427159554811·X₁₆-9.34192459275574E-02·X₁₇+3.75959418377708 E-02·X₁₈+0.140441812259277·X₁₉+0.130902438187563·X₂₀+0.3750001 35236033·X₂₁+0.407714589795041·X₂₂-0.245866506626271·X₂₃+0.129 641648557972·X₂₄-0.38961138317635·X₂₅-0.210467060864052·X₂₆-1.1 3028383697272E-02·X₂₇+0.134858971155099·X₂₈+2.92203886573371E -04·X₂₉+0.189169283541164·X₃₀-0.219903013635167·X₃₁-0.176347572 189998·X₃₂-3.82129235904802E-02·X₃₃-0.149341892845856·X₃₄+0.602 969537823784·X₃₅+0.216166784146167·X₃₆+0.143622931484571·X₃₇，

式中Y是苹果可溶性固体含量的模型预测值，X₁～X₃₇分别是苹 果样品在37个特征波长下的校正后的光谱反射率；苹果可溶性固体 含量的测算程序化步骤，将数学回归模型编写为MATLAB的计算程 序，对其他任何苹果，只需将该苹果在0℃下保存至少12小时，采 用高光谱成像仪测出苹果在光波长400～1000nm范围内的光谱反射 率，提取上述37个特征波长下的光谱反射率，并输出为txt文件，在 MATLAB中读入txt文件，并运行计算程序，即可计算得到苹果的可 溶性固体含量。

本发明所提供的检测苹果可溶性固体含量的方法，还可以具有这 样的特征：其中，各个苹果在0℃下保存相同的时间。

发明的作用与效果

根据本发明所提供的检测苹果可溶性固体含量的方法，因为取若 干苹果样品，采用高光谱成像仪测得苹果样品的光谱反射率信息，取 至少37个特征光波长下的光谱反射率值，与国标方法测得的苹果可 溶性固体含量平均值进行关联，得到线性回归模型，所以在以后的检 测中，只需采用高光谱成像仪测得苹果的光谱反射率信息，提取特征 波长下的光谱反射率值，带入线性回归模型中，即可计算得出该苹果 的可溶性固体含量，因此该方法无需破坏苹果就可以快速准确地测算 出苹果的可溶性固体含量，非常适合应用于苹果以及苹果制品质量等 级的判定。

附图说明

图1是高光谱成像仪的结构示意图；

图2是苹果样品的光波长-光谱反射率平均值曲线；以及

图3是苹果样品可溶性固体含量的模型预测值与国标测试值的 对比图。

具体实施方式

以下结合附图，对本发明所提供的检测苹果可溶性固体含量的方 法作详细阐述。

<实施例>

图1是高光谱成像仪的结构示意图。

如图1所示，高光谱成像仪10通过线扫描对物体进行高光谱成 像，包括水平位移台11、支架12、高光谱相机13、矩形聚光透镜14、 矩形聚光透镜15、玻璃光纤16、玻璃光纤17以及光源18。

水平位移台11放置在水平台(图中未示出)上。支架12是十字 形支架，支架12的竖直部分的底端与水平位移台11固定相连，支架 12的水平部分位于竖直部分的中部。高光谱相机13可移动地连接在 支架12的竖直部分的顶部。矩形聚光透镜14和矩形聚光透镜15的 中轴分别连接在支架12的水平部分的左右两端，并可以绕各自的中 轴转动。玻璃光纤16和玻璃光纤17的一端分别与矩形聚光透镜14 和矩形聚光透镜15相连，玻璃光纤16和玻璃光纤17的另一端与光 源18相连。

本实施例所使用的高光谱成像仪10中高光谱相机13的型号为 VNIR Imspector V10E，光源18是功率可调的卤素灯。测试过程中， 测试者首先将苹果样品放置在水平位移台11上，然后打开光源18和 高光谱相机13，光源18发出的光被玻璃光纤16和玻璃光纤17分成 两束光，这两束光分别被矩形聚光透镜14和矩形聚光透镜15汇聚到 苹果样品表面，高光谱相机13接收被苹果样品反射的光线，得到苹 果样品的高光谱图像和光谱反射率信息。

本发明所提供的检测苹果可溶性固体含量的方法是基于苹果样 品的光谱反射率信息来测定苹果的可溶性固体含量，具体步骤如下：

准备苹果样品步骤，从当地超市购买山东乔纳金、山东红富士和 新疆阿克苏三种苹果各99个，作为苹果样品，将所有的苹果样品清 洗干净，在0℃下保存12个小时。

调节高光谱成像仪10步骤，取一个苹果样品，放置在高光谱成 像仪10的水平位移台11上，调节矩形聚光透镜14和矩形聚光透镜 15的角度，使苹果样品的表面处在透镜焦平面处，调节两个矩形聚 光透镜的景深，使高光谱相机13能够拍摄到整个苹果样品的清晰像， 调节高光谱相机13的曝光时间使相机的拍摄效果最佳。

黑板报校正步骤，移开苹果样品，保持光源18的强度不变，将 白板放置在水平位移台11上，记录白板在光波长为400～1000nm范 围内的光谱反射率r₁，然后将矩形聚光透镜14和矩形聚光透镜15均 用黑板封闭，记录黑板在光波长为400～1000nm范围内的光谱反射率 r₂。

测量苹果样品的光谱反射率步骤，移开黑板和白板，将苹果样品 放置在水平位移台11上，用黑板封闭矩形聚光透镜14和矩形聚光透 镜15，记录样品黑板光谱反射率r₀，然后移开黑板，通过移动水平位 移平台11，扫描苹果样品的检测面，得到苹果样品在光波长为 400～1000nm范围内的光谱反射率r_s，按如下公式校正苹果样品的光 谱反射率r_s，得到校正后苹果样品的光谱反射率R：

测量其他苹果样品的光谱反射率步骤，重复上述测量苹果样品的 光谱反射率步骤，分别测量每种99个苹果样品在光波长为400～1000 nm范围内的光谱反射率，计算相应的校正后的光谱反射率，然后计 算出每种99个苹果样品在光波长为400～1000nm范围内校正后的光 谱反射率平均值。

图2是苹果样品的光波长-光谱反射率平均值曲线。

如图2所示，求得三种苹果样品的光谱反射率平均值后，分别将 光谱反射率平均值对光波长作图，即得到三种苹果样品的光波长-光 谱反射率平均值曲线。

用国标方法测定苹果样品可溶性固体含量步骤，根据国标ISO 2173规定的水果和蔬菜制品可溶性固体含量的测定方法，将苹果样 品切碎，再放到高速组织匀浆机中将苹果碎块捣成浆，然后提取浆液， 采用糖度计测量苹果样品的可溶性固体含量，采用该方法测得所有苹 果样品的可溶性固体含量后，分别计算每种苹果的可溶性固体含量平 均值，计算结果如表1所示；

表1

建立数学回归模型步骤，根据苹果样品的光波长-光谱反射率平 均值曲线，在有吸收峰等位置选取37个特征波长，在本实施例中， 37个特征波长分别为：628.738nm、629.997nm、629.368nm、628.108 nm、630.627nm、631.258nm、626.849nm、627.478nm、631.888nm、 632.518nm、626.219nm、625.589nm、633.148nm、624.960nm、633.779 nm、624.331nm、635.040nm、635.670nm、623.701nm、618.670nm、 640.718nm、641.349nm、1017.012nm、1013.845nm、1018.912nm、 1016.379nm、1012.444nm、999.255nm、1006.871nm、1014.478nm、 1019.545nm、1017.646nm、1020.678nm、1018.279nm、997.985nm、 1011.944nm、618.042nm。

选择198个苹果样品的校正后的光谱反射率信息，将198个苹果 样品在37个特征波长下的校正后光谱反射率与采用国标方法测得的 相应的各个苹果样品的可溶性固体含量进行关联，运用偏最小二乘回 归法建立数学回归模型，得到数学回归模型如下：

Y＝11.9607537294788-0.143642520620794·X₁-0.105863600481352 ·X₂+0.109246444282623·X₃-0.254882848332586·X₄-2.78691407831185 E-02·X₅-7.86541002189555E-02·X₆-7.24996914024644E-02·X₇-0.1443 04439911315·X₈-5.07006523779758E-02·X₉+5.95223777462032E-02·X ₁₀-5.62302521068629E-02·X₁₁-2.35601365841085E-03·X₁₂-5.211636820 57966E-02·X₁₃+1.83828562670345E-02·X₁₄-0.140737132450369·X₁₅-0. 049427159554811·X₁₆-9.34192459275574E-02·X₁₇+3.75959418377708 E-02·X₁₈+0.140441812259277·X₁₉+0.130902438187563·X₂₀+0.3750001 35236033·X₂₁+0.407714589795041·X₂₂-0.245866506626271·X₂₃+0.129 641648557972·X₂₄-0.38961138317635·X₂₅-0.210467060864052·X₂₆-1.1 3028383697272E-02·X₂₇+0.134858971155099·X₂₈+2.92203886573371E -04·X₂₉+0.189169283541164·X₃₀-0.219903013635167·X₃₁-0.176347572 189998·X₃₂-3.82129235904802E-02·X₃₃-0.149341892845856·X₃₄+0.602 969537823784·X₃₅+0.216166784146167·X₃₆+0.143622931484571·X₃₇，

式中Y是苹果可溶性固体含量的模型预测值，X₁～X₃₇分别是苹 果样品在400～1000nm范围内37个特征波长下的校正后光谱反射率。

图3是苹果样品可溶性固体含量的模型预测值与国标测试值的 对比图。

如图3所示，三种苹果样品按照上述数学回归模型计算出的可溶 性固体含量的模型预测值与国标方法测定的国标测试值非常接近，该 数学回归模型的相关性达到0.9349。

苹果可溶性固体含量的测算程序化步骤，将上述数学回归模型编 写为MATLAB的计算程序，在以后需要检测苹果的可溶性固体含量 时，只需将该苹果洗净后在0℃下保存12小时，然后重复上述测量 苹果样品的光谱反射率步骤，得到苹果在光波长400～1000nm范围内 的光谱反射率信息，提取37个特征波长下的光谱反射率值，并输出 为txt文件，在MATLAB中读入该txt文件，运行上述计算程序，即 可计算得到该苹果的可溶性固体含量值。

实施例的作用与效果

根据本实施例所提供的检测苹果可溶性固体含量的方法，因为采 用高光谱成像仪测得每种99个苹果样品的光谱反射率信息，取37个 特征光波长下的校正后光谱反射率，与国标方法测得的相应的苹果样 品的可溶性固体含量进行关联，得到线性回归模型，所以在以后的检 测中，只需采用高光谱成像仪测得苹果的光谱反射率信息，提取37 个特征波长下的光谱反射率带入线性回归模型中，即可计算得出该苹 果的可溶性固体含量，因此该方法无需破坏苹果就可以快速准确地测 算出苹果的可溶性固体含量。

另外，本实施例中选取37个特征光波长下的光谱反射率来作偏 最小二乘回归分析，得到的数学回归模型相关性超过0.93，既保证了 数学回归模型的可靠性，又避免了选取过多参数导致运算过于复杂。

当然，本发明所提供的检测苹果可溶性固体含量的方法并不仅仅 限定于以上实施例所述的内容。以上仅为本发明构思下的基本说明， 而依据本发明的技术方案所作的任何等效变换，均应属于本发明的保 护范围。

另外，在实施例中选择37个特征波长所对应的光谱反射率与国 标测试值进行偏最小二乘回归分析，得到数学回归模型，还可以选择 更多的特征波长下的光谱反射率用于最小二乘回归分析，所得到的数 学回归模型相关性更高，但是数据处理过程比较复杂。