GEE支持下的河南省冬小麦面积提取及长势监测

2021-06-16周珂柳乐张俨娜苗茹杨阳

中国农业科学 2021年11期

周珂，柳乐，张俨娜，苗茹，杨阳

1河南大学计算机与信息工程学院，河南开封 475004；2河南大学实验室与设备管理处，河南开封 475004；3河南省大数据分析与处理重点实验室/河南大学，河南开封 475004

【】使用遥感技术对2017—2020年河南省冬小麦的空间分布信息进行高精度的提取，然后对2020年冬小麦的长势进行高频度的监测并结合气象条件进行分析。本文基于谷歌地球引擎（Google Earth Engine，GEE）云平台，对选取的Landsat 8影像数据根据NDVI最大值进行合成，然后进行特征构建，添加地形特征、纹理特征、NDVI以及一个新特征NDVI增幅，使用随机森林分类方法对样本数据按照构建的特征进行训练提取河南省2017—2020年冬小麦的播种面积信息；经过精度验证后对提取的河南省2020年的冬小麦种植区域生成掩膜，对掩膜区域（冬小麦种植区域）结合MODIS高时间分辨率影像数据，使用NDVI同期差值法对2020年2—4月份的冬小麦进行高频度的长势监测。使用GEE云平台能够对河南省冬小麦种植区域的空间分布信息进行快速制图；使用随机森林方法加入地形特征、纹理特征、NDVI后再加入新特征NDVI增幅，能够有效提高冬小麦的提取精度以及降低与统计数据的相对误差，基于混淆矩阵计算的平均总体分类精度为95.2%、平均kappa系数为0.909、冬小麦的平均分类精度为95.3%，与河南省统计年鉴数据相比，本文方法提取的2017—2019年河南省冬小麦播种面积相对误差均低于3%，河南省冬小麦主要种植区域的冬小麦播种面积的平均相对误差低于6%；使用MODIS影像数据结合NDVI差值模型能够对河南省2020年的冬小麦进行高频度的长势监测，河南省冬小麦在返青初期长势较往年及2019年好，到生育后期大部分区域长势与往年及2019年持平，总体上2020年冬小麦的长势较往年及2019年好。本文提出的方法能够对河南省冬小麦进行高精度的提取以及高频度的长势监测，且能够为地方政府或者一些农业部门在安排指导农事活动上提供科学依据。

冬小麦；长势监测；谷歌地球引擎；随机森林；归一化植被指数；Landsat；MODIS

0 引言

【研究意义】冬小麦作为我国三大主粮作物之一，在国民经济中占据着重要地位[1]。河南省作为我国农业大省之一，以冬小麦种植为主，其播种面积占据全省粮食作物的54%，且其产量多年一直保持全国第一，占全国的20%以上[2]，及时获知河南省冬小麦的分布及长势情况对确保粮食产量具有重要意义。【前人研究进展】卫星遥感数据具有覆盖面广、光谱信息丰富、周期性强等特点，目前已成为监测冬小麦分布及长势的重要手段之一[3-10]。常用的提取冬小麦播种面积的手段主要包括：（1）结合遥感影像，利用NDVI等植被指数的时序变化来设置合适阈值，对冬小麦的空间分布信息进行提取[11-14]；（2）使用机器学习的方法，对遥感影像各类型地物进行分类，进而提取冬小麦的播种面积[15-19]。在机器学习方法中随机森林（RF）自动化程度高、较容易实现，且其计算开销小[20]，在遥感影像分类上被广泛应用。在冬小麦的空间分布信息提取上，李旭青等[17]、Liu等[19]使用RF分类算法对研究区域冬小麦的种植面积进行提取，均取得了较好的效果。在使用机器学习方法对遥感影像进行分类时，You等[21]、杨蕙宇等[22]、何昭欣等[23]的研究数据表明，使用机器学习方法分类时加入地形、纹理特征能够提高分类精度。近年来，随着信息技术的不断发展，云计算也快速进入卫星遥感数据处理领域，并取得了很好的效果。Google Earth Engine（GEE）是谷歌公司开发的一款用于遥感数据分析的工具，有着强大的数据处理、分析、存储和可视化能力[24]。国内遥感领域的科学研究者基于GEE云平台已进行了一些在遥感科学领域的研究。何昭欣等[23]使用GEE云平台对江苏省夏收作物进行较好的遥感提取，且使用GEE处理较本地处理有明显优势；朱德海等[25]使用GEE对山东省近30年来的农业大棚进行较高精度的提取。Liu等[26]使用GEE中大量Landsat影像进行全球城市土地分类，并绘制了多时相的全球城市土地利用类型图。郝斌飞等[27]对GEE进行分析，数据表明GEE是一个高效的科研工具，在相关领域的研究上粗略统计其综合效率（成本、效益等）能提升90%以上。【本研究切入点】目前结合云平台以及机器学习的方法在省域范围内提取冬小麦的研究还较少。在冬小麦长势监测方面[28-29]，较多研究使用MODIS影像数据，基于植被指数时序变化进行面积提取并构建长势检测模型[28-29]，而MODIS影像数据存在空间分辨率较低会产生大量混合像元的问题。【拟解决的关键问题】本文以河南省为研究区，基于Landsat 8 OLI、MODIS两种影像各自的优势，在GEE云平台上首先结合Landsat 8 OLI影像数据对训练特征进行构建，以机器学习的方法提取2017—2020年河南省冬小麦播种面积及其空间分布，然后使用MODIS MO9GQ产品数据对2020年2、3月份以及4月上中旬河南省的冬小麦进行长势监测及分析，为实现高频度的中高分辨率影像协同支持的农作物动态监测，进而构建基于云平台和机器学习的农作物信息提取系统以及农业生产提供参考依据和客观数据支撑。

1 材料与方法

1.1 研究区域与数据源

1.1.1 研究区域概况河南省位于中原腹地，介于31°23′—36°22′N、110°21′—116°39′E之间（图1）。河南省在地势上西高东低，中东部地区多为平原地带，西部多为山区丘陵地带。气候四季分明，属于亚热带和暖温带季风气候，夏季高温且多雨，冬季寒冷且干燥。河南省全年无霜期为201—285 d，全省由南向北的年平均气温为10.5—16.7℃，降水以6—8月份最多，全省年平均降水量为407.7—1 295.8 mm，年平均日照为1 285.7—2 292.9 h，适宜多种农作物的生长[30]。

图1 研究区域地形图

1.1.2 影像数据集及预处理本研究选用Landsat 8 OLI影像数据以及MODIS MO9GQ影像数据提取河南省冬小麦种植面积并进行长势监测。Landsat 8搭载OLI陆地成像仪，包括9个波段，Landsat 8卫星大约每2周对整个地球表面成像，包括多光谱和热数据。MODIS MO9GQ为地表反射率产品，空间分辨率为250 m，时间分辨率为1 d。Landsat 8影像数据具有较高的空间分辨率，为30 m，能够减少混合像元对冬小麦播种面积提取的影响；MODIS MO9GQ影像数据有很高的时间分辨率，能够对冬小麦长势进行持续跟踪。

河南省的冬小麦一般在秋末冬初进行播种，下一年的夏季进行收割，生育周期如表1所示。冬小麦在播种期到起身拔节期的NDVI变化较大，所以选取此区间的影像进行面积提取，考虑到受云量及其他因素的影响，同一个时间区间获取的影像可能不完整，而相邻年份同时间段的冬小麦种植区域差异不是太大，所以选取相邻年份的影像数据作为补充，具体选取的Landsat 8影像数据详情如表2所示。2月至4月上中旬，冬小麦从返青期到孕穗—抽穗期，在这个时间区段冬小麦的颜色是在不断变绿的，结合NDVI能够有效对河南省的冬小麦长势进行监测，因此选取2017—2020年2—3月份以及4月上中旬的MODIS MO9GQ影像数据，对研究区域2019—2020年冬小麦长势进行监测。选取影像后对其进行裁剪、镶嵌等数据预处理，所有操作均在GEE中进行。

1.1.3 样本数据及验证数据根据2019—2020年12—3月份合成影像的各类地物不同的纹理信息，结合Google Earth中的高分辨率影像，在研究区域内共选取3 889个样本点。样本数据分布均匀，覆盖全省范围，其中冬小麦样本点2 554个，非冬小麦样本点1 335个，非冬小麦样本包括水体、不透水面以及其他植被。样本选取依据如表3所示。验证数据来自河南省统计年鉴（http://www.henan.gov.cn/zwgk/zfxxgk/fdzdgknr/ tjxx/tjnj/）发布的冬小麦播种面积数据。

表1 河南省冬小麦生育周期

表2 Landsat 8选取影像详情

表3 样本选取依据

1.2 研究方法

1.2.1 技术流程本文基于GEE云平台Landsat影像数据，使用机器学习分类方法提取研究区2017—2020年的冬小麦播种面积，结合MODIS影像数据，利用NDVI差值模型对研究区2020年2、3月份以及4月上中旬的冬小麦长势进行高频度监测。步骤如下：（1）从GEE云平台获取Landsat影像数据集，进行筛选、影像镶嵌等数据预处理操作；（2）按照NDVI最大值镶嵌影像；（3）特征构建；（4）对遥感影像分类提取研究区冬小麦播种面积；（5）精度验证；（6）对监测时段的MODIS影像数据计算NDVI，并对每旬影像进行NDVI最大值合成，然后进行差值计算；利用同期对比法对2020年2、3月份的冬小麦进行长势监测。流程图如图2所示。

1.2.2 冬小麦种植面积信息获取本研究使用GE高分辨率影像选取4类地物的样本点，结合Landsat 8影像绘制各类地物在9月份到次年4月份的NDVI变化（图3）。从图3中可以看出，研究区冬小麦从12月到次年4月的NDVI值在不断增长且明显高于其他地类，此时间区间为冬小麦的返青及起身期，冬小麦在不断变绿。为了在影像上突出冬小麦的特征，本次研究在影像集合成上做了一些改变。本研究选择12月份到次年3月份冬小麦返青到起身期的影像数据，计算每景影像的NDVI值，在GEE中对影像集的每个像元按照NDVI值从小到大进行排序，提取根据NDVI最大值合成的影像（包含原始光谱信息）。NDVI计算公式如下：

式中，为红波段，为近红外波段。

河南省自西向东地势从高变低，在西部地区由于地势高山体丘陵较多，冬小麦种植区域较少且较难识别。为了提高冬小麦种植区域的分类精度，在训练特征中加入了地形特征和纹理特征。本文研究利用航天飞机雷达地形任务[31]（shuttle radar topography mission，SRTM）的SRTMG数据产品获取研究区的地形特征，其分辨率为30 m。

纹理信息是图像的一种重要特征，刻画了图像像素邻域灰度空间分布的规律，不同的地物类型在遥感影像上有不同的纹理信息。冬小麦在遥感影像上具有连续的、规则的纹理，为了更好地提取研究区域的冬小麦空间分布信息，本研究在GEE上使用灰度共生矩阵（gray-level co-occurrence matrix，GLCM）构建纹理特征。由于NDVI值能够很好地反映冬小麦的特征，所以本研究在影像合成后，使用NDVI计算影像的纹理特征，并选择4种相关性较低的纹理信息[32]。分别是角二阶矩（ASM）、逆差距（IDM）、对比度（CON）、相关性（CORR），计算公式如下：

图2 冬小麦长势监测流程

图3 研究区各类地物时序NDVI变化

式中，为原始图像的灰度级数，(,)为基于原始图像生成的灰度共生矩阵，为均值，为方差。

除上述特征外，本文加入了一个新的特征NDVIincrease，NDVIincrease由王九中等[33]提出。结合图3和冬小麦的生育周期中可以发现，冬小麦播种初期到冬小麦起身拔节期NDVI有明显变化，增长幅度较其他3类地物高，基于这种变化计算冬小麦从播种期到起身拔节期NDVI增长的幅度，计算公式如下：

式中，NDVImin为播种时节9月份到11月份影像集合的最小NDVI合成影像；NDVImax为12月份到次年3月份影像集合的最大NDVI合成影像。

在分类方法上，本研究使用机器学习中的随机森林分类（random forest，RF）算法对影像进行分类，提取研究区2017— 2020年冬小麦的播种面积。RF是由BREIMAN等[34]提出的包含多棵决策树的分类方法，其随机主要体现在两方面：一方面是子模型的训练样本是有放回的随机抽取；另一方面是子模型的特征变量也是随机抽取。随机森林的每一棵决策树都是一个分类器，当输入一个样本时，随机森林的任一棵树都会产生相对应的分类结果，而随机森林收集了所有子树的分类结果，将涵盖子树最多的类别作为最终的输出结果。本研究在GEE平台对构建的特征以独立的光谱波段添加到影像以及样本数据中，然后使用RF分类方法提取研究区冬小麦的种植面积以及空间分布情况。参数设置上，综合考虑冬小麦提取的精度、与统计数据的误差以及计算效率，将决策树的数量设置为500，最大叶节点数设为无限制，分割节点的特征数量设为输入特征总数的平方根。

1.2.4 精度验证获得研究区冬小麦播种面积后，对2017—2019年遥感提取的冬小麦播种面积与河南省统计年鉴统计的冬小麦播种面积数据进行比较，计算绝对误差以及相对误差。相对误差能够客观描述提取精度且能够更好地反映提取的冬小麦播种面积的可信度，相对误差越小表示面积提取精度越高。绝对误差和相对误差的计算公式如下：

式中，为绝对误差，I为遥感提取面积，I为统计的播种面积，为相对误差。

除相对误差验证外，本研究结合研究区域范围及算法本身的特点，将选取的样本数据按照4﹕1的比例进行分配，在样本数据中随机选取80%的样本数据作为训练数据，剩下的20%的样本数据作为验证数据用来验证分类的精度。本研究在对影像分类后使用混淆矩阵对验证数据的识别结果进行检验，计算总体分类精度（）、系数以及冬小麦的分类精度（），计算公式如下：

式中，为样本种类的个数，混淆矩阵，主对角线位置为验证样本数据中各个类别正确分类的个数。

2 结果

2.1 面积监测结果

本文通过GEE云平台对选取的Landsat 8影像进行特征构建，使用RF分类对2017—2020年河南省冬小麦的种植面积进行计算，并与河南省统计局统计的冬小麦种植面积进行比较计算相对误差（表4），发现加入NDVI增幅后，2017—2019年的冬小麦的提取相对误差有所降低，且相对误差均在3%以内，在省域范围内提取误差较小。为了更好地检验此方法提取的精确性，对河南省冬小麦主要种植区域的冬小麦种植面积进行提取并计算相对误差（表5），结果表明加入NDVI增幅这一特征后，提取2017—2019年河南省冬小麦主要种植区域的冬小麦的相对误差有很大降低，2017—2018年冬小麦的平均相对误差从12.10%降低到5.69%，2018—2019年冬小麦的平均相对误差从10.84%降低到5.49%，在市域范围上提取效果也较好。除上述验证外，本文基于混淆矩阵对验证样本计算分类后的总体分类精度、kappa系数以及冬小麦分类精度（表6），发现加入NDVI增幅后，总体精度、kappa系数、冬小麦分类精度均有提高。通过以上分析可以看出本研究所使用的方法提取河南省冬小麦的精度较高、效果较好。

本研究使用GEE对河南省2017—2020年的冬小麦播种面积区域进行提取，其能够快速对遥感影像的数据预处理进行操作，并能够对河南省冬小麦种植区域的空间分布信息进行快速制图。河南省2017—2020年冬小麦的播种面积的空间信息如图4所示，从中可以看出冬小麦主要分布在河南中东部平原和南阳盆地，由于西部山体丘陵较多以及信阳地区靠南，小麦种植较少。从表4可以看出遥感所测的2017—2019年河南省冬小麦播种面积处于增长状态，与河南省统计局统计的冬小麦播种面积变化一致，且2020年遥感所测的冬小麦面积也在增加，因此预计河南省2020年冬小麦的播种面积要高于往年。

表4 河南省冬小麦遥感提取的播种面积与统计的播种面积比较

表5 河南省冬小麦主要种植区域的遥感提取的播种面积与统计的播种面积比较

表6 研究区冬小麦基于混淆矩阵提取精度比较

2.2 长势监测结果

本研究所使用的方法对河南省冬小麦的播种面积提取精度较高，提取结果可靠，所以将河南省2020年的冬小麦种植面积分布应用式（10）的NDVI差值模型与往年同期进行对比，统计NDVI对应增量的像元个数（图5），发现与往年同期相比，从2月上旬到3月中旬，NDVI增量大于0的像元个数明显多于NDVI增量小于0的像元个数，说明冬小麦种植区域绝大数区域的NDVI值大于往年，在2月上旬到3月中旬冬小麦长势情况较往年好；3月下旬到4月中旬的NDVI增量大于0与小于0的像元个数相差不大，为了更好地观察此时间段的冬小麦长势，本文应用式（11）对NDVI差值结果进行分类，发现从3月下旬到4月中旬，曲线变化幅度较小，NDVI增量在-0.1到0.1区间内的像元个数最多。王利民等[35]、孙丽等[28]、黄青等[29]研究数据表明在冬小麦生育后期大部分长势与往年持平，因此本文将式（11）中的n1设为-0.1，n2设为0.1，根据设定的值对2020年冬小麦与往年及2019年相比各类长势的面积占比进行统计（图6—7），并与往年同期相比冬小麦各类长势情况分布区域进行制图（图8）。结果发现，与2019年相比，2020年2月上旬到3月上旬冬小麦长势好于2019年同期；从3月中旬到4月中旬，长势较好的区域渐渐减少，较差的区域增多，但总体上依然是长势好的区域高于长势差的区域。与往年相比，2020年2月上旬到3月中旬大部分区域冬小麦的长势好于往年，这个结果与图5统计的结果一致，在3月下旬到4月中旬，长势较好的区域大幅减少，90%以上的冬小麦种植区域长势与往年持平，但总体上长势好的区域多于长势差的区域。

图5 2020年冬小麦种植区域与往年相比的NDVI变化

图6 与2019年相比研究区2020年2—4月的冬小麦的长势情况

图7 与往年相比研究区2020年2—4月的冬小麦的长势情况

图8 与往年相比研究区2020年2—4月的冬小麦的长势分布

气候变化对冬小麦的长势有很大影响，因此本文对河南省气象局现代农业网（http://www.hnnw. net/）发布的2017—2020年2—4月份的农业气象周报数据进行统计（图9），然后根据统计数据，对河南省冬小麦的长势做出分析。从图5—8可以看出，在2月上旬到2月下旬，无论是与2019同期相比还是与往年同期相比，2020年冬小麦的长势要远好于2019年及往年同期，通过图9对比可以看出2020年2月份气温较2019年及往年同期要高，且有充足的光照时间，有利于冬小麦生长，因此冬小麦进入返青期要快于往年及2019年；3月下旬到4月中旬期间，2020年与往年及2019年同期相比气温在降低，光照时间在不断变化，但光照时间总体较往年及2019年短，因此冬小麦长势放缓。虽然冬小麦长势受到气候的影响，但总体上2020年冬小麦的长势依旧要好于往年及2019年。

通过2020年2—4月份的冬小麦长势与往年及2019年对比分析，发现河南省2020年的冬小麦长势相比2019年要好。在冬小麦产量方面，盛磊等[2]研究表明在冬小麦从返青期到抽穗期植被指数与产量的相关系数最高，千怀遂[36]研究认为小麦单产的最佳时相为拔节后期到抽穗期，而本文的研究数据表明无论是从返青期到抽穗期，还是拔节后期到抽穗期，NDVI增量大于0的区域要大于NDVI增量小于0的区域，冬小麦长势好的区域要多于冬小麦长势差的区域，因此判断河南省2020年冬小麦的总体产量要好于往年及2019年。

图9 2017—2020年2—4月份的平均气温及周平均日照

3 讨论

本研究所使用的方法能够很好地对河南省2017—2020年的冬小麦的空间分布信息进行快速提取并制图，并对2020年的冬小麦进行高频度的长势监测，但在平台以及方法上还存在一些不足之处。

在冬小麦种植信息提取方面，本研究使用Landsat 8影像数据构建训练特征，通过RF分类方法提取河南省2017—2020年的冬小麦种植空间分布信息，与李方杰等[37]、邓荣鑫等[38]的研究方法相比，本文使用中高分辨率Landsat 8 影像数据在影像合成上突出了冬小麦的特征，且在模型训练时加入纹理特征、地形特征、NDVI以及一个新特征NDVI增幅，提取的冬小麦种植区域存在的混合像元较少，基于混淆矩阵计算平均总体精度有所提高，且与统计数据相比总体相对误差也较小。目前研究虽然整体提取效果较好，但在冬小麦种植区域较多的地区（如驻马店、周口）提取的面积还存在较大误差，Landsat 8影像为中高分辨率影像，也不可避免地存在一些混合像元，此外，原始光谱波段间还存在较大的相关性，下一步研究方向将结合更高空间分辨率影像与高时间分辨率影像进行融合，然后对原始光谱特征进行优化，使用深度学习方法对更大范围内冬小麦的种植面积进行提取。

在冬小麦长势监测方面，孙丽等[28]、黄青等[29]基于MODIS影像数据，使用NDVI阈值法提取冬小麦的空间分布信息，在此基础上构建模型对研究区域的冬小麦进行长势监测，此方法提取的冬小麦种植区域存在较多的混合像元以及自动化程度不高。为了更好地提取冬小麦的种植区域，从而对冬小麦进行更加准确地监测，本研究在Landsat 8合成影像上，使用随机森林方法对河南省冬小麦进行较高精度地提取，对提取的区域生成掩膜，对掩膜区域结合MODIS影像数据，利用NDVI差值法进行长势监测。试验表明，本研究使用的方法既能够获取较高精度的冬小麦种植分布信息，又能进行高频度的长势监测。对于冬小麦长势监测的方法，本文还存在一定的不足，本文利用每旬MODIS合成影像的NDVI与次年以及往年进行比较，在部分地区偶尔会出现NDVI的误差情况，下一步将考虑对NDVI时序曲线进行修正，进一步提高监测的精度；关于苗情分级，较多的文献是根据农学和地面实测数据进行分级分类，部分借助于遥感手段的方式方法比较复杂，增加了苗情判读的复杂度，本文尝试在尽可能提高精度的情况下简化苗情分类判读方式，计划在下一步的研究中，通过试验对比结合数据分析寻求较为实用的长势评估手段。

GEE能够高效地处理海量遥感数据，使用GEE云平台能够对河南省2017—2020年的冬小麦进行提取并对冬小麦种植区域进行快速制图。但由于使用的是GEE平台免费资源，所以提供的计算能力有限，当研究区域过大时，在计算过程中偶尔会出现“计算超时”“计算超出限制”或者浏览器崩溃等错误，随着后续的改进，会逐渐完善。在应用价值上，本文以GEE云平台的海量遥感数据为支撑，开展了小麦长势的监测，具有一定的实用性，在后续工作中可以以本文工作为基础，基于GEE云平台开发供农业统计遥感的业务化运行系统，以利于提升监测分析的效率且为地方政府或一些农业部门在农事活动的安排指导上提供科学依据。

4 结论

本研究基于GEE云平台，使用Landsat 8影像数据对特征进行构建，使用RF分类法提取河南省2017—2020年冬小麦面积，然后对河南省2020年的冬小麦的播种面积区域进行掩膜处理，在MODIS影像上对掩膜区域进行NDVI差值计算，对冬小麦进行高频度的长势监测，并对2020年河南省2、3月份以及4月上中旬的冬小麦长势结合气候因素进行监分析，得到结论如下：

（1）GEE提供的高性能计算能力以及Landsat 8、MODIS影像数据，可以为冬小麦的种植面积提取以及长势的动态监测提供数据保障和计算能力保障。本文所使用的方法能够对河南省2017—2020年的冬小麦播种面积进行有效提取，基于混淆矩阵计算的平均总体分类精度为95.2%，平均kappa系数为0.909，平均冬小麦分类精度为95.3%，在省域范围内与河南省统计的数据相比相对误差均低于3%。

（2）本文结合Landsat8以及MODIS 2种影像各自的优势，能够很好地对河南省2020年的冬小麦进行高频度的长势监测。受到气候等因素的影响，河南省绝大部分区域2020年冬小麦在2月份返青初期长势较2019年及往年好，到了3月上旬长势放缓，从3月下旬开始，冬小麦各类长势变化波动不大，趋于稳定状态，绝大部分种植区域冬小麦长势与往年持平。总体上，河南省2020年冬小麦的播种面积要高于往年，长势较往年及2019年要好。

[1] 郑文倩. 我国小麦价格形成机制及波动特征分析[D]. 北京: 中国农业科学院, 2016.

ZHENG W Q. Study on wheat price formation mechanism and price volatility characteristics in China[D]. Beijing: Chinese Academy of Agricultural Sciences, 2016. (in Chinese)

[2] 盛磊, 何亚娟, 吴全, 王飞. 河南省冬小麦产量遥感监测精度比较研究. 中国农业信息, 2018, 30(2): 95-102.

SHENG L, HE Y J, WU Q, WANG F. Comparative study on accuracy of winter wheat production by remote sensing monitoring in Henan

province. China Agricultural Information, 2018, 30(2): 95-102. (in Chinese)

[3] HAN J C, ZHANG Z, CAO J, LUO Y C, ZHANG L L, LI Z Y, ZHANG J. Prediction of winter wheat yield based on multi-source data and machine learning in China. Remote Sensing, 2020, 12(2): 236.

[4] 刘新杰, 魏云霞, 焦全军, 孙奇, 刘良云. 基于时序定量遥感的冬小麦长势监测与估产研究. 遥感技术与应用, 2019, 34(4): 756-765.

LIU X J, WEI Y X, JIAO Q J, SUN Q, LIU L Y. Growth monitoring and yield prediction of winter wheat based on time-series quantitative remote sensing data. Remote Sensing Technology and Application, 2019, 34(4): 756-765. (in Chinese)

[5] 杨珺雯, 张锦水, 潘耀忠, 孙佩军, 朱爽. 基于遥感识别误差校正面积的农作物种植面积抽样高效分层指标研究——以冬小麦为例. 中国农业科学, 2018, 51(4): 675-687.

YANG J W, ZHANG J S, PAN Y Z, SUN P J, ZHU S. An efficient hierarchical indicator based on the correction area of remote sensing identification error for planting acreage sampling-A case study of winter wheat. Scientia Agricultura Sinica, 2018, 51(4): 675-687. (in Chinese)

[6] FRANCH B, VERMOTE E F, SKAKUN S, ROGER J C, BECKER- RESHEF I, MURPHY E, JUSTICE C. Remote sensing based yield monitoring: Application to winter wheat in United States and Ukraine. International Journal of Applied Earth Observation and Geoinformation, 2019, 76: 112-127.

[7] HAO Z, ZHAO H L, ZHANG C, WANG H, JIANG Y, YI Z. Estimating winter wheat area based on an SVM and the variable fuzzy set method. Remote Sensing Letters, 2019, 10(4): 343-352.

[8] ZHANG X W, LIU J F, QIN Z Y, QIN F. Winter wheat identification by integrating spectral and temporal information derived from multi- resolution remote sensing data. Journal of Integrative Agriculture, 2019, 18(11): 2628-2643.

[9] MENG S Y, ZHONG Y F, LUO C, HU X, WANG X Y, HUANG S X. Optimal temporal window selection for winter wheat and wapeseed mapping with sentinel-2 images: A case study of Zhongxiang in China. Remote Sensing, 2020, 12(2): 226.

[10] 谭昌伟, 罗明, 杨昕, 马昌, 严翔, 周健, 杜颖, 王雅楠. 用PLS算法由HJ-1A/1B遥感影像估测区域冬小麦理论产量. 中国农业科学, 2015, 48(20): 4033-4041.

TAN C W, LUO M, YANG X, MA C, YAN X, ZHOU J, DU Y, WANG Y N. Remote sensing estimation of winter wheat theoretical yield on regional scale using partial least squares regression algorithm based on HJ-1A/1B images. Scientia Agricultura Sinica, 2015, 48(20): 4033-4041. (in Chinese)

[11] DONG C, ZHAO G X, QIN Y W, WAN H. Area extraction and spatiotemporal characteristics of winter wheat–summer maize in Shandong province using NDVI time series. PloS one, 2019, 14(12): 226508.

[12] 王利民, 刘佳, 姚保民, 季富华, 杨福刚. 基于GF-1影像NDVI年度间相关分析的冬小麦面积变化监测. 农业工程学报, 2018, 34(8): 184-191.

WANG L M, LIU J, YAO B M, JI F H, YANG F G. Area change monitoring of winter wheat based on relationship analysis of GF-1 NDVI among different years. Transactions of the Chinese Society of Agricultural Engineering, 2018, 34(8): 184-191. (in Chinese)

[13] 申健, 常庆瑞, 李粉玲, 王力. 基于时序NDVI的关中地区冬小麦种植信息遥感提取. 农业机械学报, 2017, 48(3): 215-220, 260.

SHEN J, CHANG Q R, LI F L, WANG L. Extraction of winter wheat information based on time-series NDVI in Guanzhong area. Transactions of the Chinese Society for Agricultural Machinery, 2017, 48(3): 215-220, 260. (in Chinese)

[14] 安塞, 沈彦俊, 赵彦茜, 郭英, 郭硕, 肖捷颖. 基于NDVI时间序列数据的冬小麦种植面积提取. 江苏农业科学, 2019, 47(15): 236-240.

AN S, SHEN Y J, ZHAO Y Q, GUO Y, GUO S, XIAO J Y. Extraction of winter wheat planting area based on NDVI time series data. Jiangsu Agricultural Science, 2019, 47(15): 236-240. (in Chinese)

[15] SHAN J, WANG Z M, SUN L, QIU L, YU K, WNAG J G. Study on extraction methods of winter wheat area based on GF-1 satellite images//2019 8th International Conference on Agro-Geoinformatics (Agro-Geoinformatics). IEEE, 2019: 1-5.

[16] 罗桓, 李卫国, 景元书, 徐向华, 陈华. 基于SVM的县域冬小麦种植面积遥感提取. 麦类作物学报, 2019, 39(4): 81-88.

LUO H, LI W G, JING Y S, XU X H, CHEN H. Remote sensing extraction of winter wheat planting area based on SVM. Journal of Triticeae Crops, 2019, 39(4): 81-88. (in Chinese)

[17] 李旭青, 刘世盟, 李龙, 金永涛, 范文磊, 吴怜. 基于RF算法优选多时相特征的冬小麦空间分布自动解译. 农业机械学报, 2019, 50(6): 218-225.

LI X Q, LIU S M, LI L, JIN Y T, FAN W L, WU L. Automatic interpretation of spatial distribution of winter wheat based on random forest algorithm to optimize multi-temporal features. Transactions of the Chinese Society for Agricultural Machinery, 2019, 50(6): 218-225. (in Chinese)

[18] 贺原惠子, 王长林, 贾慧聪, 陈方. 基于随机森林算法的冬小麦提取研究. 遥感技术与应用, 2018, 33(6): 1132-1140.

HE Y H Z, WANG C L, JIA H C, CHEN F. Research on extraction of winter wheat based on random forest. Remote Sensing Technology and Application, 2018, 33(6): 1132-1140. (in Chinese)

[19] LIU J T, FENG Q L, GONG J H, ZHOU J P, LIANG J M, LI Y. Winter wheat mapping using a random forest classifier combined with multi-temporal and multi-sensor data. International Journal of Digital Earth, 2018, 11(8): 783-802.

[20] 周志华. 机器学习. 北京: 清华大学出版社, 2016: 179-180.

ZHOU Z H. Machine Learning. Beijing: Tsinghua University Press, 2016: 179-180. (in Chinese)

[21] YOU J, PEI Z Y, WANG F, WU Q, GUO L. Area extraction of winter wheat at county scale based on modified multivariate texture and G satellite images. Transactions of the Chinese Society of Agricultural Engineering, 2016, 32(13): 131-139.

[22] 杨蕙宇, 王征强, 白建军, 韩红珠. 基于多特征提取与优选的冬小麦面积提取. 陕西师范大学学报(自然科学版), 2020, 48(1): 40-49.

YANG H Y, WANG Z Q, BAI J J, HAN H Z. Winter wheat area extraction based on multi-feature extraction and feature selection. Journal of Shaanxi Normal University(Natural Science Edition), 2020, 48(1): 40-49. (in Chinese)

[23] 何昭欣, 张淼, 吴炳方, 邢强. Google Earth Engine支持下的江苏省夏收作物遥感提取. 地球信息科学学报, 2019, 21(5): 752-766.

HE Z X, ZHANG M, WU B F, XING Q. Extraction of summer crop in Jiangsu based on Google Earth Engine. Journal of Geo-information Science, 2019, 21(5): 752-766. (in Chinese)

[24] GORELICK N, HANCHER M, DIXON M, ILYUSHCHENKO S, THAU D, MOORE R. Google Earth Engine: Planetary-scale geospatial analysis for everyone. Remote Sensing of Environment, 2017, 202: 18-27.

[25] 朱德海, 刘逸铭, 冯权泷, 欧聪, 郭浩, 刘建涛. 基于GEE的山东省近30年农业大棚时空动态变化研究. 农业机械学报, 2020, 51(1): 168-175.

ZHU D H, LIU Y M, FENG Q L, OU C, GUO H, LIU J T. Spatial- temporal dynamic changes of agricultural greenhouses in Shandong province in recent 30 years based on Google Earth Engine. Transactions of the Chinese Society for Agricultural Machinery, 2020, 51(1): 168-175. (in Chinese)

[26] LIU X P, HU G H, CHEN Y M, LI X, XU X C, LI S Y, PEI F S, WANG S J. High-resolution multi-temporal mapping of global urban land using Landsat images based on the Google Earth Engine Platform. Remote sensing of environment, 2018, 209: 227-239.

[27] 郝斌飞, 韩旭军, 马明国, 刘一韬, 李世卫. Google Earth Engine在地球科学与环境科学中的应用研究进展. 遥感技术与应用, 201 8, 33(4): 600-611.

HAO B F, HAN X J, MA M G, LIU Y T, LI S W. Research progress on the application of Google Earth Engine in geoscience and environmental sciences. Remote Sensing Technology and Application, 2018, 33(4): 600-611. (in Chinese)

[28] 孙丽, 王蔚丹, 陈媛媛, 董沫. 2019年美国冬小麦长势遥感监测分析. 安徽农业科学, 2020, 48(1): 241-244.

SUN L, WANG W D, CHEN Y Y, DONG M. Analysis of winter wheat growth of united states with remote sensing data in 2019. Journal of Anhui Agricultural Sciences, 2020, 48(1): 241-244. (in Chinese)

[29] 黄青, 李丹丹, 陈仲新, 刘佳, 王利民. 基于MODIS数据的冬小麦种植面积快速提取与长势监测. 农业机械学报, 2012,43(7): 163-167.

HUNAG Q, LI D D, CHEN Z X, LIU J, WANG L M. Monitoring of planting area and growth condition of winter wheat in China based on MODIS data. Transactions of the Chinese Society for Agricultural Machinery, 2012, 43(7): 163-167. (in Chinese)

[30] 何亚萍. 河南省循环农业发展研究[D]. 北京: 中国农业科学院, 2018.

HE Y P. Study on the development of circular agriculture in Henan province[D]. Beijing: Chinese Academy of Agricultural Sciences, 2018. (in Chinese)

[31] FARR T G, ROSEN P A, CARO E. The shuttle radar topography mission. Reviews of geophysics, 2007, 45(2):1-13.

[32] ULABY F T, KOUYATE F, BRISCO B, WILLIAMS T L. Textural Infornation in SAR Images. Geoscience and Remote Sensing, IEEE Transactions on GE-24, 1986, 24(2): 235-245.

[33] 王九中, 田海峰, 邬明权, 王力, 王长耀. 河南省冬小麦快速遥感制图. 地球信息科学学报, 2017, 19(6): 846-853.

WANG J Z, TIAN H F, WU M Q, WANG L, WANG C Y. Rapid mapping of winter wheat in Henan province. Journal of Geo- information Science, 2017, 19(6): 846-853. (in Chinese)

[34] BREIMAN L. Random forests. Machine Learning, 2001, 45(1): 5-32.

[35] 王利民, 刘佳, 唐鹏钦, 姚保民, 刘荣高. 农作物长势遥感监测需求、系统框架及业务应用. 中国农业信息, 2019, 31(2): 1-10.

WANG L M, LIU J, TANG P Q, YAO B M, LIU R G. Demands, system framework, and operational application of crop growth remote sensing monitoring. China Agricultural Information, 2019, 31(2): 1-10. (in Chinese)

[36] 千怀遂. 农作物遥感估产最佳时相的选择研究——以中国主要粮食作物为例. 生态学报, 1998,18(1):48-55 .

Qian H S. Selection of the optimum temporal for crop estimation using remote sensing data-Main food crops in China. Acta Ecologica Sinica, 1998,18(1):48-55. (in Chinese)

[37] 李方杰, 任建强, 吴尚蓉, 陈仲新, 张宁丹. 河南省冬小麦种植频率时空变化及影响因素分析. 中国农业科学, 2020, 53(9): 1773-1794.

LI F J, REN J Q, WU S R, CHEN Z X, ZHANG N D. Spatio-temporal variations of winter wheat planting frequency and their analysis of influencing factors in Henan province. Scientia Agricultura Sinica, 2020, 53(9): 1773-1794. (in Chinese)

[38] 邓荣鑫, 王文娟, 魏义长, 张富, 李春静, 刘文玉. 河南省冬小麦种植面积遥感监测及其时空特征研究. 灌溉排水学报, 2019, 38(9): 49-54.

DENG R X, WANG W Q, WEI Y C, ZHANG F, LI C J, LIU W Y. Remote estimation of winter wheat area and its spatio-temporal characteristics in Henan province. Journal of Irrigation and Drainage, 2019, 38(9): 49-54. (in Chinese)

Area extraction and growth Monitoring of Winter Wheat in Henan Province supported by Google Earth Engine

1School of Computer and Information Engineering, Henan University, Kaifeng 475004, Henan;2Department of Laboratory and Equipment Management, Henan University, Kaifeng 475004, Henan;3Henan Key Laboratory of Big Data Analysis and Processing/ Henan University, Kaifeng 475004, Henan

【】The aim of this study was to use remote sensing technology to extract the spatial distribution information of winter wheat in Henan province from 2017 to 2020, and then to monitor the growth of winter wheat in 2020 with high frequency and to analyze the meteorological conditions. 【】Based on the cloud platform of Google Earth engine (GEE), the selected Landsat 8 image data were synthesized according to the maximum value of NDVI, and then the features were constructed to add terrain features, texture features, NDVI and a new feature NDVI amplification. Random forest classification method was used to train the sample data according to the constructed features to extract the winter wheat planting area in Henan province from 2017 to 2020. The accuracy of the extracted winter wheat sown area was verified by confusion matrix and Henan statistical yearbook data. After accuracy verification, a mask was generated for the extracted winter wheat planting area in Henan province in 2020. In the mask area (winter wheat planting area) combined with MODIS high time resolution image data, the NDVI synchronization difference method was used to monitor the winter wheat growth from February to April in 2020. 【】The GEE cloud platform could be used to quickly map the spatial distribution information of winter wheat planting areas in Henan province. Using random forest method to add terrain feature, texture feature, NDVI and new feature NDVI could effectively improve the extraction accuracy of winter wheat and reduce the relative error with statistical data. Based on confusion matrix, the average overall classification accuracy was 95.2%, the average kappa coefficient was 0.909, and the average classification accuracy of winter wheat was 95.3%. Compared with the statistical yearbook data of Henan province, the relative errors of winter wheat sown area extracted by this method in Henan province from 2017 to 2019 were all less than 3%. The average relative error of winter wheat sown area in the main planting areas of winter wheat in Henan province was less than 6%. MODIS image data combined with NDVI difference model could be used to monitor the growth of winter wheat in Henan province in 2020. The growth of winter wheat in Henan province was better than that of previous years and 2019 during the return to green period. In the later growth stage of winter wheat, the growth of most areas was the same as that of previous years and 2019. On the whole, the growth of winter wheat in 2020 was better than that of previous years and 2019. 【】The method proposed in this paper could carry out high-precision extraction and high-frequency growth monitoring of winter wheat in Henan province, and could provide a scientific basis for local governments or some agricultural departments in arranging and guiding agricultural activities.

winter wheat; growth monitoring; Google Earth Engine; random forests; NDVI; Landsat; MODIS