基于激光超声融蚀机制的粮食品质检测系统及方法

1.基于激光超声融蚀机制的粮食品质检测系统，其特征在于：包括激光器、分光镜、衰减片、透镜、光电二极管、放大器、计算机和超声探测器，该超声探测器的输出端连入放大器的输入端，放大器的输出连入计算机的输入接口，所述分光镜、衰减片和透镜顺次设于激光器的光轴上，所述分光镜水平设置，激光器的发射光束平行于分光镜入射角的入射光线，所述衰减片和透镜垂直于激光器的光轴设置，所述光电二极管置于分光镜反射光束的光轴上，光电二极管与计算机有线连接；所述超声探测器为带宽电容式超声探测器。
2.根据权利要求1所述的基于激光超声融蚀机制的粮食品质检测系统，其特征在于：
所述分光镜的入射角为45°，激光器的发射光束与分光镜呈45°角。
3.利用权利要求1所述的检测系统进行粮食品质检测的方法，其特征在于，包括以下步骤：
（1）调节检测系统的透镜与粮食籽粒表面的距离，使粮食籽粒表面处于透镜焦平面上，激光器发出的激光脉冲通过透镜汇聚于粮食籽粒表面，在粮食籽粒表层产生热应力区，从而在粮食籽粒内部产生应力波即超声波；
（2）超声探测器接收超声波信号，该超声波信号经采样放大后，存储于计算机中；
（3）取M粒同一种类相同品质粮食籽粒，重复步骤（2），收集M粒粮食籽粒的超声波信号，其中M≥300；（4）计算机对接收到的超声波信号进行时域分析，提取与粮食品质有关的时域特征参数，该时域特征参数包括幅值总和、幅度方差、幅度均方根值、波形指标、峰值因子、脉冲因子和过零率参数；
（5）使用小波包变换技术提取超声波信号的频域特征参数；
（6）取同一种类粮食的N个不同品质样本，其中N≥10，重复步骤(3)~(5)，得到同一种类不同品质粮食的时、频特征参数；
（7）建立粮食品质参数模型，以得到的时域特征参数和频域特征参数作为自变量，品质指标为因变量，反演出粮食的品质参数；
（8）在得到粮食的品质参数后，根据国家相关分类标准，对粮食品质进行分级。
4.根据权利要求3所述的方法，其特征在于：所述步骤（1）中的激光器通过调节激光能量以增加激发效率，使粮食籽粒表层发生融蚀效应，以激励出更强的超声纵波。
5.根据权利要求3所述的方法，其特征在于：所述步骤（5）是使用Daubechies4（db4）小波基函数对超声波信号进行小波包分解，其消失矩为4，支撑区间宽度为7，通过对超声探测器所采集到的信号进行3层小波包分解，得到8个小波包重构信号分量即8个频域特征参数。
6.根据权利要求3所述的方法，其特征在于：所述步骤（7）中建立粮食品质参数的反演模型用到的特征参数，是通过计算时域特征参数与品质指标间的相关性，频域特征参数与品质指标间的相关性，挑选相关系数大于0.8的特征参数，以这些特征参数作为自变量，建立粮食品质参数的模型。
7.根据权利要求6所述的方法，其特征在于：所述步骤（7）中以挑选的特征参数为自变量，品质指标为因变量，建立多元回归方程：Y＝a0+a1X1+a2X2+…+aLXL,其中L为挑选的特征参数的个数，Xi(i＝1,2,…,L,L≤15)为挑选的特征参数，Y为品质指标，通过对N个品质样本的测量，得到N个数据集为(X1i,X2i，…,XLi;Yi)，i＝1,2,…,N,则可以写成矩阵形式Y＝XA，即为：
T -1 T
使用最小二乘法求解系数矩阵：A＝(XX) XY，由此可以通过超声波时频信号的特征参数反演出粮食品质的指标参数。

基于激光超声融蚀机制的粮食品质检测系统及方法\n技术领域\n[0001] 本发明涉及粮食品质在线检测系统及方法，特别是一种基于激光超声融蚀机制的粮食品质检测系统及方法，属于粮食品质检测技术领域。\n背景技术\n[0002] 粮食的生产、贸易及利用是以粮食的品质为基础的。在粮食生产过程中，粮食品质是育种和指导种植的重要目标；在收购环节，粮食品质是定等作价的依据；在粮食储藏过程中，随着储藏时间的延长，粮食的品质会发生很大变化，甚至发生劣变，粮食品质是指导国家储备粮合理轮换的最重要指标。必须不断监测粮食品质的变化，避免给国家造成重大损失，危及国家粮食安全；在粮食加工及贸易过程中，粮食品质决定粮食的最佳用途与市场价值，是有关食品工业质量控制的重要基础；“民以食为天”，粮食品质是人们身体健康的基本保证。任何粮食的数量都是以粮食的基本品质为前提的，撇开粮食的品质只谈粮食的数量是没有意义的。确保粮食品质优质、安全是贯穿我国粮食工作始终的生命线。客观、快速、准确地测定粮食的品质是维护我国粮食贸易基本利益、加强流通、合理利用粮食资源、引导种植结构调整，提高农民收入的技术保证，是粮食工作最重要的内容之一。\n[0003] 在粮食检验和品质控制中，目前有许多常规方法用于检测粮食的物理特性、化学成分和功能特征。根据工作原理，现行的这些分析方法大致可以分为计算机视觉检测法、化学分析方法、近红外光谱分析法和声学法（可听声频段）等。上述几种方法在粮食品质检测方面都获得了广泛应用并取得较好的结果，到目前为止，这些方法仍然是其它方法不可替代，今后仍将被广泛应用的。然而，这些常规方法均存在一些难以克服的缺陷。例如，基于计算机视觉的粮食品质检测方法主要对粮食的外观品质参数检测较为适用，并且这种方法属于静态检测，限制了检测的实时性；化学分析方法不同程度地存在着检测时间长、成本高、污染环境和影响工作人员健康的问题；近红外光谱分析法的模型维护较为困难，建立的模型并不能一劳永逸，并且品质参数测量的精度很大程度上依赖于已知样品的化学分析准确度；而声学方法（可听声频段）对环境噪声屏蔽的要求较为苛刻。随着社会对粮食安全和生态环境的关注，人们越来越欢迎准确、快速、无损、成本低，不需要化学药剂，不污染环境、不会对工作人员健康造成伤害、环境适应能力强，并且能够在线检测的新型检验技术。粮食品质检测，特别是准确、快速、在线检测是目前我国粮食生产、收购、储藏、加工、流通和消费过程中最突出的技术“瓶颈”之一。\n发明内容\n[0004] 本发明的目的是提供一种基于激光超声融蚀机制的粮食品质检测系统及方法，以解决现有粮食检测存在的检测时间长、成本高、污染环境、影响工作人员身体健康、检测环境苛刻的问题。\n[0005] 为实现上述目的，本发明的基于激光超声融蚀机制的粮食品质检测系统，包括激光器、分光镜、衰减片、透镜、光电二极管、放大器、计算机和超声探测器，该超声探测器的输出端连入放大器的输入端，放大器的输出连入计算机的输入接口，所述分光镜、衰减片和透镜顺次设于激光器的光轴上，所述分光镜水平设置，激光器的发射光束平行于分光镜入射角的入射光线，所述衰减片和透镜垂直于激光器的光轴设置，所述光电二极管置于分光镜反射光束的光轴上，光电二极管与计算机有线连接。 \n[0006] 进一步的，所述超声探测器为带宽电容式超声探测器。\n[0007] 进一步的，所述分光镜的入射角为45°，激光器的发射光束与分光镜呈45°角。\n[0008] 利用上述检测系统进行粮食品质检测的方法包括以下步骤：\n[0009] （1）调节检测系统的透镜与粮食籽粒表面的距离，使粮食籽粒表面处于透镜焦平面上，激光器发出的激光脉冲通过透镜汇聚于粮食籽粒表面，在粮食籽粒表层产生热应力区，从而在粮食籽粒内部产生应力波即超声波；\n[0010] （2）超声探测器接收超声波信号，该超声波信号经采样放大后，存储于计算机中；\n[0011] （3）取M粒同一种类相同品质粮食籽粒，重复步骤（2），收集M粒粮食籽粒的超声波信号，其中 ；\n[0012] （4）计算机对接收到的超声波信号进行时域分析，提取与粮食品质有关的时域特征参数；\n[0013] （5）使用小波包变换技术提取超声波信号的频域特征参数；\n[0014] （6）取同一种类粮食的N(N≥10)个不同品质样本，重复步骤(3)~(5)，得到同一种类不同品质粮食的时、频特征参数；\n[0015] （7）建立粮食品质参数模型，以得到的时域特征参数和频域特征参数作为自变量，品质指标为因变量，反演出粮食的品质参数；\n[0016] （8）在得到粮食的品质参数后，根据国家相关分类标准，对粮食品质进行分级。\n[0017] 进一步的，所述步骤（1）中调节激光能量以增加激发效率，使粮食籽粒表层发生融蚀效应，以激励出较强的超声纵波。\n[0018] 进一步的，所述步骤（4）中时域特征参数包括：幅值总和、幅度方差、幅度均方根值、波形指标、峰值因子、脉冲因子、过零率参数。\n[0019] 进一步的，所述步骤（5）是使用Daubechies4（db4）小波基函数对超声波信号进行小波包分解，其消失矩为4，支撑区间宽度为7，通过对超声探测器所采集到的信号进行3层小波包分解，得到8个小波包重构信号分量即8个频域特征参数。\n[0020] 进一步的，所述步骤（7）中建立粮食品质参数的反演模型用到的特征参数，是通过计算时域特征参数与品质指标间的相关性，频域特征参数与品质指标间的相关性，挑选相关系数大于0.8的特征参数，以这些特征参数作为自变量，建立粮食品质参数的模型。\n[0021] 进一步的，所述步骤（7）中以挑选的特征参数为自变量，品质指标为因变量，建立多元回归方程： ,其中L为挑选的特征参数的个数，\n为挑选的特征参数，Y为品质指标，通过对N个品质样本的测量，得\n到N个数据集为 ， ,则可以写成矩阵形式 ，即为：\n[0022] ，\n[0023] 使用最小二乘法求解系数矩阵： ，由此可以通过超声波时频信\n号的特征参数反演出粮食品质的指标参数。\n[0024] 本发明的基于激光超声融蚀机制的粮食品质检测系统及方法具有如下优点：\n[0025] （1）由于使用脉冲激光激励超声信号，使用宽带电容超声换能器接收超声纵波信号，因此本检测系统和方法具有非接触的特点；\n[0026] （2）由于以超声纵波信号作为超声探测器的接收信号，因此该系统和方法对环境噪声的抗干扰能力较强，易于实现实时在线检测；\n[0027] （3）由于采用了激光超声的融蚀机制激发超声波，因此增加了超声波的激发效率，改善了接收信号的信噪比，可有效提高测量精度，且粮食籽粒的表层融蚀深度仅几个微米，接近于无损检测；\n[0028] （4）由于该系统和方法在检测过程中，无需使用任何化学制剂、无放射性物质产生，数据的采集与处理由计算机自动完成，因此该系统和方法具有无污染、无需人工干预、检测效率高等优点。\n[0029] 综上所述，本发明由于以脉冲激光作为激励源，以超声波作为检测信号，因此该方法具有快速、准无损、非接触、抗干扰能力强、效率高、无污染等许多优点，为粮食品质检测提供了一种新型的实时在线检测系统和方法。\n附图说明\n[0030] 图1 基于激光超声融蚀机制的粮食品质检测系统原理图；\n[0031] 图2 基于激光超声融蚀机制的粮食品质检测方法流程图;\n[0032] 图3 超声信号的三层小波包分解二叉树。\n具体实施方式\n[0033] 如图1所示，本发明的基于激光超声融蚀机制的粮食品质检测系统包括激光器1、分光镜2、衰减片3、透镜4、光电二极管5、放大器6、计算机8和超声探测器7，该超声探测器7的输出端连入放大器7的输入端，放大器7的输出连入计算机7的输入端，放大器7、光电二极管5分别与计算机8有线连接，分光镜2、衰减片3和透镜4顺次设于激光器1的光轴上，分光镜2水平设置，激光器1发射光束的光轴与分光镜2呈45°的入射夹角，衰减片\n3和透镜4垂直于激光器1的光轴，光电二极管5置于分光镜2反射光束的光轴上，光电二极管5与计算机8有线连接，用于产生计算机8的同步接收信号，超声探测器6为带宽电容式超声探测器。主要器件的性能描述如下表：\n[0034] \n[0035] 本发明的基于激光超声融蚀机制的粮食品质检测方法，是以粮食籽粒激光超声融蚀机制的机理研究为出发点；以脉冲激光器为激光源、电容超声换能器为接收器件建立检测系统，获取粮食籽粒的超声纵波信号；利用小波包变换技术提取信号的时频特征参数；\n进而基于最小二乘法建立粮食品质参数的反演模型；以国标为依据，利用品质参数最终实现粮食品质的分级与评价。\n[0036] 以检测小麦籽粒为例，流程图如图2所示，具体实施方案如下：\n[0037] （1）调节检测系统的透镜与小麦籽粒表面的距离，使小麦籽粒表面处于透镜焦平面上，使激光器发出的激光脉冲通过透镜汇聚于小麦籽粒表面，利用高能激光脉冲与小麦籽粒表层的瞬时热作用，在小麦籽粒表层产生热应力区，从而在小麦籽粒内部产生应力波即超声波；脉冲激光作用于小麦籽粒表面，激发超声波位移场，该过程可划分成两个相互耦合的子过程：一是脉冲激光在小麦籽粒内部产生温度场，二是由温度场引起的超声位移场：\n[0038] ， \n[0039] ， \n[0040] 这里， 是小麦籽粒内部t时刻温度的瞬态空间分布， 表示t 时\n刻引起的超声波位移场， 、 、k分别表示小麦籽粒的密度、定容比热容和热传导率，为热弹耦合系数，其中 为线膨胀系数， 和 为粮食的Lame常数。\n[0041] 在极坐标系中，激光辐照区域的热流边界条件可表示为:\n[0042] ，\n[0043] 是样品表面光吸收率，这是一个随温度而改变的量， 为入射激光的功率密度。 和 分别为激光脉冲的空间和时间分布， 和 由脉冲激光器生产厂\n家提供，可表示为:\n[0044] ，\n[0045] ，\n[0046] 是脉冲激光线源的辐照半宽度， 为激光脉冲的上升时间，采用有限元方法得到超声波场 的数值模拟解。小麦激光超声的有限元数值模拟分析可划分成两个过程：一是脉冲激光在小麦内部产生的温度场 的数值分析，二是由温度场引起的超声位移场 的数值分析。\n[0047] （2）带宽电容式超声探测器距离小麦籽粒表面1-2mm左右以接收超声波信号，该超声波信号经采样放大后，存储于计算机中；超声探测器将接收的超声波信号转换成与强度呈正比关系的电信号 ，即： 。\n[0048] （3）对于同一种类相同品质的粮食，例如豫麦49，使用同样的处理对M粒小麦籽粒的超声波信号，得到一系列电信号 ，一般M取300~500个籽粒。\n[0049] （4）计算机采集到这些信号后，提取其时域特征参数，主要包括：幅值总和、幅度方差、幅度均方根值、波形指标、峰值因子、脉冲因子、过零率。各参数的定义为：\n[0050] ① 幅度总和： ，\n[0051] ② 幅度方差： ，\n[0052] ③ 幅度均方根值： ，\n[0053] ④ 波形指标： ，\n[0054] ⑤ 峰值因子： ，\n[0055] ⑥ 脉冲因子： ，\n[0056] ⑦ 过零率： 。\n[0057] （5）对采集到的超声电信号 进行小波包变换，这里采用的小波基函数是Daubechies4（db4），其消失矩为4，支撑区间宽度为7，通过对电容超声探测器所采集到的信号进行3层小波包分解，得到8个小波包重构信号分量，如图3所示。经三层小波包分解重构后各频段的能量: 。\n[0058] 其中， 是重构信号 的各点的幅值。\n[0059] （6）使用步骤（4）、（5）对同一种类不同品质粮食的N个样本(例如同属于小麦种类的豫麦34、豫麦49、豫麦50、内乡188等，N≥10，N值越大，拟合效果越好)分别提取特征参数 、 。对每一样本，首\n先以特征参数 为自变量，品质指标 为因变量（通过标准方法测量得到），对于N个样本点， ，…， ，则一元线性回归方程可表示为：\n[0060] ；\n[0061] 其中， ，\n[0062] ，\n[0063] ，\n[0064] 。\n[0065] （7）计算特征参数 和品质指标 的相关性，计算公式为：\n[0066] ，\n[0067] ，当 时 和 完全相关，当 时则 和 之间没有相关性， 越\n接近1， 和 的相关性越好。\n[0068] 对于其他特征参数（ ），重复使用步骤（6）、（7）分别计算一元线性回归方程的系数 以及相应特征参数和品质指标之间的相关性。\n[0069] （8）通过前面的计算，比较特征参数与品质指标间的相关系数，挑选相关系数大于\n0.8的特征参数，假设为 ,以这些特征参数为自变量，品质指标\n为因变量，建立多元回归方程： ，通过对N个样本的测量，\n得到N个数据集为 ， ,则可以写成矩阵形式 ，即\n为：\n[0070] ，\n[0071] 使用最小二乘法求解系数矩阵： ，由此可以通过超声波时频信\n号的特征参数 ,反演出粮食品质的指标参数 。\n[0072] （9）在得到粮食品质的指标参数后，根据国家相关分类标准，完成粮食品质的分级。