花生作为我国主要的油料作物和经济作物之一,因其具有较高的产量和出油率,在油料作物中占据着独特地位。与其他油料作物如大豆相比,花生果实不仅具有较高的油脂含量,较强的耐旱性和适应性,而且具有极高经济价值
[1]。准确地获取花生种植区域的生物量信息对于农田生产管理和保障农户经济利益至关重要。因此,及时、准确且无损地获取花生生物量信息对于监控植株的生长健康状况、优化农田管理措施,实现花生精细化种植管理具有重要的实际意义
[2-4]。
传统的作物生物量估算方法,如统计调查和气象估产等,虽然被广泛使用,但往往耗时费力,且估算结果存在时间和空间精度上的误差
[5],难以满足精细化种植管理的需求。近年来,无人机遥感技术因其操作简单、机动性强和高时空分辨率的特点,在农业领域得到了广泛应用。特别是在作物产量、生物量监测方面,使用无人机遥感技术已经成为现代农业研究和实践中的重要趋势
[6-7]。其中,高光谱能够探究农作物光谱特征上的微小差异,常用于建立作物的各种反演模型,能够提供快速、准确的生物量数据,极大地增强了对农业生产的监控和管理
[8]。陶惠林等
[9]利用无人机搭载的高光谱相机获取冬小麦各个时期的遥感影像,选取了9种植被指数和5种红边参数,构建了不同生育期的小麦产量估算模型;Ashapure 等
[10]成功地应用高分辨率无人机影像和机器学习方法构建了棉花产量估算模型,该模型综合了时序、非时序和定性特征,在多次测试的平均情况下,模型展示了出色的预测能力;向友珍等
[11]以西北地区的大豆为研究对象,从无人机多光谱影像中筛选出最具相关性的植被指数,基于支持向量机(Support Vector Machine, SVM)、随机森林(Random Forest, RF)和反向神经网络(Back Propagation Neural Network, BPNN)分别构建了大豆LAI、地上部生物量和产量的估计模型。为了提高模型的实用性并增强其稳定性,一些研究在原始光谱数据的基础上采用了连续小波变换和光谱微分等技术,能够充分挖掘更深层次的和潜在的有效光谱信息,从而有效提高对产量和生物量等参数的估算精度。比如:刘杨等
[12]结合使用多个植被指数和连续小波转换技术,从中提取小波系数并进行相关性分析,筛选出最具相关性的前10个波段,并使用偏最小二乘回归(Partial Least Squares Regression, PLSR)方法构建了马铃薯地上生物量的估算模型;李长春等
[13]将冬小麦在4个生育期的原始光谱与小波能量系数进行相关性分析,后进行筛选,并基于支持向量回归、随机森林和高斯过程回归3种算法构建了冬小麦生物量估算模型。
已有研究主要集中在玉米、小麦、大豆等作物的地上生物量估算,而利用无人机高光谱技术探究花生整体生物量估算的研究仍较为缺乏
[14]。随着近年来作物秸秆在发酵制氢
[15]和制肥
[16]等农业废弃物高值化利用技术中的应用日益受到关注,本研究综合考虑了花生地上部分和地下部分的生物量,全面评估了花生的整体生物量,涵盖其副产品的应用价值。因此,本研究拟以河南荥阳花生种植试验田为研究对象,通过对高光谱影像使用小波变换得到小波能量系数,计算植被指数,结合SVR、RF、CNN及使用粒子群优化这些模型PSO-RF、PSO-SVR以及PSO-CNN的方法对花生生物量进行建模和评估,分析不同模型对于生物量的预测效果,以期寻找构建花生生物量反演的最佳模型,为农田地块尺度的花生生物量监测提供一种有效的参考。
1 材料与方法
1.1 研究区概况
研究区域荥阳市位于河南省中部,地处黄河的中下游平原地区,研究区的地理位置以及试验田样本点示意图见
图1。该地区地势相对平坦,土壤肥沃,土地资源丰富,适宜农业种植,特别是粮食作物的生产。荥阳市的气候类型属于温带大陆性季风气候,四季分明,历年平均气温14.8 ℃,降水集中在6—9月,历年平均降水量608.8 mm,该地区以农作物种植为主,常见的作物包括花生、小麦、玉米等。
1.2 数据获取及处理
1.2.1 光谱数据获取和处理
本研究于2022年10月10日花生收获当天,利用大疆M600 Pro六旋翼无人机搭载四川双利合谱有限公司的GaiaSky-mini2-VN高光谱成像仪,获取了研究区域内花生田块的高光谱影像数据。飞行高度设置为50 m,旁向和航向重叠率均为80%,空间分辨率约为1.8 cm/像素。该高光谱成像仪具有176个光谱反射率通道,覆盖的光谱范围为400~1 000
,光谱分辨率为3.5 nm。采集的影像数据通过高光谱系统处理软件SpecView进行了一系列预处理,包括黑白帧校准、反射率校准、大气校正以及镜头变换等。在光谱数据的采集过程中,高光谱传感器可能受到外部环境因素或传感器自身的影响,通常表现为光谱曲线上的不平滑或毛刺。为解决上述问题,本研究采用Savitzky-Golay(以下简称S-G)滤波器进行卷积平滑处理,以提高光谱的信噪比并减少外部干扰的影响;具体而言,使用了窗口大小为13、阶数为5的S-G滤波参数,从而显著改善了影像的平滑效果。平滑处理前后的影像反射率效果如
图2所示。
1.2.2 花生生物量获取
本研究选取多品种花生试验田进行数据采集,该试验田种植:豫科花、龙乡红、驻花、周科花、星花、濮花、宛花、周花、远杂、豫花等在内11个大类花生品种,种植试验田共66个小区,每个小区的长为6.7 m,宽为2 m,种5行空一行,每个品种重复3次。在花生饱果成熟期,对采样点0.5 m×0.5 m范围内的所有花生植株进行采摘,剔除根部泥土后装袋密封,将采集的花生植株带回实验室进行清洗;放置在通风良好处,自然风干后,将采样点植株放入烤箱中65~75 ℃烘至恒重,进行称重作为各采样点的生物量干重实测值。为了进行模型训练和验证,从70个采样点中随机选取80%样本点作为训练集,剩余的20%样本点作为验证集。样本点生物量的统计结果如
表1所示。
1.3 数据处理及模型构建
1.3.1 小波变换
小波变换是一种常见的数据处理方法
[17-18]。该方法基于傅里叶变换增加了时间尺度,使得信号可以在时间域和频域进行分解,从而得到原始信号中的深层次隐含的信息。在高光谱数据处理中,基于小波母函数的平移和缩放对高光谱数据进行卷积运算,在不改变原始波段范围和位置的情况下将高光谱信息分解成不同尺度下和不同波段长度的小波能量系数。其计算公式如下:
式中 为小波能量系数;为地上花生叶高光谱反射率;为光谱波段(400~1 000 nm);为选取Gaussian4作为小波基函数;为尺度因子;为平移因子。本研究将经过SG滤波平滑处理的高光谱数据进一步通过使用Matlab 2023a软件进行连续小波变换分析。在1~10尺度(21、22、23、24、25、26、27、28、29、210)进行分解,使用Gaussian4作为小波分析的母小波函数,以提取高光谱数据中的隐含有效信息。
1.3.2 植被指数构建方法
利用小波变换处理后的光谱反射率数据来构建植被指数。植被指数通过线性或非线性计算获得,不仅减少了干扰信息,还突出了植被的关键信息,从而能够反映出植被的深层次特征。基于已有研究成果
[19-21],选用了53个与作物产量或生物量相关的经验植被指数,分别与实测花生生物量进行皮尔逊相关性进行分析。考虑到模型的简洁性以及特征冗余易造成的多重共线性的问题,本研究仅筛选出相关性最高的前1%,即前5个植被指数作为特征向量。
所筛选的植被指数其计算公式如
表2所示。其中,PSRI对于监测植被冠层胁迫性变化十分敏感,常用于作物生产力的预测模型;DCNIⅠ则被视为预测植物氮浓度的最佳光谱指标,对作物生物量预测具有极大的潜力;MMTCI作为经过改性的陆地叶绿素指数(MTCI),展现了比原始MTCI更高的提取精度和抗干扰能力,常用于估算作物的叶绿素含量;SAVI是一个广泛使用于作物估产的植被指数;VARI因其在预测作物水分胁迫方面的性能而被广泛应用于小麦和苎麻等作物的估产模型。
1.3.3 模型构建方法
1)粒子群优化算法。粒子群优化算法(PSO)是一种启发式优化算法,每个粒子被初始化速度和位置之后在解空间中迭代寻找全局最优解,并通过适应性评价来记录和比较解的品质,具有收敛速度快且易于实现等优点,PSO算法被广泛应用于模型优化等问题
[22-23]。为了提升RF、SVR和CNN模型的拟合能力,并减少参数设置的误差影响,本研究采用了PSO来优化这些模型的参数,具体算法流程如
图3所示。具体而言,设置学习因子c1和c2均为1.5,惯性权重设置为0.73。在初始化种群和模型超参数后,并将其应用于各种模型中,对每次迭代的适应度进行记录,并判断是否达到终止条件;如果未达到终止条件,则继续更新粒子的位置和速度参数,直至满足终止条件。
2)随机森林。随机森林(RF)是一种应用于分类和回归问题的集成学习方法,基于多个决策树构建。此方法通过构建并组合多个决策树的预测结果来提高整体预测性能,每棵决策树都是独立地从原始数据集中随机选择样本和特征进行训练。这种随机选择增加了模型的多样性,使得RF在回归预测中显示出优异的鲁棒性
[24]。在本研究中,利用PSO算法对RF模型的关键参数进行优化,具体包括决策树的数量和叶节点的最小样本数。参数的初始范围详见
表3。
3)支持向量回归。支持向量回归(SVR)是从支持向量机算法衍生而来,可用于解决回归拟合问题。SVR通过使用基核函数,如高斯核和多项式核,以及通过梯度下降算法优化权重系数,来构建最优的决策边界;这一过程包括构建一个超平面,将输入样本映射到高维空间中,从而在该空间中最小化数据点的距离,使得数据更加线性可分。这样构建出来的回归模型不仅预测性能出色,而且具有很高的鲁棒性
[25]。在本研究中,进行优化的参数包括惩罚因子和径向基函数参数。参数的初始范围详见
表3。
4)卷积神经网络。卷积神经网络(CNN)是一种深度学习架构,主要由输入层、卷积层、池化层和全连接层组成,本研究使用的CNN网络结构如
图4所示。其中,卷积层负责提取输入数据的特征,池化层用于减小数据的维度,而全连接层则利用Softmax激活函数整合这些特征以产生最终输出
[26]。CNN的性能和准确性在很大程度上受到其超参数的影响,如初始化学习率、批处理大小和L2正则化系数等
[27]。因此,优化这些超参数是提高模型训练效率和预测精度的关键。本研究采用PSO算法对CNN的超参数进行优化,以加速收敛过程并提高模型的总体性能。超参数的初始优化范围详见
表3。这种方法通过系统地调整关键参数,旨在实现更优的网络训练结果和提升模型的泛化能力。
1.3.4 模型精度评价
为了评估花生生物量估算模型的精度,将随机选取的20%样本点的样本数据带入已训练好的模型中进行测试,采用决定系数()、均方根误差(RMSE)、均方误差(MSE)、平均绝对误差(MAE)对模型的性能进行评价。具体计算公式如下:
式中:为样本集的个数,;为花生生物量实际测量值;为花生生物量的预测值;为花生生物量的平均值。其中,接近于1, RMSE、MSE和MAE接近于0,说明模型的拟合效果越好。
2 结果与分析
2.1 小波分析结合植被指数特征构建
图5展示了基于高光谱原始反射率直接计算的植被指数与花生生物量之间的相关性。包括相关系数以及通过小波变换分解为10个尺度后计算的相关系数。可见:不同植被指数与花生生物量的相关系数范围为-0.49~0.45,部分植被指数与生物量呈正相关;另一些则呈负相关。这种广泛的相关性范围反映了不同植被指数在捕捉植被生长状况及其生理反应方面的差异,尤其是在不同小波尺度下,揭示了不同尺度的小波分解在捕捉与生物量相关的重要光谱信息方面的有效性。
由
图5可见,相关性较高的植被指数主要集中在第3、4、5、7、8和9尺度,表明这些尺度的小波变换能够更有效地捕捉与花生生物量密切相关的重要光谱特征。具体而言,本研究使用第3尺度下的土壤调节植被指数(SAVI)、可见光大气修正植被指数(VARI)、第4尺度下的植被衰老反射率指数(PSRI)、修正地面叶绿素指数(MMTCI)和第5尺度下的双峰冠层氮素指数Ⅰ(DCNIⅠ)作为模型输入,其相关性值分别为0.448、-0.486、0.377、0.389和0.370。可以看出,不同尺度的植被指数在捕捉与花生生物量相关的光谱特征方面具有显著差异,且某些特定尺度下的指数与生物量之间存在较强的相关性。因此,小波变换有助于解析植被指数与生物量之间的复杂关系,经过小波变换分解后的植被指数能够有效保留植被指数在捕捉生物量信息方面的敏感性,为多尺度分析方法提供了有力依据。
2.2 不同模型构建及精度评价
利用筛选出的植被指数作为特征变量,以构建花生生物量反演模型。具体而言:使用Matlab 2023a软件与多种机器学习模型建模,包括SVR、RF、CNN及其通过粒子群优化算法(PSO)增强的版本,PSO-SVR、PSO-RF和PSO-CNN,共计获得6种模型。
表4中展示了每个模型的 4 种性能评价指标的检测结果,包括:决定系数(
R²)、均方根误差(RMSE)、均方误差(MSE)和平均绝对误差(MAE)。由
表4可见,PSO 优化后的模型普遍表现更优,尤其是在
R² 和 RMSE 方面,显示出更高的预测准确性和较低的误差,这表明使用PSO优化后的超参数有助于提高模型拟合效果。
6种模型在训练集和测试集上的检验结果如
图6所示。可见在这6个模型中,大部分数据点紧邻1∶1虚线,表明所建模型能够有效地预测训练集和测试集的花生生物量,反映了模型的预测准确性和泛化能力。由
图6(a)、(c)和(e)可见:在机器学习模型中,RF在测试集上的决定系数(
R²)为0.674,均方根误差(RMSE)为0.429 kg/m
2,均方误差(MSE)和平均绝对误差(MAE)分别为0.155和0.28 kg/m
2; SVR测试集的
R²为0.651,RMSE为0.407 kg/m
2,MSE和MAE分别为0.165和0.346 kg/m
2。在深度学习模型中,使用默认参数的CNN测试集的
R²为0.710,RMSE为0.371 kg/m
2,MSE和MAE分别为0.138和0.329 kg/m
2,表明深度学习CNN模型的预测性能相比机器学习模型(RF和SVR)有一定提升,即深度学习算法在对花生生物量的建模中的拟合效果要优于机器学习模型。
由
图6(b)、(d)和(f)可见:PSO-RF和PSO-SVR在测试集中的
R²、RMSE、MSE和MAE分别为0.689、0.384、0.147、0.274和0.674、0.393、0.154、0.334 kg/m
2。PSO-CNN在测试集上表现最佳,
R²为0.768,RMSE为0.331 kg/m
2,MSE和MAE分别为0.11和0.283 kg/m
2。与PSO-RF 和PSO-SVM模型相比,PSO-CNN模型的拟合效果达到最优,说明经过PSO优化后的深度学习模型优于PSO优化后的机器学习模型。将模型预测精度从高到低排序,未经PSO优化的模型为CNN、RF、SVR;而经过PSO优化的模型排序为PSO-CNN、PSO-RF、PSO-SVR。因此经过PSO优化的深度学习模型(PSO-CNN)在预测花生生物量方面表现最佳,显著优于其他未经优化和经过优化的机器学习模型,进一步验证了PSO方法在提升模型预测精度方面的有效性。
3 讨 论
相比以往研究利用原始光谱反射率构建植被指数估算作物生物量,本研究利用小波变换对相关性较强的植被指数进行分解,并利用分解后的特征构建生物量估算模型。结果发现在未分解的53种植被指数中,双峰冠层氮素指数Ⅰ(DCNI Ⅰ)与实测生物量的相关性最高,相关系数为0.403,这与裴浩杰等
[20]研究结论相似,但仍低于本研究中使用小波变换后的植被指数SAVI(0.448),本研究使用相同的方法,挑选相关性最强的5个未经小波变换的植被指数,并构建6种花生生物估算模型(SVR、RF、CNN、PSO-SVR、PSO-RF、PSO-CNN),结果表明PSO-CNN模型精度最高,决定系数
R²为0.65。
利用PSO优化的模型精度普遍高于未优化之前,PSO-CNN模型在使用高光谱影像预测花生生物量方面表现出最高的精度,试验结论与王晗等
[22]和刘兰军等
[23]研究结果一致。其中,PSO对于RF和SVR的提升效果有限,对于CNN提升效果较为明显,训练集的决定系数从0.710提升到了0.768,提升了约为8.2%。这一结果也突显了PSO在优化深度学习模型参数中的有效性,为精确估算花生生物量提供了一个强有力的工具。
本研究展示了无人机高光谱遥感技术在精准农业中的应用潜力,特别是在花生生物量估算方面,所提出的高效、无损的花生生物量估测方法,不仅增强了现代化作物管理的科学性,为农业监测、农业保险等应用提供了可靠的依据,有助于在农业估产领域实现更高水平的推广应用,为发展低空经济和数字化农业提供重要参考。
4 结 论
为提高田块尺度花生生物量估算精度、实现精准农业管理,本研究基于无人机高光谱成像技术,探索了结合小波变换和植被指数的特征构建方法,从53种经验植被指数中筛选出相关性较强的特征,并通过机器学习模型、深度学习模型以及粒子群优化算法(PSO)构建花生生物量估算模型。主要结论如下:
1)在小波变换第3尺度下的SAVI 和 VARI 、第4尺度下的PSRI 和 MMTCI 和第5尺度下的DCNIⅠ 与实测花生生物量具有较强的相关性。相较于直接利用原始波段构建的植被指数,本研究构建的最优模型在决定系数(R²)上提升了约16.8%;
2)深度学习模型CNN在花生生物量预测精度上优于传统的机器学习模型,如RF和SVR,在花生生物量预测模型精度上表现出更优的效果;
3)通过粒子群算法PSO进行参数优化后,RF、SVR、CNN模型的预测精度都有提升,其中CNN的提升较为明显,模型精度提升约为8.2%。使用PSO-CNN模型对于花生生物量的估算能够获得最优结果,测试集R²为0.768。
河南省科技研发计划联合基金青年科学家项目(225200810089)
河南省高校科技创新人才支持计划(23HASTIT028)
河南省青年人才托举工程项目(2024HYTP018)
河南省博士后科研项目(202101041)