无人机高光谱的玉米冠层大斑病监测

梁 辉，何 敬*，雷俊杰

1.成都理工大学地球科学学院，四川 成都 610059 2.国土资源部地学空间信息技术重点实验室，四川 成都 610059

引 言

玉米种植面积约占我国全部粮食作物种植面积五分之一，玉米大斑病，是一种危害较重的常见玉米病害，该病害严重引起玉米的产量与品质下降，因此，如何通过遥感技术精准快速，高效可靠的获取玉米大斑病病情，已成为急需解决的问题。

在农业遥感中，植物在病虫害胁迫下会在光谱响应上表现出吸收和反射特性的改变，常出现可见光区域的反射率增加，红边发生“蓝移”，短波红外区域呈现较低反射率等现象[1]。Jones等[2]研究了感染叶斑病番茄叶的病情与反射光谱在紫外、可见光和近红外处的定量关系，表明使用光谱检测田间番茄细菌叶斑是可行的。诸多研究表明不同的光谱特征总能与病害间表现出积极的相关性，这些研究为植物病虫害光学遥感监测奠定了基础。

目前应用较多的是通过人工识别和地面高光谱结合的方法来研究植物病虫害，存在作业效率低，不利于在农业遥感中大面积推广使用等问题。无人机应用系统具有高效、低成本、能够获取高时空分辨率影像等特点。国内外在无人机遥感农情监测领域已开展诸多研究与探索。杨贵军等[3]利用农用无人机多传感器遥感辅助小麦育种信息获取。Juliane Bendig等[4]将基于无人机的遥感测量与植被光谱指数相结合，用于大麦的生物量监测。裴浩杰等[5]运用无人机高光谱影像对2015年小麦多生育期的长势监测。

而鲜有学者利用无人机应用系统对玉米冠层病害进行监测研究，因此设计了基于无人机高光谱的玉米冠层大斑病监测实验，利用无人机搭载高光谱仪对感病玉米冠层进行光谱分析，提取光谱响应敏感波段，研究以无人机高光谱仪为平台对玉米冠层大斑病监测的可行性与可靠性，为精准农业发展提供技术支撑。

1 实验部分

1.1 研究区概况

田间试验区位于四川省绵阳市马角镇，地处四川省东北部，地理坐标为东经105°05′08″，北纬32°05′35″。该地区以平坝和丘陵地貌为主，棕红土壤，农作物种植条件优越。为构建大斑病病情梯度，共设3块试验田(B0，B2，B4)，每块试验田总面积约300 m2，图1为B0研究区。

图1 B0研究区Fig.1 The study area B0

1.2 供试品种及接种

为保证对照实验变量的单一性，本试验选取抗大斑病等级为中感的云甜玉8号玉米作为试验对象。种植距离为行距0.8 m，株距0.3 m，并设置隔离带确保试验田不互相影响。内设3个大斑病接种水平，7月17日小喇叭口期至7月28日抽雄期共接种次数分别为B0(0次)、B2(2次)和B4(4次)，每次采用喷雾法接种浓度为1×106个·mL-1的大斑病分生孢子悬浮液。

1.3 监测系统的组成

监测系统由无人机平台和光谱仪组成。无人机平台为大疆M600 Pro型六轴无人机，有效载荷6 kg,续航时间25 min。采用D-RTK GNSS定位系统，地面平台通过GJI GS pro软件设置无人机飞行参数和航道，实现对无人机的飞行控制。高光谱遥感影像获取设备为芬兰生产的SENOP RIKOLA型高光谱成像仪。该光谱系统提供了VIS-VNIR光谱范围内的快照图像，在高度100 m时分辨率为6.5 cm，光谱范围为500～900 nm，光谱分辨率最高可达1 nm，有效波段数为380个。

1.4 影像获取与数据处理

(1)高光谱影像获取

高光谱影像于2018年8月5日11:00—13:00在试验区第一次采集，天空晴朗无风。此时玉米处于抽雄期，冠层刚出现大斑病病斑特征，随后分别在灌浆期和完熟期再次采集影像。无人机起飞之前，在地面对高光谱相机进行暗电流校正，并使用参考板进行辐射标定。无人机飞行高度为120 m，设定航速5 m·s-1。光谱仪镜头垂直向下，焦距为9 mm，曝光时间为7 ms，图2为不同病症下的玉米冠层。

(2)病情指数(DI)值调查

无人机高光谱影像采集的同时，在试验区以5株玉米为一个样点，在每个试验区均匀取样120个样点，根据GB/T23391.1—2009《玉米大、小斑病和玉米螟防治技术规范第1部分：玉米大斑病》[6]指导，结合玉米冠层实际病害程度，利用病情指数计算公式得出样点平均DI值。

(1)

式(1)中，x为各梯度病情的级值，n是最高梯度值9，f为各梯度的叶片数。

(3)高光谱影像处理

将RIKOLA采集的影像用自带软件Rikola HSI进行辐射校正和处理，在ENVI 5.3软件中对处理后的影像进行解译，根据地面采集样点在影像中对应位置构建兴趣区(ROI)，以ROI平均光谱反射率作为该样点光谱反射率，取各试验区90个样点作为建模样本，剩下30个样点作为检验样本。

1.5 光谱特征分析与参数选择

B0由于未接种且隔离较好，所以几乎未感病，将其设为正常组，B2和B4根据感病程度不同，将其分为中度组与严重组。如图3所示，光谱曲线呈现明显的“峰和谷”特征。在可见光波段(500～721 nm)内，冠层光谱反射率主要受叶片叶绿素含量的影响，健康玉米冠层单位面积上叶绿素含量高，对光吸收的多，因而反射率较低。大约在553 nm绿波段处有一反射峰-绿峰，在红光(600～680 nm)处有一个吸收带，这是因为叶绿素主要吸收红光反射绿光；在680 nm附近，有一红光吸收谷-红谷，并具有很低的反射率；在近红外区(770～900 nm)，冠层叶片光谱反射率主要受叶片内部物质及细胞结构影响，反射率较高，是近红外波段的高反射区。其中在680～756 nm波段有一反射的“陡坡”，因为植被中的叶绿素会吸收大部分的可见光，而细胞的反射机制类似小型的角反射器，对于波长700 nm以上电磁波是几乎透明的，此时植被细胞结构是影响电磁波反射的主要因素，因此植被在680～756 nm的反射率可从5%快速变化为50%。

图3 玉米冠层不同生育期光谱反射曲线(a)：抽雄期；(b)：灌浆期；(c)：完熟期Fig.3 Spectral reflection curves of corn canopy at different growth stages(a)：Tasseling stage;(b)：Filling stage; (c)：Full ripeness period

玉米冠层3个不同生育期的一阶微分光谱见图4。从图4可见，红边区域(680～756 nm)内的一阶微分变化最大，总体趋势是随病情指数增大，红边内一阶微分值逐渐减小，特别是在728 nm峰值附近表现出明显规律性，同样在绿峰532 nm附近也有相似现象。不同病情的玉米冠层一阶微分光谱形状相似，但呈梯度显现，随病情增加，红边位置逐渐向左(短波方向)移动，即红边斜率逐渐减小，红边位置发生“蓝移”，病害越重，蓝移越明显。

图4 玉米冠层不同生育期一阶微分光谱反射曲线(a)：抽雄期；(b)：灌浆期；(c)：完熟期Fig.4 First-order differential spectral reflection curves of corn canopy at different growth stages(a)：Tasseling stage;(b)：Filling stage; (c)：Full ripeness period

RIKOLA数据中包含丰富的光谱参数，将样点各波段的反射率在SPSS软件中进行线性或非线性组合计算，可以得到相应的光谱参数。结合玉米冠层的光谱特征，所选择的植被指数包括:仅利用绿光波段、绿光和近红外波段、红光和近红外波段、红边构造的光谱参数。同时考虑到病虫害探测植被指数应用的广泛性和实用性，以遥感研究中13种常见植被病害光谱参数为基础，根据玉米冠层光谱特征对参数进行改进，用于玉米冠层大斑病监测模型的构建，表1为改进后所用光谱参数。

表1 病害相关指数Table 1 Disease related index

注：RNir，RRe，RG及Rλi分别表示光谱在近红外、红、绿、及波长λ处波段反射率的均值；SDG，SDR及SDλi分别为绿光、红边及波长λ处波段一阶微分光谱值

Note:RNir,RRe,RGandRλirepresent the mean reflectance of the spectrum at near-infrared,red,green and wavelengthλ；SDG,SDRand SDλiare the first-order differential spectral values of green,red edge and wavelengthλ

2 结果与讨论

2.1 玉米冠层DI与光谱参数相关性分析

将试验区各生育期建模样点DI值与对应影像的光谱参数值进行相关性分析，可了解不同波段对玉米冠层DI的敏感性。结果表明，上述13个光谱参数除PRI(Rλi,Rλj)与OSAVI外，其余在所有生育期均与玉米冠层DI呈显著相关关系，R2达到0.6以上。其中，抽雄期光谱参数NDVI(SDλi,SDλj)，RERDVI，灌浆期光谱参数GNDVI(Rλi,Rλj)，NDVI(SDλi,SDλj)，SDR/SDG，RERDVI和完熟期光谱参数(NDVI(SDλi,SDλj)，SDR/SDG，RERVI(RNir,RRe))与玉米冠层DI的R2达到0.8以上，呈极显著相关关系，结果如表2所示。

根据表2相关性分析得知，R2达到0.8以上的光谱参数在灌浆期最多，其次是完熟期，最后是抽雄期。原因可能是玉米植株在抽雄期刚开始出现大斑病病症，叶片上病斑范围与数量较小，此时无人机高光谱仪区分正常叶片与感病叶片能力较弱，而灌浆期是玉米植株发育旺盛期，此时大斑病病斑扩散较快，病症表现明显，易于无人机光谱仪识别，发育至完熟期时，正常不感病植株叶片逐渐干枯，其生理状态与感病植株相似，因此无人机光谱仪对感病冠层识别能力下降。通过对13组光谱参数分析结果对比发现，近红外波段(770～818 nm)参与的光谱参数比绿光波段(500～556 nm)参与的光谱参数有更高的相关性，说明近红外波段对玉米冠层病情变化更敏感，解释了仅绿光波段构成的光谱参数PRI(Rλi,Rλj)在不同生育期相关性均较低的原因，在削弱土壤和环境影响方面，一阶微分类参数(SDR，SDG，NDVI(SDλi,SDλj)，SDR/SDG)比土壤调节类参数(OSAVI，GSAVI)体现了更高的相关性，说明一阶微分类参数更适用于监测玉米冠层DI。

表2 玉米冠层DI与光谱参数相关性Table 1 Correlation between DI and spectral parameters of corn canopy

2.2 监测模型精度分析

综上所述，选取各生育期R2达到0.8以上的光谱参数用于玉米冠层大斑病DI值监测模型的构建，将各生育期模型(抽雄期：DI-NDVI(SDλi,SDλj)，DI-RERDVI；灌浆期：DI-GNDVI(Rλi,Rλj)，DI-NDVI(SDλi,SDλj)，DI- SDR/SDG，DI-RERDVI；完熟期：DI-NDVI(SDλi,SDλj)，DI- SDR/SDG，DI-RERVI(RNir,RRe))设为自变量，玉米冠层DI值设为因变量，将监测模型应用于3块试验田各生育期90个检验样点，并在SPSS软件中分别进行DI实测值和预测值相关性分析，对其预测的可靠性和稳定性进行检验，采用预测精度(R2)、均方根差(RMSE)和平均相对误差(RE)3项统计指标，R2越接近1，RMSE和RE越小，表明模型精度越高，分析结果如图5、表3所示。

根据检验结果得知，抽雄期模型DI-NDVI(SDλi,SDλj)的回归斜率(0.829 3)和决定系数(R2=0.842 7)都最接近1，均方根误差(RMSE=4.59)和相对误差(RE=12.3)更小，说明模型DI-NDVI(SDλi,SDλj)的预测能力和精度更高；在灌浆期和完熟期，模型DI-RERDVI和DI-NDVI(SDλi,SDλj)有更高的预测能力和精度，可分别作为各自生育期的监测模型。

2.3 玉米冠层DI高光谱影像反演

使用各生育期对应监测模型对研究区不同生育期的高光谱影像像元进行解算，得到试验地不同生育期的玉米冠层大斑病DI值反演估算图，图6为B2研究区部分反演估算图。为了对估算图结果进行检验，根据影像在已有样点中随机选取30个，将DI实测值和估算值进行拟合分析(图7)，在RMSE最小、R2最接近1的情况下，拟合结果越接近1∶1线(图中虚线)，表明估测结果越精确。结果可以看出，3个监测模型各项检验指标均体现出较高相关性，证明基于无人机高光谱对玉米冠层大斑病监测是切实可行的。

图5 玉米冠层DI的实测值与估测值相关性分析(a)：抽雄期DI-NDVI(SDλi,SDλj)；(b)：抽雄期DI-RERDVI；(c)：灌浆期GNDVI(Rλi,Rλj)； (d)：灌浆期DI-NDVI(SDλi,SDλj)；(e):灌浆期DI-SDR/SDG;(f):灌浆期DI-RERDVI; (g):完熟期DI-NDVI(SDλi,SDλj);(h):完熟期DI-SDR/SDG;(i):完熟期DI-RERVI(RNir,RRe)Fig.5 Correlation analysis of DI measured value and estimated value of corn canopy(a)：Tasseling stage DI-NDVI(SDλi,SDλj)；(b)：Tasseling stage DI-RERDVI；(c)：Filling stage GNDVI(Rλi,Rλj)；(d)：Filling stage DI-NDVI(SDλi,SDλj)；(e):Filling stage DI-SDR/SDG;(f):Filling stage DI-RERDVI;(g):Full ripeness period DI-NDVI(SDλi,SDλj);(h):Full ripeness period DI-SDR/SDG;(i):Full ripeness period DI-RERVI(RNir,RRe)

表3 玉米冠层DI值估算模型精度检验Table 3 Accuracy test of inversion models for corn canopy DI values

3 结 论

对无人机高光谱监测玉米冠层大斑病病情指数的方法进行了研究，利用其实效性强的特点分析了玉米冠层多生育期的病情，结合高光谱仪的高分辨率从不同梯度病害胁迫入手，展现了由此导致的光谱反射率、一阶微分光谱及病害参数在反映DI值方面的差异，对比了具有代表性的13种病害参数在预测DI值方面的能力，并在小尺度上分析了监测模型对监测结果的影响。通过精度检验和分析比较，结果表明：无人机高光谱遥感技术在监测玉米冠层大斑病病情方面具有准确度高，快速便捷，时效性强等诸多优势。

图6 不同生育期玉米冠层DI值反演估测图(a)：抽雄期；(b)：灌浆期；(c)：完熟期Fig.6 Distribution maps of different corn canopy DI values at different growth stages(a):Tasseling stage;(b) Filling stage;(c) Full ripeness period

图7 不同模型DI预测值与实测值拟合结果(a)：DI-NDVI(SDλi,SDλj)；(b)：DI-RERDVI；(c)：DI-NDVI(SDλi,SDλj)Fig.7 Regression results of predicted and measured corn canopy DI values(a):DI-NDVI(SDλi,SDλj);(b):DI-RERDVI;(c):DI-NDVI(SDλi,SDλj)

研究结果表明，将无人机遥感监测玉米冠层大斑病病情的技术运用到实际生产或科学研究，能为小尺度空间范围提供高效快捷、精确可靠的玉米冠层大斑病监测服务；与前人的工作研究相比，本试验进行了多生育期下不同梯度病害胁迫的研究，得出了不同生育期最适用的病情监测指数，证明了利用无人机高光谱技术监测玉米冠层大斑病病情是切实可行的。虽然在各生育期均取得较高监测精度，但是玉米冠层大斑病一般从抽雄期开始发病，为了尽早控制病情，应在抽雄期监测。建议在实际生产应用或科学研究时，可以考虑使用红光及近红外的多光谱仪或更高分辨率的高光谱仪来提高监测精度和效率。