检索
E-mail
RSS
高级检索 图表检索

当期目录 我要投稿 过刊浏览
高级检索 图表检索

作物杂志, 2021, 37(2): 183-190 doi: 10.16035/j.issn.1001-7283.2021.02.027

农业信息技术

结合植被指数与纹理特征的玉米冠层FAPAR遥感估算研究

王思宇,1,2, 聂臣巍2, 余汛2,3, 邵明超2,4, 王梓旭2,5, 努热曼古丽·托乎提2,6, 刘亚东2, 程明瀚2,7, 官云兰,1, 金秀良,2

1东华理工大学测绘工程学院,330013,江西南昌

2中国农业科学院作物科学研究所,100081,北京

3河南理工大学,454003,河南焦作

4河北工程大学地球科学与工程学院,056006,河北邯郸

5长安大学地球科学与资源学院,710061,陕西西安

6中国地质大学（武汉）地球科学学院,430074,湖北武汉

7河海大学农业工程学院,210098,江苏南京

Maize Canopy FAPAR Remote Sensing Estimation Combining Vegetation Indexes and Texture Characteristics

Wang Siyu,1,2, Nie Chenwei2, Yu Xun2,3, Shao Mingchao2,4, Wang Zixu2,5, Nuremanguli· Tuohuti2,6, Liu Yadong2, Cheng Minghan2,7, Guan Yunlan,1, Jin Xiuliang,2

1Faculty of Geomatics, East China University of Technology, Nanchang 330013, Jiangxi, China

2Institute of Crop Sciences, Chinese Academy of Agricultural Sciences, Beijing 100081, China

3Henan Polytechnic University, Jiaozuo 454003, Henan, China

4School of Earth Science and Engineering, Hebei University of Engineering, Handan 056006, Hebei, China

5School of Earth Sciences and Resources, Chang'an University, Xi’an 710061, Shaanxi, China

6School of Earth Sciences, China University of Geosciences, Wuhan 430074, Hubei, China

7College of Agricultural Engineering, Hohai University, Nanjing 210098, Jiangsu, China

通讯作者: 金秀良,主要从事作物表型鉴定与精准农业研究,E-mail： jinxiuliang@caas.cn; 官云兰,主要从事摄影测量与遥感技术应用、三维激光扫描等研究,E-mail： 275332432@qq.com

收稿日期: 2020-11-15   修回日期: 2021-01-9   网络出版日期: 2021-04-15

基金资助: 国家自然科学基金-面上项目(42071426)
中国农业科学院基本科研业务费专项院级统筹项目(Y2020YJ07)
中国农业科学院科技创新工程和基本科研业务费专项项目(ICS2020YJ01BX)

Received: 2020-11-15   Revised: 2021-01-9   Online: 2021-04-15

作者简介 About authors

王思宇,主要从事地理信息系统和农业遥感研究,E-mail： wsy_427@163.com

摘要

光合有效辐射吸收比率（fraction of absorbed photosynthetically active radiation,FAPAR）是反映作物产量的重要参数之一。无人机遥感能够快速无损地获取高分辨率植被冠层光谱信息,已成为进行物理化参数反演的重要手段。以不同播期玉米为研究对象,基于无人机搭载多光谱传感器,提取植被指数与植被纹理特征,使用偏最小二乘（partial least squares regression,PLSR）方法将二者结合反演玉米FAPAR,并与传统单独使用植被指数或植被纹理特征反演植被FAPAR的方法进行比较。结果表明：使用传统方法单独利用植被指数反演FAPAR（验证RMSE最低为7.33×10-2,rRMSE最低为8.66%）的效果比单独利用纹理特征反演FAPAR（验证RMSE最低为9.50×10-2,rRMSE最低为11.23%）的精度更高;使用PLSR方法单独利用植被指数或纹理特征估算FAPAR的效果比传统方法精度更高（植被指数与纹理特征的验证RMSE最低分别为6.77×10-2和5.24×10-2,rRMSE最低分别为8.01%和6.19%）;使用PLSR方法将植被指数与纹理特征相结合估算FAPAR（验证RMSE最低为4.72×10-2,rRMSE最低为5.57%）的效果比单独使用植被指数或纹理特征估算FAPAR的精度更高。综上,使用PLSR方法将植被指数和植被纹理特征相结合来反演玉米冠层FAPAR可行,为作物FAPAR遥感反演研究提供了新的思路。

关键词： FAPAR ; 多光谱影像 ; 植被指数 ; 纹理特征 ; PLSR

Abstract

The fraction of absorbed photosynthetically active radiation (FAPAR) is one of the most important parameters reflecting the crop yields. The unmanned aerial vehicle (UAV) remote sensing, which can obtain high-resolution vegetation canopy spectral information quickly and nondestructively, has become an important method of inverting the physicochemical parameters of crops. Taking maize with different sowing dates as the research object, this research extracted vegetation indexes (VIs) and texture features based on the UAV multispectral images. After that, a partial least squares regression (PLSR) method was used to invert maize FAPAR combining the two indexes, and compared with the traditional methods using VIs or texture features alone. The results showed that the accuracy of FAPAR inversion by using VIs alone (the validated RMSE is as low as 0.0733, and the validated rRMSE is as low as 8.66%) was higher than that by using texture feature alone (the validated RMSE is as low as 0.0950, and the validated rRMSE is as low as 11.23%). Moreover, the PLSR method was more accurate than the traditional methods in estimating FAPAR by VIs or texture features alone (the validated RMSE of VIs and texture features is as low as 0.0677 and 0.0524, and the validated rRMSE is as low as 8.01% and 6.19%). The accuracy of combining VIs and texture features using the PLSR method (the validated RMSE is as low as 0.0472, and the rRMSE is as low as 5.57%) was higher than that using VIs or texture features alone. The research indicated that it is feasible to invert FAPAR of maize canopy by combining VIs and texture features using the PLSR method, which provides a new idea for crop FAPAR remote sensing inversion.

Keywords： FAPAR ; Multispectral images ; Vegetation index ; Texture features ; PLSR

PDF (842KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

王思宇, 聂臣巍, 余汛, 邵明超, 王梓旭, 努热曼古丽·托乎提, 刘亚东, 程明瀚, 官云兰, 金秀良. 结合植被指数与纹理特征的玉米冠层FAPAR遥感估算研究[J]. 作物杂志, 2021, 37(2): 183-190 doi:10.16035/j.issn.1001-7283.2021.02.027

Wang Siyu, Nie Chenwei, Yu Xun, Shao Mingchao, Wang Zixu, Nuremanguli· Tuohuti, Liu Yadong, Cheng Minghan, Guan Yunlan, Jin Xiuliang. Maize Canopy FAPAR Remote Sensing Estimation Combining Vegetation Indexes and Texture Characteristics[J]. Crops, 2021, 37(2): 183-190 doi:10.16035/j.issn.1001-7283.2021.02.027

光合有效辐射吸收比率（fraction of absorbed photosynthetically active radiation,FAPAR）是植被冠层绿色部分所吸收的光合有效辐射（photosynthetically active radiation,PAR）占到达冠层顶部PAR的比例,是直接反映植被冠层对光能的截获能力与吸收能力的重要参数[1],对作物产量形成和籽粒品质都有着直接影响[2],同时也是基于光能利用效率模型估算作物籽粒产量的重要参数之一[3,4]。因此,实时且精准监测玉米FAPAR对筛选玉米品种和指导田间栽培管理,进而实现玉米优质高产具有重要意义。

目前,基于遥感数据估算FAPAR的方法大多是通过建立植被指数与FAPAR的数学模型进行回归统计[5],这种方法操作简便、参数少且效率高,在FAPAR的估算中应用最为广泛[6,7,8,9,10]。如Dong 等[11]基于模拟Sentinel-2数据选择叶绿素相关植被指数,检验其对小麦和玉米的FAPAR估算能力,结果表明使用近红外和红边反射率的植被指数对FAPAR的估算结果最好。刘桂鹏等[12]基于高光谱影像使用几种植被指数和高光谱反射率及其导数等进行了春玉米FAPAR估算,结果表明使用归一化植被指数（normalized difference vegetation index,NDVI）构建的二次多项式模型对FAPAR的估算效果最好,但是当FAPAR较大时其估算结果会受到影响。田春燕等[13]基于植被指数构建棉花FAPAR估算模型,用以监测棉花长势并预测产量,结果表明使用增强植被指数（enhanced vegetation index,EVI）反演FAPAR的精度最好。Qin等[14]基于机载激光雷达和高光谱数据建立了多元线性回归（multiple linear regression,MLR）模型来估算玉米FAPAR,结果表明,将激光雷达和高光谱指标相结合,比单独使用激光雷达或高光谱指标估算玉米FAPAR的精度更好。Zhao等[15]使用13种植被指数提出了基于二次函数的回归模型来估算玉米FAPAR,结果表明,结合使用冠层结构相关的植被指数与土壤调节、叶绿素和生理相关的植被指数可以更好地估算玉米FAPAR。以上研究表明,采用植被指数反演FAPAR的研究已经取得较好的结果,但该方法在冠层覆盖度较大时使用会存在一定的饱和现象,导致反演精度降低。

植被的纹理特征也可以用来表征植被冠层结构[16],因此也有一些研究通过纹理特征来反演生物量[17]、叶绿素含量[18]和冠层高度[19]等植被结构或冠层信息。如谢士琴等[20]基于SPOT5影像通过使用植被纹理等多个特征反演森林结构参数,结果表明使用纹理等特征可以较好地对森林结构参数进行估算。蔡文婷等[21]将Sentinel-2影像与冬小麦纹理信息相结合估算茬覆盖度,结果表明该方法提高了估算冬小麦茬覆盖度的精度。贾丹等[22]利用无人机影像结合光谱与纹理信息来反演冬小麦氮浓度,结果表明该方法比单一使用光谱或纹理信息估算的精度有所提高。蒙诗栎等[23]基于WorldView-2影像将植被指数与纹理信息相结合反演森林地上生物量,结果表明将二者相结合的反演方法精度最高。杨俊等[24]通过无人机影像将小麦植被指数与纹理特征相结合估算其生物量与产量,结果表明结合纹理特征与植被指数可以提高生物量与产量的估算精度。纹理特征可以反映植被冠层覆盖情况,冠层覆盖度大时到达冠层底部的光较少,FAPAR也会相应变大。因此,当光谱特征不足以反映FAPAR情况时,纹理特征也可以提供植被冠层的结构信息[25]。以上研究表明,可以使用植被纹理特征反演植被结构及冠层参数,且将植被纹理特征与光谱特征相结合可以在一定程度上改善植被光谱特征的饱和问题[26,27],但目前利用植被纹理特征反演植被FAPAR的研究较少。

因此本研究使用PLSR方法,基于无人机多光谱影像提取出的植被指数与植被纹理特征估算玉米FAPAR,并与传统直接使用植被指数或纹理特征估算FAPAR的方法进行比较。

1 材料与方法

1.1 研究区概况

试验于2020年7月-9月在中国农业科学院新乡综合基地（35.2°N,113.8°E）进行。河南新乡的气候为温带大陆性气候,四季分明,冬冷夏热。年均气温14.1℃,1月最冷,7月最热。年均降水量548.3mm,多集中在7-8月。年均蒸发量1908.7mm。年均日照2407.7h,年均无霜期200.5d。试验基地土壤为黏质壤土,耕层土壤pH 8.21,含有机质12.6g/kg、速效氮61.2mg/kg、速效磷16.2mg/kg、速效钾110.0mg/kg。

1.2 试验设计

研究对象为夏玉米,共设置8个播期,获取关键时期的玉米分层FAPAR数据,每个播期间隔10d左右,播种时间从4月下旬至7月上旬,播期1的播种时间最早,播期8的播种时间最晚,其中播期6为当地适播期。每个播期包含3个黄淮海地区主栽玉米品种,分别为丰垦139、京农科728和郑单958,每个品种设置3个重复。每个播期相同品种面积之和为64.8m2,田间水肥与病虫害管理与实际大田管理一致。

1.3 地面FAPAR数据获取

使用SunScan专业版植物冠层分析仪进行采集地面分层FAPAR数据,在2020年7月-9月有稳定太阳光直射时进行采集。获取数据前架设漫射系数传感器（beam fraction sensor,BFS）,使其水平以减小误差,BFS上方不能有任何遮挡,与探测器通过无线方式进行连接,探测器与掌上电脑（PAD）相连;之后打开SunData软件将BFS与探测器保持在同一水平面上对仪器进行校正和数据采集。获取数据时选取长势均匀的区域,将探测器水平放在冠层最底部进行数据采集。每个小区采用米字型方式测4个重复,然后将4次测量结果的均值作为该小区的代表值,使用PAD记录数据。获取数据后将数据导出进行后续处理。实测的地面FAPAR数据详细信息如表1所示。

表1   地面测量FAPAR数据统计信息

Table 1  Statistics of in situ FAPAR

日期
Date		品种
Variety		样本数
Sample number		取值范围(×10-2)
Value range
2020-07-23丰垦1392151.62~92.49
京农科7282150.07~93.82
郑单9582161.20~94.15
2020-08-02丰垦1392438.25~94.19
京农科7282444.10~97.42
郑单9582446.05~97.43

新窗口打开| 下载CSV


1.4 多光谱影像数据的获取与处理

多光谱数据使用大疆M600-PRO六旋翼无人机搭载的RedEdge-M多光谱相机进行获取。多光谱相机获取的波段有蓝、绿、红、红边和近红外,视场角47.2°,数据采样间隔1.5s。无人机飞行高度30m,设置航向重叠和旁向重叠为80%。无人机数据采集选择在晴朗无风的10∶00-15∶00之间进行,获取时间与地面FAPAR测量日期相同。

1.4.1 多光谱影像预处理 多光谱影像的预处理主要为影像的拼接与辐射校正两部分。拼接部分使用AgisoftPhtotscan Pro软件进行处理,首先将需要拼接的照片添加到软件中之后对影像进行筛选,剔除冗余照片以缩短处理时间;然后选择高精度对齐照片并构建密集点云;之后进行颜色校正并构建纹理与正射影像;最后将构建的正射影像以TIFF形式导出。辐射校正是将原始影像转换为反射率影像,首先使用地面ASD高光谱仪测得彩色定标布的高光谱反射率Rcolor_i,通过多光谱各通道的波段响应函数Fband_i计算地面定标布对应多光谱各波段的反射率（公式1）;然后提取无人机多光谱影像中彩色定标布的像元值（digital number,DN）,建立影像每波段DN值与反射率的回归关系,之后将各波段的回归关系应用到整个试验地块的无人机影像,最后得到校正后的反射率影像。

Rcolor_band_n=∑(Rcolor_i×Fband_i)/∑Fband_i

1.4.2 植被指数的选择与提取 使用ENVI软件对处理好的影像提取所需要播期单个小区的感兴趣区（region of interest,ROI）,并进行统计,最后得出小区中每个波段的平均反射率。

根据前人[28,29,30,31,32,33,34,35,36,37,38,39,40,41]研究,从文献中选取了14个用于估算FAPAR较为常用的植被指数（表2）,根据公式利用每个波段的反射率计算出14种植被指数。

表2   选取的14种植被指数及其计算公式

Table 2  Fourteen vegetation indices and their formulas

新窗口打开| 下载CSV


1.4.3 纹理信息的提取 根据传统纹理特征提取方法,选择多光谱影像的R、G和B三个波段使用基于二阶概率统计的滤波（co-occurrence measures）方法提取植被纹理特征,根据前人研究结果[42]并进行比较分析,发现当滤波窗口大小设为7×7时,处理速度快且能够保证提取纹理的精度。因此使用7×7窗口对纹理特征进行提取,空间相关矩阵偏移X和Y默认为1,最后输出得到的8种纹理特征如表3所示。

表3   选取的8种纹理特征

Table 3  Features of the eight textures

新窗口打开| 下载CSV


表3中,Pi,j表示灰度共生矩阵中第i行j列对应的两个灰度值同时出现的概率;N为图像的灰度级;μi和μj为灰度共生矩阵的均值;σi和σj为标准差[43]。

1.5 模型精度评定

随机选取样本的2/3作为训练集,剩余1/3样本作为验证集。选取植被指数或纹理特征分别与实测FAPAR的决定系数R2、估算FAPAR与实测FAPAR验证的均方根误差RMSE和相对均方根误差rRMSE作为评价这些植被指数或纹理特征性能的指标。R2越接近1,则该模型的拟合精度越高;RMSE与rRMSE越小,表示估算值与实测值差异越小,模型的估算能力越好,其计算公式如下：

$R^{2}=\frac{\sum_{i=1}^{n}\left(x_{i}-\bar{x}\right)^{2} \times\left(y_{i}-\bar{y}\right)^{2}}{\sum_{i=1}^{n}\left(x_{i}-\bar{x}\right)^{2} \times \sum_{i=1}^{n}\left(y_{i}-\bar{y}\right)^{2}}$
$R M S E=\sqrt{\frac{\sum_{i=1}^{n}\left(y_{i}-x_{i}\right)^{2}}{n}}$
$rRMSE=\frac{RMSE}{\bar{X}}$

其中,xi和yi分别为实测值和预测值;x和y分别为实测值和预测值的均值;n为样本个数。

2 结果与分析

2.1 利用植被指数估算FAPAR

使用表2中选取的14种植被指数分别与FAPAR建立指数函数、线性函数、二次多项式函数、对数函数和幂函数5种统计回归模型,得到每个植被指数对FAPAR估算的最佳回归检验结果,如表4所示。

表4   各植被指数估算玉米FAPAR的最佳回归检验结果

Table 4  The best regression test results of maize FAPAR estimating using each vegetation indice

植被指数
Vegetation indice		回归模型
Regression model		建模Modeling验证Validating
R2(×10-2)RMSE(×10-2)rRMSE(%)R2(×10-2)RMSE(×10-2)rRMSE(%)
DVIy=-8.8039x2+6.3709x-0.247666.137.989.4266.359.1010.75
GNDVIy=-1.5792x2+2.4437x-0.039178.526.357.49278.007.338.66
MTCIy=-0.1751x2+0.4515x+0.624861.138.5410.0870.637.909.34
TVIy=0.5125ln(x)-0.693456.489.0410.6757.3910.0511.88
MTVI2y=-1.2913x2+2.3547x-0.149371.747.288.5966.939.1710.84
OSAVIy=1.169x1.006972.187.248.5472.187.809.22
NDVIy=0.9303x1.244973.747.038.3074.977.598.97
RVI1y=0.0497ln(x)+0.63633.9511.1413.1437.5311.1413.16
RVI2y=-0.0024x2+0.062x+0.565144.5510.2012.0452.069.8011.58
MNLIy=-5.4349x2+2.9606x+0.507177.836.457.6175.598.5910.15
SAVIy=1.3688x0.820866.587.939.3667.878.219.70
MSRy=0.0908ln(x)+0.67240.9810.5312.4246.4510.3112.18
NLIy=0.3779x+0.551574.546.918.1574.597.428.77
RDVIy=1.487x0.905268.167.749.1369.258.089.55

新窗口打开| 下载CSV


从表4综合比较可以看出,14种植被指数中有4种建模R2达到0.7以上,验证RMSE小于0.08,rRMSE小于10%;其中GNDVI的二次多项式模型对FAPAR的估算结果最好,R2为78.52×10-2,RMSE

为7.33×10-2,rRMSE为8.66%。原因可能是由于GNDVI用绿波段代替了NDVI中的红波段,而绿波段对叶绿素浓度更为敏感,可以在很大范围内精确评估色素浓度的变化[29],使用绿波段建立的植被指数加强了对FAPAR估算结果的稳定性。图1为GNDVI与FAPAR的回归模型。

图1

新窗口打开| 下载原图ZIP| 生成PPT

图1   玉米FAPAR与GNDVI的回归关系

Fig.1   Regression relationship of the FAPAR to GNDVI in maize


2.2 利用纹理信息估算FAPAR

提取多光谱影像R、G和B波段的8种纹理特征之后,分别与FAPAR建立指数函数、线性函数、二次多项式函数、对数函数和幂函数5种统计回归模型,综合比较之后发现绿波段对纹理特征的反演效果最好,因此使用绿波段中的纹理特征对玉米FAPAR进行估算,每种纹理特征的最佳回归检验结果如表5所示。

表5   各纹理信息估算玉米FAPAR的最佳回归检验结果

Table 5  The best regression test results of maize FAPAR estimating using each texture

纹理特征
Texture		回归模型
Regression model		建模Modeling验证Validating
R2(×10-2)RMSE(×10-2)rRMSE(%)R2(×10-2)RMSE(×10-2)rRMSE(%)
均值Meany=-0.0031x2+0.0493x+0.716958.738.8010.3855.539.5011.23
方差Variancey=-0.0003x2+0.0192x+0.64418.4412.3814.6125.0012.2814.51
协同性Homogeneityy=-3.1473x2+1.8839x+0.60075.1213.3515.7519.2713.0615.43
对比度Contrasty=-0.0005x2+0.0235x+0.654216.4212.5314.7822.5612.5014.77
相异性Dissimilarityy=-0.0199x2+0.1895x+0.518411.8512.8715.1919.6512.7915.11
信息熵Entropyy=-2.1224x2+14.464x-23.75612.5312.8215.1336.2311.4513.53
二阶矩Second momenty=-399x2+38.016x-0.02886.7513.2315.6139.6412.2614.49
相关性Correlationy=-14.067x2+19.182x-5.625421.6812.1314.3129.3911.9414.11

新窗口打开| 下载CSV


综合比较可以看出,Mean的二次多项式模型对FAPAR的估算结果最好,R2为58.73×10-2,RMSE为9.50×10-2,rRMSE为11.23%。其原因可能是由于试验田纹理较规则,规则的纹理特征能够更好地体现植被冠层覆盖度情况,覆盖度均匀则光的透过率也均匀,与FAPAR的相关性较高,进而能够更好地反映FAPAR,而Mean特征反映了植被纹理的规则程度;同时提取纹理信息的窗口大小和研究区植被类型也会影响模型的精度[23],因此,综合试验区作物类型与纹理提取情况,Mean特征最适用于对本研究区域中纹理信息的表达,能够更好地反演玉米FAPAR。图2为纹理特征Mean与FAPAR的回归模型。

图2

新窗口打开| 下载原图ZIP| 生成PPT

图2   玉米FAPAR与纹理特征Mean的回归关系

Fig.2   Regression relationship of the FAPAR to mean of texture feature in maize


2.3 利用PLSR方法估算FAPAR

使用PLSR方法分别对14种植被指数、8种纹理特征以及将二者相结合估算玉米FAPAR,并进行验证,结果如表6所示。

表6   利用PLSR方法估算玉米FAPAR的验证结果

Table 6  Validation results of FAPAR estimation using PLSR method

指标
Indice		建模Modeling验证Validating
R2
(×10-2)		RMSE
(×10-2)		rRMSE
(%)		R2
(×10-2)		RMSE
(×10-2)		rRMSE
(%)
植被指数
Vegetation indice		82.158.199.6681.406.778.01
纹理Texture89.126.397.5487.055.246.19
结合Combination94.394.595.4190.814.725.57

新窗口打开| 下载CSV


可以看出,对14种植被指数和8种纹理特征分别使用PLSR方法估算玉米FAPAR之后,其R2均有所提升,RMSE和rRMSE均有所下降,将所有指标相结合的PLSR估算结果最好,R2可达到94.39×10-2,RMSE和rRMSE分别为4.72×10-2和5.57%,相比传统回归统计方法有所提高,该结果说明使用PLSR方法可以提高玉米FAPAR估算的精度,将植被指数与植被纹理特征相结合的效果更好。这主要是由于植被的光谱信息和纹理特征都可以反映植被冠层的信息,而对于多元回归模型来说,加入有效的自变量可以提高构建模型的精度,因此将二者相结合构建多元回归模型,可以将植被指数和纹理特征的优势相融合,从多个角度充分对玉米冠层FAPAR进行反演,能够提高其估算精度,这与前人对植物其他指标进行估算研究的结果具有一致性[21-24,44]。图3为使用PLSR方法结合植被指数与植被纹理特征估算FAPAR的验证结果。

图3

新窗口打开| 下载原图ZIP| 生成PPT

图3   使用PLSR方法结合植被指数与纹理特征估算的FAPAR验证结果

Fig.3   Validation results of FAPAR estimation using PLSR method combined with vegetation index and texture features


3 讨论

FAPAR不仅是植被冠层的重要结构参数,也是作物生长和作物估产等模型的重要参数,因此快速准确地获取作物FAPAR对于估测作物产量及评价作物生长状况有重要意义。由于单独使用光谱特征对FAPAR进行估算的精度相对较差,而纹理特征提供了更多的植被冠层结构信息,能够提高FAPAR的估算精度。因此,本研究将光谱特征与纹理特征相结合,给FAPAR的估算模型增加了更多信息,在一定程度上提高了单独使用植被指数或纹理特征估算FAPAR的精度,为以后估算作物的FAPAR研究提供了新的思路。然而,本研究只选择了1个试验站点采集的FAPAR数据,没有在其他站点进行试验,所使用的试验数据较少,建立的回归模型结果可能不够精确;此外,由于错过了一些播期玉米的生长状况,使得处在生育初期的FAPAR数据相对较少,对建模结果也有一定的影响。因此,在接下来的研究中,将开展多站点试验,同时获取更多的田间及无人机影像数据,以更好地阐明光谱信息与FAPAR数据之间的光谱机理与回归关系。

4 结论

对玉米冠层FAPAR进行测量并结合基于无人机多光谱影像提取的植被指数和纹理特征分别对其进行估算,之后使用PLSR方法对植被指数、纹理特征和将二者结合分别对FAPAR进行估算。将3种估算FAPAR的方法进行比较验证。结果如下：

（1）使用传统方法直接利用植被指数或纹理特征估算FAPAR时,植被指数的精度更高,R2比纹理特征提高约0.2,RMSE和rRMSE分别降低约0.02和2.60%;

（2）单独利用植被指数或纹理特征估算FAPAR时,使用PLSR方法比传统的统计回归方法精度更高。其中植被指数的R2提高约0.04,RMSE和rRMSE分别降低约0.06×10-2和0.70%,纹理特征的RMSE和rRMSE分别降低约4.80×10-2和5.04%;

（3）将植被指数和纹理特征相结合并使用PLSR估算FAPAR的精度最高,R2达到94.39×10-2,RMSE和rRMSE分别达到4.72×10-2和5.57%。

参考文献

董泰锋, 蒙继华, 吴炳方.

基于遥感的光合有效辐射吸收比率(FAPAR)估算方法综述

生态学报, 2012,32(22):7190-7201.

[本文引用: 1]

周元刚.

不同受旱条件下冬小麦和夏玉米冠层光截获特征研究

杨凌:西北农林科技大学, 2016.

[本文引用: 1]

李贺丽, 罗毅, 薛晓萍, 等.

冬小麦冠层对入射光合有效辐射吸收比例的估算方法评价

农业工程学报, 2011,27(4):201-206.

[本文引用: 1]

Moriondo M, Maselli F, Bindi M.

A simple model of regional wheat yield based on NDVI data

European Journal of Agronomy, 2007,26(3):266-274.

[本文引用: 1]

张宇.

基于光学和雷达遥感的水稻FAPAR协同反演

成都:电子科技大学, 2019.

[本文引用: 1]

Zhang Y, Li S, He Z, et al.

Synergistic Inversion of Rice FPAR Based on Optical and Radar Remote Sensing Data. IGARSS 2019-2019 IEEE International Geoscience and Remote Sensing Symposium

IEEE, 2019.

[本文引用: 1]

王方永, 韩焕勇, 陈兵, 等.

数字图像估测棉花光合有效辐射吸收比例

西北农业学报, 2015,24(10):74-79.

[本文引用: 1]

Majasalmi T, Stenberg P, Rautiainen M.

Comparison of ground and satellite-based methods for estimating stand-level fPAR in a boreal forest

Agricultural and Forest Meteorology, 2017,232:422-432.

[本文引用: 1]

Liang S, Hou X, Sui X, et al.

Deciduous broadleaf forests green FPAR and its relationship with spectral vegetation indices. 2017 10th International Congress on Image and Signal Processing,BioMedical Engineering and Informatics (CISP-BMEI)

IEEE, 2017.

[本文引用: 1]

Tan C W, Wang D L, Zhou J, et al.

Remotely assessing fraction of photosynthetically active radiation (FAPAR) for wheat canopies based on hyperspectral vegetation indexes

Frontiers in Plant Science, 2018,9:776.

[本文引用: 1]

Dong T, Meng J, Shang J, et al.

Evaluation of chlorophyll-related vegetation indices using simulated sentinel-2 data for estimation of crop fraction of absorbed photosynthetically active radiation

IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2015,8(8):4049-4059.

[本文引用: 1]

刘桂鹏, 贺婷, 王国骄, 等.

玉米LAI和FAPAR的高光谱遥感预测模型研究

玉米科学, 2016,24(2):115-119,128.

[本文引用: 1]

田春燕, 黄春燕, 郭晓飞, 等.

基于植被指数的棉花冠层光合有效辐射截获量和叶片净光合速率估算研究

新疆农业科学, 2017,54(6):981-987.

[本文引用: 1]

Qin H M, Wang C, Zhao K G, et al.

Estimation of the fraction of absorbed photosynthetically active radiation (FAPAR) in maize canopies using LiDAR data and hyperspectral imagery

PLoS ONE, 2018,13(5):e0197510.

[本文引用: 1]

Zhao L, Liu Z G, Xu S, et al.

Retrieving the diurnal FAPAR of a maize canopy from the jointing stage to the tasseling stage with vegetation indices under different water stresses and light conditions

Sensors, 2018,18(11):3965.

[本文引用: 1]

万亮, 岑海燕, 朱姜蓬, 等.

基于纹理特征与植被指数融合的水稻含水量无人机遥感监测

智慧农业(中英文), 2020,2(1):58-67.

[本文引用: 1]

张华, 张玉红, 张改改.

民勤绿洲青土湖植被优势种地上生物量估算

干旱区地理, 2020,43(1):201-210.

[本文引用: 1]

陈鹏, 冯海宽, 李长春, 等.

无人机影像光谱和纹理融合信息估算马铃薯叶片叶绿素含量

农业工程学报, 2019,35(11):63-74.

[本文引用: 1]

高凯旋, 焦海明, 王新闯.

融合影像纹理、光谱与地形特征的森林冠顶高反演模型

国土资源遥感, 2020,32(3):63-70.

[本文引用: 1]

谢士琴, 赵天忠, 王威, 等.

结合影像纹理、光谱与地形特征的森林结构参数反演

农业机械学报, 2017,48(4):125-134.

[本文引用: 1]

蔡文婷, 赵书河, 王亚梅, 等.

结合Sentinel-2光谱与纹理信息的冬小麦作物茬覆盖度估算

遥感学报, 2020,24(9):1108-1119.

[本文引用: 2]

贾丹, 陈鹏飞.

低空无人机影像分辨率对冬小麦氮浓度反演的影响

农业机械学报, 2020,51(7):164-169.

[本文引用: 1]

蒙诗栎, 庞勇, 张钟军, 等.

WorldView-2纹理的森林地上生物量反演

遥感学报, 2017,21(5):812-824.

[本文引用: 2]

杨俊, 丁峰, 陈晨, 等.

小麦生物量及产量与无人机图像特征参数的相关性

农业工程学报, 2019,35(23):104-110.

[本文引用: 2]

Wulder M A, LeDrew.

Aerial image texture information in the estimation of northern deciduous and mixed wood forest leaf area index (LAI)

Remote Sensing of Environment, 1998,64(1):64-76.

[本文引用: 1]

Cen H, Wan L, Zhu J, et al.

Dynamic monitoring of biomass of rice under different nitrogen treatments using a lightweight UAV with dual image-frame snapshot cameras

Plant Methods, 2019,15(1):32.

[本文引用: 1]

Wan L, Li Y, Cen H, et al.

Combining UAV-based vegetation indices and image classification to estimate flower number in oilseed rape

Remote Sensing, 2018,10(9):1484.

[本文引用: 1]

Tucker C J.

Red and photographic infrared linear combinations for monitoring vegetation

Remote Sensing of Environment, 1979,8(2):127-150.

[本文引用: 1]

Anatoly A G, Merzlyak M N.

Remote sensing of chlorophyll concentration in higher plant leaves

Advances in Space Research, 1998,22(5):689-692.

[本文引用: 2]

Dash J, Curran P J.

The MERIS terrestrial chlorophyll index

International Journal of Remote Sensing, 2004,25(23):5403-5413.

[本文引用: 1]

Broge N H, Leblanc E.

Comparing prediction power and stability of broadband and hyperspectral vegetation indices for estimation of green leaf area index and canopy chlorophyll density

Remote Sensing of Environment, 2001,76(2):156-172.

[本文引用: 1]

Haboudane D, Miller J R, Pattey E, et al.

Hyperspectral vegetation indices and novel algorithms for predicting green LAI of crop canopies: modeling and validation in the context of precision agriculture

Remote Sensing of Environment, 2004,90(3):337-352.

[本文引用: 1]

Rondeaux G, Steven M, Baret F.

Optimization of soil-adjusted vegetation indices

Remote Sensing of Environment, 1996,55(2):95-107.

[本文引用: 1]

蔡薇.

基于NDVI时间序列的荷泽市东明县植被动态变化研究

经济研究导刊, 2009(18):136-137.

[本文引用: 1]

陈拉, 黄敬峰, 王秀珍.

不同传感器的模拟植被指数对水稻叶面积指数的估测精度和敏感性分析

遥感学报, 2008(1):143-151.

[本文引用: 1]

Xue L, Cao W, Luo W, et al.

Monitoring leaf nitrogen status in rice with canopy spectral reflectance

Agronomy Journal, 2004,96(1):135-142.

[本文引用: 1]

Gong P, Pu R, Biging G S, et al.

Estimation of forest leaf area index using vegetation indices derived from hyperion hyperspectral data

IEEE Transactions on Geoscience and Remote Sensing, 2003,41(6):1355-1362.

[本文引用: 1]

湛青青, 王辉源.

基于综合指数分析的城市建筑用地信息提取研究

测绘与空间地理信息, 2014,37(2):62-65,69.

[本文引用: 1]

Chen J M.

Evaluation of vegetation indices and a modified simple ratio for boreal applications

Canadian Journal of Remote Sensing, 1996,22(3):229-242.

[本文引用: 1]

Goel N S, Qin W H.

Influences of canopy architecture on relationships between various vegetation indices and LAI and FPAR: A computer simulation

Remote Sensing Reviews, 1994,10(4):309-347.

[本文引用: 1]

Roujean J L, Breon F M.

Estimating PAR absorbed by vegetation from bidirectional reflectance measurements

Remote Sensing of Environment, 1995,51(3):375-384.

[本文引用: 1]

Sun Q, Sun L, Shu M, et al.

Monitoring maize lodging grades via unmanned aerial vehicle multispectral image

Plant Phenomics, 2019(2):1-16.

[本文引用: 1]

刘思涵.

基于高分辨率影像纹理特征的落叶松人工林提取研究

西安:西安科技大学, 2020.

[本文引用: 1]

刘畅, 杨贵军, 李振海, 等.

融合无人机光谱信息与纹理信息的冬小麦生物量估测

中国农业科学, 2018,51(16):3060-3073.

[本文引用: 1]

/

〈 〉