基于BP神经网络的夏玉米多生育期叶面积指数反演研究
1.
2.
3.
4.
Leaf Area Index Inversion of Summer Maize at Multiple Growth Stages based on BP Neural Network
1.
2.
3.
4.
通讯作者:
收稿日期: 2018-08-23 修回日期: 2019-11-15 网络出版日期: 2020-03-31
基金资助: |
|
Received: 2018-08-23 Revised: 2019-11-15 Online: 2020-03-31
作者简介 About authors
刘俊(1993-),男,湖南衡阳人,硕士研究生,主要从事定量遥感以及农业遥感研究E⁃mail:
关键词:
Keywords:
本文引用格式
刘俊, 孟庆岩, 葛小三, 刘顺喜, 陈旭, 孙云晓.
Liu Jun, Meng Qingyan, Ge Xiaosan, Liu Shunxi, Chen Xu, Sun Yunxiao.
1 引 言
反演模型、反演策略及遥感数据对作物不同生育期LAI定量反演十分重要[10,11,12]。目前,遥感影像LAI反演模型主要分两类:一类是基于遥感影像光谱数据与实测LAI值之间建立数学关系式的统计模型,常用的光谱数据包括植被指数、光谱反射率、导数光谱等[13];另一类是基于非各向同性辐射传输模型及植被—土壤波谱特性的物理模型,利用辐射传输模型的反向过程反演获得输入参数从而达到反演地表参量的目的[14,15]。高林等[16]基于高光谱数据、无人机影像及GF-1 WFV影像,分别采用5种植被指数构建植被指数(LAI)经验回归模型,对3类遥感数据反演大豆叶面积指数进行差异分析。武佳丽等[17]基于辐射传输模型和人工神经网络算法,针对环境一号卫星CCD相机特征,设计出新的一种植被指数HJVI,提高了精度且有效避免数据饱和现象。刘洋等[18]将物理模型和神经网络模型结合,利用MODIS地表反射率和4-SCALE模型反演LAI。统计模型简单高效,但忽视了光子在植被冠层内的传输过程,其普适性较差,缺乏可移植性[19,20]。物理模型机理性强,物理意义明确,不需要大量实测数据,但模型参数多而且较难获得,同时,存在“病态反演”问题[21,22]。此外,随着遥感数据对于非线性物理模型的优化要求越来越高,常规反演手段处理复杂的非线性问题有点束手无策[23,24],而神经网络在应对这些问题的拟合及模式识别方面有着无可比拟的优势,因此成为混合反演方法中的一种常用方法[25,26,27]。
植物不同生育期的叶面积指数是变化的,地表覆盖度不同且土壤背景的影响也在逐步变化,因此不同反演方法反演精度也不相同。各生育期反演方法的选取和使用及不同方法的反演特点是一个亟待解决的问题[28]。贺佳等[29]通过不同的统计方法,构建了冬小麦不同生育期LAI监测模型。在不同氮磷情况下,冬小麦LAI随施肥量增加呈现递增的趋势,随生育期的改变呈抛物线趋势变化。赵娟等[30]根据地表作物覆盖度和反射率的变化,在冬小麦的不同生育期,选择不同植被指数建立冬小麦LAI反演模型所得到的反演精度均高于在整个生育期使用NDVI模型得到的反演结果。但上述研究大都针对地面实测光谱数据进行LAI多生育期研究,较少针对遥感数据进行反演分析研究。
以河北省石家庄市栾城区为研究区域,以国产高分卫星(GF-1 WFV影像)和LAI实测数据为基础,通过建立基于不同隐含层的BP神经网络LAI反演模型,与统计模型进行不同方法的应用和综合比对,分析反演精度和误差,并对夏玉米不同生育期进行BP1模型和BP2模型反演分析,建立夏玉米随时间LAI动态变化图,为进一步使用GF卫星多光谱数据定量反演作物多生育期LAI提供参考。
2 研究区
栾城区位于冀中平原西部,河北省西南部,省会石家庄东南方。地理坐标37°59′20″~37°47′34″ N,140°28′36″~114°47′35″ E。栾城位于东部季风区,处于暖温带半湿润地区,属温带大陆性季风气候,气候温和,光照充足,降水适中,四季分明,春季干燥多风,夏季炎热多雨,秋季凉爽多雾,冬季寒冷少雪。年平均气温12.8℃,年平均降水量474 mm,年平均无霜期205 d,年日照总时数2 521.9 h,年平均太阳辐射总量125.438,年平均风速为2.6 m/s。当地农作物以玉米,小麦为主,小麦播种面积1.77万hm2,玉米播种面积1.65万hm2。研究区位置及采样点分布如图1所示。
图1
图1
研究区位置及采样点分布
Fig.1
The location of the study area and the distribution of sampling points
3 数据与方法
3.1 数据获取与处理
3.1.1 遥感影像
3.1.2 LAI实测数据
研究区采样点LAI数据的测定采用美国LI-COR LAI-2200C植物冠层分析仪,地面测量时间均采用为与遥感影像成像时间相同步的日期。在研究区选定42个样方,每个样方点选择3处测量并记录LAI值,取3次测量平均值作为该样点的LAI值,同时记录GPS经纬度值。其中,所有测量时间都选择在06:30~09:30之间或者16:00~19:00之间,尽量避免因为太阳光线直射的原因而带来的测量误差。
3.1.3 植被指数提取
其中:
3.2 BP神经网络算法
BP(Back Propagation)神经网络是一种按误差反向传播(简称为误差反传)训练的多层前馈网络,其算法称为BP算法,它的基本思想是梯度下降法,利用梯度搜索技术,以期使网络的实际输出值和期望输出值的误差均方差达到最小。BP网络是在输入层与输出层之间增加若干层(一层或多层)神经元,这些神经元被称为隐单元,它们与外界没有直接联系,但其状态的改变,能影响输入与输出之间的关系,每一层可以有若干个节点。其结构框图如图2所示。
图2
BP神经网络的计算过程由正向计算过程和反向计算过程组成。正向传播过程,输入模式从输入层经隐单元层逐层处理,并转向输出层,每一层神经元的状态只影响下一层神经元的状态。如果在输出层不能得到期望的输出,则转入反向传播,将误差信号沿原来的连接通路返回,通过修改各神经元的权值,使得误差信号最小。
BP神经网络相对其他非线性方法最重要的优势是神经网络能够全局近似,并且具有高度准确性。同时,BP神经网络的建模不需要事先假设,它很大程度上由数据的特性确定。BP神经网络能描述复杂、非线性和不确定的系统[39]。
激活函数的选择是构建神经网络过程中的重要环节,以下是两种常用的激活函数。S型函数(Log-Sigmoid Function)(式(7))和双曲正切S型函数(Tan-Sigmoid Function)(式(8)),曲线如图3所示,两者在输入很大或很小的时候,输出都几乎平滑,梯度很小。两者的主要区别在于函数的值域,双曲正切S型函数值域是(-1,1),而S型函数值域是(0,1)。
图3
3 精度评价方法
将42个样本点按照2∶1的比例随机分为3组,一组28个样本用来构建模型,训练样本组;另一组14个样本用于校验模型,测试样本组。模型精度用决定系数(R2)、均方根误差(Root-Mean-Square Error,RMSE)和相对误差(Relative Error)评价,R和RMSE计算方法如式(9)~(11)所示:
其中:
4 结果与讨论
4.1 基于BP神经网络的LAI模型构建
以4个多光谱波段值和6种植被指数作为输入变量,实测LAI值作为输出变量,隐含层采用Tan-Sigmoid函数作为传递函数,输出层采用Purelin线性传递函数,运用Matlab神经网络工具箱,训练迭代次数设置为5 000次,学习率设置为0.001,目标误差设置为0.001,采用Traindx训练函数,根据以往学者给出的隐含层节点数经验公式[40],确定隐含层节点数区间为[3,25],在此区间内进行不断地迭代循环,直至网络性能最佳时停止循环,输出隐含层节点数及最佳结果,应用此方法分别训练一层隐含层(BP1)和两层隐含层(BP2)的神经网络模型。研究采用Matlab R2016a实现BP神经网络模型。
4.2 基于BP神经网络的LAI反演
图4
图5
图5
BP神经网络模型反演值与实测值拟合图
Fig.5
Fitting diagram of inversion value and measured value of BP neural network model
图6
4.3 BP神经网络方法分析与精度评价
图7
图7
BP1和BP2反演值 - LAI实测值拟合图
(注:BP1为一层隐含层的BP神经网络模型,BP2为两层隐含层的BP神经模型;S1为训练样本点,S2为测试样本点,S3为所有样本点)
Fig.7
BP1 and BP2 inversion values - LAI measured values fitting diagram
表1 BP1和BP2 - LAI反演精度
Table 1
评价指标 | 训练样本 | 测试样本 | 所有样本 | |||
---|---|---|---|---|---|---|
BP1 | BP2 | BP1 | BP2 | BP1 | BP2 | |
R2 | 0.77 | 0.85 | 0.52 | 0.56 | 0.70 | 0.76 |
RMSE | 0.23 | 0.23 | 0.34 | 0.31 | 0.27 | 0.26 |
同时,基于GF-1 WFV影像的6种植被指数(NDVI、RVI、DVI、EVI、SAVI、ARVI)与实测LAI值(训练样本)拟合,建立线性和指数(非线性)两种关系式,对遥感影像进行LAI反演,并与实测LAI值(测试样本)进行验证分析,得到各决定系数R2和均方根误差RMSE值,结果如表2所示。
表2 统计模型各植被指数反演精度
Table 2
植被指数 | 拟合公式 | 决定系数(R2) | 均方根误差(RMSE) |
---|---|---|---|
NDVI | 0.78 | 0.42 | |
0.83 | 0.39 | ||
RVI | 0.71 | 0.41 | |
0.72 | 0.41 | ||
DVI | 0.68 | 0.51 | |
0.70 | 0.47 | ||
EVI | 0.72 | 0.47 | |
0.75 | 0.45 | ||
SAVI | 0.74 | 0.43 | |
0.76 | 0.42 | ||
ARVI | 0.74 | 0.45 | |
0.76 | 0.42 |
由表2可知,在6种植被指数LAI反演中,LAI与各植被指数模型均有良好相关性。其中,以NDVI的指数方程式回归模型拟合度最优;而其与BP神经网络模型相比,BP1和BP2的R2普遍都比各植被指数的反演结果小,相关性较之统计模型反演结果低,然而RMSE则均比各植被指数的反演结果小,拟合差异程度低于统计模型反演结果。在此次研究中,统计模型反演聚焦于小区域研究,所以其相关性更为优异,但是BP神经网络反演结果的误差更小,与实测值之间的吻合程度更好。所以,BP神经网络模型和统计模型各有优劣之处。
同时,BP神经网络模型中,基于不同隐含层的反演情况中,BP2模型在训练样本、测试样本、所有样本3组中的决定系数R2,均方根误差RMSE均优于BP1模型。由图7可见,在不同样本组中,BP2模型和BP1模型线性回归线基本分布在1∶1关系线两侧,表明BP1模型与BP2模型拟合效果较好,但是BP2模型在各样本组中回归线均更为接近1∶1关系线,预测精度更高。
4.4 夏玉米多生育期LAI反演分析
将已有的夏玉米实测数据分别按照要求代入已经建立好的BP神经网络模型(BP1和BP2)中,分别进行1层和2层隐含层反演得到夏玉米4个月不同生育期反演值,结果如图8所示。
图8
图8
BP1和BP2玉米6~9月反演结果
Fig.8
BP1 and BP2 varying curve of maize from June to September
图9
由以上分析结果可知,BP神经网络模型的非线性映射能力、自学习和自适应能力、泛化能力、容错能力对解决LAI反演复杂关系的模型拥有良好的反演结果。BP1和BP2均能够快速准确地估测夏玉米LAI值,同时BP2的LAI反演结果略优于BP1的LAI反演结果。
5 结 语
本研究以GF-1 WFV影像结合夏玉米实测数据,分别建立以4个多光谱波段和6种植被指数作为输入层,基于不同隐含层的BP神经网络模型(BP1和BP2),将反演结果与统计模型(6种植被指数)的反演结果进行精度比对分析,同时对夏玉米四个月多生育期进行反演分析,结果表明:
(1)BP1模型和BP2模型反演结果与基于GF-1 WFV影像6种植被指数(NDVI、RVI、DVI、EVI、SAVI、ARVI)分别建立线性和指数关系式得到LAI反演结果相比,其R2相对较小,相关性不如统计模型的反演结果,但RMSE更小,BP1和BP2与实测值差异性更小,与实测值的吻合程度更优。
(2)BP1模型与BP2模型的R2和RMSE值情况均较好,拟合精度较高,符合实际情况,基于BP神经网络模型拟合情况良好,能够很好适用于GF-1 WFV影像数据的LAI反演,表现出了较强的自学习和自适应能力。
(3)BP2模型反演结果在R2和RMSE上优于BP1模型,反演结果与实测LAI值的线性回归线更接近1∶1关系线,预测精度更高,反演结果更精确。
(4)基于BP1和BP2两种模型对栾城区夏玉米进行反演,绘制了LAI随时间变化的拟合图,6~9月夏玉米生长过程中LAI呈现缓慢增长—快速增长—逐渐减小的变化过程,曲线反映当地夏玉米叶面积指数增长动态变化全过程,基本吻合夏玉米实际生长情况。
但是,BP1模型和BP2模型反演结果均存在决定系数R2偏小、相关性偏低的缺点,这是因为训练样本点过少所导致,后续研究需要增加足够多的样本点个数,开展多年度夏玉米LAI动态变化曲线图,提高BP模型的适用性。
本研究使用一层和两层隐含层的BP神经网络模型2套方案,基于GF-1 WFV影像的4个波段值和6个植被指数数据来反演4个生育期的夏玉米LAI,通过决定系数(R2)和均方根误差(RMSE)来比较和分析其反演精度和拟合情况。结果表明,基于BP神经网络模型的夏玉米LAI反演精度整体较好,符合实际情况,基本吻合夏玉米实际生长情况,而隐含层个数不同对反演精度有影响,两层隐含层比一层隐含层拟合精度更高。基于BP神经网络模型的夏玉米LAI反演精度受样本参数个数、反演波段选择与模型参数三方面因素影响。当样本参数个数足够,波段选择恰当时,模型参数好的数据表现出更优的LAI反演精度和稳定性。反之,样本参数个数不够,模型参数的改变难以显著提升LAI反演精度,存在偶然性的反演情况发生。
参考文献
Leaf Area Index Retrieval based on Landsat 8 OLI Multi-spectral Image Data and BP Neural Network
[J].基于Landsat 8 OLI多光谱影像数据和BP神经网络的叶面积指数反演
[J].
Recent Progress in Land Surface Quantitative Remote Sensing
[J].陆表定量遥感反演方法的发展新动态
[J].
Comparative Physiological Studies on the Growth of Field Crops: I. Variation in Net Assimilation Rate and Leaf Area between Species and Varieties, and within and between Years
[J].
The Study of Inversion of Chlorophyll A in Taihu based on GF-1 WFV Image and BP Neural Network
[J].基于GF-1WFV影像和BP神经网络的太湖叶绿素a反演
[J].
Leaf Area Index Estimation Using Chinese GF-1 Wide Field View Data in an Agriculture Region
[J].
Study of Environmental Vegetation Index based on Environment Satellite CCD Data and LAI Inversion
[J].基于环境星CCD数据的环境植被指数及叶面积指数反演研究
[J].
The Neural Network Algorithm for Estimating Plantation Forest Leaf Area Index based on ENVISAT/ASAR
[J].高帅, 牛铮, 邬明权
.基于ENVISAT/ASAR的神经网络反演人工林叶面积指数研究
[J].
Retrieving Leaf Area Index of Corn Canopy based on Sentinel-2 Remote Sensing Image
[J].基于Sentinel-2遥感影像的玉米冠层叶面积指数反演
[J].
Modeling and Predicting of MODIS Leaf Area Index Time Series based on a Hybrid SARIMA and BP Neural Network Method
[J],基于SARIMA-BP神经网络组合方法的MODIS叶面积指数时间序列建模与预测
[J].
Comparison of Four Radiative Transfer Models to Simulate Plant Canopies Reflectance : Direct and Inverse Mode
[J].
Comparison between Multispectral and Hyperspectral Remote Sensing for LAI Estimation
[J].基于多光谱与高光谱遥感数据的冬小麦叶面积指数反演比较
[J].
Evaluate Quantitatively Effect of Understory on Leaf Area Index (LAI) Estimation Combining Active and Passive Remote Sensing Data
[J],联合森林主被动遥感数据定量评价林下植被对叶面积指数估算的影响
[J].
Comparison of Precision in Retrieving Soybean Leaf Area Index based on Multi-source Remote Sensing Data
[J].基于多源遥感数据的大豆叶面积指数估测精度对比
[J].
Analysis of Canopy Reflectance Spectra between Winter Wheat Cultivars with Different Gluten Strength
[J].不同筋力冬小麦冠层反射光谱的差异分析
[J].
Inversion of Leaf Area Index based on A Simple Physical Model
[J].基于一种简单物理模型的叶面积指数反演
[J].
The Algorithm Research of LAI Inversion of HJ-1 CCD Camera Data base on SAIL Model
[J].基于SAIL模型的HJ-1卫星LAI反演算法研究
[J].
Estimation of Crop LAI based on the Neural Network Trained from Physical Model
[J].基于物理模型训练神经网络的作物叶面积指数遥感反演研究
[J].
Inversion of LAI and Chlorophyll Content of Winter Wheat at Different Stages Using Two-layer Canopy Reflectance Model
[J].李宗南, 陈仲新, 王利民
.基于ACRM模型不同时期冬小麦LAI和叶绿素反演研究
[J].
A Feasibility Stuudy of Leaf Area Index Inversion Using Radiative Transfer Model based on TM Data
[J].陈艳华, 张万昌, 雍斌
.基于TM的辐射传输模型反演叶面积指数可行性研究
[J].
Effects of Plastic Film Residue on Evapotranspiration and Soil Evaporation in Cotton Field of Xinjiang
[J].地膜残留量对新疆棉田蒸散及棵间蒸发的影响
[J].
Retrieving Leaf Area Index Using A Neural Network based on Classification Knowledge
[J].基于分类知识利用神经网络反演叶面积指数
[J].
Estimation of Maize Leaf Area Index based on GF-1 WFV Image and Machine Learning Random Algorithm
[J].基于GF-1 WFV影像和机器学习算法的玉米叶面积指数估算
[J].
Leaf Area Index Estimates Using Remotely Sensed Data and BRDF Models in A Semiarid Region
[J].
Leaf Area Index Retrieval of Winter Wheat Using Artificial Neural Network
[J].马茵驰
,基于人工神经网络方法的冬小麦叶面积指数反演
[J].
Neural Network Estimation of LAI, fAPAR, fCover, and LAI × Cab, from Top of Canopy MERIS Reflectance Data: Principles and Validation
[J].
Neuroet : An Easy-to-Use Artificial Neural Network for Ecological and Biological Modeling
[J].
Soybean LAI Estimation with In-situ Collected Hyperspectral Data based on BP-Neural Networks
[J].
Accuracy Comparison of the Leaf Area Index Inversion Method of Winter Wheat at Different Growth Stages based on Landsat_TM Data
[D].苗乃哲. 基于Landsat_TM数据的冬小麦不同生育期叶面积指数反演方法精度比较
[D].
Monitoring Model of Leaf Area Index of Winter Wheat based on Hyperspectral Reflectance at Different Growth Stages
[J].贺佳, 刘冰锋, 李军
.不同生育时期冬小麦叶面积指数高光谱遥感监测模型
[J].
Inversion of Leaf Area Index during Different Growth Stages in Winter Wheat
[J].冬小麦不同生育时期叶面积指数反演方法
[J].
Effects of Atmospheric Correction on Remote Sensing Estimation of LAI of Broadleaved Forest
[J].大气校正模型对阔叶林叶面积指数遥感估算的影响
[J].
Comparative Study of Different Atmospheric Correction Model on Estimation of Leaf Area Index Using Remote Sensing
[J],不同大气校正模型对叶面积指数遥感反演精度的影响比较
[J].
Monitoring Vegetation Systems in the Great Plains with ERTS
[C]∥
Monitoring Vegetation Systems in the Great Plains with ERTS
[C]∥
A Soil-adjusted Vegetation Index (SAVI)
.[J].
Derivation of Leaf Area Index from Quality of Light on the Forest Floor
[J].
Atmospherically Resistant Vegetation Index (ARVI) for EOS-MODIS
[J].
A Feedback based Modification of the NDVI to Minimize Canopy Background and Atmospheric Noise
[J].
Modeling Canopy Transpiration of Young Poplar Trees ( Populus × Euramericana cv. N3016 ) based on Back Propagation Artificial Neural Network
[J].BP人工神经网络模拟杨树林冠蒸腾
[J].
Some Problems to be Considered in Using BP Algorithm
[J].郭海涛, 张殿伦, 马国芳
.使用BP算法时应考虑的若干问题
[J].
/
〈 |
|
〉 |
