基于Google Earth Engine和NDVI时序差异指数的作物种植区提取

doi:10.12082/dqxxkx.2021.200291

地球信息科学学报 ›› 2021, Vol. 23 ›› Issue (5): 938-947.doi: 10.12082/dqxxkx.2021.200291

基于Google Earth Engine和NDVI时序差异指数的作物种植区提取

姜伊兰¹^,²(), 陈保旺¹, 黄玉芳¹, 崔佳琪¹^,², 郭宇龙¹^,²^,^*()

1.河南农业大学资源与环境学院,郑州 450002
2.河南省土地整治与生态重建工程技术研究中心,郑州 450002

收稿日期:2020-06-08 修回日期:2020-08-24 出版日期:2021-05-25 发布日期:2021-07-25
通讯作者: 郭宇龙
作者简介:姜伊兰（1997— ）,女,河南洛阳人,硕士生,主要从事土地资源管理研究。E-mail:934121453@qq.com
基金资助:
国家自然科学基金项目(41701422);河南省重点研发与推广专项(科技攻关)(192102310251)

Crop Planting Area Extraction based on Google Earth Engine and NDVI Time Series Difference Index

JIANG Yilan¹^,²(), CHEN Baowang¹, HUANG Yufang¹, CUI Jiaqi¹^,², GUO Yulong¹^,²^,^*()

1. College of Resources and Environment , Henan Agricultural University, Zhengzhou 450002, China
2. Henan Engineering Research Center of Land Consolidation and Ecological Restoration, Zhengzhou 450002, China

Received:2020-06-08 Revised:2020-08-24 Online:2021-05-25 Published:2021-07-25
Contact: GUO Yulong
Supported by:
National Natural Science Foundation of China(41701422);Key Research and Development and Promotion Projects in Henan Province(192102310251)

摘要/Abstract

摘要：

为提高农作物种植信息遥感监测的效率,扩展数据适用范围,本文提出了一种基于时间序列NDVI差异指数的作物种植区提取方法。随着海量遥感与云计算的发展,Google Earth Engine作为一个全球尺度地理空间分析云平台,弥补了单机计算耗时长的不足,为快速遥感分类带来了新机遇。基于Google Earth Engine平台,以河南省开封市杞县为研究区,以2019—2020年杞县地区多时相Sentinel-2影像为数据源,结合物候信息,根据不同作物在时间序列NDVI曲线上的差异构建NDVI时序差异指数,从而提取作物种植区,区分不同作物类型,并与其他方法进行了精度验证和对比。结果表明：① NDVI时序差异指数法以作物物候信息为基础,与GEE高性能的计算能力相结合,形成了作物种植信息快速提取框架,可以方便快捷地进行作物种植区提取,较本地处理具有明显优势;② 杞县冬小麦和大蒜种植区有明显的空间分异性,冬小麦种植区主要集中在研究区西北部以及南部的农村居民点周围,而杞县大蒜则由于产品流通需要,主要集中在研究区中部以及东北部,居民点较为密集,交通便利的城市周边;③ 与时间序列支持向量机法和最大似然法相比较, NDVI时序差异指数进行作物种植区提取的总体精度达到83.72%, Kappa系数为0.67,分别比最大似然法提高了10.02%和0.21,比支持向量机法提高了4.18%和0.09,表明该方法能更高效率,更高精度地提取作物种植信息,实现区域作物种植信息的高效准确监测。总体来看,该方法在一定程度上可拓展遥感数据在农业领域的应用范围,具有推广价值。

关键词: Sentinel-2, NDVI, 时间序列, 物候特征, Google Earth Engine, 作物提取, 精度对比, 杞县

Abstract:

In order to improve the efficiency of remote sensing monitoring of crop planting and expand applications of remote sensing data, a method of crop planting area extraction based on NDVI time series difference index is proposed. With the development of remote sensing and cloud computing technologies, Google Earth Engine, as a global-scale geospatial analysis cloud platform, overcomes the disadvantages of traditional single-machine computing and brings new opportunities for rapid remote sensing classification. In this study, taking Qi County in Henan province as the study area, the NDVI time series difference index of different crops is constructed according to the characteristics of time series NDVI curve of each crop to extract crop planting information and distinguish different crop types using multi-temporal Sentinel-2 images in 2019-2020 based on the Google Earth Engine platform. The extraction accuracy is verified and compared with other existing methods. The results show that the NDVI time series difference index is based on crop phenology information and developed using GEE's high-performance computing capability, which forms a framework for rapid crop planting information extraction and has obvious advantages over traditional local computing. The winter wheat and garlic planting areas in Qi County have obvious spatial variation. The winter wheat planting areas are mainly concentrated in the northwest and southern rural residential areas of the study area. While the garlic in Qi County is mainly concentrated in the central and northeastern part of the study area due to the needs of transportation. Compared with other methods using support vector machine and maximum likelihood, the overall accuracy of crop planting area extraction using the NDVI time series difference index reaches 83.72%, and the Kappa coefficient is 0.67. The overall accuracy and the Kappa coefficient are 10.02% and 0.21 respectively higher than the maximum likelihood method, and are 4.18% and 0.09 respectively higher than the support vector machine method, which indicates that our method can extract crop planting information with high efficiency and high accuracy. We develop an efficient and accurate monitoring method for regional crop planting information extraction and expand the application of remote sensing data in the agricultural field, which has significant value for future agricultural applications.

Key words: Sentinel-2, NDVI, time-series, phenological characteristics, Google Earth Engine, crop extraction, accuracy comparison, Qi County

姜伊兰, 陈保旺, 黄玉芳, 崔佳琪, 郭宇龙. 基于Google Earth Engine和NDVI时序差异指数的作物种植区提取[J]. 地球信息科学学报, 2021, 23(5): 938-947.DOI:10.12082/dqxxkx.2021.200291

JIANG Yilan, CHEN Baowang, HUANG Yufang, CUI Jiaqi, GUO Yulong. Crop Planting Area Extraction based on Google Earth Engine and NDVI Time Series Difference Index[J]. Journal of Geo-information Science, 2021, 23(5): 938-947.DOI:10.12082/dqxxkx.2021.200291

图/表 9

图1

表1

表2

图2

图3

图4

图5

图6

表3

参考文献 30

[1]	白燕英, 高聚林, 张宝林. 基于Landsat8影像时间序列NDVI的作物种植结构提取[J]. 干旱区地理, 2019,42(4):893-901.
	[ Bai Y Y, Gao J L, Zhang B L. Extraction of crop planting structure based on time-series NDVI of Landsat8 images[J]. Arid Land Geography, 2019,42(4):893-901. ]
[2]	Jia K, Wu B F, Li Q Z. Crop classification using HJ satellite multispectral data in the North China Plain[J]. Journal of Applied Remote Sensing, 2013,7(1):073576. doi: 10.1117/1.JRS.7.073576
[3]	贾坤, 李强子, 田亦陈, 等. 微波后向散射数据改进农作物光谱分类精度研究[J]. 光谱学与光谱分析, 2011,31(2):483-487.
	[ Jia K, Li Q Z, Tian Y C, et al. Accuracy improvement of spectral classification of crop using microwave backscatter data[J]. Spectroscopy and Spectral Analysis, 2011,31(2):483-487. ]
[4]	胡琼, 吴文斌, 宋茜, 等. 农作物种植结构遥感提取研究进展[J]. 中国农业科学, 2015,48(10):1900-1914.
	[ Hu Q, Wu W B, Song Q, et al. Recent progresses in research of crop patterns mapping by using remote sensing[J]. Scientia Agricultura Sinica, 2015,48(10):1900-1914. ]
[5]	马丽, 徐新刚, 刘良云, 等. 基于多时相NDVI及特征波段的作物分类研究[J]. 遥感技术与应用, 2008,23(5):520-524.
	[ Ma L, Xu X G, Liu L Y, et al. Study on crops classification based on multi-temporal NDVI and characteristic bands[J]. Remote Sensing Technology and Application, 2008,23(5):520-524. ]
[6]	康峻, 侯学会, 牛铮, 等. 基于拟合物候参数的植被遥感决策树分类[J]. 农业工程学报, 2014,30(9):148-156.
	[ Kang J, Hou X H, Niu Z, et al. Decision tree classification based on fitted phenology parameters from remotely sensed vegetation data[J]. Transactions of the Chinese Society of Agricultural Engineering, 2014,30(9):148-156. ]
[7]	张焕雪, 曹新, 李强子, 等. 基于多时相环境星NDVI时间序列的农作物分类研究[J]. 遥感技术与应用, 2015,30(2):304-311.
	[ Zhang H X, Cao X, Li Q Z, et al. Research on crop identification using multi-temporal NDVI HJ images[J]. Remote Sensing Technology and Application, 2015,30(2):304-311. ]
[8]	李晓东, 姜琦刚. 基于多时相遥感数据的农田分类提取[J]. 农业工程学报, 2015,31(7):145-150.
	[ Li X D, Jiang Q G. Extraction of farmland classification based on multi-temporal remote sensing data[J]. Transactions of the Chinese Society of Agricultural Engineering, 2015,31(7):145-150. ]
[9]	陈君颖, 田庆久. 高分辨率遥感植被分类研究[J]. 遥感学报, 2007,11(2):221-227.
	[ Chen J Y, Tian Q J. Vegetation classification based on high-resolution satellite image[J]. Journal of Remote Sensing, 2007,11(2):221-227. ]
[10]	Bolton, Douglas K, Josh M Gray, et al. Continental-scale land surface phenology from harmonized Landsat 8 and Sentinel-2 imagery[J]. Remote Sensing of Environment, 2020,240:111685. doi: 10.1016/j.rse.2020.111685
[11]	Arvor D, Jonathan M, Meirelles M S P, et al. Classification of MODIS EVI time series for crop mapping in the state of Mato Grosso, Brazil[J]. International Journal of Remote Sensing, 2011,32(22):7847-7871. doi: 10.1080/01431161.2010.531783
[12]	郝鹏宇, 牛铮, 王力, 等. 基于历史时序植被指数库的多源数据作物面积自动提取方法[J]. 农业工程学报, 2012,28(23):123-131,297.
	[ Hao P Y, Niu Z, Wang L, et al. Multi-source automatic crop pattern mapping based on historical vegetation index profiles[J]. Transactions of the Chinese Society of Agricultural Engineering, 2012,28(23):123-131,297. ]
[13]	潘耀忠, 李乐, 张锦水, 等. 基于典型物候特征的MODIS-EVI时间序列数据农作物种植面积提取方法——小区域冬小麦实验研究[J]. 遥感学报, 2011,15(3):578-594.
	[ Pan Y Z, Li L, Zhang J S, et al. Crop area estimation based on MODIS-EVI time series according to distinct characteristics of key phenology phases: A case study of winter wheat area estimation in small-scale area[J]. Journal of Remote Sensing, 2011,15(3):578-594. ]
[14]	陈安旭, 李月臣. 基于Sentinel-2影像的西南山区不同生长期水稻识别[J]. 农业工程学报, 2020,36(7):192-199.
	[ Chen A X, Li Y C. Rice recognition of different growth stages based on Sentinel-2 images in mountainous areas of Southwest China[J]. Transactions of the Chinese Society of Agricultural Engineering, 2020,36(7):192-199. ]
[15]	刘佳, 王利民, 姚保民, 等. 基于多时相OLI数据的宁夏大尺度水稻面积遥感估算[J]. 农业工程学报, 2017,33(15):200-209.
	[ Liu J, Wang L M, Yao B M, et al. Ningxia rice area remote sensing estimation on large scale based on multi-temporal OLI data[J]. Transactions of the Chinese Society of Agricultural Engineering, 2017,33(15):200-209. ]
[16]	黄健熙, 贾世灵, 武洪峰, 等. 基于GF-1 WFV影像的作物面积提取方法研究[J]. 农业机械学报, 2015,46(S1):253-259.
	[ Huang J X, Jia S L, Wu H F, et al. Extraction method of crop planted area based on GF-1 WFV image[J]. Transactions of the Chinese Society for Agricultural Machinery, 2015,46(S1):253-259. ]
[17]	刘佳, 王利民, 杨福刚, 等. 基于HJ时间序列数据的农作物种植面积估算[J]. 农业工程学报, 2015,31(3):199-206.
	[ Liu J, Wang L M, Yang F G, et al. Remote sensing estimation of crop planting area based on HJ time-series images[J]. Transactions of the Chinese Society of Agricultural Engineering, 2015,31(3):199-206. ]
[18]	周静平, 李存军, 史磊刚, 等. 基于决策树和面向对象的作物分布信息遥感提取[J]. 农业机械学报, 2016,47(9):318-326,333.
	[ Zhou J P, Li C J, Shi L G, et al. Crops distribution information extracted by remote sensing based on decision tree and object-oriented method[J]. Transactions of the Chinese Society of Agricultural Machinery, 2016,47(9):318-326,333. ]
[19]	康峻, 侯学会, 牛铮, 等. 基于拟合物候参数的植被遥感决策树分类[J]. 农业工程学报, 2014,30(9):148-156.
	[ Kang J, Hou X H, Niu Z, et al. Decision tree classification based on fitted phenology parameters from remotely sensed vegetation data[J]. Transactions of the Chinese Society of Agricultural Engineering, 2014,30(9):148-156. ]
[20]	李鑫川, 徐新刚, 王纪华, 等. 基于时间序列环境卫星影像的作物分类识别[J]. 农业工程学报, 2013,29(2):169-176,298.
	[ Li X C, Xu X G, Wang J H, et al. Crop classification recognition based on time-series images from HJ satellite[J]. Transactions of the Chinese Society of Agricultural Engineering, 2013,29(2):169-176,298. ]
[21]	张荣群, 王盛安, 高万林, 等. 基于时序植被指数的县域作物遥感分类方法研究[J]. 农业机械学报, 2015,46(S1):246-252.
	[ Zhang R Q, Wang S A, Gao W L, et al. Remote-sensing classification method of County-level agricultural crops using time-series NDVI[J]. Transactions of the Chinese Society for Agricultural Machinery, 2015,46(S1):246-252. ]
[22]	杜保佳, 张晶, 王宗明, 等. 应用Sentinel-2A NDVI时间序列和面向对象决策树方法的农作物分类[J]. 地球信息科学学报, 2019,21(5):740-751. doi: 10.12082/dqxxkx.2019.180412
	[ Du B J, Zhang J, Wang Z M, et al. Crop mapping based on sentinel-2A NDVI time series using object-oriented classification and decision tree model[J]. Geo-Information Science, 2019,21(5):740-751. ]
[23]	杨闫君, 占玉林, 田庆久, 等. 基于GF-1/WFVNDVI时间序列数据的作物分类[J]. 农业工程学报, 2015,31(24):155-161.
	[ Yang Y J, Zhan Y L, Tian Q J, et al. Crop classification based on GF-1/WFV NDVI time series[J]. Transactions of the Chinese Society of Agricultural Engineering, 2015,31(24):155-161. ]
[24]	张焕雪, 曹新, 李强子, 等. 基于多时相环境星NDVI时间序列的农作物分类研究[J]. 遥感技术与应用, 2015,30(2):304-311.
	[ Zhang H X, Cao X, Li Q Z, et al. Research on crop identification using multi-temporal NDVI HJ images[J]. Remote Sensing Technology and Application, 2015,30(2):304-311. ]
[25]	刘焕军, 闫岩, 张新乐, 等. 面向农业区划的作物种植结构遥感提取[J]. 中国农业资源与区划, 2017,38(8):43-54.
	[ Liu H J, Yan Y, Zhang X L, et al. Remotesensing extraction of crop planting structure oriented to agricultural regionalization[J]. Chinese Journal of Agricultural Resources and Regional Planning, 2017,38(8):43-54. ]
[26]	王林江, 吴炳方, 张淼, 等. 关键生育期冬小麦和油菜遥感分类方法[J]. 地球信息科学学报, 2019,21(7):1121-1131. doi: 10.12082/dqxxkx.2019.180421
	[ Wang L J, Wu B F, Zhang M, et al. Winter wheat and rapeseed classification during key growth period by integrating multi-source remote sensing data[J]. Journal of Geo-information Science, 2019,21(7):1121-1131. ]
[27]	徐涵秋. 一种快速提取不透水面的新型遥感指数[J]. 武汉大学学报·信息科学版, 2008,33(11):1150-1153,1211.
	[ Xu H Q. A new remote sensing index for fastly extracting impervious surface information[J]. Geomatics and Information Science of Wuhan University, 2008,33(11):1150-1153,1211. ]
[28]	赵英时. 遥感应用分析原理与方法[M]. 北京: 科学出版社, 2003.
	[ Zhao Y S. Principles and methods of remote sensing application analysis[M]. Beijing: Science Press, 2003. ]
[29]	王利军, 郭燕, 贺佳, 等. 基于决策树和SVM的Sentinel-2A影像作物提取方法[J]. 农业机械学报, 2018,49(9):146-153.
	[ Wang L J, Guo Y, He J, et al. Classification method by fusion of decision tree and SVM based on Sentinel-2A image[J]. Transactions of the Chinese Society for Agricultural Machinery, 2018,49(9):146-153. ]
[30]	Congalton, Russell G. A review of assessing the accuracy of classifications of remotely sensed data[J]. Remote Sensing of Environment, 1991,37(1):35-46. doi: 10.1016/0034-4257(91)90048-B

	属性	中心波长/μm	分辨率/m
Band 1	气溶胶	0.443	60
Band 2	蓝	0.490	10
Band 3	绿	0.560	10
Band 4	红	0.665	10
Band 5	植被红边1	0.705	20
Band 6	植被红边2	0.740	20
Band 7	植被红边3	0.783	20
Band 8	近红外	0.842	10
Band 8A	植被红边4	0.865	20
Band 9	水蒸气	0.945	60
Band 10	短波红外-卷云	1.375	60
Band 11	短波红外1	1.610	20
Band 12	短波红外2	2.190	20
QA60	云掩膜	-	-

	作物类型
	小麦	大蒜
1月	越冬	幼苗期
2月	越冬	幼苗期
3月	返青	幼苗期
4月	拔节孕穗	花芽和鳞茎分化期
5月	抽穗灌浆	鳞茎肥大期
6月	成熟	休眠期
7月	-	休眠期
8月	-	休眠期
9月	-
10月	播种出苗	播种发芽
11月	-	幼苗期
12月	分蘖	幼苗期

分类方法	大蒜		小麦		总体精度/%	Kappa 系数
分类方法	用户精度/%	制图精度/%	用户精度/%	制图精度/%	总体精度/%	Kappa 系数
NTDI指数法	88.94	82.59	77.44	85.33	83.72	0.67
最大似然法	77.78	77.40	67.94	68.43	73.71	0.46
支持向量机法	91.74	61.14	75.19	95.28	79.54	0.58

基于Google Earth Engine和NDVI时序差异指数的作物种植区提取

Crop Planting Area Extraction based on Google Earth Engine and NDVI Time Series Difference Index

RichHTML

PDF (PC)

赞

可视化

摘要/Abstract

引用本文

使用本文

图/表 9

参考文献 30

相关文章 15

Metrics

本文评价

推荐阅读 0

[1]	葛中曦, 黄静, 赖佩玉, 郝斌飞, 赵银军, 马明国. 耕地复种指数遥感监测研究进展[J]. 地球信息科学学报, 2021, 23(7): 1169-1184.
[2]	冯子钰, 施润和. 中国近地面PM_2.5浓度与排放的时空分布及其关联分析[J]. 地球信息科学学报, 2021, 23(7): 1221-1230.
[3]	熊皓丽, 周小成, 汪小钦, 崔雅君. 基于GEE云平台的福建省10 m分辨率茶园专题空间分布制图[J]. 地球信息科学学报, 2021, 23(7): 1325-1337.
[4]	陈赛楠, 蒋弥. Sentinel-1 SAR在洪水范围提取与极化分析中的应用研究[J]. 地球信息科学学报, 2021, 23(6): 1063-1070.
[5]	管琪卉, 丁明军, 张华, 王鹏. ESTARFM算法在长江中下游平原地区的适用性研究[J]. 地球信息科学学报, 2021, 23(6): 1118-1130.
[6]	彭妍菲, 李忠勤, 姚晓军, 牟建新, 韩伟孝, 王盼盼. 基于多源遥感数据和GEE平台的博斯腾湖面积变化及影响因素分析[J]. 地球信息科学学报, 2021, 23(6): 1131-1153.
[7]	李宇君, 张磊. 基于沙地指数模型的沙地监测方法[J]. 地球信息科学学报, 2021, 23(4): 680-691.
[8]	谢聪慧, 吴世新, 张晨, 孙文涛, 何海芳, 裴韬, 罗格平. 基于谱系聚类的全球各国新冠疫情时间序列特征分析[J]. 地球信息科学学报, 2021, 23(2): 236-245.
[9]	陈如如, 胡中民, 李胜功, 郭群. 不同数据源归一化植被指数在中国北方草原区的应用比较[J]. 地球信息科学学报, 2020, 22(9): 1910-1919.
[10]	刘新, 赵宁, 郭金运, 郭斌. 基于LSTM神经网络的青藏高原月降水量预测[J]. 地球信息科学学报, 2020, 22(8): 1617-1629.
[11]	李培林, 刘小平, 黄应淮, 张鸿辉. 基于GEE平台的广州市主城区不透水面时间序列提取[J]. 地球信息科学学报, 2020, 22(3): 638-648.
[12]	杨丹, 周亚男, 杨先增, 郜丽静, 冯莉. LSTM支持下时序Sentinel-1A数据的太白山区植被制图[J]. 地球信息科学学报, 2020, 22(12): 2445-2455.
[13]	郑诗晨, 盛业华, 吕海洋. 基于粒子滤波的行车轨迹路网匹配方法[J]. 地球信息科学学报, 2020, 22(11): 2109-2117.
[14]	刘怡媛, 李鹏, 肖池伟, 刘影, 谢正磊. 老挝北部Sentinel-2 A/B云覆盖适宜阈值界定及其影像获取概率时空差异[J]. 地球信息科学学报, 2020, 22(11): 2267-2276.
[15]	周玲, 李云梅, 赵焕, 房晟忠, 颜翔, 吕恒. 滇池湖滨带植被时空变化及其对滇池总磷浓度的影响[J]. 地球信息科学学报, 2020, 22(10): 1946-1958.