本发明涉及一种稻种发芽率无损检测方法，特别是涉及一种基于荧光光谱法的稻种发芽率检测方法，属于稻种无损检测领域。
背景技术：
我国是水稻生产大国，种植面积居世界第一。水稻种子的质量直接影响水稻产量，因此，加强稻种质量的检验对确保水稻产量具有重要意义。其中，发芽率是水稻种子检测的重要指标之一。传统的稻种发芽率检测如四唑(2，3，5-triphenyltetrazoliumchloride，ttc)染色法(氯化三苯基四氮唑染色法)，浸种时间长且ttc溶液易受光、温度等环境因素干扰，存在效率低、精度差、有损等问题；为降低人为检测误差并提高检测效率，采用机器视觉方法结合四唑法对染色后的水稻种子进行图像采集和识别，但图像处理速度慢且专业化要求高；近红外光谱法检测灵敏度差，且需采取不同的预处理方法才能获取不同稻种样品特异的吸收峰。热红外、偏振光及光声光谱等方法检测效果较好，但特征光谱较宽，有用信息较为分散，基于此技术进行产品化开发时光学部分较为复杂。
荧光光谱法通过检测特定波长激发光下物质产生的荧光强度对物质含量进行定性或定量的分析，因其扫描速度快、精度高且使用方便，已广泛应用于无损检测中，如采用新荧光扫描法对大败毒胶囊中小檗碱含量进行测定，检测精度约为紫外扫描法的8倍；应用日立f-4600型荧光分光光度计对貉体内砷含量进行测定，测量误差小于0.5μg；通过建立最佳激发波长的选择方法对八子补肾胶囊中蛇床子素含量进行测定，检测灵敏度较原波长的提高4倍；采用荧光分光光度计在547nm的激发波长下对水样品中罗丹明b的含量进行测定，在免去提取过程中有机溶剂使用的同时提高了检测效率；基于荧光共振能力转移(fret)系统对乳制品中三聚氰胺(mb)含量进行测定，可准确得到荧光强度与mb间的相关系数为0.01亩摩尔/升。因此亟需研究一种高精度、高效率、无损检测稻种发芽率新方法。
技术实现要素：
为克服现有技术的缺陷，本发明提出一种基于荧光光谱的稻种发芽率无损检测方法，具有高精度、高效率等优势。
为实现上述目的，本发明采用以下技术方案：
本发明所述的一种基于荧光光谱的稻种发芽率无损检测方法：
步骤1：将建模所需的不同发芽率稻种样品置入装有蒸馏水的试管中浸泡；
步骤2：取出试管中的浸泡液置于荧光分光光度计的样品池中，扫描浸泡液得到波长范围为360～650nm的荧光光谱；
步骤3：利用小波去噪对步骤2中得到的荧光光谱进行平滑处理；
步骤4：利用主成分分析(pca)提取步骤3平滑处理后的荧光光谱中的特征波段；
步骤5：对建模所需的不同发芽率的稻种样品按照步骤1至步骤4进行操作，得到特征波段荧光光强pi(为第i个样品的特征波段荧光光强，i＝1，2…n，n为建模样本数量)；对建模用的不同发芽率的稻种样品按照gb/t3543.4表1规定的条件(20℃)和时间(5d)进行发芽试验，得到不同稻种样品的发芽率gi(为第i个样品的发芽率，i＝1，2…n，n为建模样本数量)；
步骤6：基于深度神经网络(dnn)建立以步骤5中建模所需的不同发芽率的稻种样品特征波段荧光强度pi为输入、稻种发芽率gi为输出的稻种发芽率预测模型；
步骤7：对待测稻种按照步骤1至步骤4进行处理，将步骤4中提取的待测稻种特征波段荧光光谱输入到步骤6中建立的稻种发芽率预测模型，得到发芽率。
本发明所述的步骤3中小波去噪的步骤为：
步骤①：先对含噪光谱f(k)进行小波分解，采用‘sym8’小波基构造小波，分解层数为9，实现该步骤的matlab代码为：
xd＝wden(x，′sqtwolog′，′s′，′sln′，9，′sym8′)
其中x为原始信号，sqtwolog为阈值信号处理，s为软阈值，sln为根据第一层小波分解的噪声水平估计进行调整的参数，“9”代表分解为9层，sym8代表采用‘sym8’小波基构造小波；
步骤②：对小波分解得到的噪声部分进行sqtwolog阈值处理。获取域值的matlab代码如下：
[thr，sorh，keepapp]＝ddencmp(′den′，′wv′，nx)
其中x为信号向量，wv表示使用小波分解，den表示去噪声；返回参数thr为阈值，sorh决定硬阈值或软阈值，keepapp为判断是否对近似分量进行阈值处理的参数。
步骤③：根据小波分解的第9层低频系数和各层高频系数对步骤①用于分解含噪光谱的小波进行重构，利用重构得到的小波对f(k)进行去噪处理。
本发明步骤4中所述的pca提取特征波段的过程为：
步骤(1)：将步骤3得到的荧光光谱设为数据矩阵a，计算矩阵a的协方差矩阵b，b＝(aij)p×p，其中，
步骤(2)：求出b的特征值λi及相应的正交变化单位特征向量ai，b的前m个较大的特征值λ1≥λ2≥…λm＞0，为前m个主成分的方差，λi对应的单位特征向量ai就是主成分fi的关于数据矩阵b的系数，则数据矩阵b的第i个主成分fi为
fi＝ai·b
步骤(3)：通过方差(信息)累计贡献率g(m)进行主成分提取，g(m)的计算公式为：
当累计贡献率达到95％时，则提取对应的前m个主成分，否则继续执行步骤(3)。
本发明步骤6中所述的稻种发芽率预测模型，其特征为，基于深度神经网络(dnn)的单输入单输出的稻种发芽率预测模型，dnn建模包括以下步骤：
步骤a：将步骤4所得到特征波段的荧光光谱作为样本数据c，对所述的样本数据c进行归一化，使样本数据范围在[0，1]之间，具体代码如下
fori＝1：m
c(i，：)＝c(i，：)/norm(c(i，：))；
步骤b：定义深度神经网络为三层结构，其中第一层与第三层均是由一个输入单元，一个隐含单元和一个输出单元组成的径向基函数rbf神经网络，第二层是由一个隐含单元和一个可视单元组成的受限玻尔兹曼机rbm神经网络；
步骤c：通过学习样本数据c的特征，训练基于深度神经网络的单输入单输出的稻种发芽率预测模型：
i.随机选取包含稻种老化时间，稻种浸泡时间两个特征的样本数据c的2/3数据作为第一个实验对象，选取样本数据c其余1/3数据作为第二个实验对象，根据权利要求1中步骤4提取两个实验对象训练样本的老化时间和浸泡时间对发芽率贡献的特征，即实验对象训练样本的低级特征；
ii.将实验对象训练样本的低级特征作为深度神经网络中第一层输入单元的输入特征，通过径向基函数rbf神经网络对低级特征进行训练，得到第一个rbf神经网络中隐含单元的输出，作为深度神经网络的第一层输出，即实验对象的高级特征，完成深度神经网络的第一层训练；
iii.通过受限玻尔兹曼机rbm神经网络对实验对象的高级特征进行训练，得到第二个rbm神经网络中隐含单元的输出，作为深度神经网络的第二层输出，即实验对象的更高级特征，完成深度神经网络的第二层训练；
iv.通过径向基函数rbf神经网络对实验对象的更高级特征进行训练，得到第三个rbf神经网络输出单元的输出，作为深度神经网络的第三层输出，即实验对象训练样本的纹理分类特征，完成深度神经网络第三层的训练；
v.若rc≤0.9，返回步骤i；否则结束dnn模型训练，得到训练好的基于荧光光谱的稻种发芽率预测模型。
附图说明
图1是一种基于荧光光谱的稻种发芽率无损检测方法流程图
具体实施方式
下面结合附图及具体实施例对本发明作进一步说明。
实施例1
如图1所示，一种基于荧光光谱的稻种发芽率无损检测方法，其按照以下步实施：
步骤1：将建模所需的不同发芽率稻种样品置入装有蒸馏水的试管中浸泡；
步骤2：取出试管中的浸泡液置于荧光分光光度计的样品池中，扫描浸泡液得到波长范围为360～650nm的荧光光谱；
步骤3：利用小波去噪对步骤2中得到的荧光光谱进行平滑处理；
步骤4：利用主成分分析(pca)提取步骤3平滑处理后的荧光光谱中的特征波段；
步骤5：对建模所需的不同发芽率的稻种样品按照步骤1至步骤4进行操作，得到特征波段荧光光强pi(为第i个样品的特征波段荧光光强，i＝1，2…n，n为建模样本数量)；对建模用的不同发芽率的稻种样品按照gb/t3543.4表1规定的条件(20℃)和时间(5d)进行发芽试验，得到不同稻种样品的发芽率gi(为第i个样品的发芽率，i＝1，2…n，n为建模样本数量)；
步骤6：基于深度神经网络(dnn)建立以步骤5中建模所需的不同发芽率的稻种样品特征波段荧光强度pi为输入、稻种发芽率gi为输出的稻种发芽率预测模型；
步骤7：对待测稻种按照步骤1至步骤4进行处理，将步骤4中提取的待测稻种特征波段荧光光谱输入到步骤6中建立的稻种发芽率预测模型，得到发芽率。