群体小麦条锈病发病动态无人机遥感监测方法

    Analysis for stripe rust dynamics in wheat population using UAV remote sensing

    • 摘要: 针对当前育种群体小麦条锈病表型分析手段单一、效率低下等问题,该研究提出了一种基于无人机低空遥感和多光谱成像技术的群体小麦田间条锈病高通量表型动态分析方法。该方法利用无人机采集自然发病的育种群体小麦(共600个样本,516个基因型)冠层多时相的光谱图像,并提取22个植被指数作为后续分析的表型,同时按照发病后的时间顺序与传统条锈病人工鉴定标准记录条锈病发病阶段和发病严重度数据;使用随机森林算法建立22个光谱植被指数同条锈病发病阶段与病害严重度的分类模型,并筛选出对上述两个分类问题敏感的植被指数;同时,使用随机蛙跳算法对特征进行筛选以降低仅使用随机森林算法对特征筛选的偶然性,并将随机蛙跳算法给出被选择概率排名在约前1/3的特征作为SVM算法的输入,构建发病阶段与病害严重度模型以验证随机蛙跳算法对特征筛选的有效性;综合两次特征选择的结果,分别筛选出对发病阶段和病害严重度敏感的3个植被指数,并基于这些指数响应的时间序列分析了群体中6个参考品种发病动态的差异。对条锈病发病阶段的分类模型构建中,随机森林和SVM模型测试集的F1分数分别为0.970和0.985;对条锈病严重度等级分类中,二者的F1分数为0.741和0.780,表明通过所建立的模型可以实现对群体小麦发病阶段和病害严重度等级的分类,且随机森林算法和随机蛙跳算法都能够筛选出对条锈病发病阶段和病害严重度敏感的特征。筛选出的差分植被红边指数(Difference Vegetation Index - Rededge,DVIRE)的响应对病害胁迫较为敏感,可用于同时描述田间条锈病发病阶段和严重度。该研究提出的高通量表型分析方法,基于无人机成像光谱提取的植被指数对群体小麦田间条锈病进行时间序列动态分析,能够精准量化群体小麦受条锈病胁迫状态,并可为其他作物抗病育种的表型分析提供一定的参考。

       

      Abstract: Breeding for wheat stripe rust-resistance could effectively reduce the harm caused by stripe rust. Accurate and efficient phenotyping is an essential step in crop breeding. Furthermore, UAV multispectral imaging had been used for crops growth monitoring, pathological phenotyping, biochemical index estimation and yield prediction widely for it could overcome the limitation of traditional phenotyping methods which has low efficiency and accuracy. In order to solve the problems of poor means and low efficiency of phenotyping of stripe rust in wheat, a dynamic high-throughput phenotyping framework for field stripe rust in population wheat using the UAV multispectral imaging was proposed. The experiment was carried out in Caoxinzhuang experimental farm, Yangling, Shaanxi Province, China. The breeding wheat population with 516 genotypes was used as experimental materials and planted into 600 plots according to the augmented design. Using UAV low altitude remote sensing and multispectral imaging, canopy multispectral images of experimental population wheat infected with stripe rust naturally were collected. At the same time, the data of disease stage and disease severity of stripe rust were recorded according to the disease time and the traditional manual identification standard of stripe rust. After image mosaic, radiometric correction, geometric correction, segmentation of region of interest and background, image clipping and other operations are carried out on the UAV remote sensing raw data, 22 vegetation indices were calculated to analyze. Using 22 spectral vegetation indices as classification features, we firstly established classification models of disease stage and disease severity in stripe rust, and screen out the vegetation indices sensitive to stripe rust disease stage and disease severity through random forest algorithm. At the same time, the random frog algorithm is also used to filter the original features to reduce the contingency caused by using the random forest algorithm to realize the feature screening merely. Using vegetation index with the selected probability ranking of about the top 1/3 according to the random frog algorithm as the input feature, the Support Vector Machine (SVM) algorithm was used to establish the classification model of disease stage and disease severity level in stripe rust, so as to verify the effectiveness of the features screened by the random frog algorithm. Combining the results of two feature selections, three vegetation indices sensitive to disease stage and disease severity were selected, and the differences of disease dynamics of six reference varieties in the population were analyzed based on the time series of response of these indices. The results showed that 1) for disease stages classification of stripe rust, the F1-score of the random forest classifier based on 22 original features and the SVM classifier based on the features sensitive to disease stage selected by random frog were 0.970 and 0.985 respectively, which proved the sensitive features selected by both random forest and random frog are reasonable and the two machine learning classifiers can be used in the classification of the disease stage of stripe rust in wheat population; 2) for disease severity classification of stripe rust, the F1-score of the random forest classifier based on 22 original features and the SVM classifier based on the sensitive features selected by random frog were 0.741 and 0.780 respectively. Although they are lower than the classifier for the disease stage of stripe rust, they have also achieved relatively satisfactory results, which also proved the features sensitive to disease severity selected by both random forest and random frog are reasonable and the two machine learning classifiers can be used in the classification of the disease severity of stripe rust in wheat population; 3) the response of Difference Vegetation Index RedEdge (DVIRE), Normalized Difference Vegetation Index (NDVI) and Normalized Difference Vegetation Index rededge (NDVIrededge) were more sensitive to both disease stage and disease severity in classification models. Especially, DVIRE was selected as the most sensitive feature to disease stress, and can describe both disease stage and severity of stripe rust in field meanwhile. So, dynamic phenotyping by using multitemporal vegetation indices extracted from UAV multispectral images may accurately and finely quantify stripe rust stress state in wheat population and could be considered as a reliable high-throughput method. This study proposed a new method for the field phenotype investigation for wheat stripe rust-resistance breeding, and also provided a reference for the phenotyping of stress resistance breeding of other crops.

       

    /

    返回文章
    返回