留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

多时相NDVI与丰度综合分析的油菜无人机遥感长势监测

刘雅婷 龚龑 段博 方圣辉 彭漪

刘雅婷, 龚龑, 段博, 方圣辉, 彭漪. 多时相NDVI与丰度综合分析的油菜无人机遥感长势监测[J]. 武汉大学学报 ● 信息科学版, 2020, 45(2): 265-272. doi: 10.13203/j.whugis20180161
引用本文: 刘雅婷, 龚龑, 段博, 方圣辉, 彭漪. 多时相NDVI与丰度综合分析的油菜无人机遥感长势监测[J]. 武汉大学学报 ● 信息科学版, 2020, 45(2): 265-272. doi: 10.13203/j.whugis20180161
LIU Yating, GONG Yan, DUAN Bo, FANG Shenghui, PENG Yi. Combining Multi-temporal NDVI and Abundance from UAV Remote Sensing Data for Oilseed Rape Growth Monitoring[J]. Geomatics and Information Science of Wuhan University, 2020, 45(2): 265-272. doi: 10.13203/j.whugis20180161
Citation: LIU Yating, GONG Yan, DUAN Bo, FANG Shenghui, PENG Yi. Combining Multi-temporal NDVI and Abundance from UAV Remote Sensing Data for Oilseed Rape Growth Monitoring[J]. Geomatics and Information Science of Wuhan University, 2020, 45(2): 265-272. doi: 10.13203/j.whugis20180161

多时相NDVI与丰度综合分析的油菜无人机遥感长势监测

doi: 10.13203/j.whugis20180161
基金项目: 

国家高技术研究发展计划 2013AA102401

中央高校基本科研业务费专项资金 2042017kf0236

详细信息

Combining Multi-temporal NDVI and Abundance from UAV Remote Sensing Data for Oilseed Rape Growth Monitoring

Funds: 

The National High-tech Research and Development Program of China 2013AA102401

the Fundamental Research Funds for the Central Universities 2042017kf0236

More Information
  • 摘要: 从无人机油菜遥感数据精准分析的角度,提出了一种新型的基于多时相混合光谱分析的油菜长势监测方法。根据油菜冠层结构与田间特征确定不同生长期端元组合,并通过线性解混模型获取丰度。针对油菜产量积累过程,分析3个生长期端元丰度与产量的相关性。在此基础上,对比分析典型归一化植被指数(normalized difference vegetation index,NDVI)、端元丰度以及两者综合3种自变量方案下建立的多时相线性估产模型精度。实验结果表明,相比广泛使用的NDVI,选取合适的丰度组合数据结合NDVI进行多时相估产具有更好的效果,且不同种植方式下的最优丰度组合相同,适用于产量的提前预测。由此说明,丰度数据能从植株发育与产量构成角度对NDVI进行亚像元层次的修正,具有较强的稳定性。
  • 图  1  多时相相关分析流程图

    Figure  1.  Flowchart of Multi-temporal Correlation Analysis

    图  2  叶期移栽部分区域丰度图

    Figure  2.  Abundance Images of the Transplanting Plots in Leaves Stage

    图  3  叶期直播部分区域丰度图

    Figure  3.  Abundance Images of the Direct Seeding Plots in Leaves Stage

    图  4  花期地物光谱反射曲线

    Figure  4.  Spectral Reflectance Curves of Features in Flowering Period

    图  5  花期移栽部分区域丰度图

    Figure  5.  Abundance Images of the Transplanting Plots in Flowering Stage

    图  6  花期直播部分区域丰度图

    Figure  6.  Abundance Images of the Direct Seeding Plots in Flowering Stage

    图  7  角果期移栽部分区域丰度图

    Figure  7.  Abundance Images of the Transplanting Plots in Podding Stage

    图  8  角果期直播部分区域丰度图

    Figure  8.  Abundance Images of the Direct Seeding Plots in Podding Stage

    表  1  皮尔逊相关系数与相关强度

    Table  1.   Pearson Correlation Coefficients and Correlation Intensities of Variables

    相关系数绝对值取值范围 相关强度
    0.8~1.0 极强相关
    0.6~0.8 强相关
    0.4~0.6 中等程度相关
    0.2~0.4 弱相关
    0.0~0.2 极弱相关或无相关
    下载: 导出CSV

    表  2  花期丰度与产量相关系数

    Table  2.   Correlation Coefficients of Abundance and Yield in Flowering Period

    端元丰度 移栽 直播
    叶丰度 0.885 -0.475
    花丰度 0.473 0.837
    下载: 导出CSV

    表  3  角果期丰度与产量相关系数

    Table  3.   Correlation Coefficients of Abundance and Yield in Podding Stage

    端元丰度 移栽 直播
    叶丰度 0.880 0.531
    角果丰度 0.256 0.766
    下载: 导出CSV

    表  4  方案1模型评价指标

    Table  4.   Evaluation Indexes of Regression Models of the First Scheme

    种植方式 R2 RMSE
    移栽 0.847 284.982
    直播 0.849 264.521
    下载: 导出CSV

    表  5  方案2自变量组合

    Table  5.   Combinations of Independent Variables of the Second Scheme

    种植方式 组合序号 叶期 花期 角果期
    移栽 A 叶丰度 叶丰度 叶丰度
    直播 A 叶丰度 叶丰度 叶丰度
    B 叶丰度 叶丰度 角果丰度
    C 叶丰度 花丰度 叶丰度
    D 叶丰度 花丰度 角果丰度
    下载: 导出CSV

    表  6  方案2模型评价指标

    Table  6.   Evaluation Indexes of Regression Models of the Second Scheme

    种植方式 丰度组合序号 R2 RMSE
    移栽 A 0.890 241.153
    直播 A 0.811 295.873
    B 0.797 306.500
    C 0.794 308.677
    D 0.749 340.746
    下载: 导出CSV

    表  7  方案3模型评价指标

    Table  7.   Evaluation Indexes of Regression Models of the Third Scheme

    种植方式 丰度组合序号 R2 RMSE
    移栽 A 0.909 219.278
    直播 A 0.862 252.676
    B 0.847 266.452
    C 0.822 286.948
    D 0.810 296.859
    下载: 导出CSV

    表  8  方案3 A丰度组合估产模型公式

    Table  8.   Formula of Regression Models of A Abundance Combination of the Third Scheme

    种植方式 估产模型
    移栽 y = 1 236.118x1 + 2 961.941x2 + 2 918.376x3 + 291.675
    直播 y = 2 137.388x1 - 5 134.258x2 + 2 046.832x3 + 1 754.69
    下载: 导出CSV

    表  9  3种方案提前估产模型评价指标

    Table  9.   Evaluation Indexes of Advanced Regression Models of the Three Schemes

    种植方式 评价指标 方案1 方案2 方案3
    移栽 R2 0.838 0.814 0.867
    RMSE 293.379 313.762 265.565
    直播 R2 0.846 0.791 0.853
    RMSE 267.299 311.359 260.744
    下载: 导出CSV

    表  10  方案3提前估产模型公式

    Table  10.   Formula of Advanced Regression Models of the Third Scheme

    种植方式 估产模型
    移栽 y = 2 499.708x1 + 2 457.636x2 + 198.883
    直播 y = 2 692.617x1 - 4 567.276x2 + 1 494.461
    下载: 导出CSV
  • [1] 王汉中.我国油菜产需形势分析及产业发展对策[J].中国油料作物学报, 2007, 29(1):101-105 doi:  10.3321/j.issn:1007-9084.2007.01.022

    Wang Hanzhong. Strategy for Rapeseed Industry Development Based on the Analysis of Rapeseed Production and Demand in China[J]. Chinese Journal of Oil Crop Sciences, 2007, 29(1):101-105 doi:  10.3321/j.issn:1007-9084.2007.01.022
    [2] 王人潮, 王珂, 沈掌泉, 等.水稻单产遥感估测建模研究[J].遥感学报, 1998, 2(2):119-124

    Wang Renchao, Wang Ke, Shen Zhangquan, et al. A Study on Modeling of Rice Yield Estimation by Remote Sensing[J]. Journal of Remote Sensing, 1998, 2(2):119-124
    [3] Gitelson A A, Kaufman Y J, Stark R, et al. Novel Algorithms for Remote Estimation of Vegetation Fraction[J]. Remote Sensing of Environment, 2002, 80(1):76-87 doi:  10.1016/S0034-4257(01)00289-9
    [4] Viña A, Gitelson A A, Nguy-Robertson A L, et al. Comparison of Different Vegetation Indices for the Remote Assessment of Green Leaf Area Index of Crops[J]. Remote Sensing of Environment, 2011, 115(12):3468-3478 doi:  10.1016/j.rse.2011.08.010
    [5] Becker-Reshef I, Vermote E, Lindeman M, et al. A Generalized Regression-Based Model for Forecasting Winter Wheat Yields in Kansas and Ukraine Using MODIS Data[J]. Remote Sensing of Environment, 2010, 114(6):1312-1323 doi:  10.1016/j.rse.2010.01.010
    [6] 秦占飞, 申健, 谢宝妮, 等.引黄灌区水稻叶面积指数的高光谱估测模型[J].武汉大学学报·信息科学版, 2017, 42(8):1159-1166 http://ch.whu.edu.cn/CN/abstract/abstract5811.shtml

    Qin Zhanfei, Shen Jian, Xie Baoni, et al. Hyperspectral Estimation Model for Predicting LAI of Rice in Ningxia Irrigation Zone[J]. Geomatics and Information Science of Wuhan University, 2017, 42(8):1159-1166 http://ch.whu.edu.cn/CN/abstract/abstract5811.shtml
    [7] 王蕾, 王鹏新, 李俐, 等.应用条件植被温度指数预测县域尺度小麦单产[J].武汉大学学报·信息科学版, 2018, 43(10):1566-1573 http://ch.whu.edu.cn/CN/abstract/abstract6228.shtml

    Wang Lei, Wang Pengxin, Li Li, et al. Wheat Yield Forecasting at County Scale Based on Time Series Vegetation Temperature Condition Index[J]. Geomatics and Information Science of Wuhan University, 2018, 43(10):1566-1573 http://ch.whu.edu.cn/CN/abstract/abstract6228.shtml
    [8] 江东, 王乃斌, 杨小唤, 等.NDVI曲线与农作物长势的时序互动规律[J].生态学报, 2002, 22(2):247-252 http://d.old.wanfangdata.com.cn/Periodical/stxb200202015

    Jiang Dong, Wang Naibin, Yang Xiaohuan, et al. Principles of the Interaction Between NDVI Profile and the Growing Situation of Crops[J]. Acta Ecologica Sinica, 2002, 22(2):247-252 http://d.old.wanfangdata.com.cn/Periodical/stxb200202015
    [9] Bartholome E. Radiometric Measurements and Crop Yield Forecasting Some Observations over Millet and Sorghum Experimental Plots in Mali[J]. International Journal of Remote Sensing, 2007, 9(10-11):1539-1552 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=10.1080/01431168808954959
    [10] 张良培.高光谱遥感[M].武汉:武汉大学出版社, 2005

    Zhang Liangpei. Hyperspectral Remote Sensing[M].Wuhan:Wuhan University Press, 2005
    [11] 吕长春.混合像元分解模型综述[J].遥感信息, 2003(3):55-58, 60 doi:  10.3969/j.issn.1000-3177.2003.03.016

    Lv Changchun. A Review of Pixel Unmixing Models[J]. Remote Sensing Information, 2003(3):55-58, 60 doi:  10.3969/j.issn.1000-3177.2003.03.016
    [12] Rao M, Mendham N J, Buzza G C. Effect of the Apetalous Flower Character on Radiation Distribution in the Crop Canopy, Yield and Its Components in Oilseed Rape (Brassica Napus)[J]. Journal of Agricultural Science, 1991, 117(2):189-196 doi:  10.1017/S0021859600065278
    [13] Fielding A. A Review of Methods for the Assessment of Prediction Errors in Conservation Presence/Absence Models[J]. Environ Conser, 1997, 24(1):38-49 doi:  10.1017/S0376892997000088
    [14] 赵懿.油菜光合功能衰退研究[D].北京: 中国农业科学院, 2006

    Zhao Yi. Study on Photosynthetic Functional Decline for Oilseed Rape[D]. Beijing: Chinese Academy of Agricultural Sciences, 2006
    [15] Heinz D C, Chang C I. Fully Constrained Least Squares Linear Spectral Mixture Analysis Method for Material Quantification in Hyperspectral Imagery[J]. IEEE Transactions on Geoscience and Remote Sensing, 2001, 39(3):529-545 doi:  10.1109/36.911111
    [16] 刘后利.实用油菜栽培学[M].上海:上海科学技术出版社, 1987

    Liu Houli. Practical Rapeseed Cultivation[M]. Shanghai:Shanghai Science and Technology Press, 1987
    [17] 杨光.油菜高效结角层结构的研究[D].扬州: 扬州大学, 2002

    Yang Guang. The Study of High-Efficient Pod Canopy of Oilseed Rape[D]. Yangzhou: Yangzhou University, 2002
    [18] 赵明, 李少昆, 王志敏, 等.论作物源的数量、质量关系及其类型划分[J].中国农业大学学报, 1998, 3(3):53-58 doi:  10.3321/j.issn:1007-4333.1998.03.010

    Zhao Ming, Li Shaokun, Wang Zhimin, et al. Exposition on Genotypes Classification and Relationship Between Quality and Quantity of Crop Source[J]. Journal of China Agricultural University, 1998, 3(3):53-58 doi:  10.3321/j.issn:1007-4333.1998.03.010
    [19] 宋亚斌, 林辉.基于多元回归模型的叶面积指数遥感反演[J].中南林业调查规划, 2012, 31(2):24-27 doi:  10.3969/j.issn.1003-6075.2012.02.007

    Song Yabin, Lin Hui.The Inversion of Leaf Area Index Multivariable Through Remote Sensing Based on the Regression Model[J]. Central South Forest Inventory and Planning, 2012, 31(2):24-27 doi:  10.3969/j.issn.1003-6075.2012.02.007
    [20] 任建强, 陈仲新, 周清波, 等. MODIS植被指数的美国玉米单产遥感估测[J].遥感学报, 2015, 19(4):568-577 http://d.old.wanfangdata.com.cn/Periodical/ygxb201504004

    Ren Jianqiang, Chen Zhongxin, Zhou Qingbo, et al. MODIS Vegetation Index Data Used for Estimating Corn Yield in USA[J]. Journal of Remote Sensing, 2015, 19(4):568-577 http://d.old.wanfangdata.com.cn/Periodical/ygxb201504004
  • [1] 张王菲, 文哲, 张亚红, 张庭苇, 李云.  Stokes参数在油菜长势监测中的可行性分析 . 武汉大学学报 ● 信息科学版, 2020, 45(2): 242-249. doi: 10.13203/j.whugis20180375
    [2] 李维炼, 朱军, 朱秀丽, 付林, 张昀昊, 胡亚, 何秋玲, 黄鹏诚.  无人机遥感数据支持下滑坡VR场景探索分析方法 . 武汉大学学报 ● 信息科学版, 2019, 44(7): 1065-1072. doi: 10.13203/j.whugis20190047
    [3] 吴柯, 何坦, 杨叶涛.  基于混合像元分解与EM算法的中低分辨率遥感影像变化检测 . 武汉大学学报 ● 信息科学版, 2019, 44(4): 555-562. doi: 10.13203/j.whugis20160546
    [4] 李鹏飞, 孙开敏, 李德仁, 王玮.  无人机影像应急并行处理负载均衡方法 . 武汉大学学报 ● 信息科学版, 2018, 43(2): 268-274. doi: 10.13203/j.whugis20130095
    [5] 李鹏飞, 孙开敏, 李德仁, 王玮.  无人机影像应急并行处理负载均衡方法 . 武汉大学学报 ● 信息科学版, 2018, 43(2): 268-274. doi: 10.13203/j.whugis20170022
    [6] 眭海刚, 冯文卿, 李文卓, 孙开敏, 徐川.  多时相遥感影像变化检测方法综述 . 武汉大学学报 ● 信息科学版, 2018, 43(12): 1885-1898. doi: 10.13203/j.whugis20180251
    [7] 樊邦奎, 张瑞雨.  无人机系统与人工智能 . 武汉大学学报 ● 信息科学版, 2017, 42(11): 1523-1529. doi: 10.13203/j.whugis20170177
    [8] 杨楠, 邵振峰, 郭丙轩, 彭哲, 黄蕾.  基于非固定初始面元的无人机影像点云优化算法 . 武汉大学学报 ● 信息科学版, 2016, 41(8): 1013-1020. doi: 10.13203/j.whugis20140750
    [9] 李粉玲, 常庆瑞, 刘佳岐, 刘京.  基于多纹理和支持向量机的ZY-102C星HR数据分类 . 武汉大学学报 ● 信息科学版, 2016, 41(4): 455-461,486. doi: 10.13203/j.whugis20140356
    [10] 汪炯骅, 许涛, 李响.  低空无人机交叉口悬停影像抖动纠正 . 武汉大学学报 ● 信息科学版, 2015, 40(6): 738-743. doi: 10.13203/j.whugis20130834
    [11] 甘文霞, 沈焕锋, 张良培, 龚威.  采用6S模型的多时相MODIS植被指数NDVI归一化方法 . 武汉大学学报 ● 信息科学版, 2014, 39(3): 300-304. doi: 10.13203/j.whugis20120731
    [12] 李德仁, 李明.  无人机遥感系统的研究进展与应用前景 . 武汉大学学报 ● 信息科学版, 2014, 39(5): 505-513. doi: 10.13203/j.whugis20140045
    [13] 王楠, 张良培, 杜博.  最小光谱相关约束NMF的高光谱遥感图像混合像元分解 . 武汉大学学报 ● 信息科学版, 2014, 39(1): 22-26.
    [14] 张春森, 李辉.  顾及混合像元分解的遥感图像光谱模拟 . 武汉大学学报 ● 信息科学版, 2013, 38(9): 1052-1056.
    [15] 刘汉丽, 裴韬, 周成虎, 朱阿兴.  结合MNF变换与灰值形态学的三江平原多光谱、多时相MODIS遥感影像分类 . 武汉大学学报 ● 信息科学版, 2011, 36(2): 153-156.
    [16] 刘晋苏, 任沁源, 韩波, 李平.  基于微型无人机平台的航拍动态阴影检测 . 武汉大学学报 ● 信息科学版, 2010, 35(11): 1275-1278.
    [17] 宋妍, 田玉刚, 贾小霞.  基于Zernike矩的多时相遥感影像匹配方法研究 . 武汉大学学报 ● 信息科学版, 2009, 34(10): 1226-1230.
    [18] 吴波, 张良培, 李平湘.  基于支撑向量机概率输出的高光谱影像混合像元分解 . 武汉大学学报 ● 信息科学版, 2006, 31(1): 51-54.
    [19] 邓文胜, 关泽群, 秦昆.  命题逻辑演算方法在多时相遥感影像专题解译中的应用 . 武汉大学学报 ● 信息科学版, 2005, 30(2): 170-173.
    [20] 沙宗尧, 边馥苓.  远程智能诊断与决策支持系统的设计与开发 . 武汉大学学报 ● 信息科学版, 2004, 29(12): 1055-1058.
  • 加载中
图(8) / 表(10)
计量
  • 文章访问数:  412
  • HTML全文浏览量:  48
  • PDF下载量:  98
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-07-11
  • 刊出日期:  2020-02-05

多时相NDVI与丰度综合分析的油菜无人机遥感长势监测

doi: 10.13203/j.whugis20180161
    基金项目:

    国家高技术研究发展计划 2013AA102401

    中央高校基本科研业务费专项资金 2042017kf0236

    作者简介:

    刘雅婷, 硕士生, 主要研究方向为定量遥感。liuyating@whu.edu.cn

    通讯作者: 龚龑, 博士, 教授。gongyan@whu.edu.cn
  • 中图分类号: P237

摘要: 从无人机油菜遥感数据精准分析的角度,提出了一种新型的基于多时相混合光谱分析的油菜长势监测方法。根据油菜冠层结构与田间特征确定不同生长期端元组合,并通过线性解混模型获取丰度。针对油菜产量积累过程,分析3个生长期端元丰度与产量的相关性。在此基础上,对比分析典型归一化植被指数(normalized difference vegetation index,NDVI)、端元丰度以及两者综合3种自变量方案下建立的多时相线性估产模型精度。实验结果表明,相比广泛使用的NDVI,选取合适的丰度组合数据结合NDVI进行多时相估产具有更好的效果,且不同种植方式下的最优丰度组合相同,适用于产量的提前预测。由此说明,丰度数据能从植株发育与产量构成角度对NDVI进行亚像元层次的修正,具有较强的稳定性。

English Abstract

刘雅婷, 龚龑, 段博, 方圣辉, 彭漪. 多时相NDVI与丰度综合分析的油菜无人机遥感长势监测[J]. 武汉大学学报 ● 信息科学版, 2020, 45(2): 265-272. doi: 10.13203/j.whugis20180161
引用本文: 刘雅婷, 龚龑, 段博, 方圣辉, 彭漪. 多时相NDVI与丰度综合分析的油菜无人机遥感长势监测[J]. 武汉大学学报 ● 信息科学版, 2020, 45(2): 265-272. doi: 10.13203/j.whugis20180161
LIU Yating, GONG Yan, DUAN Bo, FANG Shenghui, PENG Yi. Combining Multi-temporal NDVI and Abundance from UAV Remote Sensing Data for Oilseed Rape Growth Monitoring[J]. Geomatics and Information Science of Wuhan University, 2020, 45(2): 265-272. doi: 10.13203/j.whugis20180161
Citation: LIU Yating, GONG Yan, DUAN Bo, FANG Shenghui, PENG Yi. Combining Multi-temporal NDVI and Abundance from UAV Remote Sensing Data for Oilseed Rape Growth Monitoring[J]. Geomatics and Information Science of Wuhan University, 2020, 45(2): 265-272. doi: 10.13203/j.whugis20180161
  • 油菜是中国传统优势油料作物,总种植面积约占全国油料作物总面积的50%,是重要的食用油原料和油粕原料[1]。对油菜进行快速准确的监测,不仅关系到保障食用油稳定供应,也关系到中国油料产品贸易与油料产业的发展。

    传统的农作物产量预测方法主要有抽样、气象估产、农学估产和统计模型估产等,但在监测过程中,这些估产方式都表现出明显的局限性[2]。迅速发展的遥感技术可提供对地宏观、动态、快速的观测,为作物长势的动态监测和产量估算提供一种全新的科学手段。作物遥感监测研究多使用卫星遥感数据,借助光谱植被指数等对作物长势进行监测。Gitelson等[3]基于小麦冠层光谱特性,研究了不同植被指数与植被覆盖率的关系;Viña等[4]利用植被指数对玉米与大豆的叶面积指数进行遥感估算;Becker-Reshef等[5]通过时间序列归一化植被指数(normalized difference vegetation index, NDVI)数据对堪萨斯州小麦产量进行估算,并在乌克兰地区实现提前预测;秦占飞等[6]分析了不同波段组合,构建了4种植被指数与水稻叶面积指数反演模型;王蕾等[7]利用条件植被温度指数预测了县域尺度的小麦单产;江东等[8]以中国河南省为例,分析了不同地区冬小麦高、中、低产县NDVI时间线曲线和不同生长阶段的关系。有研究[9]认为,NDVI与产量相关,进行遥感估产最普遍的方式是在NDVI与作物产量之间建立经验模型。广泛应用的产量与遥感光谱指数的单纯统计模式操作简单,便于推广。然而,传统植被遥感中的NDVI无法体现农作物的具体植株器官形态及发育,忽略了不同器官在光合作用、淀粉累计等产量形成环节中的差异,无法满足农作物精准分析的客观需求。另外,从产量构成角度来看,物理意义不够明确也使植被指数模型在跨地区与跨年的应用效果不佳。因此,从作物成像和产量形成机理上对油菜进行定量研究尤为迫切。

    遥感影像中的光谱信息可以反映作物生长情况,不同时期各端元丰度数据可反映作物的物候发育、物质生产与分配,与作物最终产量紧密相联。由于混合像元的广泛存在,丰度数据可以改善植被指数单纯基于像元反射率的估产方式,能从亚像元层次与植株器官发育的角度对NDVI估产模型进行修正。本文基于无人机遥感的方式,从单一时期丰度数据与产量构成关系分析出发,紧密结合油菜自身生长的农学特性,提出了新型的利用多时相丰度信息与NDVI相结合的分析方法,建立了适合油菜种植实际特点的估产模型。

    • 本文实验数据包括无人机影像数据、地物端元光谱数据和油菜产量数据。实验田内共有48块实验小区,其中24块小区直播种植,另24块小区移栽种植。无人机影像数据包括40 m航高的油菜叶期影像、50 m航高的花期影像以及50 m航高的角果期影像,使用的Mini-MCA(multiple camera array)6波段相机的波段范围分别为485~495 nm、545~555 nm、665~675 nm、715~725 nm、780~820 nm、890~910 nm。对3个时期无人机影像的6个波段分别进行几何纠正和辐射定标。几何纠正时,以第5波段作为基准波段进行配准。经过几何纠正后,使用与Mini-MCA传感器配套的标准定标毯进行场地定标,得到地表反射率数据。地物端元光谱使用ASD(analytical spectral devices)光谱仪采集,探测波段范围为350~2 500 nm。采集时,光谱仪探头垂直向下,距地物0.2 m左右,每次记录5个光谱作为一采样光谱,选定位置分布均匀的3~5个点,将其光谱平均值作为该地物端元的光谱反射值。由于采集到的ASD光谱数据在对应MCA相机6个波段范围内的波动较小,对每个端元参照无人机影像的6个波段范围,分别将3个时期的地物光谱数据进行均值处理,得到平均反射率数据。

    • 油菜作物结构复杂,传感器成像时,单一像元内往往包含不同的地物覆盖成分,混合光谱现象广泛存在。为准确获取目标地物的光谱信息,提高遥感影像的应用精度,必须进行混合光谱分解[10]。具体步骤分为端元光谱获取和混合像元分解。端元获取的常见方式有图像端元提取法、物理端元提取法、光谱库获取法等。本文使用物理端元提取法,通过野外采集地物光谱曲线的方式获取端元光谱。常见的混合光谱分解模型有线性模型、概率模型、随机几何模型等。其中,线性模型使用最为广泛,原理也较为简单,一般形式如下:

      $$\mathit{\boldsymbol{P}} = \sum\limits_{i = 1}^N {{c_i}{e_i} + e = } Ec + n$$ (1)

      式中,P为像元的波谱向量;N为地物的端元数量;E为在每个波段范围内端元的反射向量;c为端元所占比例系数向量;ci为第i个端元所占像元的面积比;ei为第i个端元的反射率;en为误差项。在实际光谱解混中,线性模型会受到两个约束条件限制——非负性约束和归一化约束[11],即:

      $$\sum\limits_{i = 1}^N {{c_i} = 1,0 \le {c_i} \le 1} $$ (2)
    • 油菜产量的形成需要经过漫长的生育期,不同时期植株器官发育均影响着作物最终经济产量。本文通过探讨不同时期各端元丰度与产量的相关性,将NDVI估产模式置于产量形成过程的具体生物分析基础上,建立结合油菜生育期特征的估产模型。

    • 为了探究油菜在不同生长期植株器官发育对产量的影响,分别将叶期、花期、角果期经解混得到的端元丰度数据与产量进行相关分析,探寻对产量最为敏感的要素。在此基础上,进行多时相相关分析并建立多时相估产模型,如图 1所示。首先对无人机影像进行几何纠正和辐射定标处理,得到作物地表反射率影像;然后根据无人机波段范围,对采集到的地物端元光谱数据进行均值处理,得到端元光谱矩阵;接着基于全约束混合光谱线性解混模型对地表反射率影像进行解混,得到端元丰度数据;依据丰度与产量的相关性分析,得到产量最为敏感的端元;最后对比结合NDVI估产模式,通过多时相相关性分析,建立多时相估产模型。

      图  1  多时相相关分析流程图

      Figure 1.  Flowchart of Multi-temporal Correlation Analysis

      决定油菜籽粒产量的主要因素有单位面积角果数、角粒数和千粒重。最终产量是干物质的积累,约90%~95%产量来源于光合产物。产量的形成主要分为3个阶段:从出苗到花芽分化前,叶片同化作用产生的干物质供植株建成和后续花芽分化,该阶段为角粒数和单位面积角果数的积淀过程;初花期到终花期是角粒数和单位面积角果数的定型过程;在角果期,角果层的成熟过程决定了千粒重大小[12]。由于作物每一生长时段周围环境的可变性,只通过作物某一生长阶段的瞬时信息预测产量,往往考虑得不够全面。为较为准确地表征这种积累过程,使用多元线性回归的方式结合3个生长期的数据建立多时相估产模型。多元线性回归模型的一般形式为:

      $$ y = {\beta _0} + {\beta _1}{x_1} + {\beta _2}{x_2} + \ldots + {\beta _i}{x_i} $$ (3)

      式中,y为因变量;x为自变量;i为自变量个数;β0βi分别是常数项与回归系数。建立多时相线性模型时,自变量的选取需考虑油菜不同生长期中有助于最终产量形成的部分。

    • 皮尔逊相关系数常用于度量两个变量XY之间的相关性(线性相关),对单一时期敏感要素与产量进行相关分析时,本文使用皮尔逊相关系数(ρXY)进行衡量:

      $$ {\rho _{X,Y}} = \frac{{ \sum \nolimits^ ({X_i} - \bar X)\left( {{Y_i} - \bar Y} \right)}}{{\sqrt[{}]{{ \sum \nolimits^ {{({X_i} - \bar X)}^2}}}\sqrt[{}]{{ \sum \nolimits^ {{({Y_i} - \bar Y)}^2}}}}} $$ (4)

      式中,XiYi分别为第i个样本的实际值;XY分别为XY的平均值。通常情况下,通过表 1所示的取值范围来判断相关强度。

      表 1  皮尔逊相关系数与相关强度

      Table 1.  Pearson Correlation Coefficients and Correlation Intensities of Variables

      相关系数绝对值取值范围 相关强度
      0.8~1.0 极强相关
      0.6~0.8 强相关
      0.4~0.6 中等程度相关
      0.2~0.4 弱相关
      0.0~0.2 极弱相关或无相关

      由于实际种植条件的限制,一般难以获取大量的产量数据。针对本文实验样本数量较少的情况(一组24个),使用留一交叉验证的方式来建立最终估产模型。留一交叉验证是最为常见的交叉验证方法之一,广泛应用于模型建立和验证,能够充分利用实验数据[13]。其基本思想是每次从n个样本中选取1个作为验证集,剩下的n-1个样本作为训练集进行训练。训练过程重复n次,将n次训练过程的结果取平均值作为泛化误差的估计,得到最终模型的系数和精度,具体计算公式如下:

      $$C = \frac{1}{n} \sum\limits_{i = 1} ^n {C_i}$$ (5)
      $${\overline {R{\rm{}}} ^2} = 1 - \left( {1 - {R^2}} \right)\frac{{n - 1}}{{n - p - 1}}$$ (6)
      $${\rm{RMSE}} = \sqrt[{}]{{\frac{1}{n} \sum\limits_{i = 1} ^n {{({X_{{\rm{obs}},j}} - {X_{{\rm{model}},j}})}^2}}}$$ (7)

      式中,n为样本个数;C与$C_i$分别为最终模型系数与第i次训练集训练得到的模型系数;$R^2$与R2分别为判定系数与调整后的判定系数;$p$为自变量个数;RMSE(root mean squared error)表示均方根误差;Xobs, jXmodel, j分别为第j个样本的样本实际值和模型预测值。

    • 由于移栽油菜排列较整齐,植株间距较大,而直播油菜生长较为密集,故两种种植方式下油菜区域结构有所差别,需分别进行敏感要素分析与建模,在此基础上分析二者的共性规律,进而得到一般性估产方法。

    • 在油菜叶期,实验田内主要为叶子和土壤。从油菜植株的生长发育角度来看,叶片作为油菜光合作用和蒸腾作用的主要场所,其同化作用过程中干物质的积累对油菜植株、角果以及籽粒的生长发育有很大影响。叶期叶片分为长柄叶和短柄叶,其中长柄叶主要供应根系养分;短柄叶养分供应兼顾上下,是春壮稳长所调控的一组叶片[14]。随着油菜的生长,不同类型叶片净光合速率的变化存在差异。

      考虑到本文实验土壤环境与叶片的发育情况,选取5种地物端元进行光谱采集,分别为长柄叶、短柄叶、干土、半干土、湿土。基于全约束混合像元线性解混的方法计算得到叶期移栽和叶期直播各端元的丰度图(见图 2图 3),最终叶丰度为长柄叶与短柄叶两种端元丰度之和。对叶端元丰度与产量进行相关分析可知,移栽油菜与产量相关系数为0.889,直播油菜与产量相关系数为0.818。两种种植方式下油菜叶丰度与产量均为极强相关,相关系数均高于0.8,能很好地反映叶丰度与产量之间的联系。

      图  2  叶期移栽部分区域丰度图

      Figure 2.  Abundance Images of the Transplanting Plots in Leaves Stage

      图  3  叶期直播部分区域丰度图

      Figure 3.  Abundance Images of the Direct Seeding Plots in Leaves Stage

    • 在开花后,油菜花序会逐渐遮挡植株茎秆,叶片光合作用会逐渐减弱,光合产物主要提供给茎干生长、花序分化。从产量积累过程来看,叶片仍然是光合作用的主要场所。花期叶片包括无柄叶、短柄叶和长柄叶,其中短柄叶的光合贡献更显著。由于干物质逐渐向花转移,油菜叶子会开始变黄,此时,实验田内地物组成较为复杂,不同地物的光谱反射曲线如图 4所示。

      图  4  花期地物光谱反射曲线

      Figure 4.  Spectral Reflectance Curves of Features in Flowering Period

      图 4可以看到,花端元的光谱曲线走势较为接近;不同类型的叶片在绿波段的反射率有着较大差异;土壤和油菜植株的反射率差别很大。由于无人机相机只有6个波段,满足线性解混模型条件下求解丰度时最多只能选取6种端元的要求[15]。选取上层花、无柄叶、短柄叶、干土、半干土、湿土6种端元进行解混,得到各端元丰度图(见图 5图 6)。分别将花丰度和叶片总丰度与产量进行皮尔逊相关分析,统计相关系数,见表 2

      图  5  花期移栽部分区域丰度图

      Figure 5.  Abundance Images of the Transplanting Plots in Flowering Stage

      图  6  花期直播部分区域丰度图

      Figure 6.  Abundance Images of the Direct Seeding Plots in Flowering Stage

      表 2  花期丰度与产量相关系数

      Table 2.  Correlation Coefficients of Abundance and Yield in Flowering Period

      端元丰度 移栽 直播
      叶丰度 0.885 -0.475
      花丰度 0.473 0.837

      表 2可以看出,不同种植方式下,不同端元与产量相关性不同。对移栽油菜而言,叶丰度与产量的相关性极强,花丰度与产量的相关性中等,而直播油菜则相反。值得注意的是,直播油菜叶丰度与产量是负相关的,这是由于直播油菜植株的单株产能不高,但最终亩产密度大。

    • 角果期是油菜生殖生长阶段进行养分积累和物质合成的关键时期。终花后,角果和叶片一同进行光合作用;角果层形成后,由于对叶片的严重遮盖,叶片开始衰老和掉落,短柄叶和无柄叶表面积减少。角果皮产生的同化物转移分配给籽粒,供籽粒充分灌浆[16]。随着角果鲜重逐渐加大,果序轴被压弯,导致角果下坠或重叠在一起[17]。根据角果期实验田地物组成,选取无柄叶、短柄叶、上层角果、下层角果、干土、湿土6种端元,并解算丰度图(见图 7图 8)。分别将角果丰度与叶片丰度与产量进行回归分析,统计相关系数,如表 3所示。可以看到,移栽油菜的叶丰度与产量仍保持极强的相关性,角果丰度与产量相关性较弱;而直播油菜的情况较为复杂,叶丰度和产量相关系数为0.531,角果丰度和产量相关系数为0.766,两者差异不明显。

      图  7  角果期移栽部分区域丰度图

      Figure 7.  Abundance Images of the Transplanting Plots in Podding Stage

      图  8  角果期直播部分区域丰度图

      Figure 8.  Abundance Images of the Direct Seeding Plots in Podding Stage

      表 3  角果期丰度与产量相关系数

      Table 3.  Correlation Coefficients of Abundance and Yield in Podding Stage

      端元丰度 移栽 直播
      叶丰度 0.880 0.531
      角果丰度 0.256 0.766
    • 油菜籽粒产量的高低和品质好坏主要受油菜植株光合作用的调控,大约90%~95%的产量来自光合产物[18]。根据产量形成过程与单一时期丰度数据分析可知,在油菜生殖生长过程中,不同器官发育状况都影响着最终产量。相比仅用一个自变量估测产量,采用多个自变量的最优组合共同估测因变量的方式更符合实际,也更有效[19]。大量估产研究表明,NDVI是反映植被生长状态及植被覆盖度的最佳因子[20]。为更好地评价多时相丰度数据综合估产的效果,分别选取以下两种不同的自变量方案进行综合回归分析。

      1)方案1:将3个生长期NDVI数据作为自变量建立估产模型,模型精度如表 4所示。

      表 4  方案1模型评价指标

      Table 4.  Evaluation Indexes of Regression Models of the First Scheme

      种植方式 R2 RMSE
      移栽 0.847 284.982
      直播 0.849 264.521

      2)方案2:分别在3个时期内各选取一种端元丰度数据作为自变量建立估产模型。根据单时期端元丰度与产量相关分析可知,对于移栽油菜,3个生长期的叶丰度与产量均极强相关,且明显强于其他端元丰度;而直播油菜情况较为复杂,没有在多个时期保持与产量相关性突出的端元丰度。因此,针对移栽油菜,方案2分别选择3个时期的叶丰度作为自变量;而对于直播油菜,尝试4种自变量组合(见表 5)。使用不同的自变量组合时,模型精度如表 6所示。

      表 5  方案2自变量组合

      Table 5.  Combinations of Independent Variables of the Second Scheme

      种植方式 组合序号 叶期 花期 角果期
      移栽 A 叶丰度 叶丰度 叶丰度
      直播 A 叶丰度 叶丰度 叶丰度
      B 叶丰度 叶丰度 角果丰度
      C 叶丰度 花丰度 叶丰度
      D 叶丰度 花丰度 角果丰度

      表 6  方案2模型评价指标

      Table 6.  Evaluation Indexes of Regression Models of the Second Scheme

      种植方式 丰度组合序号 R2 RMSE
      移栽 A 0.890 241.153
      直播 A 0.811 295.873
      B 0.797 306.500
      C 0.794 308.677
      D 0.749 340.746

      表 6可以看到,方案2中不同丰度组合下模型调整后测定的系数均在0.74以上(A组合在0.8以上),说明从植株发育与干物质的积累方面预测最终产量是可行的。移栽油菜的方案2相比方案1拥有更好的拟合效果;而直播油菜4种端元丰度组合方式的效果均略逊于方案1,其中A组合表现最佳,这是由于直播植株生长密集、冠层结构复杂造成的。方案1主要通过增强近红外与红光通道反射率的对比度来监测植株的生长情况;而方案2从植株器官发育的角度出发,利用不同端元丰度数据反映干物质的积累与转移过程,进而估测产量,解混得到的丰度数据也可在亚像元层次反映植株的生长情况。

      3)为结合光谱指数和丰度信息的优势,本文提出方案3,结合方案2中提出的不同端元丰度组合方式,将3个生长期的3种端元丰度与对应时期的NDVI相乘,乘积作为自变量与产量进行建模。模型最终评价指标统计结果如表 7所示。

      表 7  方案3模型评价指标

      Table 7.  Evaluation Indexes of Regression Models of the Third Scheme

      种植方式 丰度组合序号 R2 RMSE
      移栽 A 0.909 219.278
      直播 A 0.862 252.676
      B 0.847 266.452
      C 0.822 286.948
      D 0.810 296.859

      对比方案1中的R2和RMSE,方案3中的A丰度组合方式可明显提高估产精度,这说明丰度数据可以改善植被指数单纯基于像元反射率的估产效果。由于移栽油菜冠层结构较直播油菜简单,方案3的估产精度提高效果更为明显,虽然不同种植方式的精度提高程度有所差异,但两种种植条件下最优丰度自变量组合选取一致,均为A组合。这说明在选取合适的端元丰度组合时,丰度对NDVI均有一定的修正效果,且在不同种植方式下均适用。在实际生产中, 中国油菜多采用移栽种植,而国外部分地区则采用直播种植。表 8列出了选取A组合时不同种植方式下方案3估产模型的回归公式。

      表 8  方案3 A丰度组合估产模型公式

      Table 8.  Formula of Regression Models of A Abundance Combination of the Third Scheme

      种植方式 估产模型
      移栽 y = 1 236.118x1 + 2 961.941x2 + 2 918.376x3 + 291.675
      直播 y = 2 137.388x1 - 5 134.258x2 + 2 046.832x3 + 1 754.69

      在实际应用中,产量往往需要提前预测,以便采取及时有效的增产措施。仅选取叶期和花期数据分别按照上述3种自变量方案(其中方案2、方案3的自变量均使用A丰度组合)建立多时相估产模型,模型回归统计结果如表 9所示。

      表 9  3种方案提前估产模型评价指标

      Table 9.  Evaluation Indexes of Advanced Regression Models of the Three Schemes

      种植方式 评价指标 方案1 方案2 方案3
      移栽 R2 0.838 0.814 0.867
      RMSE 293.379 313.762 265.565
      直播 R2 0.846 0.791 0.853
      RMSE 267.299 311.359 260.744

      表 9可知,仅用前两个生育期数据(叶期和花期)的模型效果略逊于三时期综合估产,但方案3的R2仍在0.85以上,能较准确地反映最终产量。对比方案1,不同种植方式下A丰度组合的方案3估产效果仍保持提高,说明丰度对NDVI的修正效果是稳定且可靠的。表 10给出了提前估产时方案3的估产模型公式(A丰度组合)。

      表 10  方案3提前估产模型公式

      Table 10.  Formula of Advanced Regression Models of the Third Scheme

      种植方式 估产模型
      移栽 y = 2 499.708x1 + 2 457.636x2 + 198.883
      直播 y = 2 692.617x1 - 4 567.276x2 + 1 494.461
    • 本文从油菜产量形成的农学意义出发,对油菜3个典型生长期进行长势监测分析,综合考虑反映油菜整体生长状态的NDVI和体现植株器官构成的丰度信息,得到了能进行产量预估的无人机遥感模型。实验结果表明:

      1)在不同的种植方式下,不同器官端元丰度与产量的相关规律有所差异,移栽油菜叶丰度在3个生长期与产量均有极强的相关性,而直播油菜由于冠层结构复杂,植株密度大,没有明显突出的强相关端元丰度。

      2)从干物质的生产分配的角度出发,利用多时相丰度数据建立的估产模型(方案2)也是可行和有效的,但精度略低于NDVI建立的估产模型。

      3)在选取合适的端元丰度组合时,丰度数据与NDVI结合进行综合估产的模型精度最高,说明丰度可从产量构成及亚像元层次改善NDVI的估产效果。虽然不同种植方式的提高程度有所差异,但两种种植条件下最优丰度自变量组合的选取是一致的(均为A组合),说明丰度对NDVI的修正效果是稳定且可靠的,适用于产量的提前预测。

      由于油菜植株是立体结构,使用线性混合光谱模型进行分解时不能完全描述地物情况。后续研究还需考虑不同种植条件下植株的生长差异,从更严密的角度进行混合光谱分析,建立更为一般化的估产模型。

参考文献 (20)

目录

    /

    返回文章
    返回