Prediction and study of strawberry hardness based on hyperspectral index segmentation
-
摘要: 为了快速无损检测草莓硬度, 连续5天采集了草莓高光谱数据和硬度信息, 提出一种基于高光谱多指数阈值逐层分割的硬度预测方法。首先分析各组分(果肉、霉变果肉、草莓籽和萼片)的光谱反射率差异并确定特征波段, 利用特征波段构建新的归一化特征指数, 完成分割阈值的确定, 采用逐层分割的方法以排除无关部分的干扰; 通过连续投影算法、主成分分析法及2次组合降维来降低光谱信息冗余度并提取特征, 利用随机森林与偏最小二乘法分别对原始光谱及降维后特征建立回归模型, 并确立最佳预测模型; 最后利用最佳预测模型对草莓果肉部分进行硬度拟合, 得到硬度分布图像, 实现了草莓硬度预测结果的直观显示。结果表明, 基于2次降维建立的偏最小二乘模型效果最好, 测试集和预测集的相关系数分别为0.9101和0.9099, 测试集均方根误差为0.1344。该研究为草莓硬度的无损检测和显示提供了参考。Abstract: To achieve rapid and non-destructive detection of strawberry hardness, strawberry hyperspectral data, and hardness information were collected for five consecutive days, and a hardness prediction method based on high spectral multi-index threshold layer-by-layer segmentation was proposed. Firstly, the spectral reflectance differences of different components (pulp, moldy pulp, strawberry seeds, and sepals) were analyzed, and the characteristic bands were identified. Subsequently, new normalized feature indices were constructed based on the characteristic bands, which were selected based on spectral reflectance differences, and the segmentation thresholds were determined. The layer-by-layer segmentation method was used to eliminate the interference of irrelevant parts. Three methods(successive projections algorithm, principal component analysis, and quadratic combination dimensionality reduction) were used to reduce the spectral information redundancy and extract features. The regression models were established for the original spectral data and the reduced feature data by random forest and partial least squares regression, respectively. The best prediction model was determined to fit the hardness of the strawberry pulp. The hardness distribution image was obtained for the intuitive display of the strawberry hardness prediction result. The result shows that the partial least squares model based on quadratic dimensionality reduction yielded the best performance, with correlation coefficients of 0.9101 and 0.9099 for the test set and prediction set, respectively, and with a root-mean-square error of 0.1344 for the test set. This study provides a reference for non-destructive detection and display of strawberry hardness.
-
0. 引言
激光半自动制导是一种典型的激光制导模式[1],其基本原理是通过位于载机或地面的激光目标指示器发射激光束照射目标,由位于弹头的探测器接收目标漫反射的激光信号, 探测器通过对接收波门内多个脉冲回波的处理以实现对目标的捕获, 实现目标捕获后,通过对目标持续确认以实现对目标的跟踪。
其它激光目标指示器的散射信号和激光干扰机发射的信号都可能影响探测器的性能。前者是非恶意干扰,具有低重频干扰效果;后者通常工作在高重频工作模式,以保证干扰效果。干扰信号对探测器的捕获和跟踪环节都可能产生影响[2-3],具体的干扰效果与干扰机布设方位[4-5]、干扰信号重频[6-7]、干扰信号编码方式[8-9]有关,也与激光目标指示器信号波形[10]、探测器波门设计[11]、脉冲锁定技术[12]有关。
激光制导抗干扰的技术途径主要可以分为两种: 一是对指示器发射信号进行编码,提升目标回波信号信息熵,降低干扰信号与目标回波的相似性[13];二是在目标捕获和跟踪处理环节充分挖掘目标回波信号的特征信息,剔除不满足目标特征的干扰信号[14]。研究者深入分析了编码波形抗干扰技术,并分别提出了性能良好的随机编码波形和抗干扰算法[15-17]。WU等人[18]提出了基于接收波门临近回波的高重频干扰排查方法,改善了抗干扰效果。CHEN等人[19]提出了基于干扰重频估计的高重频干扰剔除技术,可有效对抗固定重频的高重频干扰。CAO等人[20]提出了一种自适应扩展实时波门技术,有效避免了干扰脉冲进入波门导致目标脉冲漏检的情况。
从原理上看,在满足激光目标指示器、探测器工程现实基础上,充分考虑感兴趣目标的行为特征和干扰机工作特点,开展发射波形精细化编码设计、脉间接收波门设计和抗干扰信号处理算法研究,可有效提升抗干扰性能。基于此,本文作者设计了一种随机重频发射信号,在收发准同步条件下,结合激光目标指示器、目标和探测器的运动特性,提出了端到端的基于自适应交叠波门的高重频干扰鉴别、目标探测方法, 仿真分析了该方法抗固定重频和随机重频干扰的性能。
1. 随机重频激光信号抗干扰模型
1.1 高重频干扰激光原理分析
探测器通过对特定距离上目标脉冲多个回波的连续探测实现目标捕获,并通过连续接收目标脉冲回波的连续观测实现目标跟踪。受激光目标指示器、目标和探测器运动影响,不同脉冲回波通常在时延上存在差异。为了有效实现脉间回波关联,通常设定与激光目标指示器、目标和探测器先验信息相关的具有一定宽度的距离波门。
当特定波门出现连续脉冲时,探测器进入脉冲锁定状态,并通过对锁定波门内的回波连续探测的关联实现目标捕获。干扰机通过发射高重频干扰信号,在探测器上产生大量密集虚假目标信号充斥于探测器接收波门之内,以较高的概率打断激光雷达目标捕获跟踪链路,严重影响激光制导效能。图 1是高重频干扰的原理图。图 1a中给出了无干扰时探测器接收目标回波信号的情况,此时目标回波信号被设置的波门较好地接收;图 1b中给出了高重频干扰信号的特征,即重复频率越高,在时间轴上分布越密集;图 1c中给出了在高重频干扰场景下探测器设置波门录取回波信号时,部分干扰信号进入了接收波门内。
在特定漏射率、检测概率条件下,探测器连续观测到特定距离上N个发射信号的回波脉冲时,可认为成功捕获到目标,其中N与指示器漏射率、探测器目标检测概率等有关,在当前工艺下N为3或4,本文中取干扰方案中,通常将脉冲编码技术、时间波门技术、脉冲锁定技术这3种技术相结合来得到更好的抗干扰效果,如图 2所示。假设探测器接收到的首信号是激光目标指示器发射的第l个发射脉冲的回波,分析第l个~第l+3个发射脉冲的回波,目标回波序列为[1, 3, 6, 9]T,其中上标T表示向量的转置。由于高重频干扰信号超前于目标回波进入接收波门,接收波门内的首脉冲被干扰信号占据,接收下一个脉冲回波信号的实时波门向前偏移,导致在接收第l+3个脉冲时无法捕获目标回波7。根据首脉冲锁定技术,此时系统判定的目标回波序列[1, 2, 4, 6]T,而实际目标回波序列为[1, 3, 5, 7]T,此时表明高重频干扰信号已经成功干扰了制导系统。
在高重频情况下,干扰信号可能在目标回波信号之前就进入探测器接收波门,导致首个接收信号的可靠性下降。此时,以首脉冲为代表的传统脉冲锁定方法存在误判干扰信号的概率高以及易被牵引的问题。因此,在高重频干扰环境中,需要在波形设计、时间波门和信号处理技术方面开展端到端的优化。
1.2 随机脉冲重复频率发射信号模型
激光目标指示器以随机重频方式发射信号,即脉冲周期为\boldsymbol{T}=\left[T_1, T_2, \cdots, T_{k-1}\right]^{\mathrm{T}},如图 3所示。激光目标指示器在工作时间内发射k个信号,可划分k-1个重复周期,定义Tk为第k-1个和第k个脉冲的重复周期,k取值范围为1≤k≤k-1,应有Tmin < Tk < Tmax,Tmin和Tmax分别为激光目标指示器发射信号的最小周期和最大周期,则激光指示器总的工作时间T_{\text {life }}=\sum\limits_{k=1}^{k-1} T_k。受传播损耗、目标散射特征、发射功率、探测器灵敏度等工程限制,通常情况下可探测回波最大时延τmax=2R/c≤Tmin,其中c为光速,R为有效探测距离。发射信号、脉冲重复周期、最大探测距离等参数如图 3所示。
假设发射脉冲信号波形为x(t),则发射波形为:
s_0(t)=\sum\limits_{k=0}^{k-1} x\left(t-t_k\right) (1) 式中: t_k=\sum\limits_{m=0}^k T_m,表示第k个脉冲相对第0个脉冲的发射迟延。第0个脉冲的发射时间为起始时间,T0=0。
1.3 干扰信号模型
当干扰机侦察到照射源发射信号时,干扰机发射高重频干扰信号以影响探测器效能,在高重频干扰场景下,探测器接收到的信号可以表示为:
s(t)=s_0{ }^{\prime}(t)+J(t)+n(t) (2) 式中: s0′(t)是发射脉冲信号对应的目标回波; J(t)为探测器接收到的干扰信号; n(t)为加性高斯白噪声。J(t)的表达式为:
J(t)=\sum\limits_{i=1}^{\infty} A_i x\left(t-{\rm i} t_{\mathrm{jam}}{ }^{\prime}\right) (3) 式中:tjam′是高重频干扰信号重复周期; Ai是干扰机发射的第i个固定高重频干扰信号的幅度。此时探测器接收到的目标回波以及干扰信号如图 4所示。
在随机重频干扰场景下,接收机接收到的干扰信号的重复周期相较于固定重频干扰信号有随机抖动,此时J(t)的表达式为:
J(t)=\sum\limits_{i=1}^{\infty} A_i x\left(t-\mathrm{i} t_{\mathrm{jam}}{ }^{\prime}-\Delta t_i\right) (4) 式中: Δti为第i个随机重频干扰脉冲的随机时间抖动,传统的反向重频脉冲信号电路[19]无法有效对抗此类干扰。
式(2)中的n(t)为加性高斯白噪声,可看作是出现概率较低、且呈现时空随机分布的背景回波,会使探测器出现虚警,在影响机理上看,与距离随机的低重频干扰相似。由于本文中考虑了随机脉冲重复周期情况下高重频J(t)干扰,可以认为噪声n(t)为特定J(t)的一部分,其影响机理可以在大样本统计实验中包含。在后续的分析中,将略去n(t)项的影响。
2. 基于交叠波门的全脉冲回波探测及抗干扰方法
传统目标捕获方法常采用固定波门和单脉冲锁定技术[11],即每个脉冲回波的波门宽度相等,而且每一个发射脉冲接收波门按前一个接收波门检测到的某一信号为时间基准设置。该技术易出现波门被高重频干扰引偏而丢失目标问题。与此不同,自适应交叠波门则通过首脉冲确定下一脉冲接收波门的开启时刻,通过末脉冲确定下一脉冲接收波门的关闭时刻[20]。本文作者在参考文献[20]的基础上,采用自适应交叠波门技术,提出基于先验信息的信号判决方法,实现目标有效捕获和抗干扰技术。
2.1 自适应交叠波门的确定
对脉冲重复周期为T=[T1, T2, …, Tk-1]T的发射激光序列,假设第l个脉冲的接收波门内接收到的首信号时间为ts,末信号的时间为te,则可以实时对第l+1个脉冲的波门进行设置:
\left\{\begin{array}{l} T_{\mathrm{g}, 1}=t_{\mathrm{s}}+\left(T_{l+1}-T_l\right)-\Delta t / 2 \\ T_{\mathrm{g}, 2}=t_{\mathrm{e}}+\left(T_{l+1}-T_l\right)+\Delta t / 2 \end{array}\right. (5) 式中: Tg, 1和Tg, 2是第l+1个脉冲实时波门的开启时刻和关闭时刻; Tl和Tl+1分别是第l个脉冲和第l+1个脉冲距离其下一个脉冲的脉冲重复周期; Δt是波门宽度,如图 5所示。
假设第l个脉冲发射之后,在第l个脉冲有效探测范围内接收到1号回波,以1号回波为基准时间,实时设置第l+1个脉冲的接收波门,并在此波门内接收到2、3号回波;以2号回波和3号回波为基准分别确定第l+2个脉冲接收波门的开启时刻和关闭时刻,并在此波门内接收到4、5、6号回波;以此递推下去,在第l+3个脉冲发射之后,探测器可以接收到7、8、9、10号回波,故将1号回波认为是疑似目标回波之后,通过自适应交叠波门得到的脉冲组合一共有4!=24种。接收机以接收4个发射脉冲的回波作为目标捕获的必要条件,在4个脉冲发射之后,在发射脉冲相对应的波门内收集到多组脉冲组合进入信号判决环节。
2.2 基于先验信息的信号判决
确定自适应交叠波门后,可以得到多个包含目标、干扰和噪声(如上文所述,噪声被看作干扰信号的一部分)的多个回波组合。由于发射脉冲时序T=[T1, T2, …, Tk-1]T已知, 通过编码先验信息和探测器运动特性,可以有效地排除固定重频干扰信号和随机重频干扰信号。
为简化问题,如第1.1节中所述,假定无干扰情况下,在4个接收波门中接收到至少3个满足时延特性的信号时, 可确认目标(具体探测脉冲数指标与漏设率、背景参数等决定,不是本文重点)。为不失一般性,假设打击目标静止不动,且激光目标指示器和探测器分别以2Ma的速度匀速飞向打击目标。假定某测试回波序列如图 5所示,此时的24种组合方式中,[1, 3, 5, 8]T组合是真实目标回波的组合,满足激光目标指示器和探测器运动特性,显然[1, 3, 6, 9]T不满足运动特性,可以直接鉴别为干扰组合。当有脉冲组合满足条件时,即可认为成功捕获到目标,否则回到目标捕获阶段,以第l+1个脉冲为起始判断脉冲,分析第l+1个脉冲至第l+4个脉冲的回波信号。算法流程图如图 6所示。
当接收到回波信号后,该算法依据激光探测器和激光目标指示器之间的时间同步关系,确定从第l个脉冲开始进行分析; 再根据发射脉冲先验信息及接收回波信号自适应生成第l个~第l+3个脉冲的交叠波门,以获得脉冲回波组合; 然后基于目标的运动特性对脉冲回波组合进行判决, 若判决结果符合运动规律,可以确定为已成功捕获目标;若未能捕获,则会将起始分析脉冲调整为第l+1个脉冲,并继续进行下一轮目标捕获。
2.3 性能评估指标
考虑到工程实际情况,设置激光目标指示器总工作时间Tlife=70 s,假定其从第l个脉冲开始检测回波信号后,判断第l个~第l+3个脉冲接收波门的回波信号是否为真实目标回波信号,为不失一般性,设置干扰信号与目标回波信号为01序列。
在仿真实验中,激光目标指示器和探测器的速率以680 m/s向目标移动,激光目标指示器的最大有效照射距离为25 km,激光漏射率为1‰,发射脉冲脉宽为10 ns,单个脉冲目标检测概率为98%,基础波门宽度为4 μs,各脉冲的自适应交叠波门由第2.1节中的算法得到。本文中随机生成了1200个重复频率均匀分布在20 Hz~30 Hz之间的发射脉冲信号,激光目标指示器按照生成的随机序列发射激光。全时段发射高重频干扰信号的干扰机位于目标的100 m处。
仿真生成与干扰信号混杂的目标回波信号,使用本文中算法对高重频干扰信号进行分析筛选,并计算出捕获到的目标的位置。当连续4个脉冲的回波信号符合运动规律时,即可认为目标捕获成功。为了验证本文中方法的性能,采用目标捕获时间(target acquisition time, TAT)指标进行性能评估,即:
T_{\mathrm{a}}=\sum\limits_{i=T_{\mathrm{s}}}^{T_{\mathrm{e}}} T_i (6) 式中: Ts为开始接收到目标回波的第1个脉冲间隔; Te为成功捕获到目标的最后一个脉冲间隔; Ta为捕获目标总消耗时间,Ta越小说明抗重频干扰算法性能越佳。
设激光目标指示器漏射率为1‰、检测概率为98%,干扰参数设置如表 1所示。共设置12种重频干扰仿真场景,其中随机重频干扰在固定重频干扰的基础上加了随机时间抖动, 仿真结果如图 7~图 9所示。图 7表示某一脉冲发射后受到200 kHz高重频干扰后,探测器接收到的回波信号,其中标注的信号为真实目标回波,对应目标距离探测器10 km;图 8表示利用本文中算法成功去除掉高重频干扰后,得到的真实目标回波脉冲组合,其起始时间为第45个脉冲的发射时间;图 9为在不同频率的固定重频干扰和随机重频干扰下各进行104次仿真实验的平均目标捕获时间。
表 1 干扰场景参数设置Table 1. Jamming scenario parameters configurationjamming types transmitted signal jamming frequency/kHz experiment fixed frequency random frequency distributed between 20 Hz~30 Hz 0.02, 0.5, 10, 100, 200, 250 10000 trials per frequency 由图 7和图 8可以看出,基于自适应交叠波门的全脉冲回波探测及抗干扰方法能够有效提取面对200 kHz的高重频干扰时目标的制导回波,并通过迟延计算确定当前目标与探测器的距离。图 9展示了不同频率干扰下目标平均捕获时间。实验结果表明,多种干扰场景下平均目标捕获时间均在0.146 s左右; 从低频到高频,目标平均捕获时间有轻微的提升; 多种干扰场景下平均目标捕获时间均小于0.2 s,这说明该算法能够在开始检测的前4个脉冲中捕获到真实目标回波的回波序列。
在多种干扰场景下的12×104次实验中,仿真结果未出现任何虚警情况,表明该算法具有良好的抗重频干扰性能。综上所述,本文中提出的抗干扰算法能够有效适应不同频率的固定重频干扰和随机重频干扰。
3. 结论
根据发射脉冲的编码规则,本文中设计了接收脉冲的自适应交叠波门,并利用探测器运动的先验信息对自适应交叠波门的接收信号进行了合理性判决。本文中的算法充分利用已知的随机脉冲规则与运动先验信息,分析连续N个脉冲接收波门的信号。由于N个脉冲对应的自适应交叠波门接收到的固定重频干扰信号以及随机重频干扰持续满足真实目标的回波规律的可能性极小,故此方法可以有效筛选去除干扰信号,对固定/随机重频干扰均具有较强的抗干扰性能。
在不同频率的固定/随机重频干扰场景下,利用仿真实验验证了本文中算法的有效性和优越性。仿真结果表明,在干扰场景下,随机重频激光制导雷达抗高重频干扰技术不仅可以成功捕获目标,而且目标捕获时间与无干扰场景下的捕获时间相当,保证了捕获目标的时效性,并且在12×104次统计试验中未出现虚警问题。此算法的优越性在于: 在高重频干扰环境下仍然具备高探测概率和低虚警概率,可用于提高激光制导雷达系统的抗干扰能力,提高其工作效率和可靠性。
-
表 1 PCA降维特征贡献度
Table 1 Feature contribution in PCA
number of features cumulative contribution rate 1 0.8669 2 0.9473 3 0.9801 4 0.9920 5 0.9970 6 0.9983 7 0.9991 8 0.9996 9 0.9997 10 0.9997 表 2 SPA-PCA降维特征贡献度
Table 2 Feature contribution in SPA-PCA
number of features cumulative contribution rate 1 0.7377 2 0.9300 3 0.9745 4 0.9902 5 0.9972 6 0.9991 7 0.9997 8 0.9997 表 3 基于PLS的模型建模结果
Table 3 PLS-based model modeling results
modeling methods CCC CCP RMSEC RMSEP RPD full-PLS 0.9084 0.8636 0.1307 0.1592 2.4973 PCA-PLS 0.9111 0.9047 0.1283 0.1352 3.1019 SPA-PLS 0.9069 0.8700 0.1323 0.1529 2.6451 SPA-PCA-PLS 0.9101 0.9099 0.1290 0.1344 3.2623 表 4 基于RF的模型建模结果
Table 4 RF-based model modeling results
modeling methods CCC CCP RMSEC RMSEP RPD full-RF 0.9611 0.8889 0.0845 0.1491 2.5344 SPA-RF 0.9321 0.8667 0.1143 0.1473 2.1535 PCA-RF 0.9522 0.8735 0.0926 0.1644 2.0531 SPA-PCA-RF 0.9620 0.8158 0.0837 0.1891 1.7907 -
[1] 白胜, 朱润华, 阳圣莹, 等. 不同草莓品种营养成分比较与品种筛选[J]. 山西农业科学, 2020, 48(1): 64-67. https://www.cnki.com.cn/Article/CJFDTOTAL-SXLX202001018.htm BAI Sh, ZHU R H, YANG Sh Y, et al. Nutrients comparison and screening of different strawberry varieties[J]. Journal of Shanxi Agricultural Sciences, 2020, 48(1): 64-67(in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-SXLX202001018.htm
[2] AFRIN S, GASPARRINI M, FORBES-HERNANDEZ T Y, et al. Promising health benefits of the strawberry: A focus on clinical studies[J]. Journal of Agricultural and Food Chemistry, 2016, 64(22): 4435-4449. DOI: 10.1021/acs.jafc.6b00857
[3] 刘凯丽, 郑洪健, 黄敏, 等. 云南香格里拉3个草莓品种的感官品质及耐储藏性对比分析[J]. 昆明学院学报, 2022, 44(3): 87-92. https://www.cnki.com.cn/Article/CJFDTOTAL-KMSG202203017.htm LIU K L, ZHENG H J, HUANG M, et al. Comparative analysis of sensory quality and storage tolerance of three strawberry varieties in Yunnan Shangri-La[J]. Journal of Kunming University, 2022, 44(3): 87-92(in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-KMSG202203017.htm
[4] 郝乾坤, 王智民, 童开林. 外源海藻糖对草莓保鲜效果的影响[J]. 陕西农业科学, 2022, 68(8): 84-87. https://www.cnki.com.cn/Article/CJFDTOTAL-SNKX202208017.htm HAO Q K, WANG Zh M, TONG K L. Effect of exogenous trehalose on preservation of strawberry[J]. Shaanxi Journal of Agricultural Sciences, 2022, 68(8): 84-87(in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-SNKX202208017.htm
[5] 陈卓, 宋俏微, 张水洞, 等. 氧化淀粉的抑菌效果及草莓涂膜保鲜应用[J]. 食品科学, 2022, 43(21): 324-331. https://www.cnki.com.cn/Article/CJFDTOTAL-SPKX202221039.htm CHEN Zh, SONG Q W, ZHANG Sh D, et al. Antimicrobial effect of oxidized starch and its application as a coating for strawberry preservation[J]. Food Science, 2022, 43(21): 324-331(in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-SPKX202221039.htm
[6] LIU Y, WU Q W, HUANG J L, et al. Comparison of apple firmness prediction models based on non-destructive acoustic signal[J]. International Journal of Food Science & Technology, 2021, 56(12): 6443-6450.
[7] 陈广大, 刘德君, 李天旭, 等. 便携式水果硬度测量仪的设计[J]. 中国农机化学报, 2016, 37(7): 85-88. https://www.cnki.com.cn/Article/CJFDTOTAL-GLJH201607019.htm CHEN G D, LIU D J, LI T X, et al. Design of the portable fruit hardness measuring instrument[J]. Journal of Chinese Agricultural Mechanization, 2016, 37(7): 85-88(in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-GLJH201607019.htm
[8] PARK B, SHIN T S, CHO J S, et al. Characterizing hyperspectral microscope imagery for classification of blueberry firmness with deep learning methods[J]. Agronomy, 2021, 12(1): 85. DOI: 10.3390/agronomy12010085
[9] QIAO M, XU Y, XIA G, et al. Determination of hardness for maize kernels based on hyperspectral imaging[J]. Food Chemistry, 2021, 366(2): 130559.
[10] 马帅帅, 于慧春, 殷勇, 等. 黄瓜水分和硬度高光谱特征波长选择与预测模型构建[J]. 食品与机械, 2021, 37(2): 145-151. https://www.cnki.com.cn/Article/CJFDTOTAL-SPJX202102028.htm MA Sh Sh, YU H Ch, YIN Y, et al. Selection of hyperspectral characteristic wavelength and construction of prediction model for cucumber hardness and moisture[J]. Food & Machinery, 2021, 37(2): 145-151(in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-SPJX202102028.htm
[11] ERKINBAEV Ch, DERKSEN K, PALIWAL J, et al. Single kernel wheat hardness estimation using near infrared hyperspectral imaging[J]. Infrared Physics and Technology, 2019, 98(6): 250-255.
[12] SU J Y, YI D W, LIU C J, et al. Dimension reduction aided hyperspectral image classification with a small-sized training dataset: Experimental comparisons[J]. Sensors, 2017, 17(12): 2726. DOI: 10.3390/s17122726
[13] XIONG J T, LIN R, BU R B, et al. A micro-damage detection method of litchi fruit using hyperspectral imaging technology[J]. Sensors, 2018, 18(3): 700. DOI: 10.3390/s18030700
[14] WU D, MENG L W, YANG L, et al. Feasibility of laser-induced breakdown spectroscopy and hyperspectral imaging for rapid detection of thiophanate-methyl residue on mulberry fruit[J]. International Journal of Molecular Sciences, 2019, 20(8): 2017. DOI: 10.3390/ijms20082017
[15] 赵凡, 闫昭如, 宋海燕. 应用高光谱鉴别黑枸杞和唐古特白刺果[J]. 光谱学与光谱分析, 2021, 41(7): 2240-2244. https://www.cnki.com.cn/Article/CJFDTOTAL-GUAN202107046.htm ZHAO F, YAN Zh R, SONG H Y. Hyperspectra used to recognize Black Goji berry and nitraria Tanggu[J]. Spectroscopy and Spectral Analysis, 2021, 41(7): 2240-2244(in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-GUAN202107046.htm
[16] 北京市农林科学院林业果树研究所; 农业部果品及苗木质量监督检验测试中心. 水果硬度的测定: NY/T 2009-2011[S]. 北京: 中华人民共和国农业部, 2011: 1-5. BEIJING ACADEMY OF FORESTRY AND POMOLOGY SCIENCES, LABORATORY OF QUALITY & SAFETY RISK ASSESSMENT FOR FRUIT. Determination of fruit firmness: NY/T 2009-2011[S]. Beijing: Ministry of Agriculture of the PRC, 2011: 1-5(in Chinese).
[17] 姜一河, 王涛, 常红伟. 高光谱图像特征提取方法研究综述[J]. 电光与控制, 2020, 27(10): 73-77. https://www.cnki.com.cn/Article/CJFDTOTAL-DGKQ202010015.htm JIANG Y H, WANG T, CHANG H W. An overview of hyperspectral image feature extraction[J]. Electronics Optics & Control, 2020, 27(10): 73-77(in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-DGKQ202010015.htm
[18] 白丽敏, 李粉玲, 常庆瑞, 等. 结合SPA和PLS法提高冬小麦冠层全氮高光谱估算的精确度[J]. 植物营养与肥料学报, 2018, 24(5): 1178-1184. https://www.cnki.com.cn/Article/CJFDTOTAL-ZWYF201805006.htm BAI L M, LI F L, CHANG Q R, et al. Increasing accuracy of hyper-spectral remote sensing for total nitrogen of winter wheat canopy by use of SPA and PLS methods[J]. Journal of Plant Nutrition and Fertilizers, 2018, 24(5): 1178-1184(in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-ZWYF201805006.htm
[19] 欧阳爱国, 万启明, 李雄, 等. 高光谱成像的水稻螟虫蛀入检测方法[J]. 光谱学与光谱分析, 2021, 41(12): 3844-3850. https://www.cnki.com.cn/Article/CJFDTOTAL-GUAN202112037.htm OUYANG A G, WAN Q M, LI X, et al. Research on rich borer detection methods based on hyperspectral imaging technology[J]. Spectroscopy and Spectral Analysis, 2021, 41(12): 3844-3850(in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-GUAN202112037.htm
[20] FAN L L, ZHAO J L, XU X G, et al. Hyperspectral-based estimation of leaf nitrogen content in corn using optimal selection of multiple spectral variables[J]. Sensors, 2019, 19(13): 2898.
[21] 刘璐, 邵慧, 孙龙, 等. 利用高光谱激光雷达检测木材的霉变与含水量[J]. 激光技术, 2023, 47(5): 620-626. DOI: 10.7510/jgjs.issn.1001-3806.2023.05.007 LIU L, SHAO H, SUN L, et al. Detection of mildew and moisture content in timber by hyperspectral LiDAR[J]. Laser Technology, 2023, 47(5): 620-626(in Chinese). DOI: 10.7510/jgjs.issn.1001-3806.2023.05.007
[22] MARSHALL M, BELGIU M, BOSCHETTI M, et al. Field-level crop yield estimation with PRISMA and Sentinel-2[J]. ISPRS Journal of Photogrammetry and Remote Sensing, 2022, 187(2): 191-210.
[23] 高升, 徐建华. 高光谱成像的红提总酸与硬度的预测及其分布可视化[J]. 食品科学, 2023, 44(2): 327-336. https://www.cnki.com.cn/Article/CJFDTOTAL-SPKX202302041.htm GAO Sh, XU J H. Hyperspectral imaging for prediction and distribution visualization of total acidity and hardness of red globe grapes[J]. Food Science, 2023, 44(2): 327-336(in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-SPKX202302041.htm
[24] SAVAŞLI E, KARADUMAN Y, ÖNDER O, et al. Estimating technological quality parameters of bread wheat using sensor-based norma-lized difference vegetation index[J]. Journal of Cereal Science, 2022, 107(8): 103535.
[25] QIN J, BAI H Y, ZHAO P, et al. Dendrochronology-based norma-lized difference vegetation index reconstruction in the Qinling Mountains, North-Central China[J]. Forests, 2022, 13(3): 443.
[26] PAVLO L. Forecasting oil crops yields on the regional scale using normalized difference vegetation index[J]. Journal of Ecological Engineering, 2021, 22(3): 53-57.
[27] PENG Y, FAN M, BAI L, et al. Identification of the best hyperspectral indices in estimating plant species richness in sandy grasslands[J]. Remote Sensing, 2019, 11(5): 558.