[发明专利]一种农田尺度作物水分利用效率遥感评价方法

权利要求书

1.一种农田尺度作物水分利用效率遥感评价方法，其特征在于，包括以下步骤：

S1、获取研究区作物全生育期蒸散量ET和全生育期总初级生产力GPP；

S2、根据ET和GPP计算得到作物水分生产率WP；

S3、对WP与ET的关系进行分析，得到作物水分生产率最大时的作物全生育期耗水量ET_min，对GPP与ET的关系进行分析，得到作物总初级生产力最大时的作物全生育期耗水量ET_max，确定ET_min～ET_max区间为作物的适宜耗水区间；

S4、根据ET与ET_min以及ET_max的关系对作物分布区进行区间划分，得到研究区作物水分利用效率评价分类结果；

所述步骤S3包括以下分步骤：

S31、对作物的WP值与ET值进行散点图分析，得到作物耗水量-水分生产率模型：

WP＝a₁×ET²+b₁×ET+c₁

其中a₁,b₁,c₁均为回归系数；

S32、根据作物耗水量-水分生产率模型，得到作物水分生产率最大时的作物全生育期耗水量ET_min：

S33、对作物的GPP值与ET值进行散点图分析，得到作物耗水量-作物总初级生产力模型：

GPP＝a₂×ET²+b₂×ET+c₂

其中a₂,b₂,c₂均为回归系数；

S34、根据作物耗水量-作物总初级生产力模型，得到作物总初级生产力最大时的作物全生育期耗水量ET_max：

S35、确定ET_min～ET_max区间为作物的适宜耗水区间；

所述步骤S4包括以下分步骤：

S41、若ET≤ET_min，则该作物分布区属于水量不足区；若ET_minET≤ET_max，则进入步骤S42；若ETET_max，则进入步骤S43；

S42、将ET作为自变量输入作物耗水量-水分生产率模型，得到作物目标水分生产率WP_max，将实际作物水分生产率WP与作物目标水分生产率WP_max进行对比，若WPWP_max，则该作物分布区属于水量适宜区，若WP≤WP_max，则该作物分布区属于灌溉不合理区；

S43、计算作物生育期内的总水分盈亏量ID并进入步骤S44：

ID＝∑(ET_i-P_ei)

其中ET_i表示作物某生育阶段的实际耗水量，P_ei表示作物某生育阶段的有效降雨量；

S44、若ID0，则该作物分布区属于过量灌溉区，若ID≤0，则该作物分布区属于水量适宜区。

2.根据权利要求1所述的农田尺度作物水分利用效率遥感评价方法，其特征在于，所述步骤S1包括以下分步骤：

S11、获取研究区作物空间分布栅格数据，确定作物生育期；

S12、下载作物生育期时段内的MODIS MOD16A2HGF蒸散量ET产品和MOD17A2HGF总初级生产力GPP产品，并根据产品导出的质量标识层，在标识有云或质量不佳的区域范围内选择3年内相同时段的同类产品进行最大值合成法处理，得到处理后的作物生育期内的ET和GPP空间分布数据；

S13、对作物生育期内的ET进行累加，并利用作物空间分布栅格数据对其进行掩模处理，得到研究区作物分布区的ET空间分布图；

S14、对作物生育期内的GPP进行累加，并利用作物空间分布栅格数据对其进行掩模处理，得到研究区作物分布区的GPP空间分布图。