基于无人机多光谱影像的棉花黄萎病监测

宋勇, 陈兵, 王琼, 王刚, 王静, 刘海建, 郑登科, 李金霞, 陈子杰, 孙乐鑫

PDF(46962 KB)
PDF(46962 KB)
棉花学报 ›› 2023, Vol. 35 ›› Issue (2) : 87-100. DOI: 10.11963/cs20220053
研究与进展

基于无人机多光谱影像的棉花黄萎病监测

作者信息 +

Monitoring of cotton Verticillium wilt based on unmanned aerial vehicle multispectral images

Author information +
History +

摘要

【目的】通过无人机多光谱影像监测棉田黄萎病发病情况,为棉花黄萎病的精准防控提供理论指导。【方法】通过分析黄萎病发病棉田冠层的光谱特征,筛选无人机多光谱影像识别棉花黄萎病的最佳植被指数、最佳波段组合及最佳时相,并基于筛选的最佳时相建立黄萎病不同发病程度的棉田影像图,利用平行六面体法、最大似然法及支持向量机径向基函数分类法对影像图进行分类对比和精度评价。【结果】结果表明,在710~760 nm波段,不同发病程度的棉株冠层光谱反射率均随着波长的增加明显增加;在760~950 nm波段,棉株冠层光谱反射率随着黄萎病的加重明显减小。随着黄萎病加重发生,棉株的叶片叶绿素含量、地上部鲜物质质量、地上部干物质质量、植株含水量以及叶面积指数均降低。无人机多光谱遥感识别棉花黄萎病的最佳植被指数和最佳波段组合分别是差值植被指数(difference vegetation index, DVI)和B3-5-8(对应的波长分别为550 nm、656 nm和800 nm)。8月中下旬为无人机多光谱遥感识别棉花黄萎病发生程度的最佳时相。支持向量机径向基函数分类法结合最佳波段组合B3-5-8与DVI综合影像对棉田黄萎病病情的分类精度最高(分类精度为96.64%,Kappa系数为95.61%)。棉田黄萎病发病程度的分类结果与棉株冠层光谱反射率和棉株农学参数的变化相对应,并与实地调查结果一致。【结论】利用支持向量机径向基函数分类法、最佳波段组合B3-5-8与DVI综合影像对棉田黄萎病发病情况进行分类是可行的,研究结果可为利用遥感技术监测作物类似病虫害提供理论依据。

Abstract

[Objective] The unmanned aerial vehicle (UAV) multi-spectral remote sensing technology was used to monitor the severity of cotton Verticillium wilt, which will provide theoretical guidance for the precise prevention and control of cotton Verticillium wilt. [Method] By analyzing the spectral characteristics of canopy in cotton field affected by Verticillium wilt, the best vegetation index, the best wavelength combination and the best time phase for multi-spectral identification of Verticillium wilt by UAV were selected. The images of cotton fields differing in the severity of Verticillium wilt were established based on the optimal time phase. Parallelepiped method, maximum likelihood method and support vector machine radial basis function classification method were used to classify and evaluate the accuracy of the images. [Result] The results showed that the canopy spectral reflectance of cotton plants differing with occurrence of Verticillium wilt increased obviously with the increase of wavelength at 710-760 nm, the spectral reflectance of cotton canopy decreased obviously with the aggravation of Verticillium wilt at the wavelength of 760-950 nm. With the aggravation of Verticillium wilt, the chlorophyll content of leaves, fresh mass of aerial tissue per plant, dry mass of aerial tissue per plant, plant water content and leaf area index of cotton plants were all decreased. The best vegetation index and best band combination for UAV multi-spectral remote sensing to identify cotton Verticillium wilt were difference vegetation index (DVI) and B3-5-8 (the corresponding wavelengths are 550 nm, 656 nm and 800 nm). Mid to late August was the best time to identify the occurrence degree of cotton Verticillium wilt by UAV multi-spectral remote sensing. Support vector machine radial basis function classification method, best band combination B3-5-8 and DVI integrated image had the highest classification accuracy for Verticillium wilt in cotton field (classification accuracy was 96.64%, and Kappa coefficient was 95.61%). The classification results of Verticillium wilt severity in cotton fields corresponded to the changes in canopy spectral reflectance and agronomic parameters of cotton plants, and were consistent with the field investigation results. [Conclusion] It is feasible to classify the severity of Verticillium wilt in cotton field by using support vector machine radial basis function classification method, optimal band combination B3-5-8 and DVI integrated images. The results can provide theoretical ground for monitoring similar crop pests and diseases by remote sensing technology.

关键词

棉花 / 黄萎病 / 无人机 / 多光谱影像 / 分类方法

Keywords

cotton / Verticillium wilt / unmanned aerial vehicle / multispectral image / classification method

引用本文

导出引用
宋勇 , 陈兵 , 王琼 , 王刚 , 王静 , 刘海建 , 郑登科 , 李金霞 , 陈子杰 , 孙乐鑫. 基于无人机多光谱影像的棉花黄萎病监测[J]. 棉花学报, 2023, 35(2): 87-100. https://doi.org/10.11963/cs20220053
Song Yong , Chen Bing , Wang Qiong , Wang Gang , Wang Jing , Liu Haijian , Zheng Dengke , Li Jinxia , Chen Zijie , Sun Lexin. Monitoring of cotton Verticillium wilt based on unmanned aerial vehicle multispectral images[J]. Cotton Science, 2023, 35(2): 87-100. https://doi.org/10.11963/cs20220053
黄萎病(Verticillium wilt)被普遍认为是棉花的“癌症”[1],具有爆发性强、危害面积广等特点,严重发生时会导致棉花纤维品质降低、大幅减产甚至绝收,严重影响棉花生产效益。传统的棉花黄萎病调查主要靠技术人员在田间进行,时效性差、费时费力,难以及时、大面积展开,且病害严重度的划分易受调查者主观因素干扰,会在一定程度上影响病害监测的准确度[2]。而无人机多光谱遥感监测技术具有快速、高效、无损等特点,为大面积棉田的病害监测提供了新方法。病虫害易导致植株光合作用减弱、外部形态和内部结构等发生明显变化[3],这些变化特征可造成特征波段范围内的光谱信息发生改变,表现出特有的诊断性光谱特征,这就是遥感技术监测作物病虫害的基本原理[4-5]。通常健康棉株受黄萎病病菌侵染后会出现“失水”、“褪绿”等特征,外部形态变化及内部结构改变造成棉株在特征波段范围内的光谱反射率发生明显改变,即在可见光波段的光谱值随病害加重而增大,而在近红外波段的光谱值随病害的加重而减小,“绿峰”和“红边”表现出“单峰”现象,且在红边位置发生“蓝移”现象[6],该光谱特征的变化为利用遥感技术监测棉花黄萎病提供了理论依据。
目前,关于利用遥感技术监测作物病害的研究大多以卫星遥感、近地遥感及数字影像等为主,以病害的光谱位置变量[7-8]、光谱特征吸收参量[9-10]、植被指数等为变量[11]进行分析并建立病情严重度的监测模型等。如陈兵[12]通过对受黄萎病危害的棉花叶片进行近地高光谱特征分析,建立线性回归诊断模型进行病害严重度估测,模型决定系数(coefficient of determination, R2)为0.72,均方根误差(root mean square error, RMSE)为0.75。竞霞等[10]基于偏最小二乘回归法和卫星遥感影像建立了棉花黄萎病病情严重度的定量估测模型,R2RMSE和相对误差(relative error,RE)分别为0.78、0.45和9.2%。石志钰等[13]利用红、绿、蓝(red, green, and blue, RGB)棉花冠层图像以及色调、亮度和饱和度(hue, intensity, and saturation, HIS)彩色空间变换特征与植被指数之间的关系,结合Relief-F算法构建了病害诊断模型,发现植被指数与棉花黄萎病病情严重度之间存在对数关系。高攀[14]以棉田监控视频为数据源获取多种病害信息相关的关键帧图像,建立的卷积神经网络对黄萎病棉花叶片的识别准确度较高(测试正确率为96.67%,预测正确率为88.75%)。陈兵等[15]基于TM(thematic mapper)卫星影像数据,利用平行六面体分类法将棉田划分为3类:健康棉田、轻病棉田和重病棉田,总体预测精度高于85%、Kappa系数高于90%。Yu等[16]将无人机高光谱成像数据与激光雷达数据进行融合,对早期受松梢甲虫危害的松林枝条的损伤率进行评估,精度达到80%以上。Dang等[17]利用无人机可见光遥感数据获取枯萎病侵染的萝卜、土壤和地膜区域的图像并采用线性谱聚类分割方法进行了分类,识别精度达到96%。Xavier等[18]选取了不同病虫害胁迫下的棉花多光谱影像,提取其光谱参数并构建分类模型,分类精度、Kappa系数分别为79%、51%。综上所述,前人大多通过获取卫星遥感影像、近地遥感数据以及监控视频数据等建立预测模型进行病害预测,其监测过程存在周期性长、易受主观因素干扰、时效性差、分辨率低等缺点,不能直观地对棉田病害的发生情况进行准确预测,而无人机多光谱遥感可以弥补卫星遥感周期性长、分辨率低的缺点,并且可以弥补近地遥感存在主观因素影响、时效性差的不足。
本研究以六旋翼无人机及高分辨率的多光谱成像光谱仪为监测平台,以棉花黄萎病为研究对象,通过筛选出识别病害严重度的最佳波段组合、最佳影像等,创新性地提出了1种影像分类方法可对棉花黄萎病病情严重度进行准确识别,更加直观地反映棉花黄萎病的发生情况,同时也可为利用无人机多光谱遥感技术识别其他作物类似病虫害的研究提供一定参考。

1 材料与方法

1.1 试验地概况及试验设计

于2019―2020年在新疆农垦科学院棉花研究所棉花黄萎病病圃(44°32′N,85°97′E)开展试验。试验地土壤肥力中等,土壤含有机质21.3 g·kg-1、全氮0.2%、有效磷55.5 mg·kg-1、速效钾664.1 mg·kg-1,土壤pH为7.8。
2年采用相同的试验设计,同一地块,棉花品种均为新陆早8号。播种期分别为2019年4月17日、2020年4月20日。采用膜上点播、膜下滴灌、宽窄行(66+10)cm的机采棉种植模式,株距为9.6 cm。试验地划分为3个小区,作为3个重复,每个小区面积为0.49 hm2(330 m×15 m)。棉花全生育期的总灌水量约为5 775 m3·hm-2。不施底肥,分别追施纯N、P2O5和K2O 420 kg·hm-2、210 kg·hm-2和1 655 kg·hm-2,追肥全部随水滴施。
按网格式布置地面监测点(共66个,每个小区22个监测点)。每个监测点均用全球定位系统(global positioning system, GPS)测量仪(合肥卓林电子科技有限公司)进行地理位置标记(图1)。
图1 试验区监测点分布示意图

Fig. 1 Distribution of monitoring points in the test region

Full size|PPT slide

1.2 无人机多光谱影像数据获取

以HY-6X电动六旋翼无人机(上海寰鹰航空技术有限公司)搭载12个通道的多光谱相机(Micro MCA12 Snap)作为主要监测平台,该相机配备1 280×1 024像素的互补金属氧化物半导体(complementary metal-oxide-semiconductor, CMOS)传感器,12通道(即12个波段)的镜头焦距均为9.6 mm,波长、波宽分别对应B1(470 nm,10 nm)、B2(515 nm,10 nm)、B3(550 nm,10 nm)、B4(610 nm,10 nm)、B5(656 nm,10 nm)、B6(710 nm,10 nm)、B7(760 nm,10 nm)、B8(800 nm,10 nm)、B9(830 nm,10 nm)、B10(860 nm,10 nm)、B11(900 nm,20 nm)、B12(950 nm,20 nm)。
选择晴朗无云、能见度高、无风等良好天气条件下进行无人机作业,时间为11:00-15:00。飞行高度为100 m,飞行速度为5 m·s-1,旁向重叠率为 80%,航向重叠率为80%,多光谱分辨率在10~20 nm,空间分辨率约5.5 cm。试验田无人机多光谱影像获取时间:2019年6月30日(蕾期)、7月17日(花铃期)、7月29日(花铃期)、8月17日(花铃期)、9月2日(吐絮期)、9月20日(吐絮期);2020年6月29日(蕾期)、7月15日(花铃期)、7月30日(花铃期)、8月15日(花铃期)、9月1日(吐絮期)、9月15日(吐絮期)。

1.3 病情调查及划分方法

待无人机多光谱影像数据采集完毕后,参考《棉花抗病虫性评价技术规范 第5部分:黄萎病》(GB/T 22101.5-2009)[19]病级分级标准立即调查监测点棉株表型进行黄萎病分级,统计病情指数。然后根据病情指数大小划分病害严重度(b0~b4,分别对应0~4病害等级)(详见表1)。同时利用相机拍照记录监测点棉花冠层发病情况。
表1 棉花黄萎病病害严重度分级标准

Table 1 Classification standard of cotton Verticillium wilt disease severity

病害等级
Disease severity level
病害严重度
Disease severity
病情指数
Disease index (DI)
b0 健康 Health 0
b1 轻度Slight 0<DI≤25
b2 中度Moderate 25<DI≤50
b3 重度Serious 50<DI≤75
b4 极严重Critical 75<DI≤100
病情指数计算公式如下:
DI=Σ(dc×nc)/(nt×4)×100
(1)
式中,DI为病情指数,dc为相应病级,nc为各病级的病株数,nt为调查总株数。

1.4 农学参数测定

农学参数测定与黄萎病发生情况调查同步进行,具体取样时间参见1.2中无人机多光谱影像获取时间。选择监测点附近的与监测点棉株长势、病害严重度一致的4株棉花。测定指标主要包括:地上部鲜物质质量、地上部干物质质量、植株含水量、叶面积指数(leaf area index, LAI)以及叶绿素含量。其中地上部鲜物质质量采用称重法进行测量;地上部干物质质量采用烘干后称重法进行测定;根据地上部鲜物质质量和地上部干物质质量计算植株含水量;采用打孔法测定LAI[9];叶绿素含量采用分光光度法[9]进行测定,测定了植株上部、中部、下部叶片的叶绿素含量取其平均值。

1.5 数据处理与分析

1.5.1 无人机多光谱影像预处理。在Pix 4D Mapper软件中进行影像拼接、叠加、辐射矫正、精度(95%以上)校正等。此外,根据传感器特点结合前人研究发现的遥感识别病害中常用的植被指数特征[20-25],本文选取了15种作为监测棉花黄萎病的植被指数(表2),利用ENVI 5.3软件提取植被指数。
表2 本研究中选取的15种植被指数

Table 2 Fifteen vegetation indexes selected in this study

植被指数
Vegetation index
公式
Formula
参考文献
Reference
归一化差值植被指数
Normalized difference vegetation index (NDVI)
(R800R656)/(R800+R656) [20]
比值植被指数
Ratio vegetation index (RVI)
R800/R656 [20]
差值植被指数
Difference vegetation index (DVI)
R800R656 [21]
重归一化差值植被指数
Re-normalized difference vegetation index (RDVI)
(R800R656)/R800+R656 [21]
绿波段归一化差值植被指数
Green band normalized difference vegetation index (GNDVI)
(R800R550)/(R800+R550) [21]
红边归一化差值植被指数
Red edge normalized difference vegetation index (RENDVI)
(R760R710)/(R760+R710) [22]
归一化差异绿度指数
Normalized difference greenness index (NDGI)
0.5×[120×(R800R550)200×(R656R550)] [22]
三角植被指数
Triangle vegetation index (TVI)
(R550R656)/(R550+R656) [23]
土壤调节植被指数
Soil adjusted vegetation index (SAVI)
1.5×(R800R656)/(R800+R656+0.5) [23]
优化土壤调节植被指数
Optimize soil adjusted vegetation index (OSAVI)
(R800R656)/(R800+R656+0.16) [24]
调整型土壤调节植被指数
Modified soil-adjusted vegetation index (MSAVI)
0.5×[(2×R800+1)(2×R800+1)28×(R800R656)] [24]
花青素反射指数
Anthocyanin reflex index (ARI)
1/R5001/R710 [25]
增强型植被指数
Enhanced vegetation index (EVI)
2.5×[(R800R656)/(R800+6×R6567.5×R550+1)] [25]
归一化差异水体指数
Normalized difference water index (NDWI)
(R950R550)/(R950+R550) [25]
水波段指数
Water band index (WBI)
R900/R950 [25]
注:R550R656R710R760R800R900R950分别表示波长550 nm、656 nm、710 nm、760 nm、800 nm、900 nm、950 nm处的棉花冠层光谱反射率。
Note: R550, R656, R710, R760, R800, R900 and R950 represent the spectral reflectance value of cotton canopy at the wavelengths of 550 nm, 656 nm, 710 nm, 760 nm, 800 nm, 900 nm, and 950 nm, respectively.
1.5.2 最佳波段组合筛选。多光谱遥感识别棉花黄萎病严重度的最佳波段组合筛选借助最佳指数因子(optimum index factor, OIF)完成。OIF[26]的计算公式如下:
OIF=S1+S2+S3R12+R13+R23
(2)
式中,S1S2S3 代表任意3个多光谱波段的灰度值标准差,R12R13R23为任意3个波段间对应的皮尔逊(Pearson)相关系数。12个波段共产生220个3波段组合,选择OIF数值较大的前15位波段组合。灰度方差分析、相关性分析、OIF提取以及最佳波段组合的筛选均在SPSS 12.0软件中完成。
1.5.3 影像分类及评价指标。在ENVI 5.3软件中,利用平行六面体分类法、最大似然分类法、支持向量机分类法(support vector machine, SVM)进行棉田黄萎病病情严重度分类。其中,利用SVM方法分类时选取径向基函数(radial basis function)分类器,模型参数设置:γ= 0. 083,惩罚系数C=100。用总体分类精度(overall accuracy)、Kappa系数对分类结果进行检验和评价。

2 结果与分析

2.1 黄萎病不同发生程度棉田冠层无人机多光谱影像特征及棉株农学参数的比较

与b0相比,b1~b4的棉株冠层光谱反射率发生变化(图2)。在可见光波段范围内(470~656 nm),470~515 nm波段内,b1~b4棉株冠层光谱反射率无明显变化,b0植株的冠层光谱反射率较低;515~550 nm波段内,b0~b4棉株的冠层光谱反射率均增加;波长为550 nm处b0~b4的棉花冠层光谱反射率均达到峰值,且光谱反射率随着发病程度的减轻而增加;550~656 nm波段内,b0~b4棉花冠层光谱反射率均呈降低趋势,且b0的冠层光谱反射率下降趋势更明显。656 nm~710 nm波段内 b0~b4棉花冠层光谱反射率均呈增大趋势,b0与b1的冠层光谱反射率变化趋势差异不明显,与b2、b3、b4的冠层光谱反射率变化趋势差异明显。在红边波段范围内(710~760 nm),b0~b4棉株的冠层光谱反射率明显增加;当波长为760 nm时,b0~b4棉株的冠层光谱反射率差异较为明显。760~800 nm波段内,b0~b4棉花冠层光谱反射率均呈增大趋势。在近红外波段范围内(800~950 nm),不同病害严重度(b0~b4)棉株冠层光谱反射率基本趋于稳定,光谱反射率表现为:b0>b1>b2>b3>b4,且b0~b4的棉花冠层光谱反射率差异较为明显。以上结果表明无人机多光谱影像特征(光谱反射率)随着棉花黄萎病发病程度的不同而发生了较大的规律性变化,可以很好地反应该病害的发生情况。
图2 黄萎病不同发生程度棉田冠层光谱反射率变化

Fig. 2 Spectral reflectance changes of cotton canopy in fields differing in the severity of Verticillium wilt

Full size|PPT slide

表3可知,随着黄萎病的加重发生,叶绿素含量、单株地上部鲜物质质量、单株地上部干物质质量、植株含水量以及LAI均呈逐渐降低的趋势。b0与b1棉株的上述5个指标值均无显著差异;b2~b4植株的叶绿素含量、单株地上部鲜物质质量、单株地上部干物质质量及LAI均显著低于b0。因此,本研究中的这5个指标也可作为棉花黄萎病发生程度的识别指标。综上,黄萎病不同发生程度(b0~b4)棉株的叶绿素含量、单株地上部鲜物质质量、单株地上部干物质质量、植株含水量以及LAI差异明显。
表3 黄萎病不同发生程度棉田棉株农学参数变化

Table 3 Changes of agronomy parameters of cotton in fields differing in the severity of Verticillium wilt

病害严重度等级
Disease severity level
叶绿素含量
Chlorophyll content/(mg·dm-2)
单株地上部鲜物质质量Fresh mass of aerial tissue per plant/g 单株地上部干物质质量
Dry mass of aerial tissue per plant/g
植株含水量
Plant water content/%
叶面积指数
Leaf area index
b0 7.95±0.02 a 323.63±29.69 a 77.70±13.36 a 76.02±1.01 a 4.79±0.35 a
b1 6.05±0.02 ab 251.74±19.81 ab 62.69±19.67 ab 74.93±1.02 a 4.09±0.28 ab
b2 4.85±0.01 b 193.20±16.57 b 51.74±14.14 b 69.26±2.98 ab 2.38±0.23 b
b3 3.30±0.01 c 111.55±30.54 bc 39.08±16.02 c 65.43±1.21 b 1.44±0.21 c
b4 2.12±0.01 cd 56.17±16.96 c 21.12±6.69 cd 61.07±1.18 bc 0.71±0.25 cd
注:数据为平均值±标准差。同列不同小写字母表示在0.05水平差异显著。
Note: Data are mean ± standard deviation. Different lowercase letters in the same column indicate significant difference at the 0.05 probability level.

2.2 无人机多光谱遥感影像识别棉花黄萎病的最佳植被指数筛选

提取监测点棉株冠层光谱反射率,构建植被指数并与黄萎病病情指数进行相关分析。结果表明,NDGI、WBI及NDWI与黄萎病病情指数存在极显著或显著正相关关系,但相关系数较小;DVI、TVI、RDVI、SAVI、MSAVI、GNDVI、ARI、OSAVI、RENDVI、NDVI、EVI和RVI与黄萎病病情指数均呈极显著负相关关系,其中DVI与黄萎病病情指数的相关性最强,相关系数达-0.86(表4)。因此,DVI可初步作为无人机多光谱影像识别棉花黄萎病发生程度的最佳植被指数。
表4 棉花黄萎病病情指数与15个植被指数的相关分析

Table 4 Correlation analysis of 15 vegetation indexes and disease index of cotton Verticillium wilt

植被指数
Vegetation index
相关系数
Correlation coefficient
植被指数
Vegetation index
相关系数
Correlation coefficient
DVI -0.86** RENDVI -0.48**
TVI -0.73** NDVI -0.37**
RDVI -0.71** EVI -0.36**
SAVI -0.70** NDGI 0.30**
MSAVI -0.68** RVI -0.27**
GNDVI -0.67** WBI 0.16**
ARI -0.56** NDWI 0.12*
OSAVI -0.55**
注:**和*分别表示极显著相关(P<0.01)和显著相关(P<0.05)。
Note: ** and * indicate extremely significant correlation (P < 0.01) and significant correlation (P < 0.05), respectively.

2.3 无人机多光谱遥感影像识别棉花黄萎病的最佳波段组合筛选

在可见光波段至近红外波段(B1~B12),12个波段下不同发病程度棉株的冠层灰度值标准差见表5。在可见光波段(B1~B5),B5的棉株冠层灰度值标准差最大,B3次之;在红边波段(B6~B7),B6的棉株冠层灰度值标准差最大;在近红外波段(B8~B12),B8的棉株冠层灰度值标准差最大。标准差越大,灰度值与均值的差值就越大,反映的图像信息量越多[26]。依据冠层灰度值的标准差分析结果结合波长划分范围,B3、B5、B6、B8分别为不同波长范围内其灰度值标准差较大的波段,其反映的图像信息量较多。
表5 12个波段下棉花冠层灰度值标准差

Table 5 Standard deviations of gray value of cotton canopy under 12 bands

波段
Band
标准差
Standard deviation
排序
Rank
波段
Band
标准差
Standard deviation
排序
Rank
B1 10.27 9 B7 18.73 3
B2 12.71 5 B8 11.41 6
B3 14.26 4 B9 7.66 11
B4 11.16 8 B10 8.56 10
B5 21.77 2 B11 11.22 7
B6 30.54 1 B12 7.56 12
基于上述波段下棉花冠层灰度值标准差的分析结果,借助OIF值进一步筛选无人机多光谱影像识别棉花黄萎病的最佳波段组合。排名前15位的波段组合的OIF值为67.83~153.44,排序前3位的波段组合分别为:B3、B5与B8(即B3-5-8),B4、B6与B8(即B4-6-8)和B4、B6与B9(即B4-6-9)(表6)。由于B3、B5、B8 这3个波段组合对应的OIF值最大,该组合的波长分别对应550 nm、656 nm、800 nm,分别代表了绿光波段、红光波段、红边波段,是植被识别的敏感波段[2],且与2.1中的分析结果一致,能较好地反映棉花黄萎病发病情况。因此,B3-5-8可作为无人机多光谱遥感监测棉花黄萎病的最佳波段组合。
表6 OIF值前15位的波段组合

Table 6 The top 15 band combinations based on OIF value

序号
Number
波段组合
Band combination
OIF
OIF value
序号
Number
波段组合
Band combination
OIF
OIF value
1 B3、B5、B8 153.44 9 B1、B6、B10 81.88
2 B4、B6、B8 132.26 10 B1、B7、B8 80.98
3 B4、B6、B9 128.28 11 B1、B6、B8 79.61
4 B4、B6、B10 109.57 12 B1、B6、B9 79.21
5 B2、B3、B8 97.48 13 B4、B5、B9 73.56
6 B4、B5、B8 91.76 14 B3、B6、B9 70.73
7 B3、B6、B8 89.27 15 B2、B8、B11 67.83
8 B3、B5、B9 83.71

2.4 基于无人机多光谱遥感监测棉花黄萎病的最佳时相选择

借鉴前人研究结果[15](本研究所采用的传感器波段信息与TM卫星遥感影像的波段信息具有相似特征),利用棉花不同生育时期的棉田多光谱原始影像结合地面监测点中黄萎病病情调查数据,采用平行六面体分类法对试验区2019年棉花黄萎病原始影像进行分类,结果如图3所示。6月30日-7月17日,试验区棉田内棉株主要表现为健康(b0)、轻度发病症状(b1);7月29日-8月17日,棉株除表现为健康(b0)、轻度症状(b1)外,中度发病棉株(b2)比例增加趋势较为明显,重度病害(b3)及极严重病害(b4)局部发生,但发病面积不明显;8月17日-9月2日,重度病害(b3)及极严重病害(b4)呈现大面积增加趋势。
图3 不同时间基于无人机多光谱影像的棉花黄萎病发病程度分类

Fig. 3 Classification of the severity of cotton Verticillium wilt based on UAV multispectral images at different times

Full size|PPT slide

根据无人机多光谱棉田原始影像分类结果及田间实际调查情况分析判断,棉田黄萎病的爆发是在8月17日以后,该阶段无人机影像特征明显,分类结果准确,验证结果一致性较好。同时,利用平行六面体分类法对试验区棉花黄萎病原始影像进行分类,发现8月17日影像分类精度(总体分类精度90.24%,Kappa系数87.29%)最高,9月2日影像的分类精度(总体分类精度86.24%,Kappa系数81.25%)次之,7月29日影像的分类精度(总体分类精度83.58%,Kappa系数78.17%)较低。9月2日后棉花陆续开始吐絮,受棉株衰老、吐絮等因素的影响,棉田内黄萎病发生程度的分类结果准确性降低,不能作为棉花黄萎病监测的最佳时相。因此,8月17日-9月2日为当地无人机多光谱遥感识别棉花黄萎病的最佳时相。

2.5 基于不同多光谱影像的棉花黄萎病发病程度分类结果及精度统计

通过筛选的最佳时相(8月17日)进行黄萎病发病程度分类,结果见图4,并进行精度评价分析,结果见表7。针对3种不同的影像类型,即12波段原始影像、最佳波段组合(B3-5-8)影像、最佳波段组合与DVI综合影像(最佳波段组合影像与DVI图组合而成的影像),利用平行六面体分类法对不同影像进行分类并评价其精度,总体分类精度及Kappa系数表现为:最佳波段组合与DVI综合影像>最佳波段组合影像>12波段原始影像。利用最大似然分类法对不同波段组合的影像进行分类及精度评价,总体分类精度及Kappa系数表现为:最佳波段组合与DVI综合影像>最佳波段组合影像>12波段原始影像。利用支持向量机径向基函数分类法对不同波段组合的影像进行分类并评价其精度,总体分类精度及Kappa系数表现为:最佳波段组合与DVI综合影像>最佳波段组合影像>12波段原始影像。
图4 基于不同光谱影像的棉花黄萎病发生程度分类
A、D、G分别为基于原始波段影像、最佳波段组合(B3-5-8)影像、最佳波段组合与DVI综合影像的平行六面体分类结果;B、E、H分别为基于原始波段影像、最佳波段组合影像、最佳波段组合与DVI综合影像的最大似然分类结果;C、F、I分别为基于原始波段影像、最佳波段组合影像、最佳波段组合与DVI综合影像的支持向量机径向基函数分类结果。

Fig. 4 Classification of cotton Verticillium wilt severity degree based on different spectral images

A, D and G are the parallelepiped classification results based on the original band images, the best band combination (B3-5-8 )images, the best combined band images and the DVI integrated images, respectively. B, E and H are the maximum likelihood classification results based on the original band images, the best band combination images, the best combined band images and the DVI integrated images, respectively. C, F and I are the support vector machine radial basis function classification results based on the original band images, the best band combination images, the best combined band images and DVI integrated images, respectively.

Full size|PPT slide

表7 不同分类模型的影像分类精度参数

Table 7 Classification accuracy parameters of images using different classification models

影像类型
Image type
平行六面体分类法
Parallelepiped classification
最大似然分类法
Maximum likehood classification
支持向量机径向基函数
分类法
Support vector machine radial basis function classification
分类精度
Classification accuracy/%
Kappa系数
Kappa coefficient/%
分类精度
Classification accuracy/%
Kappa系数
Kappa coefficient/%
分类精度
Classification accuracy/%
Kappa系数
Kappa coefficient%
原始波段影像
Original band image
90.24 87.39 88.51 85.26 61.65 53.51
最佳波段组合影像
Best band combination image
90.47 87.52 94.58 93.03 88.99 85.62
最佳波段组合与DVI综合影像
Best combined band images
and DVI integrated image
91.46 88.89 94.64 93.12 96.64 95.61
针对12波段原始影像,不同分类方法的分类精度及Kappa系数表现为:平行六面体分类法>最大似然分类法>支持向量机径向基函数分类法;针对最佳波段组合影像,不同分类方法的分类精度及Kappa系数表现为:最大似然分类法>平行六面体分类法>支持向量机径向基函数分类法;针对最佳波段组合与DVI综合影像,不同分类方法的分类精度及Kappa系数表现为:支持向量机径向基函数分类法>最大似然分类法>平行六面体分类法。3种分类方法的分类精度结果显示,基于支持向量机径向基函数分类法对最佳波段组合与DVI综合影像进行分类的总体分类精度(96.64%)及Kappa系数(95.61%)均最高。因此,基于支持向量机径向基函数分类法和最佳波段组合与DVI综合影像可以更准确地识别花铃期棉花黄萎病发病程度。

2.6 基于最佳分类影像和分类法的棉花黄萎病严重度实地验证结果

根据上述分析结果,任意选取2020年黄萎病不同发病程度的地面监测点作为验证点进行实地验证,结果如图5所示。将监测点的实拍图像与最佳分类影像的分类结果进行对比验证,结果显示基于支持向量机径向基函数分类法、最佳波段组合与DVI综合影像的分类结果与棉田实际调查选择的30个验证点的黄萎病发病情况一致,验证准确度达95%以上,符合棉田黄萎病实际发生情况。因此,支持向量机径向基函数分类法及最佳波段组合与DVI综合影像可以准确地监测试验区内棉花黄萎病的发生程度。
图5 基于最佳分类影像的棉花黄萎病发生程度分类结果实地验证

Fig. 5 Field verification of cotton Verticillium wilt occurrence degree classification results based on the optimum classification images

Full size|PPT slide

3 讨论

棉花受黄萎病病菌侵染后,会在维管束内部产生胼胝体,影响植株体内的水分运输,降低叶片蒸腾效率,且易产生有毒物质,从而使外部形态和冠层特征发生改变,典型症状表现为失绿、失水、LAI降低、冠层覆盖度下降和生物量减少等[5],该现象易引起特定波段内的光谱特征发生明显变化[3-4]。本研究通过提取不同生育时期黄萎病棉田无人机多光谱影像的光谱特征,筛选出识别黄萎病的较佳植被指数,分别是DVI、TVI、RDVI,通过对比验证后结果表明与病情指数相关性最强的植被指数——DVI的识别效果最好,其结果与陈兵等[3]利用光谱指数估测棉花黄萎病发生程度所筛选出的植被指数一致。洪帅等[27]利用高光谱数据分析比较盛蕾期与现蕾期棉花冠层叶片叶绿素含量,发现盛蕾期棉叶螨大面积发生导致棉花冠层光谱反射率降低。陈兵等[28]利用光谱红边参数建立黄萎病氮素含量诊断模型,结果显示随着病害加重,氮素含量显著减低。Wang等[29]通过获取受雹灾棉花的高光谱遥感数据并筛选相关植被指数建立了基于NDVI、RVI和DVI的产量损失预测模型,结果显示基于这3种植被指数建立的棉花产量损失预测模型效果较好,其植被指数的筛选结果与本研究的相符,可能是由于棉花受雹灾或病害侵染后其冠层光谱特征具有一定的相似性。
可见光波段(470~656 nm)的光谱反射率主要受叶绿素含量影响[9]。黄萎病发生后棉花叶片的叶绿素含量随着病害的加重逐渐降低,对蓝光和红光的吸收减弱、反射增强,同时对绿光的反射减弱。近红外波段(710~950 nm)的光谱反射率主要受叶绿素含量、植株含水量和冠层结构的影响[9]。病害引起植株内部叶绿素含量及含水量降低、生理结构遭受破坏,植株冠层结构发生相应的变化,致使近红外波段的反射率下降。本研究发现不同发病程度棉株的农学参数在光谱曲线上的光谱响应特征与前人研究结果一致。本研究采用类似方法筛选出的最佳波段组合为B3-5-8,而赵庆展等[26]运用无人机多光谱影像特征筛选出的最佳波段组合为B1-6-12,可能是由于分类目标的不同,本研究的目标主要是黄萎病棉田冠层,而赵庆展等[26]的研究目标是裸地、林地和耕地等,其光谱特征存在较大差异。在最佳时相选择方面,本研究采取的方法与陈兵等[15]的方法一致,但筛选出的最佳时相存在偏差;陈兵等[15]研究发现基于TM卫星遥感监测棉花黄萎病的最佳时相为7月中下旬至8月中旬,本研究筛选的最佳时相为8月下旬,最佳时相选择出现偏差的原因可能与棉花品种、当年的气候因素及棉田实际发病情况等有关。
Wang等[30]通过获取棉花根腐病无人机多光谱遥感影像并建立分类算法,实现了对棉花根腐病的分级分类。因此,借助无人机多光谱影像也可对棉花黄萎病发生程度进行分类。本研究利用监督分类法中的平行六面体分类法、最大似然分类法及支持向量机分类法分别对黄萎病棉田的原始波段影像、最佳波段组合影像、最佳波段组合与DVI综合影像进行了分类及精度评价。不同分类算法在3个波段和多个波段上的运行时间有差异,其中支持向量机径向基函数分类法运行时间最长,其次是最大似然分类法,最后为平行六面体分类法,但不同分类算法的运行时间差异较小,均能满足研究要求。因此,运行时间可作为病害识别的次要考虑因素。从病害识别的精度方面考虑,基于支持向量机分类法和最佳波段组合与DVI综合影像的分类精度最高,而陈兵等[15]基于TM卫星遥感影像进行棉花黄萎病严重度分类时选用的平行六面体分类法的分类精度最高。分类结果不同的原因可能是陈兵等[15]采用的影像为TM卫星原始遥感影像,而本研究采用的分类影像是最佳波段组合与DVI综合影像,2种影像在波段特征上存在一定的差异;其次,这2种影像在波段特征与分辨率方面存在较大差异,造成分类精度也存在一定的差异。朱耀辉[31]利用支持向量机分类法对预处理后的无人机高光谱小麦全蚀病影像进行分类,结果表明支持向量机的径向基函数分类法分类精度最高,可以较好地实现小麦全蚀病的实时监测,与本研究的结果一致。Xavier等[18]利用多项逻辑回归、多项增强逻辑回归、支持向量机、随机森林4种不同分类法对获取的无人机多光谱棉花枯萎病影像进行分类,结果显示基于支持向量机分类法的分类影像更符合棉田实际发病情况。因此,将支持向量机分类法、最佳波段组合与DVI综合影像相结合可以精准实现对棉田黄萎病发生程度的监测,可为病害的精准防治提供理论指导。
此外,本研究以无人机多光谱遥感为平台,以棉花黄萎病为研究对象,在棉花黄萎病病圃开展了小区试验研究,与大田示范推广相比,具有一定的局限性。由于大田试验的不确定性因素较多,未来研究过程中需要综合考虑各种因素的影响,采用多因素融合(如气象条件、土壤类型、地理特征、多种病虫胁迫并存等)试验来验证该多光谱影像分类监测模型的准确度,建立并优化棉花黄萎病发病程度分类模型。本研究筛选出了能够精准识别棉花黄萎病病情严重度的分类影像及分类方法,可为类似病害发生情况的监测提供研究思路。

4 结论

本研究条件下,差值植被指数(DVI)为无人机多光谱影像识别棉田黄萎病发病程度的最佳植被指数。B3-5-8(波长分别为550 nm、656 nm、800 nm)为无人机多光谱遥感监测棉花黄萎病的最佳波段组合。采用支持向量机径向基函数分类法基于最佳波段组合B3-5-8与DVI综合影像能够准确地对棉花黄萎病病情严重度进行分类,总体分类精度为96.64%,Kappa系数为95.61%。相关研究结果可为利用无人机多光谱影像监测作物病害发生情况提供借鉴。

参考文献

[1]
朱荷琴, 李志芳, 冯自力, 等. 我国棉花黄萎病研究十年回顾及展望[J/OL]. 棉花学报, 2017, 29(S1): 37-50[2022-09-20]. https://doi.org/10.11963/1002-7807.zhqzhq.20170825.
Zhu Heqin, Li Zhifang, Feng Zili, et al. Overview of cotton Verticillium wilt research over the past decade in China and its prospect in future[J/OL]. Cotton Science, 2017, 29(S1): 37-50[2022-09-20]. https://doi.org/10.11963/1002-7807.zhqzhq.20170825.
[2]
王姣, 李志沛, 张立福, 等. 基于棉花黄萎病多“症状”特征的植被指数构建及病情遥感监测研究[J/OL]. 地理与地理信息科学, 2019, 35(5): 46-51[2022-09-20]. https://doi.org/10.3969/j.issn.1672-0504.2019.05.008.
Wang Jiao, Li Zhipei, Zhang Lifu, et al. Study of cotton Vertici-llium wilt: construction of a vegetation index based on multiple "symptoms" characteristics and remote sensing monitoring[J/OL]. Geography and Geo-Information Science, 2019, 35(5): 46-51[2022-09-20]. https://doi.org/10.3969/j.issn.1672-0504.2019.05.008.
[3]
陈兵, 李少昆, 王克如, 等. 基于TM影像光谱指数的棉花病害严重度估测[J]. 红外与毫米波学报, 2011, 30(5): 451-457.
Chen Bing, Li Shaokun, Wang Keru, et al. Estimating severity level of cotton disease based on spectral indices of TM image[J]. Journal of Infrared and Millimeter Waves, 2011, 30(5): 451-457.
[4]
Korolev N, Perez-Artes E, Mercado-Blanco J, et al. Vegetative compatibility of cotton defoliating Verticillium dahliae in Israel and its pathogenicity to various crop plants[J/OL]. European Journal of Plant Pathology, 2008, 122(4): 603-617[2022-09-20]. https://doi.org/10.1007/s10658-008-9330-1.
[5]
陈巧玲, 范燕敏, 武红旗, 等. 棉花黄萎病冠层光谱特征分析及其严重度估测模型的建立[J/OL]. 新疆农业大学学报, 2020, 43(4): 261-269[2022-09-20]. https://doi.org/10.3969/j.issn.1007-8614.2020.04.004.
Chen Qiaoling, Fan Yanmin, Wu Hongqi, et al. Spectral characteristics analysis of cotton Verticillium wilt canopy and establishment of its severity estimation model[J/OL]. Journal of Xinjiang Agricultural University, 2020, 43(4): 261-269[2022-09-20]. https://doi.org/10.3969/j.issn.1007-8614.2020.04.004.
[6]
陈兵, 王琼, 肖春华, 等. 9种胁迫条件下棉花叶片光谱特征及其识别研究[J/OL]. 西北农业学报, 2015, 24(10): 64-73[2022-09-20]. https://doi:10.7606/j.issn.1004-1389.2015.10.010.
Chen Bing, Wang Qiong, Xiao Chunhua, et al. Exploration and recognition of spectral characteristics of cotton leaf suffered nine stresses[J/OL]. Acta Agriculturae Boreali-occidentalis Sinica, 2015, 24(10): 64-73[2022-09-20]. https://doi:10.7606/j.issn.1004-1389.2015.10.010.
[7]
陈兵, 李少昆, 王克如, 等. 棉花黄萎病病叶光谱特征与病情严重度的估测[J]. 中国农业科学, 2007, 40(12): 2709-2715.
Chen Bing, Li Shaokun, Wang Keru, et al. Spectrum characteristics of cotton single leaf infected by Vericillium wilt and estimation on severity level of disease[J]. Scientia Agricultura Sinica, 2007, 40(12): 2709-2715.
[8]
竞霞, 王纪华, 宋晓宇, 等. 棉花黄萎病病情严重度的连续统去除估测法[J/OL]. 农业工程学报, 2010, 26(1): 193-198[2022-09-20]. https://doi.org/10.3969/j.issn.1002-6819.2010.01.034.
Jing Xia, Wang Jihua, Song Xiaoyu, et al. Continuum removal method for cotton Verticillium wilt severity monitoring with hyperspectral data[J/OL]. Transactions of the Chinese Society of Agricultural Engineering, 2010, 26(1): 193-198[2022-09-20]. https://doi.org/10.3969/j.issn.1002-6819.2010.01.034.
[9]
陈兵. 基于多平台棉花黄萎病的遥感监测研究[D]. 石河子: 石河子大学, 2010.
Chen Bing. Study on monitoring cotton infected with Verticil-lium wilt based on multi-platform remote sensing[D]. Shihezi: Shihezi University, 2010.
[10]
竞霞, 黄文江, 琚存勇, 等. 基于PLS算法的棉花黄萎病高空间分辨率遥感监测[J/OL]. 农业工程学报, 2010, 26(8): 229-235[2022-09-20]. https://doi.org/10.3969/j.issn.1002-6819.2010.08.039.
Jing Xia, Huang Wenjiang, Ju Cunyong, et al. Remote sensing monitoring severity level of cotton Verticillium wilt based on partial least squares regressive analysis[J/OL]. Transactions of the Chinese Society of Agricultural Engineering, 2010, 26(8): 229-235[2022-09-20]. https://doi.org/10.3969/j.issn.1002-6819.2010.08.039.
[11]
Zheng Q, Huang W J, Cui X M, et al. Identification of wheat yellow rust using optimal three band spectral indices in diffe-rent growth stages[J/OL]. Sensors, 2018, 19(1): 35[2022-09-20]. https://doi.org/10.3390/s19010035.
[12]
陈兵. 棉花黄萎病高光谱特征及其光谱诊断模型研究[D]. 石河子: 石河子大学, 2007.
Chen Bing. Study on high-spectral characteristics of cotton infected with Verticillium wilt and diagnosing models by spectrum[D]. Shihezi: Shihezi University, 2007.
[13]
石志钰, 裴雅琨, 朱玉涛, 等. 基于数字图像棉花黄萎病诊断与防治的管理远程监测[J/OL]. 新疆农业科学, 2020, 57(6): 1166-1174[2022-09-20]. https://doi.org/10.6048/j.issn.1001-4330.2020.06.022.
Shi Zhiyu, Pei Yakun, Zhu Yutao, et al. Study on remote monitoring of cotton Verticillium wilt diagnosis and control ma-nagement based on digital image[J/OL]. Xinjiang Agricultural Sciences, 2020, 57(6): 1166-1174[2022-09-20]. https://doi.org/10.6048/j.issn.1001-4330.2020.06.022.
[14]
高攀. 基于图像识别的棉花黄萎病实时诊断方法研究[D]. 石河子: 石河子大学, 2019.
Gao Pan. Study on real-time diagnosis of Verticillium wilt of cotton based on image recognition[D]. Shihezi: Shihezi University, 2019.
[15]
陈兵, 王克如, 李少昆, 等. 棉花黄萎病疑似病田的卫星遥感监测——以TM卫星影像为例[J/OL]. 作物学报, 2012, 38(1): 129-139[2022-09-20]. https://doi.org/10.3724/SP.J.1006.2012.00129.
Chen Bing, Wang Keru, Li Shaokun, et al. Monitoring cotton field with suspected Verticillium wilt using satellite remote sensing with TM satellite image as an example[J/OL]. Acta Agronomica Sinica, 2012, 38(1): 129-139[2022-09-20]. https://doi.org/10.3724/SP.J.1006.2012.00129.
[16]
Yu R, Luo Y, Zhou Q, et al. A machine learning algorithm to detect pine wilt disease using UAV based hyperspectral imagery and LiDAR data at the tree level[J/OL]. International Journal of Applied Earth Observation and Geo-information, 2021, 101(1): 102363[2022-09-20]. https://doi.org/10.1016/j.jag.2021.102363.
[17]
Dang L, Wang H, Li Y, et al. Fusarium wilt of radish detection using RGB and near infrared images from unmanned aerial vehicles[J/OL]. Remote Sensing, 2020, 12 (17): 2863[2022-09-20]. https://doi.org/10.3390/rs12172863.
[18]
Xavier T W F, Souto R N V, Statella T, et al. Identification of Ramularia leaf blight cotton disease infection levels by multispectral, multiscale UAV imagery[J/OL]. Drones, 2019, 3(2): 33[2022-09-20]. https://doi.org/10.3390/drones3020033.
[19]
中华人民共和国农业部. 棉花抗病虫性评价技术规范第5部分: 黄萎病 GB/T 22101.5-2009[S]. 北京: 中国标准出版社, 2009.
Ministry of Agriculture of the People's Republic of China. Technical specification for evaluating resistance of cotton to disease and insect pests-Part 5: Verticillium wilt GB/T 22101.5-2009[S]. Beijing: China Standards Press, 2009.
[20]
陈兵, 李少昆, 王克如, 等. 病害胁迫下棉花叶片色素含量高光谱遥感估测研究[J/OL]. 光谱学与光谱分析, 2010, 30(2): 421-425[2022-09-20]. https://doi.org/10.3964/j.issn.1000-0593(2010)02-0421-05.
Chen Bing, Li Shaokun, Wang Keru, et al. Study on hyperspectral estimation of pigment contents in leaves of cotton under disease stress[J/OL]. Spectroscopy and Spectral Analysis, 2010, 30(2): 421-425[2022-09-20]. https://doi.org/10.3964/j.issn.1000-0593(2010)02-0421-05.
[21]
田庆久, 闵祥军. 植被指数研究进展[J/OL]. 地球科学进展, 1998(4):10-16[2022-09-20]. https://doi.org/10.3321/j.issn:1001-8166.1998.04.002.
Tian Qingjiu, Min Xiangjun. Research progress on vegetation index[J/OL]. Advances in Earth Science, 1998(4): 10-16[2022-09-20]. https://doi.org/10.3321/j.issn:1001-8166.1998.04.002.
[22]
Mirik M, Ansley R J I, Michels G J, et al. Spectral vegetation indices selected for quantifying Russian wheat aphid (Diuraphis noxia) feeding damage in wheat (Triticum aestivum L.)[J/OL]. Precision Agriculture, 2012, 13(4): 501-516[2022-09-20]. https://doi.org/10.1007/s11119-012-9264-7.
[23]
Penuelas J, Filella I. Visible and near-infrared reflectance techniques for diagnosing plant physiological status[J/OL]. Trends in Plant Science, 1998, 3(4): 151-156[2022-09-20]. https://doi.org/10.1016/S1360-1385(98)01213-8.
[24]
Phadikar S, Sil J, Das A K. Vegetative indices and edge texture based shadow elimination method for rice plant images[C]//In-ternational Conference on Radar, Communication and Compu-ting (ICRCC). New York: Institute Electrical and Electronics Engineers, 2012: 329-333.
[25]
Naidu R A, Perry E M, Pierce F J, et al. The potential of spectral reflectance technique for the detection of grapevine leaf roll associated virus-3 in two red berried wine grape cultivators[J/OL]. Computers and Electronics in Agriculture, 2009, 66(1): 38-45[2022-09-20]. https://doi.org/10.1016/j.compag.2008.11.007.
[26]
赵庆展, 刘伟, 尹小君, 等. 基于无人机多光谱影像特征的最佳波段组合研究[J/OL]. 农业机械学报, 2016, 47(3): 242-248, 291[2022-09-20]. https://doi.org/10.6041/j.issn.1000-1298.2016.03.034.
Zhao Qingzhan, Liu Wei, Yin Xiaojun, et al. Selection of optimum bands combination based on multi-spectral images of UAV[J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2016, 47(3): 242-248, 291[2022-09-20]. https://doi.org/10.6041/j.issn.1000-1298.2016.03.034.
[27]
洪帅, 张泽, 张立福, 等. 滴灌棉花不同生育时期冠层叶片叶绿素含量的高光谱估测模型[J/OL]. 棉花学报, 2019, 31(2): 138-146[2022-09-20]. https://doi.org/10.11963/1002-7807.hslx.20190319.
Hong Shuai, Zhang Ze, Zhang Lifu, et al. Hyperspectral estimation model of chlorophyll content in cotton canopy leaves under drip irrigation at different growth stages[J/OL]. Cotton Science, 2019, 31(2): 138-146[2022-09-20]. https://doi.org/10.11963/1002-7807.hslx.20190319.
[28]
陈兵, 王方永, 韩焕勇, 等. 基于光谱红边参数的棉花黄萎病叶片氮素含量诊断研究[J/OL]. 棉花学报, 2013, 25(3): 254-261[2022-09-20]. https://doi.org/10.11963/cs130310.
Chen Bing, Wang Fangyong, Han Huanyong, et al. Monitoring nitrogen contents in leaves of cotton under Verticillium wilt stress based on spectra red-edge parameters[J/OL]. Cotton Science, 2013, 25(3): 254-261[2022-09-20]. https://doi.org/10.11963/cs130310.
[29]
Wang L, Liu Y, Wen M, et al. Using field hyperspectral data to predict cotton yield reduction after hail damage[J/OL]. Computers and Electronics in Agriculture, 2020, 190(11): 1-9[2022-09-20]. https://doi.org/10.1016/J.COMPAG.2021.106400.
[30]
Wang T, Thomasson J A, Isakeit T, et al. A plant-by-plant method to identify and treat cotton root rot based on UAV remote sensing[J/OL]. Remote Sensing, 2020, 12(15): 2453[2022-09-20]. https://doi.org/10.3390/rs12152453.
[31]
朱耀辉. 基于无人机成像高光谱的冬小麦全蚀病等级监测[D]. 郑州: 河南农业大学, 2018.
Zhu Yaohui. Monitoring and classification of wheat take all in field based on UAV hyper-spectral image[D]. Zhengzhou: Henan Agricultural University, 2018.

基金

国家自然科学基金(41961054)
新疆农垦科学院智慧棉田创新团队(NCG202304)
新疆生产建设兵团英才项目(BTYC202001)
PDF(46962 KB)

622

Accesses

0

Citation

Detail

段落导航
相关文章

/