摘要
区域叶面积指数(Leaf Area Index, LAI)定量反演是开展大尺度农作物长势监测和产量估算的重要基础。针对当前区域LAI遥感定量反演存在的反演精度不理想和模型稳定性弱等问题,提出了一种基于少量训练样本进行LAI高精度反演的深度神经网络(Small Simple Learning LAI-Net, SSLLAI-Net)。该网络由2个卷积层、1个池化层和3个全连接层构成,将光谱反射率数据作为网络输入端、输出端得到LAI反演值,且该网络模型可支持小样本数据量的训练。以德国阿尔卑斯山麓高光谱遥感卫星影像Environmental Mapping and Analysis Program (EnMAP)为数据源,以该区域的谷物、玉米、油菜、其他作物为研究对象,数值实验结果表明当各作物类别的训练样本量均为50时,基于SSLLAI-Net的LAI反演精度分别为0.95、0.99、0.98、0.90;且在添加噪声的情况下,各作物类别的LAI反演精度分别为0.95、0.98、0.96、0.89。综上,提出的基于深度神经网络的区域LAI遥感定量反演方法SSLLAI-Net是鲁棒可靠的,且该模型能够支持稳定的小样本建模。
叶面积指数(Leaf Area Index, LAI)是描述陆地生态系统和植被结构特征的一个重要的植物学参数,定量获取区域时空连续的植被LAI对农作物长势监测、产量估算具有重要意
植被的LAI与光谱反射率之间通常认为是非线性关
以欧空局公开发布的EnMAP高光谱反射率和LAI产品为数据源(http://www.enmap.org/),选取了位于德国阿尔卑斯山麓(48.0514 N,111.0760 E)的一套科学数据集,数据获取时间为2006年7月22日,并以该区域内正处于中后期生长阶段的谷物、玉米、油菜和其他作物为研究对象。其中,EnMAP反射率数据已经过大气校

图1 研究区EnMAP遥感数据及产品 (a) 主要作物分布图 (b) LAI产品
Fig.1 Study area and EnMAP dataset(a) spatial distribution of main crops (b) LAI product
深度神经网络通常由卷积层、池化层、全连接层构成,其中,1)卷积层主要用于提取高光谱遥感图像中的光谱特征,其参数由多个可学习的卷积核集合构成;利用卷积核运算时,需设定其尺寸、步长、填充方式,如
, | (1) |
. | (2) |

图2 深度神经网络基本结构(a)卷积,(b)池化,(c)全连接
Fig.2 The basic structure of deep neural network(a)convolution,(b)pooling,(c)fully connected
为实现区域LAI的高精度反演并降低网络模型对训练样本数据的依赖性,本文提出了一种小样本深度学习网络SSLLAI-Net来进行LAI定量反演。SSLLAI-Net以EnMAP光谱反射率数据为输入,其大小为1×244,以LAI为输出,其大小为1×1,网络结构包括2个卷积层、1个池化层和3个全连接层,如

图3 基于SSLLAI-Net的LAI反演流程图
Fig.3 The flowchart of LAI estimation with our new proposed network
首先,利用两个卷积层提取光谱特征,设置两个卷积层的卷积核尺寸为1×3、步长为3,能够保证从光谱反射率数据中提取特征信息的同时缩小特征信息维度与输出端LAI维度的差距。其中,第一个卷积层的通道数为1,卷积核个数设为4,则第一层卷积计算后能够得到4个1×81的光谱特征;第二个卷积层的通道数为4,卷积核个数设为16,则第二层卷积计算后能够得到16个1×27的光谱特征,之后经ReLU函数激活以提升网络的非线性表达能力。
其次,将卷积层的输出作为池化层的输入进行数据降维,池化尺寸为1×3、步长为3,则池化后能够得到16个1×9的光谱特征。
最后,将池化层的输出作为全连接层的输入,全连接层神经元个数采用逐层递减的方式,第一、二、三全连接层的神经元个数分别设置为32、8、1。每个全连接层的输出均经过激活函数进行激活,其中第一和第二全连接层经ReLU函数激活,第三全连接层经Sigmoid函数激活。此外,为防止过拟合现象发生,在全连接层的输入与第一全连接层之间引入Dropout正则化,即在每次迭代中以p=0.5的概率随机让全连接层输入的一半神经元不参与网络训练。
首先从样本总体中随机选取一定数量的样本作为训练样本,包括光谱反射率数据和LAI。然后,设定网络超参数,包括初始学习率=0.01、数据批尺寸batchsize=100以及迭代次数epoch=10 000。此外,训练网络参数采用Adam算法,Adam算法中的超参数为Pytorch框架默认值;学习率衰减呈阶梯型下降,其中前1 000次迭代,每迭代250次=0.6×,之后每迭代1 000次=0.6×。最后将训练样本代入SSLLAI-Net网络,基于正向传播和反向传播过程更新网络参数,其中SSLLAI-Net网络参数包含卷积层的权重和偏置项、全连接层的权重和偏置项。在SSLLAI-Net网络中,训练样本以batchsize批量传入,通过前向传播计算得到LAI数据,并基于损失函数计算每个训练样本反演LAI与EnMAP LAI之间的差异,称为残差。这里选用均方误差损失函数,详见公式(3)。
, | (3) |
其中,W表示卷积层权重、全连接层权重,b表示卷积层偏置项、全连接层偏置项,为EnMAP LAI值,为经过SSLLAI-Net得到的LAI反演值,代表的维度,这里。
网络参数更新依据代价函数最小化,代价函数详见公式(4),其中m=batchsize,参数更新方式见公式(5)-(6)。
, | (4) |
, | (5) |
, | (6) |
因此,每次训练样本以batchsize批量传入SSLLAI-Net,其经过一次前向和一次反向传播过程,网络参数更新一次。所有训练样本均完成此过程,称为一次迭代。依据设定的迭代次数epoch,进行反复训练,不断缩小LAI反演值与EnMAP LAI的差异,最终达到收敛状态,迭代终止,完成训练。
以EnMAP高光谱数据为研究对象,其对应的LAI产品为EnMAP LAI值,开展基于SSLLAI-Net的植被叶面积指数遥感定量反演数值实验。首先,从谷物、玉米、油菜以及其他作物四类研究对象中随机选取一定数量的训练样本;然后,基于该样本训练SSLLAI-Net并更新网络参数,得到优化的网络;最后,将优化的SSLLAI-Net模型应用于该类别的验证样本数据,完成对该类作物LAI的定量反演,并依据LAI数据产品进行精度检验。
本文为测试网络模型对训练数据的依赖性,随机抽取了不同大小的训练样本,其数据量依次为:5000个、1 000个、500个、300个、100个、50个,基于训练样本的每类农作物LAI反演结果见
模型 | 谷物 总量108187 | 玉米 总量25334 | 油菜 总量28870 | 其他作物 总量27071 | ||||
---|---|---|---|---|---|---|---|---|
RMSE | RMSE | RMSE | RMSE | |||||
SSLLAI-Net5000 | 0.99 | 0.06 | 0.99 | 0.03 | 0.99 | 0.09 | 0.99 | 0.06 |
SSLLAI-Net1000 | 0.99 | 0.11 | 0.99 | 0.05 | 0.99 | 0.10 | 0.98 | 0.10 |
SSLLAI-Net500 | 0.99 | 0.15 | 0.99 | 0.06 | 0.99 | 0.13 | 0.97 | 0.14 |
SSLLAI-Net300 | 0.99 | 0.15 | 0.99 | 0.07 | 0.98 | 0.16 | 0.97 | 0.15 |
SSLLAI-Net100 | 0.98 | 0.23 | 0.99 | 0.09 | 0.98 | 0.20 | 0.95 | 0.18 |
SSLLAI-Net50 | 0.95 | 0.34 | 0.99 | 0.12 | 0.98 | 0.21 | 0.90 | 0.28 |

图4 基于SSLLAI-Net的叶面积指数定量反演散点密度图(a)谷物SSLLAI-Net5000 (b)玉米SSLLAI-Net5000 (c)油菜SSLLAI-Net5000 (d)其他作物SSLLAI-Net5000(e)谷物SSLLAI-Net300 (f)玉米SSLLAI-Net300 (g)油菜SSLLAI-Net300 (h)其他作物SSLLAI-Net300(i)谷物SSLLAI-Net50 (j)玉米SSLLAI-Net50 (k)油菜SSLLAI-Net50 (l)其他作物SSLLAI-Net50
Fig.4 Scatterplot of LAI estimation with our new proposed network(a)SSLLAI-Net5000 of cereals (b) SSLLAI-Net5000 of maize (c) SSLLAI-Net5000 of rape (d) SSLLAI-Net5000 of other crops(e)SSLLAI-Net300 of cereals (f) SSLLAI-Net300 of maize (g) SSLLAI-Net300 of rape (h) SSLLAI-Net300 of other crops(i)SSLLAI-Net50 of cereals (j) SSLLAI-Net50 of maize (k) SSLLAI-Net50 of rape (l) SSLLAI-Net50 of other crops
提出的SSLLAI-Net适用于LAI反演,用于建模的小样本数据量在50个时,每个类别LAI的反演精度均能够达到0.90以上。且训练样本数量越大,

图5 SSLLAI-Net LAI反演精度
Fig.5 The LAI estimation accuracy with our new proposed network

图6 基于SSLLAI-Net的研究区LAI反演结果(a) SSLLAI-Net5000 (b) SSLLAI-Net 50 (c) LAI差值对比(d) SSLLAI-Net5000局部 (e) SSLLAI-Net 50局部 (f) SSLLAI-Net5000局部 (g) SSLLAI-Net 50局部
Fig.6 Regional LAI mapping result(a) SSLLAI-Net5000 (b) SSLLAI-Net 50 (c) comparison of the residual error of estimated LAI based on SSLLAI-Net5000 and SSLLAI-Net50(d)subplot of estimated LAI with SSLLAI-Net5000 (e)subplot of estimated LAI with SSLLAI-Net50 (f)subplot of estimated LAI with SSLLAI-Net5000 (g)subplot of estimated LAI with SSLLAI-Net50
从
选取了两种常用的统计类植被叶面积指数遥感定量反演模型进行对比分析,分别是对数模
, | (7) |
式(7)中,为裸土的归一化植被指数值(Normalized Difference Vegetation Index,NDVI),为LAI达到无穷大时的NDVI值,为消光系数,与植被群体结构参数(特别是叶倾角分布)和叶片光学属性有关。
, | (8) |
其中,为模型待定系数,为自然常数。
基于对数模型、指数模型和本文提出的SSLLAI-Net模型的LAI遥感定量反演结果见
模型 | 谷物 | 玉米 | 油菜 | 其他作物 | ||||
---|---|---|---|---|---|---|---|---|
RMSE | RMSE | RMSE | RMSE | |||||
Model15000 | 0.63 | 1.00 | 0.99 | 0.11 | 0.54 | 0.98 | 0.77 | 0.42 |
Model11000 | 0.62 | 1.01 | 0.99 | 0.11 | 0.53 | 0.98 | 0.76 | 0.43 |
Model1500 | 0.62 | 1.01 | 0.99 | 0.11 | 0.53 | 0.98 | 0.76 | 0.43 |
Model1300 | 0.63 | 1.00 | 0.99 | 0.11 | 0.53 | 0.99 | 0.76 | 0.43 |
Model1100 | 0.62 | 1.01 | 0.99 | 0.12 | 0.54 | 1.00 | 0.76 | 0.44 |
Model150 | 0.61 | 1.02 | 0.99 | 0.12 | 0.53 | 1.00 | 0.76 | 0.44 |
Model25000 | 0.62 | 1.01 | 0.89 | 0.42 | 0.48 | 1.04 | 0.72 | 0.46 |
Model21000 | 0.62 | 1.01 | 0.89 | 0.42 | 0.48 | 1.04 | 0.71 | 0.47 |
Model2500 | 0.62 | 1.02 | 0.89 | 0.42 | 0.48 | 1.04 | 0.71 | 0.47 |
Model2300 | 0.62 | 1.02 | 0.89 | 0.42 | 0.48 | 1.05 | 0.71 | 0.47 |
Model2100 | 0.61 | 1.05 | 0.89 | 0.42 | 0.48 | 1.08 | 0.71 | 0.48 |
Model250 | 0.61 | 1.03 | 0.89 | 0.43 | 0.47 | 1.08 | 0.71 | 0.48 |
SSLLAI-Net5000 | 0.99 | 0.06 | 0.99 | 0.03 | 0.99 | 0.09 | 0.99 | 0.06 |
SSLLAI-Net1000 | 0.99 | 0.11 | 0.99 | 0.05 | 0.99 | 0.10 | 0.98 | 0.10 |
SSLLAI-Net500 | 0.99 | 0.15 | 0.99 | 0.06 | 0.99 | 0.13 | 0.97 | 0.14 |
SSLLAI-Net300 | 0.99 | 0.15 | 0.99 | 0.07 | 0.98 | 0.16 | 0.97 | 0.15 |
SSLLAI-Net100 | 0.98 | 0.23 | 0.99 | 0.09 | 0.98 | 0.20 | 0.95 | 0.18 |
SSLLAI-Net50 | 0.95 | 0.34 | 0.99 | 0.12 | 0.98 | 0.21 | 0.90 | 0.28 |
model1为对数模型,model2为指数模型
提出的SSLLAI-Net能够实现LAI定量反演的高精度,归因于网络模型结构以及参数更新方式。SSLLAI-Net优势在于采用卷积层、池化层、全连接层相结合的模式,依次实现提取光谱特征、降维、非线性回归。SSLLAI-Net模型参数包括卷积层参数和全连接层参数,卷积层参数用于特征提取,全连接层参数用于光谱特征非线性映射到LAI,由于模型参数多,以及每个网络层的参数都有各自的功能,因此可以充分逼近光谱反射率与LAI的非线性关系,而统计模型难以实现如上效果。并且小样本数据中的光谱反射率通过SSLLAI-Net正向传播输出LAI反演值,通过反向传播过程实现网络参数更新,如此反复训练,使残差不断减小最终达到收敛状态。因此,提出的SSLLAI-Net兼具支持小样本建模和反演精度高的优点。
光子采集数据通常是带有泊松噪声,为检验SSLLAI-Net在植被叶面积指数定量反演中的抗噪性,以EnMAP原始反射率数据作为均值获得了服从泊松分布的反射率数据,以模拟带有噪声数据,并对该数据进行LAI定量反演。LAI定量反演的数值实验结果如
模型 | 谷物 总量108187 | 玉米 总量25334 | 油菜 总量28870 | 其他作物 总量27071 | ||||
---|---|---|---|---|---|---|---|---|
RMSE | RMSE | RMSE | RMSE | |||||
SSLLAI-Net5000 | 0.99 | 0.11 | 0.99 | 0.05 | 0.99 | 0.17 | 0.99 | 0.07 |
SSLLAI-Net1000 | 0.99 | 0.12 | 0.99 | 0.05 | 0.99 | 0.17 | 0.98 | 0.09 |
SSLLAI-Net500 | 0.99 | 0.16 | 0.99 | 0.09 | 0.98 | 0.17 | 0.97 | 0.15 |
SSLLAI-Net300 | 0.98 | 0.17 | 0.99 | 0.09 | 0.98 | 0.18 | 0.97 | 0.15 |
SSLLAI-Net100 | 0.97 | 0.25 | 0.99 | 0.13 | 0.97 | 0.22 | 0.93 | 0.23 |
SSLLAI-Net50 | 0.95 | 0.33 | 0.98 | 0.18 | 0.96 | 0.27 | 0.89 | 0.29 |
用50个样本进行建模,实现谷物、玉米、油菜和其他作物的LAI定量反演的

图7 基于SSLLAI-Net LAI的研究区LAI反演结果(带噪声)(a) SSLLAI-Net5000 (b) SSLLAI-Net 50 (c) LAI差值对比(d) SSLLAI-Net5000局部 (e) SSLLAI-Net 50局部 (f) SSLLAI-Net5000局部 (g) SSLLAI-Net 50局部
Fig.7 LAI inversion results based on proposed network containing noise(a) SSLLAI-Net5000 (b) SSLLAI-Net 50 (c)comparison of the residual error of estimated LAI based on SSLLAI-Net5000 and SSLLAI-Net50(d)subplot of estimated LAI with SSLLAI-Net5000 (e)subplot of estimated LAI with SSLLAI-Net50 (f)subplot of estimated LAI with SSLLAI-Net5000 (g)subplot of estimated LAI with SSLLAI-Net50
为提高叶面积指数遥感定量反演精度和稳定性,本文提出了一种基于EnMAP高光谱遥感影像数据和深度神经网络模型的LAI定量反演方法。该方法在SSLLAI-Net模型构建过程中,将卷积层、池化层与全连接层相结合,完成提取光谱特征、降维、非线性回归分析过程。由于网络的前向传播和反向传播过程,使得网络可以更好地学习数据特征,完成数据预测。此外,本文还验证了网络模型的抗噪性能,结果表明在含噪声时仍能够对LAI实现高精度反演。综上,深度神经网络模型具有很好的稳定性与抗噪性,支持小样本建模,并保证LAI反演的高精度。
提出的SSLLAI-Net在LAI反演中的效果较好,具有极大的应用潜力在后续的研究中,将继续优化该方法,以期实现基于多光谱遥感数据的参数反演,以及对作物整个生育期的时序植被理化参数的定量反演应用研究。
References
Mananze S, Pôças I, Cunha M. Retrieval of Maize Leaf Area Index Using Hyperspectral and Multispectral Data[J]. Remote Sens, 2018, 10(12), 1942. [百度学术]
LU Kun, MENG Qing-Yan, SUN Yun-Xiao, et al. Estimating leaf area index of wheat at the booting stage using GF-2 data:A case study of Langfang City,Hebei Province[J]. Remote Sensing for Land & Resources, 2018, 30(1):196-202. [百度学术]
Siegmann B, Jarmer T. Comparison of different regression models and validation techniques for the assessment of wheat leaf area index from hyperspectral data[J]. International Journal of Remote Sensing, 2015, 36(18):4519-4534. [百度学术]
ZHENG Guang, LIU Rong-Gao, LIU Yang, et al. Research progress and prospect of leaf area index retrieval by remote sensing[J]. JOURNAL OF GEO-INFORMATION SCIENCE,(郑光,刘荣高,刘洋,等.叶面积指数遥感反演研究进展与展望.地球信息科学学报), 2013, 15(5). [百度学术]
HU Gu-Yue, LI Shao-Da, YANG Rong-Hao. Comparison of three regression models for remote sensing estimation of leaf area index[J]. Science of Surveying and Mapping, 2018. [百度学术]
LI Jian-Jian, ZHU Xiao-Hua, MA Liang-Liang, et al. Leaf Area Index Retrieval and Scale Effect Analysis of Multiple Crops from UAV-based Hyperspectral Data[J]. Remote Sensing Technology & Application, 2017. [百度学术]
DONG Ying-Ying, WANG Ji-Hua, LI Cun-Jun, et al. Estimating leaf area index from remote sensing data: based on data segmentation and principal component analysis[J]. J. Infrared Millim. [百度学术]
Waves,董莹莹,王纪华,李存军,等.基于数据分割与主成分分析的LAI遥感估算.红外与毫米波学报,2011,30(02):124-130. [百度学术]
Berger K, Atzberger C, Danner M, et al. Evaluation of the PROSAIL model capabilities for future hyperspectral model environments: a review study. Remote Sensing, 2018,10(1), 85. [百度学术]
LI He, CHEN Zhong-Xin, LIU Guo-Huan, et al. Improving Winter Wheat Yield Estimation from the CERES-Wheat Model to Assimilate Leaf Area Index with Different Assimilation Methods and Spatio-Temporal Scales[J]. Remote Sens, 2017, 9(3), 190. [百度学术]
DONG Ying-Ying, WANG Ji-Hua, LI Cun-Jun, et al. Comparison and analysis of data assimilation algorithms for predicting the leaf area index of crop canopies[J]. Selected Topics in Applied Earth Observations and Remote Sensing, IEEE Journal of, 2013, 6(1):188-201. [百度学术]
FENG Wei, WU Ya-Peng, HE Li, et al. An optimized non-linear vegetation index for estimating leaf area index in winter wheat. Precision Agric, 2019. [百度学术]
SUN Yuan-Heng, REN Hua-Zhong, ZHANG Tian-Yuan, et al. Crop Leaf Area Index Retrieval Based on Inverted Difference Vegetation Index and NDVI[J]. IEEE Geoscience and Remote Sensing Letters, 2018,15(11):1662-1666. [百度学术]
Sivasankar T, Kumar D, Srivastava H S, et al. Wheat leaf area index retrieval using RISAT-1 hybrid polarized SAR data[J], Geocarto International, 2019. [百度学术]
Garrigues S, Allard D, Baret F, et al. Influence of landscape spatial heterogeneity on the non-linear estimation of leaf area index from moderate spatial resolution remote sensing data[J]. Remote Sensing of Environment, 2006, 105(4):286-298. [百度学术]
YANG Min, LIN Jie, GU Zhe-yan, et al. Leaf area index retrieval based on Landsat 8 OLI multi-spectral image data and BP neural network[J],Science of Soil and Water Conservation.(杨敏,林杰,顾哲衍,等.基于Landsat 8 OLI多光谱影像数据和BP神经网络的叶面积指数反演.中国水土保持科学),2015,13(04):86-93. [百度学术]
JIANG Chun-Lei, ZHANG Shu-Qing, ZHANG Ce, et al. Modelling and predicting of MODIS leaf area index time series based on a hybrid SARIMA and BP neural network method[J]. Spectroscopy & Spectral Analysis, 2017. [百度学术]
Ercanli I, Gnl A, Senyurt M, et al. Artificial neural network models predicting the leaf area index: a case study in pure even-aged Crimean pine forests from Turkey[J]. Forest Ecosystems, 2018, 5(04):58-69. [百度学术]
Migdall S, Bach H, Kaufmann H. Simulation of enmap-like hyperspectral images based on textural information and a radiative transfer approach [C]. In Hyperspectral 2010 Workshop (p. 6)Frascati, Italy. 2010. [百度学术]
Verhoef W, Heike Bach. Coupled soil–leaf-canopy and atmosphere radiative transfer modeling to simulate hyperspectral multi-angular surface reflectance and TOA radiance data[J]. Remote Sensing of Environment, 2007, 109(2):166-182. [百度学术]
TAN Chang-Wei, WANG Ji-Hua, HUANG Yi-De, et al. Quantitative Improvement of Beer-Lambert Law with Spectral Remote Sensing Technology and Its Application[J],Scientia Agricultura Sinica.(谭昌伟, 王纪华, 黄义德, 等. 运用光谱技术改进Beer-Lambert定律的定量化及其应用研究. 中国农业科学), 2005, 38(3):498-503. [百度学术]
XIA Tian, WU Wen-Bin, ZHOU Qing-Bo, et al. Comparison of two inversion methods for winter wheat leaf area index based on hyperspectral remote sensing[J],Transactions of the Chinese Society of Agricultural Engineering.(夏天, 吴文斌, 周清波, 等. 冬小麦叶面积指数高光谱遥感反演方法对比. 农业工程学报), 2013,29(3):139-147. [百度学术]