1.本发明涉及食用菌信息育种技术领域，尤其涉及一种基于显微高光谱图像秀珍菇菌种退化无损检测方法。
背景技术：2.食用菌，即大型可食用真菌，内含丰富的营养物质和微量元素。因其味道鲜美，且具有极高的医疗保健作用，故深受国内外消费者的喜爱。我国是食用菌生产的第一大国，食用菌产量约占世界食用菌总产量的70％。近年来，随着栽培技术的不断进步、品种的不断增多、工厂化种植的日益成熟，食用菌产业已经成为我国农业的第五大种植业，对我国的经济建设起着重要作用。
3.食用菌的种子被称为食用菌菌种。菌种的质量直接影响着栽培结果，即产量的高低，是产生经济效益的关键因素之一。我国秀珍菇主要培育的品种在亲缘关系方面十分接近，菌种在多年的扩繁后容易出现退化和老化情况。菌种退化是指食用菌菌种群体经济性状发生劣变的现象，从而导致食用菌产量、品质、抗性等方面发生了背离人类需要的变化。菌种退化是众多食用菌栽培中共同存在的严重问题，在实际生产中，一些优良的食用菌品种往往用不了几代就出现菌种退化现象，一旦使用退化菌种进行栽培，往往会导致出菇晚、产量低、品质差等，给食用菌生产者带来重大的经济损失。因此，快速进行退化菌种甄别及有效消除退化菌种给食用菌生产带来的隐患，防止因退化菌种流入生产环节造成的损失,一直是食用菌生产需面对和解决的问题。
4.显微高光谱成像技术是一门新型光谱技术，融合了显微光学成像技术和光谱分析技术，且已广泛应用于医学诊断、食品安全检测等领域。由于高光谱成像技术非接触、无损、快速、精准等众多优点，及显微技术可以观察肉眼所不能分辨的微小物体形态结构及其特性，高光谱成像技术和显微技术相结合使得信息采集速度更快、探索层面更深、捕捉的信息更加持久稳定，在生物信息育种技术领域具有良好的发展前景。
技术实现要素：5.考虑目前对食用菌菌种退化的研究主要是应用化学或生物技术对退化前后菌种生理信息变化进行检测分析，传统的检测方法效率低、耗时长，且整个过程操作复杂繁琐，仍然缺乏高效、快速的食用菌菌种退化检测手段，不能很好的满足生产需求的问题。
6.为了克服现有技术中存在的不足，本发明以传代过程中易发生退化的秀珍菇菌种为对象，提出一种基于显微高光谱图像的秀珍菇菌种退化的无损检测方法。
7.本发明方案所基于的基础是：秀珍菇菌种内部生理信息是判断菌种退化与否的重要指标，而菌种的光谱图像可以有效的反应菌种内部生理信息，同时不会对菌种产生损伤。因此，利用菌种的显微高光谱图像信息建立检测模型，检测速度快，准确率高，可以快速无损的对秀珍菇菌种的退化进行辨别。
8.其方案包括：步骤s1：选取秀珍菇菌种进行活化培养；步骤s2：分析已选定秀珍菇
菌种退化情况，划定菌种退化评价等级；步骤s3：分别采集退化与未退化秀珍菇菌种的显微高光谱图像；步骤s4：针对显微高光谱图像提取感兴趣区域的平均光谱透射率；步骤s5：通过卷积神经网络提取感兴趣特征作为学习特征，建立菌种退化快速检测模型，完成分类。本发明基于一维卷积神经网络算法建立模型，有效地从高维数据中学习相应特征，避免了复杂的特征提取过程，减少网络计算时间，提高模型分类的准确率，提升了秀珍菇菌种退化的检测效率，实现了秀珍菇菌种退化的快速无损检测。
9.本发明具体采用以下技术方案：
10.一种基于显微高光谱图像秀珍菇菌种退化无损检测方法，其特征在于，包括以下步骤：步骤s1：选取秀珍菇菌种进行活化培养；步骤s2：分析已选定秀珍菇菌种退化情况，划定菌种退化评价等级；步骤s3：分别采集退化与未退化秀珍菇菌种的显微高光谱图像；步骤s4：针对显微高光谱图像提取感兴趣区域的平均光谱透射率；步骤s5：通过卷积神经网络提取感兴趣特征作为学习特征，建立菌种退化快速检测模型，完成分类。
11.进一步地，在步骤s1中，将选取的秀珍菇菌种接种至固体培养基上，于25℃培养10天，进行一级菌种活化培养；
12.在步骤s2中，将活化后的秀珍菇菌种转接至液体马铃薯葡萄糖琼脂培养基中避光震荡培养5天，得到秀珍菇二级液体菌种，经4层滤纸过滤获得秀珍菇孢子液并转接至色度培养基中，测定秀珍菇菌种在色度培养基中的脱色率，划定秀珍菇菌种退化评价等级；
13.在步骤s3中，根据步骤s2中划分的秀珍菇菌种退化等级开展扩大培养，采集退化与未退化秀珍菇菌种显微高光谱图像，设定显微镜放大倍数、波长范围、曝光时间、分辨率和采集速度。
14.进一步地，步骤s4对采集到的退化与未退化的秀珍菇菌种显微高光谱图像进行处理，创建感兴趣区域，提取平均光谱透射率，完成数据集建立；
15.步骤s5将步骤s4建立的数据集作为输入数据，分为训练集和测试集，通过卷积神经网络算法建立菌种退化检测模型，分析检测模型对菌种退化与否检测的准确率，对菌种退化检测模型进行评价。
16.进一步地，在步骤s2中未退化与已退化菌种鉴别方法具体为：设定脱色情况为t＝培养液上清液的od
615
/对照液上清液的od
615
,当t》0.7时，判定该秀珍菇菌种为退化菌种；反之为未退化菌种。
17.进一步地，在步骤s3中，通过在显微镜正上方设置高光谱成像仪，采集退化与未退化秀珍菇菌种的显微高光谱图像；采集方式为透射方式，显微镜目镜倍数为10倍，物镜倍数为10倍，光源为20w卤素灯，设置曝光时间为15ms，采集速度为0.06cm/s，波长范围为400-1046nm；同时获取参考板和暗电流的显微高光谱图像数据，用于数据处理前原始光谱图像的校正，以减少光照不均匀和暗电流对光谱图像的影响。
18.进一步地，步骤s4中创建的所述感兴趣区域为图像中任意菌丝部位，像素大小为50
×
15pixels的图像范围内的光谱值，最终建立数据集为361*200的矩阵，其中包含200条光谱信息，每条光谱信息由360个波段组成。
19.进一步地，步骤s5中所述菌种退化快速检测模型采用基于一维卷积神经网络的分类算法，模型包含6个卷积层，其中过滤器维度为3，激活函数使用双曲正切函数以提高模型的表达能力；每两层卷积层后添加一个池化层，以减少模型的计算量，提高模型的计算速
度；并在卷积神经网络和全连接层中间添加flatten层使多维数据一维化，同时设置交叉熵损失函数作为模型训练的损失函数，以描述模型特征输出与期望输出之间的差异。
20.进一步地，在步骤s3中，获取空白pda培养基和暗电流的显微高光谱图像数据，根据校正公式：
[0021][0022]
对数据处理前原始光谱图像进行校正，以减少光照不均匀和暗电流对光谱图像的影响；其中r为校正后的光谱图像，i
dark
为暗电流图像，i
raw
为原始光谱图像，i
ref
为空白pda培养基光谱图像。
[0023]
进一步地，在步骤s5中，采用的激活函数的表达式如下：
[0024][0025]
其中为在wk的某个位置(i,j)处的值，wk为第k个过滤器，x为某个位置(i,j)处的值，bk为常数；
[0026]
损失函数如下所示：
[0027][0028]
其中l为交叉熵损失函数值，为标签为1的概率，为标签为0的概率，当预测输出越接近真实样本标签1，损失函数l越小；预测输出越接近0，l越大。
[0029]
相较于现有技术，本发明及其优选方案具有以下有益效果：
[0030]
通过显微高光谱图像采集装置采集秀珍菇菌种的光谱图像，从而快速有效的获取秀珍菇菌种内部生理信息的光谱特征；
[0031]
以在短时间内获取秀珍菇菌种有效光谱信息，通过模型分析快速准确辨别秀珍菇菌种退化情况，提高秀珍菇菌种退化检测的效率和准确率。
[0032]
其综合了显微高光谱成像技术信息采集速度更快、探索层面更深、捕捉的信息更加持久稳定等众多优点，利用卷积神经网络算法，对秀珍菇菌种退化与否进行分类，提高了检测的速度和准确性。该方法通过显微高光谱成像系统，利用光学设备和传感器等将人工观察和测量的数据数字化，为食用菌菌种退化的早期鉴别提供新的思路和方法。
附图说明
[0033]
下面结合附图和具体实施方式对本发明进一步详细的说明：
[0034]
图1是本发明实施例采用试验装置的结构示意图，图中：
[0035]
1-便携式高光谱成像仪；2-目镜；3-物镜；4-秀珍菇菌种样本；5-载物台；6-光源；7-显微镜；8-计算机；
[0036]
图2是本发明实施例具体操作流程图；
[0037]
图3是本发明实施例退化与未退化秀珍菇菌种平均光谱透射率曲线图；
[0038]
图4是本发明实施例秀珍菇菌种样本训练集分类检测结果图；
[0039]
图5是本发明实施例秀珍菇菌种样本测试集分类检测结果图。
具体实施方式
[0040]
为让本专利的特征和优点能更明显易懂，下文特举实施例，并配合附图，作详细说明如下：
[0041]
如图2所示，本发明的主要设计点围绕以下步骤：
[0042]
步骤s1：选取菌种进行活化培养(a)
[0043]
将选取的秀珍菇菌种接种至固体培养基上，于25℃培养10天，进行一级菌种活化培养；
[0044]
步骤s2：划定菌种退化评价等级(b)
[0045]
将活化后的秀珍菇菌种转接至液体马铃薯葡萄糖琼脂培养基(液体pda)中避光震荡培养5天，得到秀珍菇二级液体菌种，经4层滤纸过滤获得秀珍菇孢子液并转接至色度培养基中，测定秀珍菇菌种在色度培养基中的脱色率，划定秀珍菇菌种退化评价等级；
[0046]
步骤s3：采集菌种样本显微高光谱图像(c)
[0047]
根据步骤s2中划分的秀珍菇菌种退化等级开展扩大培养，采集退化与未退化秀珍菇菌种显微高光谱图像，设定显微镜放大倍数、波长范围、曝光时间、分辨率和采集速度；
[0048]
步骤s4：提取感兴趣区域平均光谱透射率(d)
[0049]
对采集到的退化与未退化的秀珍菇菌种显微高光谱图像进行处理，创建感兴趣区域，提取平均光谱透射率，完成数据集建立；
[0050]
步骤s5：建立秀珍菇菌种退化快速检测模型，完成分类(e)
[0051]
将步骤s4建立的数据集作为输入数据，分为训练集和测试集，通过卷积神经网络算法建立菌种退化检测模型，分析检测模型对菌种退化与否检测的准确率，对菌种退化检测模型进行评价。
[0052]
以下对本发明实施例的具体设计细节分别展开叙述：
[0053]
一、装置：
[0054]
如图1本方案采用的实验装置装置包括：显微镜7和试验对象秀珍菇菌种样本4。
[0055]
并设置有便携式高光谱成像仪1、目镜2、物镜3、载物台5、光源6和计算机8；
[0056]
具体连接方式为：
[0057]
显微镜7设有目镜2、物镜3、载物台5和光源6；其中物镜3置于载物台5上方，并在载物台5上放置秀珍菇菌种样本4；光源6在载物台5正下方；
[0058]
显微镜7正上方设置便携式高光谱成像仪1；
[0059]
便携式高光谱成像仪1外接计算机8；
[0060]
上述的各功能部件均为通用件。
[0061]
其工作机理是：打开计算机8和便携式高光谱成像仪1，接通显微镜7的光源6，在便携式高光谱成像仪1自带软件上设置曝光时间、分辨率以及相机的采集速度，调节载物台6的位置，获得清晰的秀珍菇菌种样本4的图像；点击软件开始采集按钮，完成秀珍菇菌种样品显微高光谱图像的采集。
[0062]
二、方法：
[0063]
1.培养基：
[0064]
步骤a1：配置固体培养基：土豆200g，葡萄糖20g，蛋白胨10g，琼脂18g，加水至1l。
[0065]
步骤a2：配置液体马铃薯葡萄糖琼脂培养基(液体pda)：土豆200g，葡萄糖20g，磷
酸二氢钾2g，硫酸镁1g，维生素b120mg，加水至1l。
[0066]
步骤a3：配置色度培养基：酵母粉4.5g，细菌蛋白胨7.5g，溴百里酚蓝0.025g，乳糖5g，加水至1l。
[0067]
上述培养基采用121℃，30min灭菌，待冷却后备用。其中维生素b高温下会失活，试验操作时过滤加入。
[0068]
2.秀珍菇退化评价：
[0069]
步骤b1：将秀珍菇菌种接种至固体培养基上，于25℃培养10天，进行一级菌种活化培养。
[0070]
步骤b2：将活化后的秀珍菇菌种转接至液体pda培养基中，每瓶液体pda培养基中接入两块直径为5mm菌块并于25℃，160r/min避光震荡培养5天，进行二级液体培养,得到秀珍菇二级液体菌种，经4层滤纸过滤获得秀珍菇孢子液转接至色度培养基中，25℃、120r/min条件下震荡培养10天，同一水平3个重复，对照组为空白pda培养液，分别将培养液和对照液离心，12000r/min、5min离心，取上清液测定od615nm值。
[0071]
步骤b3：根据测定的吸光度计算脱色情况t＝培养液上清液的od
615
/对照液上清液的od
615
，并设定当t》0.7时,判断该秀珍菇菌种为退化菌种。
[0072]
3.数据采集：
[0073]
根据步骤s2中划分的秀珍菇菌种退化等级开展扩大培养，采用显微高光谱图像采集装置采集退化与未退化秀珍菇菌种显微高光谱图像，其中目镜倍数为10倍，物镜倍数为10倍，光源为20w卤素灯，设置曝光时间为15ms，采集速度为0.06cm/s，同时获取空白pda培养基和暗电流的显微高光谱图像数据，根据校正公式：
[0074][0075]
其中r为校正后的光谱图像，i
dark
为暗电流图像，i
raw
为原始光谱图像，i
ref
为空白pda培养基光谱图像，对数据处理前原始光谱图像进行校正，以减少光照不均匀和暗电流对光谱图像的影响。
[0076]
4.数据提取：
[0077]
将校正后的显微高光谱图像导入envi中，创建感兴趣区域为任意菌丝部分，提取感兴趣区域的平均光谱透射率并保存，获取的数据为400nm-1046nm波长下的平均光谱透射率。
[0078]
图3为提取的退化与未退化菌种平均透射率光谱曲线。
[0079]
将提取保存的平均光谱透射率导入到excel的表格中，设置每行表示一个菌种的平均光谱反射率信息，每列为不同波长下的光谱反射率，最后一列设为菌种退化真实类别，完成数据集建立。
[0080]
5.模型建立
[0081]
在matlab环境中，将步骤s4中建立的数据集作为输入数据，设置x是数据的特征信息，y表示真实类别，其中x在读取的时候矩阵增加一维，并以7：3的比例划分为训练集和验证集。
[0082]
一维卷积神经网络模型在matlab环境中，通过卷积层的卷积操作从复杂的数据集中提取出感兴趣特征，使用双曲正切函数作为激活函数提高模型的表达能力，采用的激活
函数的表达式如下：
[0083][0084]
其中为在wk的某个位置(i,j)处的值，wk为第k个过滤器，x为某个位置(i,j)处的值，bk为常数。
[0085]
同时为了保证主要特征仍旧保留以及降低计算量，每两个卷积层后添加一个最大池化层。以输入端输入数据大小360*1为例，通过16个维度为3的过滤器，经过两次卷积操作，输出数据大小为16@356*1，通过池化除去冗余的信息后，输出数据大小为16@89*1。
[0086]
在卷积层和全连接层之间加入flatten层，将多维数据一维化。
[0087]
使用交叉熵损失函数作为模型训练的损失函数，描述训练模型的特征输出与实际标定输出之间的差异，损失函数越小，两者之间差异越小，模型准确率越高。损失函数如下所示：
[0088][0089]
其中l为交叉熵损失函数值，为标签为1的概率，为标签为0的概率，当预测输出越接近真实样本标签1，损失函数l越小；预测输出越接近0，l越大。
[0090]
三、试验结果：
[0091]
将划分的测试集输入到上述基于一维卷积神经网络的分类模型中，由模型中的卷积层完成对秀珍菇菌种平均光谱信息的特征提取，池化层完成对平均光谱信息的降维，最终由全连接层完成对秀珍菇菌种样本的分类识别，并通过混淆矩阵直观的显示分类效果。
[0092]
图4、5分别为秀珍菇菌种样本训练集与测试集的分类检测结果。
[0093]
从图中可以看出，经过800次迭代训练后，秀珍菇菌种退化检测模型训练集的平均准确率为95.7％，测试集的平均准确率为93.3％，所建的检测模型得到了较好的检测结果。
[0094]
本专利不局限于上述最佳实施方式，任何人在本专利的启示下都可以得出其它各种形式的基于显微高光谱图像秀珍菇菌种退化无损检测方法，凡依本发明申请专利范围所做的均等变化与修饰，皆应属本专利的涵盖范围。
技术特征：1.一种基于显微高光谱图像秀珍菇菌种退化无损检测方法，其特征在于，包括以下步骤：步骤s1：选取秀珍菇菌种进行活化培养；步骤s2：分析已选定秀珍菇菌种退化情况，划定菌种退化评价等级；步骤s3：分别采集退化与未退化秀珍菇菌种的显微高光谱图像；步骤s4：针对显微高光谱图像提取感兴趣区域的平均光谱透射率；步骤s5：通过卷积神经网络提取感兴趣特征作为学习特征，建立菌种退化快速检测模型，完成分类。2.根据权利要求1所述的基于显微高光谱图像秀珍菇菌种退化无损检测方法，其特征在于：在步骤s1中，将选取的秀珍菇菌种接种至固体培养基上，于25℃培养10天，进行一级菌种活化培养；在步骤s2中，将活化后的秀珍菇菌种转接至液体马铃薯葡萄糖琼脂培养基中避光震荡培养5天，得到秀珍菇二级液体菌种，经4层滤纸过滤获得秀珍菇孢子液并转接至色度培养基中，测定秀珍菇菌种在色度培养基中的脱色率，划定秀珍菇菌种退化评价等级；在步骤s3中，根据步骤s2中划分的秀珍菇菌种退化等级开展扩大培养，采集退化与未退化秀珍菇菌种显微高光谱图像，设定显微镜放大倍数、波长范围、曝光时间、分辨率和采集速度。3.根据权利要求1所述的基于显微高光谱图像秀珍菇菌种退化无损检测方法，其特征在于：步骤s4对采集到的退化与未退化的秀珍菇菌种显微高光谱图像进行处理，创建感兴趣区域，提取平均光谱透射率，完成数据集建立；步骤s5将步骤s4建立的数据集作为输入数据，分为训练集和测试集，通过卷积神经网络算法建立菌种退化检测模型，分析检测模型对菌种退化与否检测的准确率，对菌种退化检测模型进行评价。4.根据权利要求1所述的基于显微高光谱图像秀珍菇菌种退化无损检测方法，其特征在于：在步骤s2中未退化与已退化菌种鉴别方法具体为：设定脱色情况为t＝培养液上清液的od
615
/对照液上清液的od
615
,当t>0.7时，判定该秀珍菇菌种为退化菌种；反之为未退化菌种。5.根据权利要求1所述的基于显微高光谱图像秀珍菇菌种退化无损检测方法，其特征在于：在步骤s3中，通过在显微镜正上方设置高光谱成像仪，采集退化与未退化秀珍菇菌种的显微高光谱图像；采集方式为透射方式，显微镜目镜倍数为10倍，物镜倍数为10倍，光源为20w卤素灯，设置曝光时间为15ms，采集速度为0.06cm/s，波长范围为400-1046nm；同时获取参考板和暗电流的显微高光谱图像数据，用于数据处理前原始光谱图像的校正，以减少光照不均匀和暗电流对光谱图像的影响。6.根据权利要求1所述的基于显微高光谱图像秀珍菇菌种退化无损检测方法，其特征在于：步骤s4中创建的所述感兴趣区域为图像中任意菌丝部位，像素大小为50
×
15pixels的图像范围内的光谱值，最终建立数据集为361*200的矩阵，其中包含200条光谱信息，每条光谱信息由360个波段组成。7.根据权利要求1所述的基于显微高光谱图像秀珍菇菌种退化无损检测方法，其特征在于：步骤s5中所述菌种退化快速检测模型采用基于一维卷积神经网络的分类算法，模型包含6个卷积层，其中过滤器维度为3，激活函数使用双曲正切函数以提高模型的表达能力；每两层卷积层后添加一个池化层，以减少模型的计算量，提高模型的计算速度；并在卷积神经网络和全连接层中间添加flatten层使多维数据一维化，同时设置交叉熵损失函数作为
模型训练的损失函数，以描述模型特征输出与期望输出之间的差异。8.根据权利要求5所述的基于显微高光谱图像秀珍菇菌种退化无损检测方法，其特征在于：在步骤s3中，获取空白pda培养基和暗电流的显微高光谱图像数据，根据校正公式：对数据处理前原始光谱图像进行校正，以减少光照不均匀和暗电流对光谱图像的影响；其中r为校正后的光谱图像，i
dark
为暗电流图像，i
raw
为原始光谱图像，i
ref
为空白pda培养基光谱图像。9.根据权利要求7所述的基于显微高光谱图像秀珍菇菌种退化无损检测方法，其特征在于：在步骤s5中，采用的激活函数的表达式如下：其中为在w
k
的某个位置(i,j)处的值，w
k
为第k个过滤器，x为某个位置(i,j)处的值，b
k
为常数；损失函数如下所示：其中l为交叉熵损失函数值，为标签为1的概率，为标签为0的概率，当预测输出越接近真实样本标签1，损失函数l越小；预测输出越接近0，l越大。
技术总结本发明提出一种基于显微高光谱图像秀珍菇菌种退化无损检测方法，包括：步骤S1：选取秀珍菇菌种进行活化培养；步骤S2：分析已选定秀珍菇菌种退化情况，划定菌种退化评价等级；步骤S3：分别采集退化与未退化秀珍菇菌种的显微高光谱图像；步骤S4：针对显微高光谱图像提取感兴趣区域的平均光谱透射率；步骤S5：通过卷积神经网络提取感兴趣特征作为学习特征，建立菌种退化快速检测模型，完成分类。本发明基于一维卷积神经网络算法建立模型，有效地从高维数据中学习相应特征，避免了复杂的特征提取过程，减少网络计算时间，提高模型分类的准确率，提升了秀珍菇菌种退化的检测效率，实现了秀珍菇菌种退化的快速无损检测。菇菌种退化的快速无损检测。菇菌种退化的快速无损检测。
技术研发人员：魏萱 刘诗洋 李思怡
受保护的技术使用者：福建农林大学
技术研发日：2022.07.25
技术公布日：2022/11/1
