基于近红外光谱分析贡梨可溶性固形物全局模型的鲁棒性
刘燕德, 廖军, 李斌, 姜小刚, 朱明旺, 姚金良, 王秋
华东交通大学机电与车辆工程学院, 智能机电装备创新研究院, 江西 南昌 330013

作者简介: 刘燕德, 女, 1967年生, 华东交通大学智能机电装备创新研究院教授 e-mail: jxliuyd@163.com

摘要

贡梨是大众喜爱的水果, 为研究不同检测方向对近红外在线检测贡梨可溶性固形物SSC的影响, 提出全局模型并分析其鲁棒性。 在贡梨六个方向上收集光谱: 茎-花萼轴垂直, 茎向上(A1)和茎向下(A5), 茎-花萼轴和水平之间45°, 茎向上倾斜(A2)和茎向下倾斜(A4), 茎-花萼轴水平, 茎朝向右侧光(A3), 茎花萼轴水平, 茎朝向带移动方向(A6)。 SSC范围为9.53~14.70的150个样品分为115个标准偏差为1.05的校准集和35个标准偏差为0.93的预测集。 采用偏最小二乘回归PLSR分别建立六个局部模型和一个全局模型, 局部模型由各方向的115个校正集数据经过Savitzky-Golay卷积平滑、 多元散射校正MSC、 高斯滤波平滑GFS三种不同的预处理方法处理后使用偏最小二乘回归PLSR建立而来; 用本方向校正集数据建立的局部模型验证本方向的35个预测集数据, 比较这三种预处理方法后所建立的PLSR模型, 结果表明经过GFS处理后建立的模型验证效果最好, 因此六个局部模型和全局模型均采用GFS处理后建立的PLSR模型。 全局模型是由A1, A2, A3, A4, A5和A6六个方向的690个校正集光谱数据经过GFS预处理后采用PLSR建立的贡梨SSC模型。 各方向的预测集分别对七个模型进行验证, 验证结果表明, 局部模型虽然在本方向的预测效果强于全局模型, 但无法验证其他方向, 鲁棒性差, 由此可知检测方向的不同对预测效果的影响很大; 全局模型能够准确预测各个检测方向的贡梨SSC, 全局模型的校正集相关系数 Rc为0.828, 校正集均方根误差RMSEC为0.424; A1, A2, A3, A4, A5和A6方向的预测集相关系数 Rp分别为0.818, 0.765, 0.799, 0.821, 0.794和0.824, 预测集均方根误差RMSEP分别为0.446, 0.525, 0.478, 0.538, 0.486和0.619; 六个方向的 Rp Rc比较接近且均在0.800左右, RMSEC与RMSEP均在0.500左右, 结果表明全局模型在检测不同方向的贡梨SSC上有着极好的鲁棒性。

关键词: 近红外; 贡梨; 可溶性固形物SSC; 全局模型; 鲁棒性
中图分类号:O657.3 文献标志码:A
Robustness of Global Model of Soluble Solids in Gongli Pear Based on Near-Infrared Spectroscopy
LIU Yan-de, LIAO Jun, LI Bin, JIANG Xiao-gang, ZHU Ming-wang, YAO Jin-liang, WANG Qiu
School of Electromechanical and Vehicle Engineering, East China Jiaotong University, Institute of Intelligent Electromechanical Equipment Innovation, Nanchang 330013, China
Abstract

Gongpear is a popular fruit. In order to study the influence of different detection directions on the online detection of soluble solid SSC in Gongpear by NIR, a global model was proposed, and its robustness was analyzed. The spectra were collected from Gongpears in six directions: stem-calyx axis vertical, stem-upward (A1) and stem-downward (A5), between the stem-calyx axis and horizontal 45°, stem-upward sloping (A2) and stem-downward sloping (A4), stem-calyx axis horizontal, stem-right light oriented (A3), stem-calyx axis horizontal, stem-band moving direction (A6). The 150 samples with SSC ranging from 9.53 to 14. 70 were divided into 115 calibration sets with a standard deviation of 1.05 and 35 prediction sets with a standard deviation of 0.93. Six local models and one global model were established by partial least-squares regression (PLSR). The local models were established by partial least-squares regression (PLSR) after 115 calibration sets of data in each direction were preprocessed by Savitzky-Golay convolution smoothing, Multiple Scattering Correction(MSC)and Gaussian Filtering Smoothing (GFS). The local model established by the local correction set was used to verify the data of 35 prediction sets in the local direction. Compared with the PLSR model established by the three pretreatment methods, the results showed that the model established by GFS processing had the best validation effect. Therefore, the PLSR model established by GFS processing was used for all the six local and global models. The global model is a Gongpears SSC model established by PLSR after GFS pretreatment from 690 calibration sets of spectral data in A1, A2, A3, A4, A5 and A6. The prediction sets in each direction verified the seven models. The verification results showed that although the prediction effect of the local model was stronger than that of the global model in the local direction, it could not be verified in other directions and the robustness was poor. Therefore, different detection directions had a great influence on the prediction effect. The global model can accurately predict the SSC of Gongpears pear in each detection direction. The global model’s correlation coefficient ( Rc)is 0.828, and the root mean square error RMSEC is 0.424. The correlation coefficients ( Rp) of A1, A2, A3, A4, A5 and A6 prediction sets were 0.818, 0.765, 0.799, 0.821, 0.794 and 0.824, and the root mean square errors RMSEP were 0.446, 0.525, 0.478, 0.538, 0.486 and 0.619, respectively. The Rp and Rc in six directions are close to each other and are around 0.800, the RMSEC and RMSEP are around 0.500. The results show that the global model has excellent robustness in detecting gongpear SSC in different directions.

Keyword: Near-infrared; Gongpear; Soluble solid SSC; Global model; Robustness
引言

贡梨是一种受人喜爱的水果, 以皮薄多汁、 味浓甘甜而闻名。 贡梨中的可溶性固形物含量是影响贡梨品质的重要内部属性之一, 也是决定水果成熟度和收获时间的关键参数之一。 因此, 开发无损、 快速的贡梨可溶性固形物含量的检测方法来满足消费者对稳定、 优质水果日益增长的需求具有重要的价值。 近红外光谱技术能够实现食品的无损检测分析, 而且几乎不需要对样品进行处理即可进行检测, 不会产生废物, 并允许同时评估几种成分[1]。 近红外在线检测已被应用于测定梨[2]、 橙子[3]、 橄榄[4]、 苹果[5]等水果的SSC。

水果的物理性质如颜色[6]、 温度[7]、 静止位置[8]、 移动速度[9]等在检测过程中可能会对光谱信息采集有很大影响, 除这些物理性质以外, 测量方向也是影响光谱信息采集的重要因素之一。 Yu等[10]研究了基于近红外光谱的红富士苹果果实方位变化对可溶性固形物在线预测的影响; 韩东海等[11]比较了3种果实方向即茎-花萼轴水平、 茎-花萼轴垂直、 茎向上和茎向下在线检测苹果中的SSC, 基于茎-花萼轴垂直、 茎向上的方向得到了较好的结果; 在透射模式下, 基于茎萼轴水平方向检测梨褐心得到了较好的结果[12], 茎萼轴垂直取向也获得了测量苹果的SSC和硬度的良好结果。 如今近红外在线检测水果内部品质的设备多采用局部模型对水果进行预测, 由于局部模型是利用单个局部方向获得的光谱数据单独开发的, 在预测其他水果方向上的SSC时, 其精度会大打折扣。 目前解决该问题主要依靠人工摆正果位或增加摆正机构, 费钱费力。 建立一个基于全方向的全局模型也许能够解决上面的问题, 也还没有人将全局模型应用于近红外在线漫透射检测装备来预测不同方向运送下的贡梨SSC。

本研究利用课题组自主研发的在线漫透射检测装备采集六个不同检测方向下的贡梨光谱, 使用Savitzky-Golay卷积平滑、 多元散射校正MSC、 高斯滤波平滑GFS对平均后的150贡梨光谱进行预处理, 对光谱进行修正和优化, 最后通过偏最小二乘回归PLSR建立了六个局部模型和一个全局模型。 六个方向的预测集分别对七个模型进行预测, 对预测结果进行比较并评判各局部模型与全局模型的鲁棒性。 为全局模型是否能够应用于近红外在线漫透射检测装备提供参考。

1 实验部分
1.1 样本

于2021年四月下旬在中国南昌某水果超市购买150个贡梨, 150个贡梨均没有表面缺陷或污染物。 贡梨的赤道直径范围为85~95 mm, 所有样品都单独清洗并编号; 实验前把样品放在实验室储存24 h, 使贡梨温度与实验室温度一致, 以减少温度对预测精度的影响。

将150个样本随机分配到校准集和预测集。 选择115个样本作为用于开发校准模型的校准集, 其余35个样本用作预测集来验证校准模型的预测性能。

1.2 光谱收集

样品的漫透射光谱由本课题组自主研发的水果内部品质近红外漫透射光谱在线检测装备(图1)测得[13]

图1 水果内部品质近红外漫透射光谱在线检测装备示意图
1: 触发开关; 2: 编码盘; 3: 链轮; 4: 链条; 5: 光纤; 6: 暗箱; 7: 光源; 8: 样品; 9: 终端PC机; 10: PLC控制柜; 11: 果杯; 12: 弹跳器; 13: 分级出口; 14: 电动机; 15: 减速器; 16: 驱动链轮
Fig.1 Schematic diagram of on-line detection equipment for fruit internal quality by near infrared diffuse transmission spectroscopy
1: Porximity switch; 2: Coding mask; 3: Chain wheel; 4: Link chain; 5: Optica fiber; 6: Camera bellows; 7: Optical source; 8: Specimen; 9: Terminal PC; 10: PLC control cabinet; 11: Fruit cup; 12: Bionic boot; 13: Grading export; 14: Electromotor; 15: Reducer; 16: Drive sprocket

使用反射率超过95%的聚四氟乙烯制成的直径ϕ =90 mm的白色球体测量参考光谱。 参考测量三次, 计算平均光谱作为最终的Rwhite。 当所有的灯都关了, 灯的入口打开时, 收集暗光谱Rdark

如图2所示[14], 茎-花萼轴垂直, 茎向上(A1)和茎向下(A5), 茎-花萼轴和水平之间45° , 茎向上倾斜(A2)和茎向下倾斜(A4), 茎-花萼轴水平, 茎朝向右侧光(A3), 茎花萼轴水平, 茎朝向带移动方向(A6)。 实验前随机选取了60个贡梨, 让它们从上一传送通道上随机传过来, 贡梨在果杯上的位置都在这六个方向附近, 这可能是因为贡梨本身的结构造成的。 且两排光源以果杯为基准为对称光源, 这六个方向相当于一个大区间, 从上一传送通道上传过来的贡梨在果杯上的位置都在这一区间内。

图2 漫透射检测机构(a)及贡梨摆放位置(b)Fig.2 Diffuse transmission testing mechanism (a), Gongli placement position (b)

分别在六个方向测量光谱数据, 光谱数据按A1, A2, A3, A4, A5和A6的顺序采集。 用手将每个样品稳定地放在水果盘上。 本课题组自主研发的水果内部品质近红外漫透射光谱在线检测装备(图1)能够保证运行平稳, 果杯上的贡梨震感极弱, 在传送贡梨的过程中可以保证贡梨的平稳传送。 仅保留346~1 125 nm的光谱数据, 由1 044个数据点即光谱变量组成, 用于进一步分析。 然后, 原始光谱Rraw根据式(1)转换为校准光谱R:

R=Rraw-RdarkRwhite-Rdark(1)

考虑到贡梨形状的不规则性, 为了减少操作人员和仪器引入的误差, 提高信噪比, 在每个水果方位采集三个光谱, 样品绕茎萼轴旋转0° , 120° 和240° , 对三个光谱进行平均得到平均光谱。 样品以30 cm· s-1的速度由传送带传送, 同时所有光谱数据存储在工业计算机中。

在光谱采集之后, 立即通过传统的破坏性试验确定SSC值。 手动挤压每个完整样品的汁液, 并将其滴到日本东京ATAGO有限公司生产的数字折光仪PAL-1上(测量精度为Brix± 0.2%, 且具有温度补偿功能, 适用温度范围10~40 ℃), 以测量SSC值, 绕茎萼轴旋转0° , 120° 和240° , 为每个样品收集三个SSC值, 并根据这三个值来获取平均值。

1.3 光谱数据预处理

由于采集的光谱可能包含无用的背景信息或噪声, 而光谱预处理往往可以提高模型的性能。 采用不同的预处理方法对原始光谱数据进行预处理, 包括Savitzky-Golay卷积平滑SGS、 多元散射校正MSC和高斯滤波平滑GFS等来优化原始光谱数据。 SGS可以有效减少光谱的噪声并提高信噪比, 多元散射校正MSC可以消除水果内部不均导致的散射效应; 高斯滤波GFS可以消除光谱中的随机噪声。

1.4 偏最小二乘回归分析

偏最小二乘回归PLSR是一种无监督的多变量统计分析方法[15], 被广泛应用于NIRS分析研究, 具有较强的抗干扰性。 PLSR的目的是通过将光谱矩阵与浓度向量投影到新的潜在空间, 使其中的潜在变量之间的协方差最大化, 并最终执行潜变量线性回归模型。

1.5 模型的建立

为了全面检验测量方向对贡梨SSC光谱分析的影响, 使用PLSR分别为每个方向和所有方向数据集建立SSC模型, 即局部模型与全局模型。 每个局部模型都是由该方向上的115个校正集光谱数据建立的。 因此, 分别建立了六个局部模型, 如图3, 校正集用长方形表示, 预测集用圆形表示; 在局部模型中, 实线所指为本方向的预测集对该模型进行验证; 虚线为其他方向的预测集对该模型进行验证, 这是为了评估局部模型的鲁棒性。 一些物理性质的全局模型已经成功地用于减少色差[8]、 温度波动[9]和生物可变性对近红外模型性能的影响。 在本研究中, 全局模型中的校正集由六个方向的校正集数据组成, 总共包括690个光谱。 为了验证全局模型的预测能力, 分别用六个预测集对全局模型进行验证。 比较七个模型的性能, 选择出鲁棒性最好的模型。

图3 局部模型、 全局模型的建立以及实验验证方案Fig.3 Establishment of local model and global model and experimental verification scheme

1.6 模型性能分析

根据校准和预测的相关系数、 校准的均方根误差RMSEC、 预测的均方根误差RMSEP来评估校准模型的性能。 Rc, Rp, RMSEC和RMSEP的计算定义如式(2)— 式(5)

Rc=i=1nc(ypi-ymi)2/i=1nc(ypi-ymeanc)2(2)

Rp=i=1np(ypi-ymi)2/i=1np(ypi-ymeanp)2(3)

RMSEC=1nci=1nc(ypi-ymi)2(4)

RMSEP=1npi=1np(ypi-ymi)2(5)

其中, ypi为第i个水果的SSC的预测值, ymi为第i个贡梨SSC的实测值, ymeancymeanp分别为校准集和预测集中SSC的平均值, ncnp分别为校准集和预测集中水果个数。 一般来说, 选择一个好的模型, 不仅要根据较高的RcRp, 较低的RMSEC和RMSEP, 还要根据RMSEC和RMSEP之间的微小差异来选择。

2 结果与讨论
2.1 光谱特征与光谱分析

在150个贡梨中随机选取第11号果, 如图4所示为11号果在A1— A6六个方向上获得的346~1 125 nm的漫透射平均原始光谱。 由图4可知, A1方向的漫透射光谱的能量明显强于其他方向, 处于其他方向的漫透射光谱上方; 其他方向的漫透射光谱的能量由强到弱依次是A5, A6, A2, A3和A4, 其中A6, A2, A3和A4的漫透射光谱的能量相差不大; 各方向的漫透射光谱的能量不同可能与贡梨所受的有效光源照射面积有关, 其中A1与A5所受的有效光源照射面积最大。 从六个方向所采集的光谱的吸收峰和波谷所处的波段大致相同, 只是光强大小不同。 其中波谷在750 nm附近, 在710与800 nm附近各有一个波峰。 750 nm附近的波谷可能由O— H三倍频伸缩振动导致。

图4 11号果六个方向的平均光谱图Fig.4 Average spectra of fruit No.11 in six directions

2.2 可溶性固形物含量分析

表1为六个方向的校正集和预测集可溶性固形物的范围, 标准偏差和平均值。 从表1可看出含115个样本的校正集的可溶性固形物含量的范围大于含35个样本的预测集的可溶性固形物含量的范围。

表1 校正集和预测集中可溶性固形物SSC 含量的范围、 标准偏差和平均值 Table 1 Range, standard deviation and average value of SSC content in calibration set and prediction set
2.3 局部模型建立及验证

以A1方向为例, 用在A1方向采集的115个校正集光谱数据建立模型, 再用另外35个预测集数据对其进行局部验证。 其他五个方向同理建立模型并进行验证。 每种预处理方法侧重的方面都有所不同, 因此对原始光谱用不同的预处理方法会有不同的效果, 为了有效减少光谱的噪声并提高信噪比, 使用SGS对六个局部模型进行预处理; 考虑到水果内部不均易导致散射效应, 采用MSC对六个局部模型进行预处理; 由于光谱中的随机噪声, 选用GFS对六个局部模型进行预处理。 表2记录了使用SGS, MSC和GFS预处理后结合PLSR建立的模型以及预测结果。 可以发现, 六个方向的原始数据使用SGS和MSC预处理后建立的PLSR模型, 它们的RcRp都会有所减小, 校正集和预测集的RMSE大多都会变大; 当使用GFS分别对六个方向原始数据预处理后, A1, A2, A4, A5和A6的Rc会变小, Rp会变大, 这就使得校正集和预测集的相关系数更接近, 拟合性更好。 因此, 结果表明通过GFS预处理后建立的六个局部模型比用其他两种预处理后建立的结果更好。

表2 六个方向的局部模型与局部预测效果 Table 2 Local model and local prediction effect in six directions

将经过GFS预处理后的光谱建立的六个局部模型及验证结果进行比较, 可以发现在A3方向的预测效果最好, 它的RcRp分别是0.913和0.872, RMSEC和RMSEP分别为0.397和0.342。 A6方向的局部模型的局部预测效果最差, 它的RcRp分别是0.831和0.794, RMSEC和RMSEP分别为0.534和0.508。 结果表明, 贡梨SSC局部模型在相同方向上的预测效果非常好。 但是在线检测时, 如果没有机构或人工摆正水果的位置, 水果在经过检测装置时不可能只是单一方向的。 因此, 要想证明局部模型该情况下的效果, 还需要看其他五个方向的预测集对其进行验证的效果。

2.4 全集模型的建立及验证

使用A1— A6六个方向上的690个校正集数据建立全局模型, 各方向的预测集分别对全局模型进行验证。 然而从表3表4中可以看到, 经过GFS预处理后建立的PLSR模型, 校正集与各个方向的预测集的拟合性变好; 校正集的Rc为0.828, A1, A2, A3, A4, A5和A6的Rp分别为0.818, 0.765, 0.799, 0.821, 0.794和0.824, 校正集的RMSEC为0.424, RMSEP分别为0.446, 0.525, 0.478, 0.538, 0.486和0.619, 可以看出全局模型对各个方向都有较好的预测效果。

为了进一步评估全局模型对全方向模型的预测效果, 把六个方向的预测集合成一个预测集, 再对全集模型进行验证, 该预测集的Rp和RMSEP分别为0.748和0.548, 结果表明依旧有较好的预测效果。 图3为贡梨SSC测量真值与预测值的散点图。

表3 全局模型以及A1, A2和A3方向预测效果 Table 3 Global model and prediction effect of A3, A4 and A5 directions
表4 全局模型以及A4, A5和A6方向预测效果 Table 4 Global model and prediction effect of A4, A5 and A6 directions

图5 全方向验证散点图Fig.5 Scatter plot of omnidirectional verification

2.5 局部模型与全集模型的对比

用六个方向的预测集数据分别对A1, A2, A3, A4, A5和A6的局部模型和全局模型进行验证。 如表5表6, 局部模型虽然对本方向的预测效果很好, 但是对于其他方向的预测效果却很差, 例如局部模型A1预测A2, A3, A4和A5的RMSEP分别为1.431, 0.816, 0.767和1.015; 局部模型A2预测A1, A4, A5和A6的Rp均在0.1以下, 不可信赖; 其他局部模型均有类似情况。 在线检测时如果没有人工对水果进行方向调整或者没有设计制造调整方向的摆正机构, 那么局部模型是不具备在线检测能力的, 局部模型在验证其他方向时所表现出来的鲁棒性太差。 六个方向的预测集对全局模型进行验证的效果虽不如局部模型预测本方向时的验证效果, 但是效果都还不错, 全局模型的校正集Rc为0.828, RMSEC为0.424; A1, A2, A3, A4, A5和A6方向的预测集分别为0.818, 0.765, 0.799, 0.821, 0.794和0.824, RMSEP分别为0.446, 0.525, 0.478, 0.538, 0.486和0.619; 结果表明六个方向的RpRc均比较接近且均在0.800左右, RMSEC与RMSEP均在0.500左右, 因此全局模型在验证六个方向时是鲁棒的。 但是总体来看R都只有0.8左右, 这可能是因为贡梨六个方向上的SSC含量分布不均匀, 而本实验所测得的SSC含量真值为贡梨整果的平均SSC含量真值。

表5 局部模型和全局模型对A1, A2, A3方向预测集的预测效果 Table 5 Prediction effects of local model and global model on A1, A2 and A3 direction validation sets
表6 局部模型和全局模型对A4, A5和A6预测集的预测效果 Table 6 Prediction effects of local model and global model on A4, A5 and A6 direction validation sets
3 结论

采用本课题组自主研发的水果内部品质近红外漫透射光谱在线检测装备采集150个贡梨六个方向的光谱, 研究了光谱测量方向的变化对贡梨SSC在线测定的影响。 使用GFS预处理后通过PLSR建立的模型效果都比较好, 局部模型能够很好的预测本方向的贡梨SSC, 全局模型可以较好的预测全方向的贡梨SSC。 通过建立局部模型和全局模型并分别用六个方向的预测集对它们进行验证, 结果表明样品的检测方向极大地影响了贡梨SSC的预测精度, 局部模型的鲁棒性差; 全局模型与局部模型相比, 基于所有方向的全局模型获得了更准确的预测结果, 六个方向的RpRc均比较接近且均在0.800左右, RMSEC与RMSEP均在0.500左右, 鲁棒性更好。 因此, 全局模型将有助于消除光谱测量中测量方向的影响, 实现贡梨SSC的整体方向在线检测, 这对近红外在线检测工业化有实际意义。

参考文献
[1] Gabrieli Alves de Oliveira, Sylvie Bureau, Catherine Marie-Geneviève Claire Renard, et al. Food Chemistry, 2014, 143: 223. [本文引用:1]
[2] QIN Shan-zhi, CHEN Bin, LU Dao-li, et al(秦善知, 陈斌, 陆道礼, ). Jiangsu Agricultural Sciences(江苏农业科学), 2014, 42(8): 284. [本文引用:1]
[3] JIANG Shui-quan, SUN Tong(江水泉, 孙通). Food and Machinery(食品与机械), 2020, 36(2): 89. [本文引用:1]
[4] Lourdes Salguero-Chaparro, Vincent Baeten, Ouissam Abbas, et al. Journal of Food Engineering, 2012, 112(3): 152. [本文引用:1]
[5] MENG Qing-long, SHANG Jing, HUANG Ren-shuai, et al(孟庆龙, 尚静, 黄人帅, ). Food and Fermentation Industry(食品与发酵工业), 2020, 46(19): 205. [本文引用:1]
[6] Guo Zhiming, Huang Wenqian, Peng Yankun, et al. Postharvest Biology and Technology, 2016, 115: 81. [本文引用:1]
[7] MO Bao-qing, ZHU Hui-juan, XIAO Li(莫宝庆, 朱慧娟, 肖黎). Chinese Journal of Health Inspection(中国卫生检验杂志), 2012, 22(1): 125. [本文引用:1]
[8] Fan Shuxiang, Zhang Baohua, Li Jiangbo, et al. Biosystems Engineering, 2016, 143: 9. [本文引用:2]
[9] Xia Yu, Huang Wenqian, Fan Shuxiang, et al. Journal of Food Process Engineering, 2018, 41(8): e12915. [本文引用:2]
[10] Xia Yu, Huang Wenqian, Fan Shuxiang, et al. Infrared Physics and Technology, 2019, 97: 467. [本文引用:1]
[11] HAN Dong-hai, ZHOU En-yang, QI Shu-ye(韩东海, 周恩洋, 戚淑叶). Transactions of the Chinese Society for Agricultural Machinery(农业机械学报), 2015, 46(11): 260. [本文引用:1]
[12] SUN Xu-dong, LIU Yan-de, LI Yi-fan, et al(孙旭东, 刘燕德, 李轶凡, ). Transactions of the Chinese Society for Agricultural Machinery(农业机械学报), 2016, 47(1): 227. [本文引用:1]
[13] LI Xiong, LIU Yan-de, SUN Xu-dong, et al(李雄, 刘燕德, 孙旭东, ). Spectroscopy and Spectral Analysis(光谱学与光谱分析), 2020, 40(10): 3265. [本文引用:1]
[14] Liu Yand e, Zhang Yu, Jiang Xiaogang, et al. Vibrational Spectroscopy, 2020, 111: 103152. [本文引用:1]
[15] HUANG Jian-jun, LI Xue-mei, TENG Hong-quan(黄建军, 李雪梅, 滕宏泉). Disaster Science(灾害学), 2021, 36(2): 60. [本文引用:1]