基于植被指数时序谱特征的大豆种植区提取方法和系统

技术特征：
1.一种基于植被指数时序谱特征的大豆种植区提取方法，其特征在于，所述方法包括：步骤s1、遥感数据预处理：对sentinel-2卫星影像进行预处理，并进行耕地掩膜；步骤s2、特征提取：提取大豆整个生育期内掩膜后的sentinel-2影像增强型植被指数时序谱特征，简称时序谱特征；提取大豆关键生育期内掩膜后的sentinel-2影像波段反射率特征，简称光谱特征；提取海拔及坡度特征，简称地形特征，并利用耕地数据对所述时序谱特征、光谱特征和地形特征分别进行掩膜；步骤s3、分类模型的构建与训练：构建分类模型并通过训练设置参数；步骤s4、以所述时序谱特征为主，将所述时序谱特征、光谱特征和地形特征进行组合得到待选特征组合，得到待选特征组合；将所述待选特征组合输入分类模型，计算待选特征组合的分类精度，根据分类精度得到最优特征组合；步骤s5、大豆种植区提取及精度评价：将最优特征组合输入随机森林分类模型提取大豆种植区分布图并进行精度评价。2.根据权利要求1所述的一种基于植被指数时序谱特征的大豆种植区提取方法，其特征在于，在所述步骤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
]形成的角度。3.根据权利要求1所述的一种基于植被指数时序谱特征的大豆种植区提取方法，其特征在于，在所述步骤s2中，所述提取大豆关键生育期内掩膜后的sentinel-2影像波段反射率特征的方法包括：利用百分位数、最大值、最小值、平均值和标准差方法提取大豆关键生育期内掩膜后的sentinel-2影像波段反射率特征。4.根据权利要求1所述的一种基于植被指数时序谱特征的大豆种植区提取方法，其特征在于，在所述步骤s2中，大豆关键生育期内掩膜后的sentinel-2影像波段包含，红边1、红边2、红边3、短波红外1和短波红外2；所述关键生育期为1-2月，或者为8-9月，具体包含大豆开花期、结荚期和灌浆期三个生长阶段。5.根据权利要求1所述的一种基于植被指数时序谱特征的大豆种植区提取方法，其特征在于，在所述步骤s2中，所述提取海拔及坡度特征的方法包括：利用航天飞机雷达地形测绘任务数字高程数据提取海拔及坡度特征。6.根据权利要求1所述的一种基于植被指数时序谱特征的大豆种植区提取方法，其特
征在于，在所述步骤s4中，所述待选特征组合的具体组合方式包括：时序谱特征、光谱特征、地形特征、时序谱特征及光谱特征、时序谱特征及地形特征、时序谱特征和光谱特征以及地形特征，共6组待选特征组合。7.根据权利要求1所述的一种基于植被指数时序谱特征的大豆种植区提取方法，其特征在于，在所述步骤s4中，计算待选特征组合的分类精度的方法包括：采用基于混淆矩阵的随机森林分类精度评价方法计算待选特征组合的分类精度。8.一种用于基于植被指数时序谱特征的大豆种植区提取系统，其特征在于，所述系统包括：第一处理模块，被配置为，遥感数据预处理：对sentinel-2卫星影像进行预处理，并进行耕地掩膜；第二处理模块，被配置为，特征提取：提取大豆整个生育期内掩膜后的sentinel-2影像增强型植被指数时序谱特征，简称时序谱特征；提取大豆关键生育期内掩膜后的sentinel-2影像波段反射率特征，简称光谱特征；提取海拔及坡度特征，简称地形特征，并利用耕地数据对所述时序谱特征、光谱特征和地形特征分别进行掩膜；第三处理模块，被配置为，分类模型的构建与训练：构建分类模型并通过训练设置参数；第四处理模块，被配置为，以所述时序谱特征为主，将所述时序谱特征、光谱特征和地形特征进行组合得到待选特征组合，得到待选特征组合；将所述待选特征组合输入分类模型，计算待选特征组合的分类精度，根据分类精度得到最优特征组合；第五处理模块，被配置为，大豆种植区提取及精度评价：将最优特征组合输入随机森林分类模型提取大豆种植区分布图并进行精度评价。9.一种电子设备，其特征在于，所述电子设备包括存储器和处理器，所述存储器存储有计算机程序，所述处理器执行所述计算机程序时，实现权利要求1至7中任一项所述的一种基于植被指数时序谱特征的大豆种植区提取方法中的步骤。10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质上存储有计算机程序，所述计算机程序被处理器执行时，实现权利要求1至7中任一项所述的一种基于植被指数时序谱特征的大豆种植区提取方法中的步骤。

技术总结
本发明提出基于植被指数时序谱特征的大豆种植区提取方法和系统。方法包括：通过构建中高时空分辨率的时间序列植被指数，利用线性谐波模型提取研究区植被指数时序谱特征。其次，基于关键生育期遥感原始影像波段反射率，通过百分位数、最大值、最小值、平均值、标准差方法提取光谱特征。另外，利用航天飞机雷达地形测绘任务数字高程数据提取地形特征。借助随机森林分类模型，分析整个生育期时序谱特征、关键生育期光谱特征及地形特征对大豆识别精度的差异，探究不同类别特征的最优组合，实现大豆种植区高精度提取。本申请提出的方案提供了大区域尺度大豆种植区快速准确的识别方法，为大豆面积、产量、病虫害监测等提供科学依据。病虫害监测等提供科学依据。病虫害监测等提供科学依据。


技术研发人员：彭代亮 罗旺 刘锦绣 陈月 楼子杭
受保护的技术使用者：中国科学院空天信息创新研究院
技术研发日：2022.05.07
技术公布日：2022/8/9