冬小麦是全球重要的粮食作物之一,其对氮元素的吸收直接影响作物的生长、产量和品质
[1]。氮肥过多或不足都会对作物产生不良影响,表现为作物晚熟、环境污染、叶绿素含量下降、产量减少等问题。在冬小麦的生长过程中,叶片氮含量随着施氮量的增减和生长时期的推进呈现出明显的变化趋势。灌浆期,冬小麦的养分吸收主要集中在籽粒上,籽粒膨大过程中需要消耗大量氮素。在此期间,大量氮素被再分配到籽粒中,植物叶片和茎秆部分的氮含量会相对降低
[2-3]。因此,通过动态监测冬小麦叶片氮含量对评估小麦生长状态进而实现合理施用氮肥有着重要意义。然而,传统的氮含量检测方法操作步骤繁琐、效率低、成本高,在现代农业的发展中难以推广应用。随着无人机和高光谱遥感技术的发展,利用无人机高光谱遥感技术评估作物氮含量成为了一种无损、高效、低耗的替代方案
[4]。无人机高光谱成像技术可以在大范围内实时获取大量的作物光谱反射率数据,包括从可见光到近红外波段的连续光谱
[5]。光谱数据通常包含多个波段,每个波段代表不同的光谱信息,波段范围的不同会显著影响模型的性能和结果
[6]。高光谱全波段数据可能存在波段冗余,不同波段之间存在高度相关性,使用冗余的波段会增加计算负担,导致模型过拟合或学习无关的特征。通过特征波段提取筛选与关键指标相关性较高的波段,可以减少不必要的复杂性,从而提高模型的准确性和泛化能力
[7]。无人机高光谱技术结合特征波段选择方法,可以提取作物在氮素含量、健康状况、光合作用效率等方面的敏感波段反射率,能够精准反映作物生长过程中的生理变化。这些技术不仅可以帮助提高农业生产效益,减少环境污染,也为精确农业的可持续发展提供了有力支持。
近年来,利用光谱遥感技术监测作物氮素含量方面的研究已经取得了丰富成果。已有研究表明,无人机高光谱遥感技术在监测大范围的作物氮素含量上具有明显优势
[8-10]。一些研究采用偏最小二乘回归(PLSR)、支持向量机(SVM)、随机森林(RF)、BP神经网络(BPNN)等方法构建了基于无人机高光谱遥感数据的氮含量评估模型。利用这些模型对小麦、水稻、玉米等作物的氮含量进行评估后,取得了良好的效果
[11-16]。ZHANG等
[14]采用连续投影算法(SPA)对光谱数据进行初始压缩,筛选出对氮含量贡献较大的波段进行建模,有效降低了模型的复杂度,并提高了氮含量评估的精确度。张黎黎等
[15]通过筛选与水稻叶片氮含量相关性较强的植被指数,结合BP神经网络模型构建了监测水稻氮素状况模型,验证集的
R2达到了0.669,该模型能够有效地估测水稻氮素变化。杨福芹等
[16]利用高光谱遥感数据分析了冬小麦各生育期的氮素营养状况,结果发现,利用无人机高光谱对作物生长状况和氮素营养动态进行监测,具有较高可靠性。上述研究表明,无人机高光谱遥感技术不仅在提高作物氮素估算的准确性方面具有潜力,而且在作物生长和氮素营养监测的动态性、时效性及可靠性上也展现出较大优势。然而,大多数研究局限于单一品种冬小麦或少量的冬小麦品种,且在同一生育期,不同氮肥处理下的冬小麦光谱特征存在明显差异
[17]。因此,研究适用于不同氮肥处理下、具有较强稳定性和泛化能力的冬小麦叶片氮含量高光谱遥感监测模型具有重要意义。
本研究以65种不同基因型的冬小麦品种为研究对象,探讨不同氮肥水平下冬小麦叶片氮含量的无人机高光谱响应特性,进而研究不同氮肥水平处理对选择冬小麦叶片氮含量监测模型的影响,实现基于无人机高光谱数据对冬小麦叶片氮含量的高通量监测,旨在为精准农业提供理论依据和技术参考。
1 材料和方法
1.1 试验地概况
试验于2022年10月至2023年6月进行。地点位于山西省太谷区山西农业大学孟家庄实验基地(37°44'N,112°59'E)。属暖温带大陆性气候,四季分明。冬季平均气温小于10 ℃,夏季平均气温大于22 ℃,春季、秋季平均气温为10~22 ℃。
1.2 试验材料
冬小麦群体试验材料为65份,包括农家品种、主要品种、新品种、核心种质和外源品种的冬小麦品种,这些小麦品种由课题组前期收集所得。
1.3 试验方法
试验共195个小区,小区长2.0 m,宽1.5 m,播种行距为20 cm。种植区域设置3个氮肥浓度梯度,分别为0(N1)、100(N2)、200(N3)kg/hm2。在施肥时,50%的氮肥用于基肥,剩余50%用于拔节期追肥,磷肥、钾肥用量均为120 kg/hm2,且均作为基肥一次性施入,其中,氮肥为尿素含N量46.1%,磷肥含P2O5量16%,钾肥含K2O量50%,其他小麦田间管理与当地小麦管理相同。
1.4 数据采集
在冬小麦生长前期,土壤中储存的氮素可以满足冬小麦早期生长的需求,不同施氮水平下的冬小麦生长差异不明显。抽穗期,冬小麦对氮素的需求增加,不同基因型的冬小麦在不同施氮水平下的表现差异明显
[17]。因此,本研究选取了抽穗期、灌浆前期和灌浆后期的冬小麦作为研究对象,以探讨不同施氮水平对冬小麦叶片氮含量的影响。利用大疆MATRICE 600 PRO无人机搭载高光谱相机,获取高光谱成像数据。数据于晴朗天气10:00—14:00进行采集。无人机飞行高度设定为30 m。高光谱波段范围为383~1 021 nm,共150个中心波段,光谱分辨率为4.1 nm。采集的数据经过影像校正和拼接,提取目标区域的反射率数据。采集试验地冬小麦叶片样品测定叶片氮含量。叶片样品于105 ℃下杀青30 min后,80 ℃下烘干至恒质量,经粉碎、消煮后,将消煮液通过全自动化学元素分析仪进行比色分析,获得叶片氮含量。
1.5 数据分析
采用连续投影算法(SPA)提取与冬小麦叶片氮含量相关的高光谱特征波段,并结合随机森林(RF)和BP神经网络(BPNN)构建冬小麦叶片氮含量监测模型。SPA通过减少矢量空间的共线性来选择特征波长,能够有效去除原始光谱数据中的冗余信息,适用于光谱特征波段的筛选
[18]。RF通过构建多个回归树,树与树之间互不干扰,最终通过综合所有决策树的输出结果来提升预测稳定性与准确性
[19]。BPNN基于误差反向传播原理,通过梯度下降法优化网络权重,以减小实际输出与期望输出之间的误差。
采用决定系数(
R²)、均方根误差(RMSE)和相对分析误差(RPD)评价模型。
R²越高、RMSE越低表示模型拟合与预测效果越好;RPD≥2.0表示模型具有较高的预测能力,1.4≤RPD<2.0表示模型预测能力一般,RPD<1.4则表明模型预测能力差,难以应用
[18]。
2 结果与分析
2.1 冬小麦叶片氮含量分析
本研究对抽穗期、灌浆前期和灌浆后期3个时期的共585份冬小麦叶片氮含量数据进行异常值剔除后,最终获取552份冬小麦叶片氮含量数据,将叶片氮含量数据按浓度梯度从小到大排序,按照3∶1的比例划分为校正集和验证集,校正集和验证集叶片氮含量的描述性统计分析如
表1所示。
总样本的氮含量为8.823~47.526 mg/g,不同处理下3个数据集的均值与变异系数(CV)均较为接近,表明不同处理下数据集的分布较为一致,且各处理下数据变化较小。从整体来看,N1、N2、N3处理下的各数据集样本量相差较小,总样本与N2处理下的均值、变异系数、最大值更为接近,不同处理下冬小麦叶片氮含量均值从大到小依次为N3>N2>N1,各处理下数据集样本的偏度小于0,数据较正态分布稍微偏向左侧,峰度小于0,峰态较为平缓,基本符合正态分布,数据的离散程度接近,各训练集与校正集CV均较相似,表明数据集划分较为合理。因此,这些数据集可以用于模型训练和验证。
2.2 不同施氮水平对冬小麦叶片氮含量的影响
由
表2可知,从抽穗期到灌浆后期,不同处理下的叶片氮含量均呈现逐渐降低的变化趋势,相同时期下的叶片氮含量随着施氮水平的增加而增加,在抽穗期和灌浆前期,N1处理下叶片氮含量显著低于N2、N3处理(
P<0.05);在灌浆后期,不同施氮水平叶片氮含量差异显著(
P<0.05),其大小顺序为N3>N2>N1。
为更直观展示冬小麦叶片氮含量在3个生长阶段不同施氮水平下的变化,绘制了冬小麦叶片氮含量箱线图如
图1所示,3个处理下的冬小麦叶片氮含量的中位线和均值均呈现抽穗期>灌浆前期>灌浆后期。随着生育时期的推进,冬小麦叶片氮含量逐渐下降。
同一时期下的冬小麦叶片氮含量在不同施氮处理下的变化趋势由小到大依次为N1< N2< N3,说明冬小麦的生长对氮肥的响应是正向的,即施氮量越大,冬小麦叶片氮含量越高。其中,抽穗期和灌浆前期在N2、N3处理下的冬小麦叶片氮含量数据中位线和均值较为接近,灌浆后期呈现较为明显的N1< N2< N3趋势,说明前2个阶段下氮肥的施用对冬小麦的氮含量影响较为均衡,冬小麦叶片对氮的吸收差异较小,氮肥的增施并未显著改变叶片的氮含量。在灌浆后期,N1处理下的叶片氮含量最低,N2适中,N3处理下的叶片氮含量最高。这表明在灌浆后期,作物对氮的需求增加,氮的积累对产量形成变得更加重要。灌浆期是作物籽粒发育的关键时期,叶片中的氮含量直接影响光合作用和营养物质的转运,进而影响籽粒的填充与发育。因此,随着灌浆的进行,冬小麦需要更多的氮来支持高强度的光合作用和营养物质合成。在这一阶段,氮肥的增施直接推动了叶片氮含量的提升。
2.3 不同施氮水平对冬小麦冠层光谱反射率的影响
选取具有代表性的冬小麦品种进行光谱反射率绘制,同一品种冬小麦在3个生长阶段不同处理下的光谱反射率如
图2所示。冬小麦冠层反射率在383~780 nm波段范围内,受不同施氮水平的影响较小,变化差异不明显。然而,在近红外光谱780~1 021 nm波段,不同施氮水平处理下冬小麦冠层反射率差异明显。这表明冬小麦冠层反射率在近红外光谱区域对施氮水平的变化更为敏感。进一步分析发现,在N3处理下,不同生长时期的冠层反射率呈现出抽穗期>灌浆前期>灌浆后期的变化趋势。
2.4 特征波段筛选
本研究采用SPA算法通过评估各波段光谱反射率数据与叶片氮含量的线性无关性,最终筛选出与叶片氮含量变化为最小共线性的波段作为不同施氮水平下冬小麦叶片氮含量的特征波段。特征波段筛选结果如
表3所示。筛选出总样本的特征波段有3个(572、926、1 016 nm),N1处理的特征波段有2个(926、1 016 nm),N2处理的特征波段有3个(720、750、944 nm),N3处理的特征波段有4个(395、703、763、944 nm)。
2.5 冬小麦叶片氮含量监测模型构建
本研究分别采用RF和BPNN构建了冬小麦叶片氮含量的监测模型,模型表现如
表4所示。当模型监测效果以验证集模型表现为准时,在不同施氮水平下,全波段下的模型表现比特征波段下的模型表现差。基于RF的叶片氮含量监测模型在校正集上的模型精度较高,而验证集的模型精度较低,表明RF模型的泛化能力较弱。3种氮肥处理和总数据集下的冬小麦叶片氮含量监测模型以SPA+BPNN模型表现最佳,N1处理下模型表现为校正集(Correction set,c)的模型精度
R2c=0.734,RMSEc=3.739 mg/g,RPDc=1.939%,验证集(Verification set,v)的模型精度
R2v=0.763,RMSEv=3.566 mg/g,RPDv=2.056%;N2处理下校正集的模型精度
R2c=0.807,RMSEc=3.173 mg/g,RPDc=2.279%,验证集的模型精度
R2v=0.802,RMSEv=3.201 mg/g,RPDv=2.246%;N3处理下校正集的模型精度
R2c=0.825,RMSEc=2.523 mg/g,RPDc=2.390%,验证集的模型精度
R2v=0.735,RMSEv=3.181 mg/g,RPDv=1.943%。总数据集校正集的模型精度
R2c=0.731,RMSEc=3.749 mg/g,RPDc=1.929%,验证集的模型精度
R2v=0.759,RMSEv=3.651 mg/g,RPDv=2.039%。综上所述,仅从验证集模型表现来看,N2处理下的SPA+BPNN模型对冬小麦叶片氮含量的监测效果最佳。
为更直观展示模型预测效果,绘制不同N处理下最优模型的散点拟合图如
图3所示,预测值与实测值分布较为均匀,靠近1∶1拟合线,SPA结合BPNN可以较好地预测冬小麦叶片氮含量。
3 结论与讨论
3.1 不同施氮水平对冬小麦叶片氮含量的影响
不同施氮水平下的冬小麦叶片氮含量变化明显。随着灌浆的进行,同一施氮水平下的冬小麦叶片氮含量呈现逐渐下降的趋势。同一时期,不同处理下的冬小麦叶片氮含量均值与中位线表现为N1<N2<N3的趋势,灌浆后期的表现尤为明显。这些变化反映出氮素在作物生长过程中的动态变化特征。灌浆后期是小麦生长的关键时期,氮素在作物体内的转移和再分配对灌浆及籽粒的发育至关重要。在这个时期,氮素的供应直接影响到籽粒的形成和产量的最终形成,这与聂胜委等
[20]对不同施肥措施下冬小麦灌浆期氮素吸收分配的研究中发现的小麦灌浆中后期氮肥供应不足的现象相符。N2、N3处理提供了足够的氮源,促进了冬小麦叶片的氮含量积累,从而有助于提高灌浆期的营养供给。而N1处理则表现出较低的叶片氮含量,可能会因为氮素不足,进而限制小麦的正常生长和灌浆过程。
从整体趋势来看,冬小麦叶片氮含量的变化与施氮水平的增加密切相关,在不同生育阶段,氮素的供给对作物生长和发育产生了显著影响。在抽穗期和灌浆前期,氮素供应对叶片氮含量的影响较小,而在灌浆后期,氮素的供应差异则对叶片氮含量产生了显著影响,进而可能影响到小麦的灌浆和最终产量,这与杨福芹等
[21]和郭燕等
[22]研究发现,冬小麦变量施肥和合理的氮素供给可以促进作物生长的研究结果一致。因此,合理的氮肥管理在不同生育阶段对于提高小麦产量和质量具有重要意义。然而,本研究的研究对象包括65种冬小麦,针对不同氮响应类型的冬小麦品种在不同施氮水平下的叶片氮含量表现未作研究,未来可从这一方面进一步研究冬小麦叶片氮含量与光谱遥感之间的关系。
3.2 特征波段对模型监测的影响
筛选包含更具辨别力的特征波段,可以使模型更加专注于有效信息,提升精度
[23]。本研究利用SPA筛选不同施氮处理下叶片氮含量特征波段,在780~1021 nm波段间,不同施氮水平下均存在与叶片氮含量相关的特征波段,但整体上的特征波段在数量及范围上存在明显差异,针对这一现象本研究尚难解释,需要进一步研究冬小麦叶片氮含量与光谱响应相关的其他因子。本研究仅采用SPA算法进行特征波段筛选,后续研究可利用不同特征选择方法进行特征波段筛选。
本研究结合RF和BPNN分别构建基于高光谱全波段和基于特征波段的冬小麦叶片氮含量监测模型,从模型构建表现结果上看,机器学习算法模型结合特征波段选择可以较好地预测冬小麦叶片氮含量,这与祝海竣
[24]利用敏感光谱作为输入量,并结合机器学习算法实现水稻叶片氮含量高精度估测结果一致。BPNN对全波段和特征波段都具有较稳定的预测能力和泛化能力,全波段光谱能够提供更多的输入特征信息,通常在建模时表现出更高的
R²和更低的RMSE。在特征选择得当的情况下,特征波段光谱也能够在模型预测中获得不错的性能,尤其是在验证集上的
RPD值有时超过全波段光谱。因此,选择合适的光谱波段并合理搭配不同的建模方法对提高监测模型的准确性至关重要,与WANG
[25]通过筛选恰当的光谱信息提高水稻茎叶氮浓度估算精度的结果相符。本研究表明,通过提取特征波段,去除无用波段可以避免噪声和无关信息影响模型表现,在提高模型精度上具有一定的应用潜力。
本研究通过无人机近红外成像高光谱技术捕捉3个关键生育时期的冬小麦反射光谱特征,结合田间实测叶片氮含量数据,探讨了不同施氮水平下的冬小麦叶片氮含量变化特征,利用SPA发掘叶片氮含量的高光谱响应特征,采用RF、BPNN等2种机器学习算法构建不同施氮水平下冬小麦叶片氮含量监测模型。结果表明,同一时期下,随着氮肥的增施,冬小麦叶片氮素含量与高光谱反射率变化明显。利用无人机高光谱遥感监测冬小麦生长状况具有较高的准确性和可靠性,在N2处理下的SPA+BPNN模型表现最佳,校正集的R2、RMSE和RPD分别为0.807、3.173 mg/g、2.279%,验证集的R2、RMSE和RPD分别为0.802、3.201 mg/g和2.246%。本研究结果为利用无人机高光谱遥感进行田块尺度的冬小麦精准施肥提供了科学依据和技术参考。