MDBPSO算法优化的全波段光谱数据协同冠层SIF监测小麦条锈病

白宗璠 竞 霞,* 张 腾 董莹莹

MDBPSO算法优化的全波段光谱数据协同冠层SIF监测小麦条锈病

白宗璠1竞 霞1,*张 腾1董莹莹2

1西安科技大学测绘科学与技术学院, 陕西西安 710054;2中国科学院遥感与数字地球研究所数字地球重点实验室, 北京 100094

为了从全波段光谱数据中提取对小麦条锈病敏感的特征参量, 提高小麦条锈病遥感探测模型的运行效率和精度, 本文首先从惯性权重和粒子更新方式两个方面对传统离散粒子群算法(discrete binary particle swarm optimization, DBPSO)进行改进, 利用改进离散粒子群算法(modified discrete binary particle swarm optimization, MDBPSO)从全波段光谱数据中优选遥感探测小麦条锈病严重度的特征变量, 然后与冠层日光诱导叶绿素荧光(solar-induced chlorophyll fluorescence, SIF)数据相结合作为自变量分别利用随机森林(random forest, RF)和后向传播(back propagation, BP)神经网络算法构建小麦条锈病遥感探测模型, 并将其与相关系数(correlation coefficient, CC)分析法和DBPSO算法提取特征参量构建模型的精度进行对比分析。结果表明: (1) MDBPSO算法比传统DBPSO算法具有更快的收敛速度和更高的寻优精度, 改进前后其迭代次数从395次减少到156次, 最优适应度函数(optimum fitness value, OFV)值从0.145减小到0.127。(2)采用MDBPSO算法选择特征变量时, RF和BP神经网络两种方法构建的模型精度均高于CC分析法和DBPSO算法, 其中RF算法预测病情指数(disease index, DI)值和实测DI值间的检验集决定系数(validation set determination coefficient,2V)比CC分析法和DBPSO算法分别提高了9%和3%, 均方根误差(validation set root mean square error, RMSEV)分别降低了28%和11%, BP神经网络算法预测DI值和实测DI值间的2V比CC分析法和DBPSO算法分别提高了13%和6%, RMSEV分别降低了21%和10%, 利用MDBPSO算法优选特征参量能够提高小麦条锈病的遥感探测精度。(3)在MDBPSO、DBPSO和CC分析法3种特征选择算法中, RF算法构建的模型精度均高于BP神经网络算法, 其中RF模型预测DI值和实测DI值间的2V比BP神经网络算法至少提高了7%, 平均提高了9%, RMSEV至少降低了15%, 平均降低了20%。以MDBPSO算法优选的特征参量为自变量利用RF方法构建的小麦条锈病遥感探测的MDBPSO-RF模型是小麦条锈病遥感探测适宜模型, 该研究结果为进一步实现作物健康状况大面积高精度遥感监测提供了新的思路。

全波段反射光谱; 改进离散粒子群; 日光诱导叶绿素荧光; 小麦条锈病; 特征波段

小麦条锈病()是一种气流传播病害, 一旦爆发会造成小麦大面积减产[1]。因此, 小麦条锈病的及时监测和防治对我国农业发展具有重要意义。与传统田间取样的病害调查方法相比, 遥感技术具有宏观、快速、无损等优势, 已被广泛应用于作物病害的探测[1-8], 尤其是高光谱遥感数据包含丰富的光谱信息, 能够很好地反映作物的生化特性[9], 提高病害监测精度[10]。但已有的研究主要集中于利用少量波段信息计算反射率光谱指数[5-7]或吸收特征[1]等算法进行作物病害监测, 这些方法在一定程度上丢失了对小麦病害遥感探测的有用信息,且未考虑作物冠层几何结构及大气等因素导致的冠层反射光谱的时空变化[11], 因此针对不同时空下的小麦条锈病, 如何利用全波段信息获取动态的冠层光谱特征变量显得尤为重要。刘琦等[10,12]在325~1075 nm全波段范围内成功建立模拟识别小麦条锈病的模型。但直接使用全波段数据, 波段之间具有较高的相关性[13], 数据冗余不仅加大了计算机的运行时间, 而且对模型精度具有一定影响, 因此从全波段数据中提取有效的特征变量成为模型构建过程中至关重要的一步。粒子群算法(particle swarm optimization, PSO)是由Eberhart和Kennedy[14]于1995年提出的一种群体智能随机搜索算法, 具有鲁棒性强、参数设置少及收敛速度快的优点[15]。已有学者将粒子群算法应用于特征变量选择中并取得了理想的结果[16-17]。但传统的粒子群算法在处理高光谱的全波段数据时易陷入局部最优, 产生“早熟”收敛[18], 为了解决这一问题, 本文从惯性权重和粒子更新方式两方面对传统离散粒子群算法进行改进, 以提高小麦条锈病遥感探测精度。

模型精度除受所选特征参量的影响外, 建模算法也是影响小麦条锈病遥感探测精度的重要因素, 选择合适的模型算法有助于提高小麦条锈病严重度的遥感估测精度。与传统的经验统计模型相比, 机器学习模型具有良好的非线性拟合能力和泛化能力, 能够通过不断训练样本数据使目标达到最优化[19]。在众多机器学习算法中BP神经网络算法具有通过学习获取知识并解决问题的能力[20], RF算法具有计算量小、能有效避免过拟合的发生以及良好的噪声容忍性等优点[21], 为此本文拟在利用MDBPSO算法优选全波段光谱特征的基础上分别采用RF和BP神经网络算法构建小麦条锈病遥感探测模型。

反射率光谱数据能够敏感地反映植被生化组分信息, 但对作物光合活性不敏感, 难以揭示植被光合生理状态[22], 而冠层SIF数据则能够直接探测到植被的光合作用状态[23], 小麦受条锈病菌侵染后, 其水分及叶绿素含量、光合速率和光能转换率等一些生理生化指标均会发生变化[13], 综合利用反射率光谱在生化参数探测方面的优势和叶绿素荧光在光合生理诊断的优势可以提高小麦条锈病的探测精度[24-25]。基于此，本文利用MDBPSO算法从全波段反射率光谱数据中选择对小麦条锈病敏感的特征因子, 并将其与冠层SIF数据结合作为模型的输入变量分别通过RF和BP神经网络2种机器学习算法构建小麦条锈病严重度的遥感估测模型, 并将其与传统的DBPSO算法及CC分析法提取特征参量所建模型的精度进行对比分析, 以确定遥感探测小麦条锈病的特征, 以及模型构建的适宜算法。

1 材料与方法

1.1 试验设计

试验区位于河北省廊坊市中国农业科学院试验站(39°30′40″N, 11°63′62″E), 小麦品种为对条锈病比较敏感的铭贤169号。试验区小麦平均种植密度为113棵 m‒2。2018年4月7日给试验田灌溉充足的水分, 并于4月9日采用浓度为9 mg 100 mL‒1的孢子溶液对小麦进行条锈病接种。将试验区小麦分为健康组(编号为A、D)和染病组(编号为B、C), 每个试验组的面积为220 m2。将每个组分为8个样方(A1~A8、B1~B8、C1~C8、D1~D8)。

1.2 数据获取

1.2.1 冠层光谱测定与处理 分别于2018年5月18日、5月24日和5月30日3个时期测定小麦条锈病不同病情指数下的冠层光谱数据, 测量所使用的仪器为ASD Field Spec 4光谱仪和QE 65 pro光谱仪。其中ASD Field Spec 4测量的光谱范围为350~2500 nm, 在350~1000 nm的波段范围内光谱分辨率为3 nm, 采样间隔为1.4 nm; 在1000~2500 nm波段范围内光谱分辨率为10 nm, 采样间隔为2 nm。QE 65 pro光谱仪的光谱分辨率为0.31 nm, 采样间隔为0.155 nm, 光谱范围为685~805 nm。冠层光谱测量时间为北京时间11:00–12:30, 探头视场角为25°, 测量高度距离地面1.3 m, 每次测量前后均用标准BaSO4参考板进行校正。然后由公式(1)计算冠层反射率值。

为了消除和减弱样本间因散射引起的光谱误差,利用标准正态校正(stand normal variate, SNV)方法对原始光谱数据进行预处理, 假设每个波段的光谱吸收值均满足正态分布, 并利用此假设对每条光谱进行校正。SNV的变换公式[26]如式(2)所示。

式中,为原始冠层反射率,target为目标辐亮度,board为参考板辐亮度,board为参考板反射率。

1.2.2 病情指数调查 冠层病情指数的调查采用5点取样法, 在每块小区内选取对称的5点, 每点约1 m2, 随机选取30株小麦分别调查其发病情况。参照国家标准“小麦条锈病测报技术规范”(GB/T15795)量化病情指数。单叶严重度分为9个梯度, 即0、1%、10%、20%、30%、45%、60%、80%和100%的叶片病斑覆盖, 分别记录不同病情严重度下的小麦叶片数, 按式(3)计算测试群体的病情指数[27]。

式中, DI为病情指数;为梯度级值;为最高梯度等级值;为各梯度的叶片数。

1.2.3 日光诱导叶绿素荧光提取 利用辐亮度基于夫琅和费暗线原理的日光诱导叶绿素荧光算法, 主要包括FLD (fraunhofer line discrimination)和3FLD (3 bands FLD)等[28]。标准FLD算法是基于吸收线内外反射率和透过率相等的假设来估测日光诱导叶绿素荧光强度, 但由于吸收线内外波段的反射率和荧光值实际上存在差异, 影响了荧光的估测精度[29]。3FLD荧光预测算法认为在吸收线波段周围叶绿素荧光和反射率光谱是线性变化的, 利用吸收线左右各一个波段的加权平均值代替标准FLD算法中的单一波段值, 从而在一定程度上减小标准FLD方法中荧光和反射率恒定假设所带来的误差, 提高了冠层SIF的预测精度[30], 而且已有研究表明3FLD算法是估测SIF最鲁棒性算法[31], 基于此本文采用3FLD算法估测叶绿素荧光强度, 计算公式如式(4)~(6)所示。

式中, λin、λleft、λright为吸收线内、左、右波段的波长; ωleft、ωright为吸收线左右2个参考波段所占的权重;left、right为吸收线左右的太阳辐照度光谱强度(μW cm‒2nm‒1);left、right为吸收线左右的植被冠层反射辐亮度光谱强度(μW cm‒2nm‒1sr‒1)。

为了减少太阳光强度对冠层SIF估测结果的影响, 将计算得到的SIF强度除以夫琅和费吸收线内参考板的太阳入射辐照度, 得到该吸收线处的冠层SIF相对强度[32]。如式(7)所示。

1.3 改进离散粒子群算法

1.3.1 离散粒子群算法 离散粒子群算法是在粒子群算法的基础上提出的一种离散群体智能随机搜索算法[14], 具有结构简单, 参数较少及易于实现等优点[33]。假定维搜索空间中第个粒子的速度和位置分别为v= (v1,v2, ...,vD),x= (x1,x2, ...,xD), 则粒子在第维的速度和位置的更新公式如式(8)~(10)所示。

1.3.2 离散粒子群改进策略 传统的离散粒子群算法一旦达到了收敛, 所有的粒子都会朝着最优解靠拢, 失去了粒子的多样性, 影响了算法的收敛速度及精度。而且传统离散粒子群算法缺乏对速度的动态调整容易发生“早熟”收敛现象, 从而导致算法的收敛精度较低, 基于此, 本文分别从惯性权重和粒子更新方式两个方面对传统粒子群算法进行改进。

1) 惯性权重: 惯性权重值越大越有利于加快收敛速度, 提高全局搜索能力, 但难以得到精确解, 反之则有利于局部搜索且能够得到更为精确的解, 但收敛速度慢, 易陷入局部极值。为了提高离散粒子群算法的全局搜索能力, 在搜索早期将设置为相对较大的值, 在搜索过程中, 随着迭代次数的增加使值逐渐减小, 从而得到精确的解[34]。基于此, 本研究对做了如式(11)所示的动态调整。

式中max和min分别为的最大值和最小值, 本文中max和min分别取0.9和0.4[35],为当前迭代次数,max为最大迭代次数。

式中各变量的含义与公式(8)~(10)一致。

MDBPSO算法具有较好的局部搜索能力, 但是全局搜索能力相对较差, 为了使算法同时具有较好的局部搜索和全局搜索能力, 本文在迭代前期(≤ 30%max)通过式(8)~(10)对粒子进行更新, 在迭代中后期(30%max<

1.3.3 适应度函数设计 适应度函数是粒子群算法中判断粒子位置优劣的重要依据, 适应度函数的设计直接影响到算法的收敛速度以及能否找到最优解。本研究基于均方根误差RMSE [式(14)]构建适应度函数, 并以此作为各粒子优胜劣汰的评判标准。将各粒子当前的适应度值与其个体历史最佳位置对应的适应度作比较, 如果当前位置的适应度值更小, 则用当前位置更新全局最佳位置。

1.4 模型构建与精度评价

以CC分析法、DBPSO和MDBPSO三种算法提取的全波段光谱特征参量协同冠层SIF数据作为模型的输入变量构建小麦条锈病严重度估测模型,并通过保留样本交叉检验的方式对模型精度进行评价, 以确定遥感探测小麦条锈病的特征变量优选算法及适宜模型。

1.4.1 模型构建方法 分别采用RF和BP神经网络两种算法构建小麦条锈病严重度估测模型, 其中RF算法的核心思想是通过bootstrap重采样的方法在原始训练集中抽取多个样本, 对每个抽取出的样本都进行决策树建模, 最后通过多数投票法得到最终的预测结果[36]。通过对训练样本的多次仿真确定RF算法中决策树的数量(ntree)为500, 内部节点随机选择属性个数(mtry)取默认值。BP神经网络算法的核心思想是通过调整各神经元之间的权值, 将误差由隐含层向输入层逐层反传, 对误差函数进行“链式求导”, 使误差逼近最小值[37], 本研究中BP神经网络算法采用3层网络标准结构, 其中隐含层设置为5个神经元, 最大训练次数为5000, 训练间隔为10, 最小均方根误差为0.001, 学习步长为0.1。

1.4.2 模型精度评价 为了能够充分利用样本中所有数据, 在有限样本容量下尽可能减少“过拟合”问题, 本文采用保留样本交叉检验的方式对模型进行构建和检验, 将52个样本数据分为训练集和检验集, 其中39个数据(35个染病样本, 4个健康样本)作为训练样本用于模型构建, 剩余的13个数据(12个染病样本, 1个健康样本)作为验证样本用以模型评价。选择病情严重度估测值与实测值之间的2和RMSE两个模型精度评价指标[式(15), (16)], 其中2越高, RMSE越低, 则表示模型的估测精度越高。

2 结果与分析

2.1 特征变量优选

2.1.1 CC分析法优选特征参量 在对小麦冠层原始反射率光谱数据进行SNV变换处理的基础上利用CC分析法对光谱反射率数据与小麦条锈病严重度进行相关性分析, 得到SNV变换后不同波长下反射率——病情严重度间相关系数的变化曲线(图1), 选择与病情严重度达到极显著相关(相关系数的绝对值高于0.443)波段范围内的极值点作为小麦条锈病遥感探测的敏感波段。由图1可以看出, 400~540 nm、742~868 nm范围内的反射率与病情指数呈极显著负相关, 561~731 nm、905~1000 nm范围内的反射率与病情指数呈极显著正相关, 通过CC法选择出50个波段作为小麦条锈病遥感探测的特征变量。

图1 SNV变换后反射率-病情指数相关系数曲线图(n = 52)

图中虚线表示0.1%显著水平,0.001[52]= 0.443。

The dotted line in the figure indicates the significant at the 0.1% propability level, and0.001[52]= 0.443.

2.1.2 MDBPSO算法优选特征参量 为了从全波段反射率光谱数据中确定小麦条锈病遥感探测的特征变量, 分别基于MDBPSO和DBPSO算法提取了小麦条锈病遥感探测的特征波段。首先通过对样本的多次仿真训练确定DBPSO和MDBPSO的最大迭代次数max分别为400和200, MDBPSO速度限定的最大值max和最小值min分别为4和–4。然后在上述参数确定的基础上对比分析了DBPSO和MDBPSO两种算法下适应度函数的收敛情况以确定OFV值(图2)。从图2可以看出, DBPSO算法迭代395次可以得到最优解, OFV最小值为0.145。MDBPSO算法得到的OFV最小值和迭代次数分别减少到了0.127和156, MDBPSO算法比DBPSO算法具有更快的收敛速度和更高的寻优精度。这是由于MDBPSO算法在迭代过程中惯性权重发生了相应的动态变化, 从而使得OFV的值变化较快, 在提高效率的同时也增加的寻优的精度。

依据已确定的OFV的最小值和迭代次数, 分别通过DBPSO和MDBPSO两种算法从全波段反射率光谱数据中提取小麦条锈病遥感探测的特征变量(图3)。从图3中可以看出, DBPSO和MDBPSO两种算法提取特征变量的个数分别为320和56, MDBPSO算法极大减少了特征参量的个数, 降低了算法运行的时间。此外, DBPSO算法选择的特征波段分布集中, 波段之间仍然存在较高的相关性, 波段信息重叠的问题依然存在, 特征波段冗余度较高的问题依然没有得到较好的解决。

图2 DBPSO和MDBPSO算法的适应度函数收敛曲线

图3 DBPSO和MDBPSO算法选择的特征波段

2.2 构建及精度评价

分别利用CC分析法、传统DBPSO和MDBPSO法从全波段反射率光谱数据中选择特征波段作为遥感探测小麦条锈病反射率数据的模型输入变量, 并将其与冠层SIF数据相结合作为模型输入的自变量, 通过RF算法和BP神经网络算法构建小麦条锈病遥感探测模型(图4), 图中横坐标为小麦条锈病严重度的估测值, 纵坐标为小麦条锈病严重度的实测值, 实线为1∶1关系线。在病情严重度实测值-估测值构成的二维特征空间中, 散点的分布越接近于1∶1关系线, 检验集和训练集的决定系数2V和2T越高, 均方根误差RMSEV和RMSET越小, 模型的估测精度越高。

由图4的模型训练集估测结果可以看出, 无论是采用RF算法还是BP神经网络算法构建模型, 以MDBPSO算法优选的特征参量为自变量构建的模型精度较CC分析法和DBPSO算法均有一定程度的提高。此外在CC、DBPSO和MDBPSO这3种特征优选算法中, RF模型的预测精度均高于对应的BP神经网络模型, RF算法构建模型的散点分布与对应的BP神经网络算法相比更加贴合于1∶1关系线, 所以RF算法更适合小麦条锈病的遥感探测。

从模型检验集的数据中可以看出, 无论采用BP神经网络算法还是RF算法构建小麦条锈病的探测模型, MDBPSO算法提取特征变量构建的模型精度均高于CC分析法和DBPSO算法, 当使用RF算法构建模型时2V与CC分析法和DBPSO算法相比分别提高了9%和3%, RMSEV分别减少了28%和11%; 当使用BP神经网络算法构建模型时2V与CC分析法和DBPSO算法相比分别提高了13%和6%, RMSEV分别减少了21%和10%。无论采用CC法、DBPSO算法还是MDBPSO算法进行特征选择, RF算法构建的模型精度均高于对应的BP神经网络算法, CC、DBPSO、MDBPSO三种特征选择算法中RF模型的2V比BP神经网络模型至少提高了7%, 平均提高了9%, RMSEV比BP神经网络模型至少减少了15%, 平均减少了20%。综上所述, 在所有的模型构建方法中, MDBPSO-RF模型精度最高, 更适合小麦条锈病病情指数的估测。

图4 小麦条锈病病情指数实测值与预测值分布

— 1:1关系线; ■ 训练集数据; △检验集数据。

— 1:1 relationship line; ■ training set; △ test set.

3 讨论

本文综合利用反射率光谱数据在作物生化参数探测方面的优势和叶绿素荧光在光合生理诊断方面的优势, 分别采用CC分析法、DBPSO和MDBPSO三种方法从全波段光谱数据中提取特征变量, 并将其与冠层SIF数据相融合构建小麦条锈病遥感探测模型, 分析不同特征变量提取方法对小麦条锈病遥感探测精度的影响。

MDBPSO算法在模型效率方面较DBPSO算法具有一定提高, 其中OFV最小值从0.145减小到0.127, 迭代次数从395次减小到156次。无论建模方法采用BP神经网络算法还是RF算法, MDBPSO算法提取特征变量构建的模型精度均高于CC分析法和DBPSO算法, 其中BP神经网络方法构建的MDBPSO-BP模型预测DI值和实测DI间的2V比CC-BP和DBPSO-BP分别提高了13%和6%, RMSEV分别减小了21%和10%, RF方法构建的MDBPSO- RF模型预测DI值和实测DI值间的V2比CC-RF和DBPSO-RF分别提高了9%和3%, RMSEV分别减小了28%和11%。这是因为CC分析法依据小麦条锈严重度与光谱反射率的相关性选取的特征变量之间存在较高的相关性和共线性影响了模型的精度。与CC法相比, DBPSO算法以模型的RMSE作为特征波段适用性的评判标准, 选择的特征变量构建的模型精度得到了一定的改善。但与MDBPSO算法相比, DBPSO算法的OFV较大, 模型运行的效率较低, 这是因为MDBPSO算法在DBPSO算法的基础上对惯性权重和粒子更新方式两个方面做出了改进, 在选择特征变量时考虑了各个波段的光谱信息对模型贡献率的问题[17], 提高了模型的反演精度。

分别以CC分析法、DBPSO和MDBPSO 3种算法筛选的特征因子为自变量构建小麦条锈病严重度估测模型时, RF算法构建的模型精度均高于BP神经网络算法, 其中CC-RF、DBPSO-RF、MDBPSO-RF模型预测DI值和实测DI值之间的2V比CC-BP、DBPSO-BP、MDBPSO-BP分别提高了10%、9%、7%, RMSEV分别减少了15%、22%、23%。RF算法比BP神经网络算法更适合于小麦条锈病的遥感探测, 这是因为RF算法对噪声有较好的容忍度, 反演精度比较高[38]。而BP神经网络算法的网络权值初始化是随机的, 而且程序的运行结果和选用的输入参数均存在差异, 致使其模型精度较低[39], 影响了模型的估测精度。

本文在基于DBPSO和MDBPSO算法利用全波段反射率光谱数据提取小麦条锈病遥感探测特征因子时, 并没有考虑因全波段反射率光谱数据量大而导致的运行速度缓慢问题, 如何在利用上述两种算法进行特征提取前对全波段数据进行有效降维处理是下一步要研究的主要内容。

冠层光谱数据不仅受到病害胁迫的影响, 也受到冠层几何结构等因素的影响, 本文在利用实验数据建立小麦条锈病遥感探测模型时并未考虑这些问题, 这是本文的不足之处。如何消除叶面积指数以及冠层几何结构等因素对冠层光谱数据的影响, 提高小麦条锈病的遥感探测精度是一个值得探讨的问题。

本研究只比较了CC、DBPSO、MDBPSO三种特征选择方法和RF、BP神经网络两种建模算法在小麦条锈病遥感探测中的适用性, 如果使用更多的特征选择算法或建模算法, MDBPSO-RF是否仍为遥感探测小麦条锈病的最适宜方法还需要进一步的比较。

4 结论

以MDBPSO算法提取的特征因子为自变量利用BP神经网络算法构建的小麦条锈病遥感探测模型预测DI值和实测DI值间的RMSEV比CC分析法和DBPSO算法分别减小了21%和10%, 利用RF算法构建模型预测DI值和实测DI值间的RMSEV比CC法和DBPSO算法分别减小了28%和11%。利用文中3种算法筛选的特征因子为自变量构建小麦条锈病严重度估测模型时, RF算法构建的模型精度均高于BP神经网络算法, 其估测DI值与实测DI值间的RMSEV比BP神经网络模型平均减少了20%, 最少减少了15%。因此在构建小麦条锈病遥感探测模型时, MDBPSO算法为特征选择的适宜算法, RF算法为模型构建的适宜算法。

[1] 董锦绘, 杨小冬, 杨贵军, 王宝山. 基于近地高光谱信息的小麦条锈病病情指数反演. 麦类作物学报, 2016, 36: 1674–1680. Dong J H, Yang X D, Yang G J, Wang B S. Inversion of wheat stripe rust disease index based on near ground hyperspectral data., 2016, 36: 1674–1680 (in Chinese with English abstract).

[2] 赵叶, 竞霞, 黄文江, 董莹莹, 李存军. 日光诱导叶绿素荧光与反射率光谱数据监测小麦条锈病严重度的对比分析. 光谱学与光谱分析, 2019, 39: 2739–2745. Zhao Y, Jing X, Huang W J, Dong Y Y, Li C J. Comparison of sun-induced chlorophyll fluorescence and reflectance data on estimating severity of wheat stripe rust., 2019, 39: 2739–2745 (in Chinese with English abstract).

[3] Liang D, Liu N, Zhang D Y, Zhao J L, Ding Y W. Discrimination of powdery mildew and yellow rust of winter wheat using high-resolution hyperspectra and imageries., 2017, 46: 50–58.

[4] Zheng Q, Huang W J, Cui X M, Shi Y, Liu L Y. New spectral index for detecting wheat yellow rust using sentinel-2 multispectral imagery.(Basel), 2018, 18: 868–886.

[5] 梁琨, 张夏夏, 丁静, 徐剑宏, 韩东燊, 沈明霞. 傅里叶中红外光谱结合稀疏表示分类方法鉴别小麦赤霉病感染等级. 光谱学与光谱分析, 2019, 39: 3251–3255. Liang K, Zhang X X, Ding J, Xu J H, Han D S, Shen M X. Discrimination of wheat scab infection level by Fourier mid-infrared technology combined with sparse representation based classification method., 2019, 39: 3251–3255 (in Chinese with English abstract).

[6] Shi Y, Huang W J, González-Moreno P, Luke B, Dong Y Y, Zheng Q, Ma H Q, Liu L Y. Wavelet-based rust spectral feature set (WRSFs): a novel spectral feature set based on continuous wavelet transformation for tracking progressive host-pathogen interaction of yellow rust on wheat., 2018, 252: 1–19.

[7] Shi Y, Huang W J, Zhou X F. Evaluation of wavelet spectral features in pathological detection and discrimination of yellow rust and powdery mildew in winter with hyperspectral reflectance data., 2017, 11: 1–18.

[8] 蒋金豹, 陈云浩, 黄文江, 李京. 冬小麦条锈病严重度高光谱遥感反演模型研究. 南京农业大学学报, 2007, 30(3): 63–67. Jiang J B, Chen Y H, Huang W J, Li J. Study on hyperspectral remote sensing retrieval models about winter wheat stripe rust severity., 2007, 30(3): 63–67 (in Chinese with English abstract).

[9] Wang L, Qu J J. Satellite remote sensing applications for surface soil moisture monitoring: a review.China, 2009, 3: 237–247.

[10] 刘琦, 谷医林, 王翠翠, 王睿, 李薇, 马占鸿. 基于偏最小二乘法的小麦条锈病潜育期冠层高光谱分析. 植物保护学报, 2018, 45: 138–145. Liu Q, Gu Y L, Wang C C, Wang R, Li W, Ma Z H. Canopy hyperspectral features analysis of latent period wheat stripe rust based on discriminant partial least squares., 2018, 45: 138–145 (in Chinese with English abstract).

[11] 黄木易, 黄文江, 刘良云, 黄义德, 王纪华, 赵春江, 万安民. 冬小麦条锈病单叶光谱特性及严重度反演. 农业工程学报, 2004, 20(1): 176–180. Huang M Y, Huang W J, Liu L Y, Huang Y D, Wang J H, Zhao C J, Wan A M. Spectral reflectance feature of winter wheat single leaf infected with stripe rust and severity level inversion., 2004, 20(1): 176–180 (in Chinese with English abstract).

[12] 刘琦, 王翠翠, 王睿, 谷医林, 李薇, 马占鸿. 潜育期小麦条锈菌的高光谱定性识别. 植物保护学报, 2018, 45: 153–160. Liu Q, Wang C C, Wang R, Gu Y L, Li W, Ma Z H. Hyperspectral qualitative identification on latent period of wheat stripe rust., 2018, 45: 153–160 (in Chinese with English abstract).

[13] Yang C, Tan Y L, Bruzzone L, Lu L J, Guan R C. Discriminative feature metric learning in the affinity propagation model for band selection in hyperspectral images., 2017, 9: 782–798.

[14] Kennedy J, Eberhart R. A discrete binary version of the particle swarm algorithm. In: IEEE, eds. International Conference on Systems, Man, and Cybernetics. Florida, USA: Computational Cybernetics and Simulation, 1997. pp 4104–4108.

[15] 沈林成, 霍霄华, 牛轶峰. 离散粒子群优化算法研究现状综述. 系统工程与电子技术, 2008, 30: 1986–1990. Shen L C, Huo X H, Niu Y F. Survey of discrete particle swarm optimization algorithm., 2008, 30: 1986–1990 (in Chinese with English abstract).

[16] Yin C, Ye Y, Zhao H, Jiang Y Z, Wang H, Shang Y Z, Wang J F. Remote sensing of water quality based on HJ-1A HSI imagery with modified discrete binary particle swarm optimization-partial least squares (MDBPSO-PLS) in inland waters: a case in Weishan Lake., 2018, 44: 21–32.

[17] 张珏, 田海清, 赵志宇, 张丽娜, 张晶, 李斐. 基于改进离散粒子群算法的青贮玉米原料含水率高光谱检测. 农业工程学报, 2019, 35(1): 285–293. Zhang J, Tian H Q, Zhao Z Y, Zhang L N, Zhang J, Li F. Moisture content detection in silage maize raw material based on hyperspectrum and improved discrete particle swarm., 2019, 35(1): 285–293 (in Chinese with English abstract).

[18] Yang J, Zhang H S, Ling Y, Cheng P. Task allocation for wireless sensor network using modified binary particle swarm optimization., 2014, 14: 882–891.

[19] 王丽爱, 周旭东, 朱新开, 郭文善. 基于HJ-CCD数据和随机森林算法的小麦叶面积指数反演. 农业工程学报, 2016, 32(3): 149‒154. Wang L A, Zhou X D, Zhu X K, Guo W S. Inverting wheat leaf area index based on HJ-CCD remote sensing data and random forest algorithm., 2016, 32(3): 149–154 (in Chinese with English abstract).

[20] 胡小平, 杨之为, 李振岐, 邓志勇, 柯长华, 汉中地区小麦条锈病的BP神经网络预测. 西北农业学报, 2000, 9(3): 28‒31. Hu X P, Yang Z W, Li Z Q, Deng Z Y, Ke C H. BP neural network prediction of wheat stripe rust in Huazhong region., 2000, 9(3): 28–31 (in Chinese with English abstract).

[21] Rodriguez-Gailano V, Mendes M P, Garcia-Soldadoo M J, Chica-Olmo M, Ribeiro L. Predictive modeling of groundwater nitratepollution using Random Forest and multisource variables related to intrinsic and specific vulnerability: a case study in an agricultural setting (Southern Spain)., 2014, 476: 189–206.

[22] Ashour loo D, Mobasheri M R, Huete A. Developing two spectral disease indices for detection for wheat leaf rust ()., 2014, 6: 4723–4740.

[23] Song L, Guanter L, Guan K, You L, Huete A, Ju W, Zhang Y. Satellite sun-induced chlorophyll fluorescence detects early response of winter wheat to heat stress in the Indian Indo-Gangetic Plains., 2018, 24: 4023–4037.

[24] 竞霞, 白宗璠, 高媛, 刘良云. 利用随机森林法协同SIF和反射率光谱监测小麦条锈病. 农业工程学报, 2019, 35(13): 154–161. Jing X, Bai Z F, Gao Y, Liu L Y. Wheat stripe rust monitoring by random forest algorithm combined with SIF and reflectance spectrum., 2019, 35(13): 154‒161 (in Chinese with English abstract).

[25] 陈思媛, 竞霞, 董莹莹, 刘良云. 基于日光诱导叶绿素荧光与反射率光谱的小麦条锈病探测研究. 遥感技术与应用, 2019, 34: 511‒520. Chen S Y, Jing X, Dong Y Y, Liu L Y. Detection of wheat stripe rust using solar-induced chlorophyll fluorescence and reflectance spectral indices., 2019, 34: 511–520 (in Chinese with English abstract).

[26] 孙红, 郑涛, 刘宁, 程明, 李民赞, Zhang Q. 高光谱图像检测马铃薯植株叶绿素含量垂直分布. 农业工程学报, 2018, 34(1): 149‒156. Sun H, Zheng T, Liu N, Cheng M, Li M Z, Zhang Q. Vertical distribution of chlorophyll in potato plants based on hyperspectral imaging., 2018, 34(1): 149–156 (in Chinese with English abstract).

[27] 鲁军景, 黄文江, 蒋金豹, 张竞成.小波特征与传统光谱特征预测冬小麦条锈病病情严重度的对比研究. 麦类作物学报,2015,35:1456–1461.Lu J J, HuangW J, Jiang J B,Zhang J C. A comparative study on the severity of stripe rust disease in winter wheat estimated by little baud sign and traditional spectral features., 2015, 35: 1456–1461 (in Chinese with English abstract).

[28] 刘新杰. 日光诱导叶绿素荧光的遥感反演研究. 中国科学院大学博士论文, 北京, 2016. pp 17–19. Liu X J. Retrieval of Solar-Induced Chlorophyll Fluorescence by Remote Sensing. PhD Dissertation of University of Chinese Academy of Sciences, Beijing, China, 2016. pp 17–19 (in Chinese with English abstract).

[29] 刘良云. 植被定量遥感原理与应用. 北京: 科学出版社, 2014. pp 146–153 Liu L Y. Principles and Applications of Quantitative Remote Sensing for Vegetation. Beijing: Science Press, 2014. pp 146–153 (in Chinese).

[30] Maier S W, Günther K P, Stellmes M. Sun-induced fluorescence: a new tool for precision farming. In: McDonald M, Schepers J, Tartly L, eds. Digital Imaging Spectral Techniques: Applications to Precision Agriculture Crop Physiology. Madison, USA: American Society of Agronomy, 2003. pp 209–222.

[31] Liu X J, Liu L Y. Improving chlorophyll fluorescence retrieval using reflectance reconstruction based on principal components analysis., 2015, 12: 1645–1649.

[32] Cordon G, Lagorio M G, Paruelo J M. Chlorophyll fluorescence, photochemical reflective index and normalized difference vegetative index during plant senescence., 2016, 199: 100–110.

[33] 曹引, 冶运涛, 赵红莉, 蒋云钟, 王浩, 严登明. 基于离散粒子群和偏最小二乘的水源地浊度高光谱反演. 农业机械学报, 2018, 49(1): 173–182. Cao Y, Ye Y T, Zhao H L, Jiang Y Z, Wang H, Yan D M. Satellite hyperspectral retrieval of turbidity for water source based on discrete particle swarm and partial least squares., 2018, 49(1): 173–182 (in Chinese with English abstract).

[34] Yang H C, Zhang S B, Deng K Z, Du P J. Research into a feature selection method for hyperspectral imagery using PSO and SVM., 2007, 17: 473–478.

[35] 张绘娟, 张达敏, 闫威, 陈忠云, 辛梓芸. 异构网络中基于吞吐量优化的资源分配机制. 计算机科学, 2019, 46(10): 109–115. Zhang H J, Zhang D M, Yan W, Chen Z Y, Xin Z Y. Throughput optimization based resource allocation mechanism in heterogeneous networks., 2019, 46(10): 109–115 (in Chinese with English abstract).

[36] 李健丽, 董莹莹, 师越, 朱溢佞, 黄文江. 基于随机森林模型的小麦白粉病遥感监测方法. 植物保护学报, 2018, 45: 395–396. Li J L, Dong Y Y, Shi Y, Zhu Y N, Huang W J. Remote sensing monitoring of wheat powdery mildew based on random forest model., 2018, 45: 395–396 (in Chinese with English abstract).

[37] 袁冰清, 程功, 郑柳刚. BP神经网络基本原理. 数字通信世界, 2008, 8(17): 28–29. Yuan B Q, Cheng G, Zheng L G. Basic principle of BP neural networks., 2008, 8(17): 28–29 (in Chinese with English abstract).

[38] 姚雄, 余坤勇, 杨玉洁, 曾琪, 陈璋昊, 刘健. 基于随机森林模型的林地叶面积指数遥感估算. 农业机械学报, 2017, 48(5): 159–166. Yao X, Yu K Y, Yang Y J, Zeng Q, Chen Z H, Liu J. Estimation of forest leaf area index based on random forest model and remote sensing data., 2017, 48(5): 159–166 (in Chinese with English abstract).

[39] 依尔夏提·阿不来提, 买买提·沙吾提, 白灯莎·买买提艾力, 安申群, 马春玥. 基于随机森林法的棉花叶片叶绿素含量估算. 作物学报, 2019, 45: 81–90.Ershat A, Mamat S, Baidengsha M, An S Q, Ma C Y. Estimation of leaf chlorophyll content in cotton based on the random forest approach.,2019, 45: 81–90 (in Chinese with English abstract).

Canopy SIF synergize with total spectral reflectance optimized by the MDBPSO algorithm to monitor wheat stripe rust

BAI Zong-Fan1, JING Xia1,*, ZHANG Teng1, and DONG Ying-Ying2

1College of Geometrics, Xi’an University of Science and Technology, Xi’an 710054, Shaanxi, China;2Institute of Remote Sensing and Digital Earth, Chinese Academy of Science, Beijing 100094, China

In order to extract the characteristic parameters sensitive to wheat stripe rust from total spectral reflectance and improve the operation efficiency and estimation accuracy of the wheat stripe rust remote sensing detection model, this paper improved the traditional discrete binary particle swarm optimization (DBPSO) algorithm from two aspects: inertia weight and particle update method. The modified discrete binary particle swarm optimization (MDBPSO) algorithm was used to select the characteristic parameters for the severity of wheat stripe rust from the total spectral reflectance. The selected characteristic variables were combined with canopy solar-induced chlorophyll fluorescence (SIF) data as independent variables to construct wheat stripe rust estimation model. Random forest (RF) algorithm and back propagation (BP) neural network algorithm were used as model construction method to compare and analyze the accuracy of the feature parameter construction model extracted by the correlation coefficient (CC) analysis method and the DBPSO algorithm. The MDBPSO algorithm had faster convergence speed and higher optimization accuracy than the DBPSO algorithm. The number of iterations before and after the improvement was reduced from 395 to 156. The optimal fitness function (OFV) value decreased from 0.145 to 0.127. When the MDBPSO algorithm was used to select feature variables, the accuracy of the models constructed by the two methods of RF and BP neural networks was higher than that by CC analysis and DBPSO. The validation set determination coefficient (2V) between the predicted disease index (DI) value and the measured DI value of the RF algorithm was 9% and 3% higher than that of the CC analysis method and the DBPSO algorithm, the validation set root mean square error (RMSEV) was reduced by 28% and 11%, respectively. The2Vbetween the predicted disease index (DI) value and the measured DI value of the BP neural network algorithm was 13% and 6% higher than that of the CC analysis method and DBPSO algorithm, respectively, and the RMSEVwas reduced by 21% and 10% respectively. The MDBPSO algorithm can improve the remote sensing detection accuracy of wheat stripe rust. Among the three feature selection algorithms of MDBPSO, DBPSO, and CC analysis, the accuracy of the model constructed by the RF algorithm was higher than that by the BP neural network algorithm. The2Vbetween the predicted DI value and the measured DI value of the RF model was at least 7% higher than that of BP neural network algorithm, with an average increase of 9%; RMSEVhad reduced by at least 15%, with an average reduction of 20%. The MDBPSO-RF model for wheat stripe rust remote sensing detection constructed by RF method using the characteristic parameters selected by the MDBPSO algorithm as independent variables is an appropriate model for wheat stripe rust remote sensing detection. The research results provide new ideas for further realizing large-area high-precision remote sensing monitoring of crop health.

total spectral reflectance; modified particle swarm optimization; solar-induced chlorophyll fluorescence; stripe rust of wheat; feature band

10.3724/SP.J.1006.2020.01004

本研究由国家自然科学基金项目(41601467)资助。

This study was supported by the National Natural Science Foundation of China (41601467).

竞霞, E-mail: jingxiaxust@163.com

E-mail: bzf1529@163.com

2020-01-09;

2020-03-24;

2020-04-17.

URL: http://kns.cnki.net/kcms/detail/11.1809.S.20200417.1057.002.html