基于水果表面不同组织热特性的水果损伤检测方法及装置的制造方法
【技术领域】
[0001] 本发明涉及一种食品检测方法及装置，尤其是涉及一种基于水果表面不同组织热 特性的水果损伤检测方法及装置。
【背景技术】
[0002] 有无损伤是水果分级的重要依据之一，在世界各国的水果评级标准中有严格的规 定。机械损伤是引起水果采后损耗的主要因素。机械损伤不仅会造成水果感官变化，降低 品质，而且增加了微生物侵染的危险性，使水果腐烂增加，货架期缩短，严重影响品质及其 经济效益（虢露葭，李萍，侯晓荣等.果蔬采后机械损伤特性研究进展.食品工业科技 [J]，2013, 34 (1) : 389-91，96.)。更为严重的是，损伤组织还会为病原菌提供滋养繁衍的场 所，进而引起正常水果的损坏，扩大经济损失。如果相关制品流入市场，还构成了食品安全 隐患。
[0003] 根据周其显等人（周其显.苹果早期机械损伤的红外热成像检测研究[D].南昌： 华东交通大学，2011.)的研究，在理想情况下，对于苹果某一特定组织的过余温度比的对数 是关于时间的一次函数。
[0004] 现有文献中出现了几种水果损伤的无损检测方法，包括机器视觉技术、光谱技术、 光谱成像技术和红外热成像技术等。但是这些无损检测方法检测损伤稳健性不强，且成本 较高，而且如何在识别损伤的同时，避免将果梗/花萼误判为损伤仍是一大难题。
【发明内容】
[0005] 本发明的目的在于提出了一种基于水果表面不同组织热特性的水果损伤检测方 法及装置，基于水果表面不同组织的热特性检测水果表面而判断损伤。
[0006] 本发明解决其技术问题所采用的技术方案是：
[0007] -、一种水果损伤检测方法，该方法的步骤如下：
[0008] 1)通过热电堆传感器采集水果样品的不同感兴趣区域的温度数据；
[0009] 2)对温度数据，依次进行压缩、傅里叶变换分析后，提取曲线峰值作为特征值；
[0010] 3)将特征值输入到最小二乘法线性分类器进行分类训练；
[0011] 4)对于待测水果，重复步骤1)~2)得到其特征值，输入到分类训练后的最小二 乘法线性分类器中进行判别得到损伤检测结果。
[0012] 所述的感兴趣区域包括水果的损伤部位、花萼部位和正常组织部位。
[0013] 所述的不同感兴趣区域的温度数据包括热电堆传感器中每个传感元件多次采集 的温度值。
[0014] 对于所述步骤2)中每次采集的温度数据具体采用以下方式处理：
[0015] 2. 1)温度场压缩：
[0016] 取热电堆传感器中每列传感元件每次采集的温度平均值形成一行矩阵，然后矩阵 转置为列矩阵，再将列矩阵按采集的时间次序排列得到温度场压缩数据矩阵，作为压缩的 温度场信号Twm表示数据矩阵的行数，i表示数据矩阵的列数，即数据采集的次数；
[0017] 2. 2)温度场时序压缩：
[0018] 对压缩的温度场信号。以每相邻的八个列矩阵为一组，计算其中每行的平均值， 将每相邻的八个列矩阵合并为一列得到温度场时序压缩数据矩阵，作为压缩的温度场时序 信号Dm_j，j表不压缩的温度场时序信号的列数；
[0019] 2. 3)特征峰提取：
[0020] 对压缩的温度场时序信号的每一行分别进行傅氏变换，得到温度场傅氏序列 ，再计算得到温度场傅氏序列的幅值，取幅值中的最高峰、次高峰、第三高峰的 值作为每行的特征元素记录到特征峰序列ΒΛ，k表示高峰种类，k= 0, 1，2分别表示最高 峰、次高峰和第三高峰；
[0021] 2. 4)特征值计算：
[0022] 取特征峰序列的每列平均值，得到特征峰均值序列ΒΜη，η= 0, 1，2,η表示特 征峰均值序列ΒΜη的列数，将特征峰均值序列ΒΜη中的所有特征峰均值相加得到第一特征值 BS1，将最高峰、次高峰的特征峰均值相加得到第二特征值BS2。
[0023] 将不同水果样品的不同感兴趣区域所有采集温度数据处理后得到的所述第一特 征值BS1和第二特征值BS2输入到最小二乘法线性分类器中进行分类训练，得到水果表面 损伤检测模型，由水果表面损伤检测模型对待测水果进行检测。
[0024] 二、一种水果损伤检测装置：
[0025] 包括托座、热电堆传感器、单片机和PC机，水果样品放置在托座顶端，托座底端中 心安装有热电堆传感器，热电堆传感器穿过托座底端朝向水果样品底部，热电堆传感器顶 部与水果样品的下部之间有间隙，热电堆传感器经单片机和PC机连接；热电堆传感器将采 集到的原始数据传输给单片机，单片机将原始数据转换为温度数据并传输给PC机，PC机处 理温度数据并显示。
[0026] 所述的热电堆传感器采用红外阵列热电堆传感器，包含有像素IR阵列。
[0027] 所述的热电堆传感器顶部与水果样品的下部之间的间隙为1~2cm。
[0028] 本发明具有的有益的效果是：
[0029] 本发明利用水果表面不同组织的热特性识别损伤部位，可靠性好，准确率高。
【附图说明】
[0030] 图1是本发明装置的结构原理示意图。
[0031] 图2是本发明托座和热电堆传感器的安装示意图。
[0032] 图3是本发明的检测流程图。 图4是本发明实施例步骤E)得到的两个特征值的离散表示图。
[0033] 图中：1、样品，2、托座，3、热电堆传感器，4、单片机，5、PC机。
【具体实施方式】
[0034] 下面结合附图和实施例对本发明作进一步说明。
[0035] 如图1所示，本发明包括托座2、热电堆传感器3、单片机4和PC机5,水果样品1 放置在托座2顶端，托座2底端中心安装有热电堆传感器3,热电堆传感器3穿过托座2底 端朝向水果样品1底部，热电堆传感器3顶部与水果样品1的下部之间有间隙，热电堆传感 器3经单片机4和PC机5连接；热电堆传感器3将采集到的原始数据传输给单片机4,单 片机4将原始数据转换为温度数据并传输给PC机5,PC机5处理温度数据并显示。
[0036] 热电堆传感器3采用红外阵列热电堆传感器，包含有像素IR阵列。如图2所示， 热电堆传感器3 -端插入托座2中，另一端通过电线连接在单片机4上。热电堆传感器3 顶部与水果样品1的下部之间的间隙为1~2cm。
[0037] 工作时，热电堆传感器采集样品原始数据，将数据传输至单片机4,单片机4接收 原始数据并将其转换为温度数据，再将温度数据传输至PC机5进行处理，PC机5处理温度 数据并显示。温度数据处理的过程是：首先对温度数据进行傅氏变换，然后提取峰值数据作 简单加和运算后作为特征量，最后将特征量代入最小二乘分类器进行分类。
[0038] 如图3所示，本发明的具体实施例及其实施过程如下：
[0039] 在本例中，热电堆传感器3采用MLX90620,MLX90620采集的样本温度场数据为4 行16列的矩阵形式，单片机4可采用ArduinoMEGA2560,水果采用黄冠梨。
[0040] 将黄冠梨样品1放置于托座2上，按每秒1次的速率分别采集损伤部位、正常组织 部位和花萼部位温度，每个部位均采集1024组温度数据。对黄冠梨样本所采集的1024点 数据作如下处理：
[0041] A)温度场压缩。取热电堆传感器中每列传感元件每次采集的温度平均值形成一行 矩阵，