基于生长盛期NMDI增减比值指数的玉米自动制图方法与流程

技术特征：

1.一种基于生长盛期NMDI增减比值指数的玉米自动制图方法，其特征在于，包括以下步骤：

步骤S01：逐像元建立研究区的植被指数时序数据集和NMDI时序数据集；

步骤S02：逐像元计算农作物在每个生长周期内的EVI最大值，获取农作物生长峰值时间；

步骤S03：依据所述农作物生长峰值时间，确定农作物生长盛期前段区间和农作物生长盛期后段区间；

步骤S04：分别建立所述农作物生长盛期前段区间的NMDI增量指标和所述农作物生长盛期后段区间的NMDI增量指标；

步骤S05：分别建立所述农作物生长盛期前段区间的NMDI减量指标和所述农作物生长盛期后段区间的NMDI减量指标；

步骤S06：基于所述NMDI增量指标和NMDI减量指标，建立NMDI增减比值指数；

步骤S07：依据所述NMDI增减比值指数进行研究区玉米制图。

2.根据权利要求1所述的基于生长盛期NMDI增减比值指数的玉米自动制图方法，其特征在于：所述步骤S03中，农作物生长盛期前段区间为农作物生长峰值时间的前20天，农作物生长盛期后段区间为农作物生长峰值时间的后20天。

3.根据权利要求1所述的基于生长盛期NMDI增减比值指数的玉米自动制图方法，其特征在于：所述步骤S04中，在农作物生长盛期前段区间，逐日计算NMDI的斜率，将所有数值为正的斜率累加，设为农作物生长盛期前段区间的NMDI增量指标，记为E1；在农作物生长盛期后段区间，逐日计算NMDI的斜率，将所有数值为正的斜率累加，设为农作物生长盛期后段区间的NMDI增量指标，记为E2。

4.根据权利要求1所述的基于生长盛期NMDI增减比值指数的玉米自动制图方法，其特征在于：所述步骤S05中，在农作物生长盛期前段区间，逐日计算NMDI的斜率，将所有数值为负的斜率累加后求绝对值，设为农作物生长盛期前段区间的NMDI减量指标，记为D1；在农作物生长盛期后段区间，逐日计算NMDI的斜率，将所有数值为负的斜率累加后求绝对值，设为农作物生长盛期后段区间的NMDI减量指标，记为D2。

5.根据权利要求1所述的基于生长盛期NMDI增减比值指数的玉米自动制图方法，其特征在于：所述步骤S6中，NMDI增减比值指数的建立方法如下：

其中，P为NMDI增减比值指数，E1为农作物生长盛期前段区间的NMDI增量指标，E2为农作物生长盛期后段区间的NMDI增量指标，D1为农作物生长盛期前段区间的NMDI减量指标，D2为农作物生长盛期后段区间的NMDI减量指标；θ为常数，取值范围为[0.001,0.004]。

6.根据权利要求1所述的基于生长盛期NMDI增减比值指数的玉米自动制图方法，其特征在于：所述步骤S7中，将所述NMDI增减比值指数与阈值ω进行比较，若NMDI增减比值指数大于阈值ω，则判断该农作物为玉米，否则判断该农作物不是玉米；然后制作研究区的玉米空间分布图。

7.根据权利要求6所述的基于生长盛期NMDI增减比值指数的玉米自动制图方法，其特征在于：所述阈值ω的取值为0.05。