基于历史植被指数时序谱曲线与产量映射模式的大豆产量预测方法技术

一种基于历史植被指数时序谱曲线与产量映射模式的大豆产量预测方法，包括以下步骤：1)历史植被指数时序谱曲线与历史产量映射模式的构建；2)基于植被指数时序谱曲线的大豆生育期划分；3)基于植被指数时序谱曲线匹配的大豆产量预测；4)基于历史植被指数时序谱曲线与产量映射模式进行时序推进的大豆产量预测；5)大豆产量预测精度验证。该方法从大豆的植被指数时序谱入手获取不同生育期的信息，依据预测当年时序谱曲线与历史时序谱曲线的相似性关系，构建基于历史植被指数时序谱曲线与历史产量映射模式的估产方法，为农作物产量预测精度改善提供了一种新的思路。改善提供了一种新的思路。改善提供了一种新的思路。

【技术实现步骤摘要】
基于历史植被指数时序谱曲线与产量映射模式的大豆产量预测方法


[0001]本专利技术涉及一种大豆产量预测方法。

技术介绍

[0002]粮食安全问题一直是我国高度重视的问题。为了满足14亿人口的温饱需求，国家正在大力发展农业及其相关产业，力求稳定这一国计民生的根基。大豆是一种重要的粮油作物，提供了人类所需的优质蛋白来源。大豆是重要基础性、战略性物资，是最具经济效益的作物。大豆是我国进口量最大的农产品，在农产品贸易方面扮演了举足轻重的角色，涉及到政治和经济等关键领域。如何保证大豆的生产和供应，是我国农业方面的一个重要课题。
[0003]目前全球大豆的主产区有美国、巴西、阿根廷和中国，这些地区集中了全世界90％以上的大豆产量。近年来随着转基因大豆商业化生产的普及，大豆种植面积和单产水平均大幅提高，进而推动了大豆产量显著增加。如何利用遥感等信息技术监测庞大的大豆种植区域，准确掌握产量信息，对确定农业发展规划、制定农产品进出口计划及调整种植业结构等方面都有着重要的战略意义。

技术实现思路

[0004]为了充分利用现有获取的海量遥感数据以及大豆产量数据，克服现有遥感数据与产量关联不确定的技术问题，本专利技术提供一种基于历史植被指数时序谱曲线与产量映射模式的大豆产量预测方法，充分利用海量的历史遥感信息和大豆产量，结合遥感快速、大面积监测能力，能够提供较为准确地区域大豆产量预测。
[0005]本专利技术解决其技术问题所采用的技术方案是：
[0006]一种基于历史植被指数时序谱曲线与产量映射模式的大豆产量预测方法，包括以下步骤：
[0007]1)历史植被指数时序谱曲线与历史产量映射模式的构建，过程如下：
[0008]首先将中高分辨率时间序列反射率数据通过云量分数筛选出云量在10％以下的数据，利用大豆高种植区数据进行掩膜，掩膜后计算时间序列植被指数，以县级行政区为基本单位，计算行政单元植被指数的平均值，构建行政单元植被指数时序谱数据，并利用Beck模型拟合时序谱曲线，利用历史植被指数时序谱曲线与历史大豆产量数据之间关系，建立由历史植被指数时序谱曲线与历史产量映射模式；
[0009][0010][0011][0012]式中表示历史植被指数时序谱曲线，Y1，Y2，Y
M
表示历史大豆产量，
VI1,VI2……
VI
n
表示植被指数时序谱中第1到第n个植被指数；
[0013]2)基于植被指数时序谱曲线的大豆生育期划分，过程如下：
[0014]利用植被指数时序谱数据，将经过曲线一阶导数最大值点且以上升最大速率为斜率的直线定义为生长线，生长线与淡季基准线相交的点定义为出苗期，生长线与高值平台相交的点定义为开花期；此外，将时间序列曲线的最高点作为结荚期，一阶导数最小值点作为转黄期，实现大豆生育期的划分；
[0015]3)基于植被指数时序谱曲线匹配的大豆产量预测，过程如下：
[0016]采用夹角θ和欧氏距离d两种匹配方法，将待预测年的大豆时序谱曲线与历史时序谱数据集中的所有单一时序谱曲线进行一对一的匹配，匹配算法如下：
[0017][0018]其中θ是两个时序谱曲线数据在弧度制下的夹角，d是两个时序谱曲线数据的空间距离，为用于预测大豆产量的时序谱数据，为历史样本时序谱数据；
[0019]分别计算θ，d，和θ*d三个变量，以三参数最小原则确定与待预测年匹配性最高的历史时序谱曲线，并将历史时序谱曲线映射的大豆产量赋值为待预测年的产量；
[0020]4)基于历史植被指数时序谱曲线与产量映射模式进行时序推进的大豆产量预测，过程如下：
[0021]在进行产量预测时，采用一种时序推进的方法，即随着时序数据获取的进度，从出苗期开始，形成出苗期
‑
开花期，出苗期
‑
结荚期，出苗期
‑
转黄期等一系列产量预测结果，开展基于历史植被指数时序谱曲线与历史产量映射模式的随时序推进的大豆产量预测，实现以数据获取时间为序从出苗期开始一系列的产量预测。
[0022]进一步，所述方法还包括以下步骤：
[0023]5)大豆产量预测精度验证
[0024]将待预测年的大豆产量Y
估测
与实测产量Y
实测
进行精度验证，
[0025][0026]式中，PRE为产量预测精度，Y
估测
为预测年大豆产量，Y
实测
为实测大豆产量。
[0027]再进一步，所述步骤1)中，植被指数包括NDVI，EVI，SAVI，DVI，MSAVI。
[0028]所述步骤3)中，考虑到历史样本中大豆种植的时间与待预测年大豆种植时间可能存在时间差，进而引起的时序谱匹配误差，采取将待匹配年时序谱前后最大滑动15天开展匹配筛选出相似度最高的曲线作预测大豆产量。
[0029]本专利技术中，首先利用遥感技术获取的农作物反射率信息提取时间序列植被指数，建立起基于历史植被指数时序谱曲线与历史大豆产量之间的映射关系，并随时序推进计算预测年植被指数时序谱曲线与历史时序谱曲线的匹配性，设计出了一种基于历史植被指数时序谱曲线与产量映射模式的大豆产量预测方法，实现县级尺度的农作物精确产量预测。
[0030]本专利技术的有益效果为：利用遥感等信息技术监测庞大的大豆种植区域，较为准确地预测产量信息。
附图说明
[0031]图1是伊利诺伊州大豆样本生育期划分示意图。
[0032]图2是伊利诺伊州开花期到结荚期以欧氏距离(a)和夹角(b)及两者结合(c)匹配产量预测结果。
[0033]图3是美国伊利诺伊州大豆使用NDVI及欧氏距离匹配的产量预测精度随生育期变化图。
[0034]图4是基于历史植被指数时序谱曲线与产量映射模式的大豆产量预测方法流程图。
具体实施方式
[0035]下面对本专利技术作进一步描述。本实施例为针对美国伊利诺伊州2020年大豆产量进行预测。
[0036]参照图1～图4，一种基于历史植被指数时序谱曲线与产量映射模式的大豆产量预测方法，包括以下步骤：
[0037]1)历史植被指数时序谱曲线与历史产量映射模式的构建，过程如下：
[0038]首先将2010年
‑
2019年共计10年的MOD09GA数据集按美国伊利诺伊州大豆生长时间过滤，之后按伊利诺伊州县级行政区裁剪，通过影像集state_1km波段计算出研究区内的云量分数，筛选出云量在10％以下的数据；然后用美国CDLs提取的当年大豆高密度种植区进行掩膜，掩膜后计算得到时间序列植被指数数据，再在研究区内统计得到植被指数的平均值，建立县级植被指数时序谱数据，并利用Beck模型拟合时序谱曲线，构建历史植被指数时序谱曲线与历史产量映射模式。对于10年美国伊利诺伊州历史数据共计构建926对植被指数时序谱曲线与产量的映射；
[0039]该步骤中植被指数包括NDVI，EVI，SAVI，DVI，MS本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于历史植被指数时序谱曲线与产量映射模式的大豆产量预测方法，其特征在于，所述方法包括以下步骤：1)历史植被指数时序谱曲线与历史产量映射模式的构建，过程如下：首先将中高分辨率时间序列反射率数据通过云量分数筛选出云量在10％以下的数据，利用大豆高种植区数据进行掩膜，掩膜后计算时间序列植被指数，以县级行政区为基本单位，计算行政单元植被指数的平均值，构建行政单元植被指数时序谱数据，并利用Beck模型拟合时序谱曲线，利用历史植被指数时序谱曲线与历史大豆产量数据之间关系，建立由历史植被指数时序谱曲线与历史产量映射模式；史植被指数时序谱曲线与历史产量映射模式；史植被指数时序谱曲线与历史产量映射模式；式中表示历史植被指数时序谱曲线，Y1，Y2，Y
M
表示历史大豆产量，VI1,VI2……
VI
n
表示植被指数时序谱中第1到第n个植被指数；2)基于植被指数时序谱曲线的大豆生育期划分，过程如下：利用植被指数时序谱数据，将经过曲线一阶导数最大值点且以上升最大速率为斜率的直线定义为生长线，生长线与淡季基准线相交的点定义为出苗期，生长线与高值平台相交的点定义为开花期；此外，将时间序列曲线的最高点作为结荚期，一阶导数最小值点作为转黄期，实现大豆生育期的划分；3)基于植被指数时序谱曲线匹配的大豆产量预测，过程如下：采用夹角θ和欧氏距离d两种匹配方法，将待预测年的大豆时序谱曲线与历史时序谱数据集中的所有单一时序谱曲线进行一对一的匹配，匹配算法如下：其中θ是两个时序谱曲线数据在弧度制下的夹角，d是两个时序谱曲线数据的空间距离，为用于预测大豆产量的时序谱数据，为历史样本时序谱数据；分别计算θ，d，...

【专利技术属性】
技术研发人员：王福民，郑崛屹，姚晓萍，徐天玥，李嘉乐，陈思婷，
申请(专利权)人：浙江大学，
类型：发明
国别省市：