作物杂志,2021, 第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. 1东华理工大学测绘工程学院,330013,江西南昌
    2中国农业科学院作物科学研究所,100081,北京
    3河南理工大学,454003,河南焦作
    4河北工程大学地球科学与工程学院,056006,河北邯郸
    5长安大学地球科学与资源学院,710061,陕西西安
    6中国地质大学(武汉)地球科学学院,430074,湖北武汉
    7河海大学农业工程学院,210098,江苏南京
  • 收稿日期:2020-11-15 修回日期:2021-01-09 出版日期:2021-04-15 发布日期:2021-04-16
  • 通讯作者: 官云兰,金秀良
  • 作者简介:王思宇,主要从事地理信息系统和农业遥感研究,E-mail: wsy_427@163.com
  • 基金资助:
    国家自然科学基金-面上项目(42071426);中国农业科学院基本科研业务费专项院级统筹项目(Y2020YJ07);中国农业科学院科技创新工程和基本科研业务费专项项目(ICS2020YJ01BX)

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

Wang Siyu1,2(), Nie Chenwei2, Yu Xun2,3, Shao Mingchao2,4, Wang Zixu2,5, Nuremanguli· Tuohuti2,6, Liu Yadong2, Cheng Minghan2,7, Guan Yunlan1(), Jin Xiuliang2()   

  1. 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
  • Received:2020-11-15 Revised:2021-01-09 Online:2021-04-15 Published:2021-04-16
  • Contact: Guan Yunlan,Jin Xiuliang

摘要:

光合有效辐射吸收比率(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.

Key words: FAPAR, Multispectral images, Vegetation index, Texture features, PLSR

表1

地面测量FAPAR数据统计信息

日期
Date
品种
Variety
样本数
Sample number
取值范围(×10-2)
Value range
2020-07-23 丰垦139 21 51.62~92.49
京农科728 21 50.07~93.82
郑单958 21 61.20~94.15
2020-08-02 丰垦139 24 38.25~94.19
京农科728 24 44.10~97.42
郑单958 24 46.05~97.43

表2

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

表3

选取的8种纹理特征

表4

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

植被指数
Vegetation indice
回归模型
Regression model
建模Modeling 验证Validating
R2(×10-2) RMSE(×10-2) rRMSE(%) R2(×10-2) RMSE(×10-2) rRMSE(%)
DVI y=-8.8039x2+6.3709x-0.2476 66.13 7.98 9.42 66.35 9.10 10.75
GNDVI y=-1.5792x2+2.4437x-0.0391 78.52 6.35 7.492 78.00 7.33 8.66
MTCI y=-0.1751x2+0.4515x+0.6248 61.13 8.54 10.08 70.63 7.90 9.34
TVI y=0.5125ln(x)-0.6934 56.48 9.04 10.67 57.39 10.05 11.88
MTVI2 y=-1.2913x2+2.3547x-0.1493 71.74 7.28 8.59 66.93 9.17 10.84
OSAVI y=1.169x1.0069 72.18 7.24 8.54 72.18 7.80 9.22
NDVI y=0.9303x1.2449 73.74 7.03 8.30 74.97 7.59 8.97
RVI1 y=0.0497ln(x)+0.636 33.95 11.14 13.14 37.53 11.14 13.16
RVI2 y=-0.0024x2+0.062x+0.5651 44.55 10.20 12.04 52.06 9.80 11.58
MNLI y=-5.4349x2+2.9606x+0.5071 77.83 6.45 7.61 75.59 8.59 10.15
SAVI y=1.3688x0.8208 66.58 7.93 9.36 67.87 8.21 9.70
MSR y=0.0908ln(x)+0.672 40.98 10.53 12.42 46.45 10.31 12.18
NLI y=0.3779x+0.5515 74.54 6.91 8.15 74.59 7.42 8.77
RDVI y=1.487x0.9052 68.16 7.74 9.13 69.25 8.08 9.55

图1

玉米FAPAR与GNDVI的回归关系

表5

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

纹理特征
Texture
回归模型
Regression model
建模Modeling 验证Validating
R2(×10-2) RMSE(×10-2) rRMSE(%) R2(×10-2) RMSE(×10-2) rRMSE(%)
均值Mean y=-0.0031x2+0.0493x+0.7169 58.73 8.80 10.38 55.53 9.50 11.23
方差Variance y=-0.0003x2+0.0192x+0.644 18.44 12.38 14.61 25.00 12.28 14.51
协同性Homogeneity y=-3.1473x2+1.8839x+0.6007 5.12 13.35 15.75 19.27 13.06 15.43
对比度Contrast y=-0.0005x2+0.0235x+0.6542 16.42 12.53 14.78 22.56 12.50 14.77
相异性Dissimilarity y=-0.0199x2+0.1895x+0.5184 11.85 12.87 15.19 19.65 12.79 15.11
信息熵Entropy y=-2.1224x2+14.464x-23.756 12.53 12.82 15.13 36.23 11.45 13.53
二阶矩Second moment y=-399x2+38.016x-0.0288 6.75 13.23 15.61 39.64 12.26 14.49
相关性Correlation y=-14.067x2+19.182x-5.6254 21.68 12.13 14.31 29.39 11.94 14.11

图2

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

表6

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

指标
Indice
建模Modeling 验证Validating
R2
(×10-2)
RMSE
(×10-2)
rRMSE
(%)
R2
(×10-2)
RMSE
(×10-2)
rRMSE
(%)
植被指数
Vegetation indice
82.15 8.19 9.66 81.40 6.77 8.01
纹理Texture 89.12 6.39 7.54 87.05 5.24 6.19
结合Combination 94.39 4.59 5.41 90.81 4.72 5.57

图3

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

[1] 董泰锋, 蒙继华, 吴炳方. 基于遥感的光合有效辐射吸收比率(FAPAR)估算方法综述. 生态学报, 2012,32(22):7190-7201.
[2] 周元刚. 不同受旱条件下冬小麦和夏玉米冠层光截获特征研究. 杨凌:西北农林科技大学, 2016.
[3] 李贺丽, 罗毅, 薛晓萍, 等. 冬小麦冠层对入射光合有效辐射吸收比例的估算方法评价. 农业工程学报, 2011,27(4):201-206.
[4] 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.
[5] 张宇. 基于光学和雷达遥感的水稻FAPAR协同反演. 成都:电子科技大学, 2019.
[6] 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.
[7] 王方永, 韩焕勇, 陈兵, 等. 数字图像估测棉花光合有效辐射吸收比例. 西北农业学报, 2015,24(10):74-79.
[8] 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.
[9] 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.
[10] 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.
[11] 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.
[12] 刘桂鹏, 贺婷, 王国骄, 等. 玉米LAI和FAPAR的高光谱遥感预测模型研究. 玉米科学, 2016,24(2):115-119,128.
[13] 田春燕, 黄春燕, 郭晓飞, 等. 基于植被指数的棉花冠层光合有效辐射截获量和叶片净光合速率估算研究. 新疆农业科学, 2017,54(6):981-987.
[14] 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.
[15] 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.
[16] 万亮, 岑海燕, 朱姜蓬, 等. 基于纹理特征与植被指数融合的水稻含水量无人机遥感监测. 智慧农业(中英文), 2020,2(1):58-67.
[17] 张华, 张玉红, 张改改. 民勤绿洲青土湖植被优势种地上生物量估算. 干旱区地理, 2020,43(1):201-210.
[18] 陈鹏, 冯海宽, 李长春, 等. 无人机影像光谱和纹理融合信息估算马铃薯叶片叶绿素含量. 农业工程学报, 2019,35(11):63-74.
[19] 高凯旋, 焦海明, 王新闯. 融合影像纹理、光谱与地形特征的森林冠顶高反演模型. 国土资源遥感, 2020,32(3):63-70.
[20] 谢士琴, 赵天忠, 王威, 等. 结合影像纹理、光谱与地形特征的森林结构参数反演. 农业机械学报, 2017,48(4):125-134.
[21] 蔡文婷, 赵书河, 王亚梅, 等. 结合Sentinel-2光谱与纹理信息的冬小麦作物茬覆盖度估算. 遥感学报, 2020,24(9):1108-1119.
[22] 贾丹, 陈鹏飞. 低空无人机影像分辨率对冬小麦氮浓度反演的影响. 农业机械学报, 2020,51(7):164-169.
[23] 蒙诗栎, 庞勇, 张钟军, 等. WorldView-2纹理的森林地上生物量反演. 遥感学报, 2017,21(5):812-824.
[24] 杨俊, 丁峰, 陈晨, 等. 小麦生物量及产量与无人机图像特征参数的相关性. 农业工程学报, 2019,35(23):104-110.
[25] 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.
[26] 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.
[27] 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.
[28] Tucker C J. Red and photographic infrared linear combinations for monitoring vegetation. Remote Sensing of Environment, 1979,8(2):127-150.
[29] Anatoly A G, Merzlyak M N. Remote sensing of chlorophyll concentration in higher plant leaves. Advances in Space Research, 1998,22(5):689-692.
[30] Dash J, Curran P J. The MERIS terrestrial chlorophyll index. International Journal of Remote Sensing, 2004,25(23):5403-5413.
[31] 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.
[32] 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.
[33] Rondeaux G, Steven M, Baret F. Optimization of soil-adjusted vegetation indices. Remote Sensing of Environment, 1996,55(2):95-107.
[34] 蔡薇. 基于NDVI时间序列的荷泽市东明县植被动态变化研究. 经济研究导刊, 2009(18):136-137.
[35] 陈拉, 黄敬峰, 王秀珍. 不同传感器的模拟植被指数对水稻叶面积指数的估测精度和敏感性分析. 遥感学报, 2008(1):143-151.
[36] 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.
[37] 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.
[38] 湛青青, 王辉源. 基于综合指数分析的城市建筑用地信息提取研究. 测绘与空间地理信息, 2014,37(2):62-65,69.
[39] 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.
[40] 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.
[41] Roujean J L, Breon F M. Estimating PAR absorbed by vegetation from bidirectional reflectance measurements. Remote Sensing of Environment, 1995,51(3):375-384.
[42] Sun Q, Sun L, Shu M, et al. Monitoring maize lodging grades via unmanned aerial vehicle multispectral image. Plant Phenomics, 2019(2):1-16.
[43] 刘思涵. 基于高分辨率影像纹理特征的落叶松人工林提取研究. 西安:西安科技大学, 2020.
[44] 刘畅, 杨贵军, 李振海, 等. 融合无人机光谱信息与纹理信息的冬小麦生物量估测. 中国农业科学, 2018,51(16):3060-3073.
[1] 王猛, 隋学艳, 梁守真, 侯学会, 梁永全. 利用无人机遥感技术提取农作物植被覆盖度方法研究[J]. 作物杂志, 2020, (3): 177–183
[2] 王连喜,刘畅,李琪,吴东丽,王清,成威. 气候变暖背景下京津冀地区冬小麦种植北界变化[J]. 作物杂志, 2017, (1): 61–67
[3] 李琪,孙晓宇,王连喜,缪淼,吴东丽. 基于不同植被指数的VSWI在河南省春季干旱监测中的应用分析[J]. 作物杂志, 2016, (1): 162–168
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
No Suggested Reading articles found!