[发明专利]基于生长前后期氮素输出累积指数的土豆制图方法

权利要求书

1.一种基于生长前后期氮素输出累积指数的土豆制图方法，其特征在于，包括如下步骤：

步骤S01：建立研究区植被指数时序数据集和氮素指数时序数据集；

步骤S02：建立氮素指数差分时序数据集；

步骤S03：建立氮素指数变化量时序数据集；

步骤S04：逐像元动态获取农作物物候期；

步骤S05：获取农作物生长前期氮素输出率指标；

步骤S06：获取农作物生长后期氮素累积率指标；

步骤S07：建立生长前后期氮素输出累积指数；

步骤S08：构建土豆制图技术流程方法；

步骤S09：获得研究区土豆空间分布图；

所述步骤S02中，针对研究区耕地区域，在研究年份内，逐日对氮素指数时序数据进行时序差分，建立氮素指数差分时序数据集；

所述步骤S03中，对氮素指数差分时序数据集中数据实施取绝对值运算，建立研究区氮素指数变化量时序数据集；

所述步骤S04中，在研究区耕地区域，逐像元动态获取包括农作物生长盛期、生长开始期和生长结束期的农作物物候期，并针对每个农作物生长期，求算植被指数时序曲线最大值，将该生长期内植被指数最大值所处时刻确定为农作物生长盛期；

所述步骤S05中，在农作物生长前期，将氮素指数差分时序数据与氮素指数变化量时序数据实施相减后取绝对值操作，获得农作物生长前期氮素输出量；在农作物生长前期，将氮素指数变化量时序数据逐日累加，获得农作物生长前期氮素变化量；将农作物生长前期氮素输出量除以农作物生长前期氮素变化量，获得农作物生长前期氮素输出率指标NORE，其计算公式为：

；

式中，NRI^t+1、NRI^t分别表示第t+1天、t天氮素指数数值；Start、Heading分别对应农作物生长开始期和生长盛期；

所述步骤S06中，在农作物生长后期，将氮素指数差分时序数据与氮素指数变化量时序数据实施相加后取绝对值操作，获得农作物生长后期氮素累积量；将农作物生长后期氮素累积量除以农作物生长后期氮素变化量，获得农作物生长后期氮素累积率指标NARL，其计算公式为：

；

式中，NRI^t+1、NRI^t分别表示第t+1天、t天氮素指数数值；Heading、End分别对应农作物生长盛期和生长结束期；

所述步骤S07中，将农作物生长前期氮素输出率指标乘以农作物生长后期氮素累积率指标，获得该农作物生长期内的生长前后期氮素输出累积指数OAI，其计算公式为：

；

式中，NORE、NARL分别表示农作物生长前期氮素输出率、农作物生长后期氮素累积率；

所述步骤S08中，利用所构建的生长前后期氮素输出累积指数OAI，建立土豆制图判别规则：如果满足OAIθ，则该像元为土豆，否则为其他农作物；阈值θ在不同区域实际应用中可作调整。