[发明专利]基于植被指数时序谱特征的大豆种植区提取方法和系统

权利要求书

1.一种基于植被指数时序谱特征的大豆种植区提取方法，其特征在于，所述方法包括：

步骤S1、遥感数据预处理：对Sentinel-2卫星影像进行预处理，并进行耕地掩膜；

步骤S2、特征提取：提取大豆整个生育期内掩膜后的Sentinel-2影像增强型植被指数时序谱特征，简称时序谱特征；提取大豆关键生育期内掩膜后的Sentinel-2影像波段反射率特征，简称光谱特征；提取海拔及坡度特征，简称地形特征，并利用耕地数据对所述时序谱特征、光谱特征和地形特征分别进行掩膜；

步骤S3、分类模型的构建与训练：构建分类模型并通过训练设置参数；

步骤S4、以所述时序谱特征为主，将所述时序谱特征、光谱特征和地形特征进行组合得到待选特征组合，得到待选特征组合；将所述待选特征组合输入分类模型，计算待选特征组合的分类精度，根据分类精度得到最优特征组合；

步骤S5、大豆种植区提取及精度评价：将最优特征组合输入随机森林分类模型提取大豆种植区分布图并进行精度评价；

在所述步骤S2中，所述提取大豆整个生育期内掩膜后的Sentinel-2影像增强型植被指数时序谱特征的方法包括：

利用线性谐波模型提取大豆整个生育期内掩膜后的Sentinel-2影像增强型植被指数时序谱特征；

所述线性谐波模型的公式为：

其中，f(t)是第t时刻的拟合的增强型植被指数的值，a_k是余弦系数，b_k是正弦系数，c是截距项，n是谐波级数的阶数，w是频率等于1.5；自变量t是一年中的某一天；n分别设置为1至5，每次增加1，依据大豆的时间序列中原始值与拟合值之间均方误差最小确定n的值，此处n取2；取相位和振幅作为时序特征输入，振幅定义为二维矢量[a_k,b_k]的长度，相位定义为二维矢量[a_k,b_k]形成的角度；

在所述步骤S2中，所述提取大豆关键生育期内掩膜后的Sentinel-2影像波段反射率特征的方法包括：

利用百分位数、最大值、最小值、平均值和标准差方法提取大豆关键生育期内掩膜后的Sentinel-2影像波段反射率特征。

2.根据权利要求1所述的一种基于植被指数时序谱特征的大豆种植区提取方法，其特征在于，在所述步骤S2中，大豆关键生育期内掩膜后的Sentinel-2影像波段包含，红边1、红边2、红边3、短波红外1和短波红外2；

所述关键生育期为1-2月，或者为8-9月，具体包含大豆开花期、结荚期和灌浆期三个生长阶段。

3.根据权利要求1所述的一种基于植被指数时序谱特征的大豆种植区提取方法，其特征在于，在所述步骤S2中，所述提取海拔及坡度特征的方法包括：

利用航天飞机雷达地形测绘任务数字高程数据提取海拔及坡度特征。

4.根据权利要求1所述的一种基于植被指数时序谱特征的大豆种植区提取方法，其特征在于，在所述步骤S4中，所述待选特征组合的具体组合方式包括：

时序谱特征、光谱特征、地形特征、时序谱特征及光谱特征、时序谱特征及地形特征、时序谱特征和光谱特征以及地形特征，共6组待选特征组合。

5.根据权利要求1所述的一种基于植被指数时序谱特征的大豆种植区提取方法，其特征在于，在所述步骤S4中，计算待选特征组合的分类精度的方法包括：

采用基于混淆矩阵的随机森林分类精度评价方法计算待选特征组合的分类精度。