基于无人机偏振遥感的水稻冠层氮素含量反演模型

时间：2024-05-24

许童羽杨佳欣白驹驰金忠煜郭忠辉于丰华

(1.沈阳农业大学信息与电气工程学院, 沈阳 110866; 2.辽宁省农业信息化工程技术研究中心, 沈阳 110866)

0 引言

氮素是植物生长必需的营养物质之一,为植物光合作用提供重要支持,对植物生长、产量和品质有着重大影响[1]。及时、准确地判断水稻氮素营养状况,进而做出合理的施肥决策,对实现少氮肥高产量的目的尤为重要[2]。目前,无人机光谱遥感技术以其覆盖面积大、速度快、无损坏等优点被广泛应用于水稻、玉米、大豆、小麦等主要农作物的生物量检测中[3-5]。通常,水稻氮素含量可以从叶片和冠层两种尺度进行估算。冠层光谱参数测量作为作物氮素含量定量检测的关键步骤,对监测作物生长发育至关重要。光谱遥感技术通过获取冠层的光谱反射特性间接反映作物叶片的营养成分,为作物营养诊断和合理施肥提供参考[6-7]。大量研究表明,水稻冠层光谱反射率与叶片含氮量呈显著正相关,且归一化差值植被指数(NDVI)、绿波段归一化植被指数(GNVI)、比值植被指数(RVI)等植被指数都是估测水稻叶片氮素含量的有效指标[8-9]。傅友强等[10]提出美香占2号(MXZ2H)的叶片氮素含量与NDVI、比值植被指数(RVI)、差值植被指数(DVI)等呈正相关;吉丰优1002(JFY1002) 的叶片氮素含量与 RVI、EVI、改进非线性植被指数(MNLI)、 DVI、NDVI呈正相关。可见,利用植被指数进行水稻冠层氮素含量的反演具有一定可行性。

上述研究中,利用无人机平台进行遥感监测时,大多都基于垂直方向进行反射强度测量,因此只能利用来自冠层顶部的光谱信息进行氮素含量的反演。而对分蘖期的水稻,垂直方向上的测量数据中包含的水稻信息较少,且由于水稻冠层结构复杂,受水体、叶倾角、太阳高度角、传感器观测角等[11]多种因素的影响,获取的高光谱影像数据中包含镜面反射分量和漫反射分量。但镜面反射信息与植物叶片内部的物质组成无关[12]。因此,进行氮素含量反演时,有必要考虑镜面反射对反演精度的影响。LI等[13]基于偏振多光谱成像系统,提出了一种新的镜面反射去除植被指数(NSRVI)来消除植物的镜面反射干扰,提高植物健康状况检测的准确性。该方法通过构建新的植被指数来消除镜面反射对结果的影响。目前,在植物应用研究方面,通常使用偏振测量技术对镜面反射分量和漫反射分量进行区分[14]。MICHIO等[15]采用不同观测天顶角对小麦冠层进行了偏振测量,发现采用偏振度表征镜面反射分量时,能够更好地描述小麦抽穗期开始时间。吕云峰[16]对抽穗前不同生长时期的玉米冠层进行了垂直方向上的高光谱偏振观测,发现偏振光在总反射光中的比例高达10%。国内外学者将偏振遥感应用于植物研究中,并取得了一定的成果,但大部分研究是利用偏振测量进行角度效应分析,缺乏偏振测量在水稻营养元素反演方面的应用。

综上,本文在进行水稻冠层氮素含量反演时,利用无人机偏振遥感技术,将镜面反射信息从高光谱中剥离,结合多角度观测,分析不同观测天顶角下,去除镜面反射信息前后的光谱信息与水稻冠层氮素含量的相关性,进而得到无人机偏振测量的最优观测角度。基于最优观测天顶角下偏振光谱数据(即去除镜面反射信息的高光谱数据),构建适用于偏振光谱数据的植被指数——偏振光谱植被指数(Polarisation spectrum vegetation index,PSVI),以期构建一种高效、精准的植被指数,为水稻冠层氮素含量快速检测提供监测方法。

1 材料与方法

1.1 研究区概况

本研究试验区位于辽宁省海城市耿庄镇沈阳农业大学精准农业航空科研基地(40°58′45.39″N,122°43′47.006 4″E),选用水稻品种为“沈农9816”。在水稻试验区(图1)共设置5个氮肥梯度,即N0(0 kg/hm2)、N1(75 kg/hm2)、N2(150 kg/hm2)、N3(225 kg/hm2)和N4(300 kg/hm2),氮肥基追比为5∶3∶2。各个试验小区之间采用田埂进行划分。每个梯度设置3个重复小区,共设3×5=15个小区。每个小区面积为5 m×8 m=40 m2,除氮肥梯度外,每小区田间管理一致,磷钾肥采用当地标准进行施用,磷肥标准施用量为144 kg/hm2,钾肥标准施用量为192 kg/hm2,基追比为1∶1,其余田间管理同常规高产管理一致。田间采样由分蘖期至抽穗期,采样间隔7 d,每次采样在各试验小区选取具有代表性的3穴水稻进行水稻氮素含量与地上干物质量的获取,最后结果取平均值作为该小区的水稻氮素含量与地上干物质量。

图1 研究区概况Fig.1 Overview of study area

1.2 数据获取

1.2.1水稻冠层偏振光谱反射率获取

采用耐司镀膜偏振镜,卡口尺寸为37 mm,厚度为3 mm。无人机是深圳大疆创新公司制造的 M600 PRO型六旋翼无人机,搭载四川双利合谱公司的GaiaSky-mini型高光谱成像系统,将偏振片置于成像系统镜头前,偏振光谱波段为400～1 000 nm,分辨率为3 nm,单幅影像记录时间为15 s,帧速率为162 f/s,设定无人机飞行高度为100 m,选取9个观测角度(-60°、-45°、-30°、-15°、0°、15°、30°、45°和60°)进行高光谱反射测量和偏振反射测量(其中“-”对应后向散射方向,“0°”对应垂直方向,“+”对应前向散射方向),如图2所示,2次测量在0.5 h内完成。正午条件下,光照效果较好且风速小,既可以确保光照的稳定,又可以实现无人机平稳飞行,故将无人机光谱数据采集时间确定在10:00—14:00。用工具软件ENVI 5.3对获取的偏振光谱影像进行提取,通过计算每小区标签附近的平均光谱反射率来表征水稻冠层的偏振光谱反射率信息。

图2 水稻冠层偏振光谱多角度观测示意图Fig.2 Multi-angle observation diagram of rice canopy polarization spectrum

1.2.2水稻冠层氮素含量测定

在每个试验小区的取样点对水稻进行整穴的破坏性取样,带回实验室后,将每穴水稻的所有新鲜叶片剪下置于干燥箱中,在120℃下杀青60 min后, 在80℃下干燥,达到恒定质量后,称量并粉碎。采用凯氏定氮法利用研磨好的粉末检测叶片的氮素含量(mg/g)。最终共得到40组水稻冠层氮素含量的有效数据。之后,使用Kennard-Stone算法[17]按照训练集与验证集比例7∶3将样本划分。氮素含量统计结果如表1所示。

表1 训练集与验证集水稻叶片氮素含量统计结果Tab.1 Nitrogen content statistics of japonica rice leaves in training and validation sets

1.3 最优观测天顶角选择

植被指数是对地表植被状况简单、有效和经验的度量,其在全球和区域土地覆盖、植被分类和环境变化等方面都有广泛应用。此外,植被指数还可以转换成叶冠生物的物理学参数。本研究结合现有文献,运用6种与水稻冠层氮素含量估算相关的植被指数,即归一化差值植被指数[18](Normalized difference vegetation index,NDVI)、比值植被指数[19](RVI)、差值植被指数[20](Difference vegetation index,DSI)、宽动态范围植被指数[21](Wide dynamic range vegetation index,WDRVI)、再归一化植被指数[22](Renormalized vegetation index,RDVI)以及转换型植被指数[23](Conversion vegetation index,TVI)。基于不同观测天顶角下的偏振光谱数据,通过Matlab平台对各植被指数针对的波段范围进行遍历,找到各植被指数的最优波长组合,在此基础上构建植被指数,计算去除偏振信息前后,各观测角度下,植被指数与对应水稻冠层氮素含量的相关系数R。根据各天顶角下相关系数R的结果,对比得出研究所需要的最优观测角度。

1.4 水稻冠层偏振光谱反射率特征波长选择

由于采集到的偏振光谱波段具有较高维度,存在大量无用和干扰波段,而构建植被指数通常是对几个特征波长进行一定的组合而形成的一种数学表达式。因此,需要基于最优观测天顶角下的偏振光谱反射率数据,提取与水稻氮素相关的特征波长,从而通过构建植被指数反演水稻氮素含量。

本研究采用连续投影法[24](Successive projections algorithm,SPA)对最优观测角度下的偏振光谱反射率信息进行特征波长的提取。连续投影算法是一种使矢量空间共线性最小化的前向波段选择算法,它是将维度较高的光谱信息进行初步压缩,把筛选出的波段按照其对氮素含量的贡献值做进一步选择,从而剔除不敏感波长,降低复杂度[25]。连续投影算法的优势在于可以提取全波段的几个特征波长,有效降低变量间的共线性,建立冗余信息最少的波段组合,实现数据降维,大大减少建模所用波段的数量。本研究采用Matlab平台实现连续投影算法提取特征波长。

1.5 水稻冠层氮素植被指数构建

目前,大多数关于水稻氮素含量植被指数的构建并未考虑镜面反射信息的影响。而冠层光谱探测过程中,受观测角度、植被冠层结构等因素的影响,传感器接收到的能量不仅包含进入植物内部多次散射后的能量(包含植物内部的生物化学信息),还包含未进入植物内部器官,在植物叶片表面发生镜面反射的能量(只与叶片表面结构的光学特性相关)[12]。

本研究利用偏振遥感手段去除光谱信息中的镜面反射信息,采用1.4节中从最优观测天顶角下去除偏振信息的光谱数据中提取得到的特征波长来构建植被指数,构建植被指数的具体步骤如下:

(1)通过连续投影法降维得到最优观测天顶角下水稻冠层偏振光谱的特征波长反射率R1、R2、…、Rn,n为光谱波段。

(3)由于氮素含量较大,均大于1 mg/g,而Qx较小,均小于1,为进一步增强二者之间的相关性,故在Qx的基础上,均加1得Bx。

(4)为提高回归模型拟合优度,在上述步骤的基础上, 采用“平方变换”方法[26],最终得到最优观测天顶角下基于偏振光谱的植被指数(PSVI)

采用线性回归法构建水稻冠层氮素含量的反演模型,以均方根误差(RMSE)和模型决定系数(R2)作为氮素反演模型的评价标准。

2 结果与分析

2.1 最优观测天顶角选择结果

基于去除偏振信息前后的各观测天顶角测量得到的光谱数据,计算6种植被指数与水稻冠层氮素含量的相关性,为了便于对比去除偏振信息前后,各观测角度下植被指数与氮素含量相关性的提升效果,本研究采用相关系数提升比率进行直观展示,如图3所示。显然,加入偏振片会对不同观测天顶角下的植被指数与氮素含量的相关性造成不同的影响。观测天顶角为-60°、0°、30°与60°时,各植被指数与水稻冠层氮素含量的相关性并没有得到提升;观测天顶角为-45°～-15°、15°与45°时,各植被指数与水稻冠层氮素含量的相关性有明显提升,且在-15°时,提升效果最好。

图3 去除偏振信息后不同观测天顶角下各植被指数与氮素含量相关性提升结果Fig.3 Results of correlation enhancement between each vegetation index and nitrogen content at different observed zenith angles after removal of polarization information

为进一步确定最优观测天顶角,基于各观测天顶角测量得到的偏振光谱数据获得的6种植被指数与水稻冠层氮素含量相关性如表2所示(表中ρ为光谱反射率,下角数字表示特征波长)。由表2可知,去除偏振信息后,不同观测天顶角下的植被指数与水稻冠层氮素含量的相关性明显存在差异。观测天顶角的变化会对偏振光谱反射率所构建的植被指数与水稻冠层氮素含量的相关性产生影响。在水稻冠层多角度偏振光谱实测数据中,当垂直观测时,植被指数NDVI、RVI、DSI、WDRVI、RDVI、TVI与水稻氮素含量的相关性最低,且相关系数最大为0.677;当观测天顶角变大时,即从垂直观测到前向观测15°和后向观测15°时,各植被指数与氮素含量的相关性都有所提升,后向观测15°时的相关性达到最高,且相关系数最小为0.802;当观测天顶角继续增大时,即从前后向15°到前后向30°时,各植被指数与氮素含量的相关性又明显下降,但仍高于垂直观测时的相关性;当观测天顶角持续增大到前后60°时,各植被指数与氮素的相关性波动不大,但较前后30°时有所提升。综上,-15°为最优观测天顶角。因此,本研究选择观测天顶角为-15°时的偏振光谱数据进行后续植被指数的构建。

表2 基于偏振光谱数据不同观测天顶角下植被指数与水稻冠层氮素含量的相关系数Tab.2 Correlation coefficients between vegetation index and nitrogen content of rice canopy at different observed zenith angles based on polarization spectral data

2.2 偏振光谱特征波长选择结果

水稻冠层偏振光谱反射率对水稻氮素的定量反演十分重要,偏振光谱的本质是去除镜面反射信息的高光谱,而连续波段的高光谱反射率信息中存在大量的冗余信息。因此,本研究应用连续投影算法对全波段的水稻冠层偏振光谱进行特征波长筛选,并利用校正集进行内部交叉验证,当内部交叉验证时,交叉验证均方根误差(Root mean square error of cross validation,RMSECV)呈现出持续下降趋势且RMSECV达到最低,此时的RMSECV即为最优波长集合包含的波长数量,也就是该时期的最优子集波长数量。根据验证结果的RMSECV筛选出6个水稻冠层氮素含量高光谱特征波长(图4):500、566、663、691、736、763 nm。

图4 水稻冠层偏振光谱特征波段SPA筛选结果Fig.4 Successive projections algorithm (SPA) screening results of hyperspectral characteristic bands of rice canopy

2.3 偏振光谱植被指数构建结果

本研究将连续投影法筛选得到的6个特征波长的反射率两两组合,分别进行PSVI指数运算,得到PSVI指数与水稻冠层氮素含量相关系数。图5中最高时的组合为500、566 nm,相关系数绝对值为0.850。

图5 PSVI指数与水稻冠层氮素含量的相关性Fig.5 Correlation between PSVI index and nitrogen content of rice canopy

2.4 水稻氮素含量反演结果

将本研究所构建的PSVI新型植被指数作为模型的输入,利用非线性关系建立水稻冠层氮素含量的反演模型,模型训练集R2为0.783 8, RMSE为 0.428 mg/g;验证集RMSE为0.662 mg/g(图6)。

图6 水稻氮素含量反演结果Fig.6 Inversion results of rice nitrogen content

2.5 反演结果分析与评价

当前,已有众多研究人员将植被指数作为模型输入,对水稻冠层氮素含量进行反演,从IndexData-Base数据库中选取6种常用于氮素反演的植被指数[17],用来与本研究所构建的新型植被指数PSVI进行对比。将6种植被指数同样采用线性回归的方式构建基于最优观测角度下偏振光谱的水稻冠层氮素含量反演模型,如图7所示,结果显示,这6种指数所构建的反演模型的决定系数R2均小于 0.783 8,且反演精度最高的R2为0.755。综上,本研究所构建的新型植被指数PSVI线性反演效果最好。

3 讨论

目前,已有较多研究人员利用无人机遥感获取的水稻冠层光谱反射率进行氮素含量的反演建模,但大多数研究仅基于垂直观测下的光谱数据,通过数据驱动的方式,构建反演模型以提升水稻冠层氮素含量的反演精度。为了进一步提升水稻冠层氮素含量的反演精度,本研究结合多角度偏振测量技术与数据驱动方式,对无人机偏振遥感的水稻冠层氮素含量反演方法开展了研究,基于水稻冠层多角度偏振光谱数据,构建了新型植被指数PSVI定量估算水稻冠层叶片氮素含量。由于不同观测天顶角下的偏振光谱存在一定的差异,本研究结合6种常用于植物氮素含量反演的植被指数,通过分析不同观测天顶角下各植被指数与水稻冠层氮素含量的相关性,得到最优观测天顶角为后向观测15°,这一结果与HE等[27]研究结果一致。这是由于后向观测主要来自具有较高反射率的水稻冠层上方阳光照射的叶片,而前向观测包含更多来自遮阴叶片的信息,较少来自阳光照射冠层水稻叶片的信息,因此后向观测15°在去除偏振信息的同时能够更多地收集到来自水稻冠层的光谱信息,从而具有更高的相关性。此外,偏振光谱受植被冠层的几何结构、平均叶倾角等影响较大。不同观测天顶角下去除镜面反射信息后,相关性的提升效果存在一定的差异,后续还需进一步研究。

在得到最优观测天顶角下的偏振光谱数据后,为尽可能减少该偏振光谱信息中的冗余信息,减小变量间的共线性,本研究采用连续投影法得到与水稻冠层氮素含量相关性最好的6个特征波长,分别为500、566、663、691、736、763 nm。由于该算法可以在确保不同波段之间共线性最小的同时,保留原始数据中的大部分特征,因此筛选到的波长分别处于与叶片氮素含量较高的绿光和红光波段,与前期的研究结果相吻合[28-29]。本研究在构建植被指数的过程中,基于表2的相关性结果可知,比值植被指数(RVI)的整体相关性较好,因此本研究将比值植被指数的比值形式作为指数构建的基础形式,为进一步提升指数与氮素含量的相关性,利用数学变换的方法,构建了新型植被指数PSVI。该指数基于后向观测15°下的水稻冠层偏振光谱数据进行计算,将计算结果作为模型输入,采用线性回归方式进行反演,得到的PSVI植被指数具有较好的反演精度。然而在对PSVI植被指数的反演效果进行验证时,虽然该指数的反演效果优于其它6个植被指数,但该指数用于其它数据集时可能出现反演精度下降的情况,因此存在普适性不足的问题,这也是数据驱动模型的弊端,缺乏农学机理上关于氮素含量与波段之间关系的解释。未来的科研中需要对上述内容进行更深入探讨。但由于植被指数计算简单,检测装置易于开发与集成,基于植被指数的水稻冠层氮素反演方法仍具有一定的研究价值[30]。本文采用线性回归方法构建了水稻叶片氮含量的反演模型。与其它神经网络算法相比,线性回归建模方法思想简单,易于实现,建模迅速,对于小数据量、简单的关系很有效,能够实现对水稻冠层氮素含量的反演。