[发明专利]基于D-fG

权利要求书

1.一种基于D-f_G参数遥感获取的冬小麦动态收获指数遥感估算方法，其特征在于，包括以下步骤：

A1、根据地面实测动态生物量数据，构建作物开花期-成熟期期间不同时期累积的地上生物量与对应时期地上生物量间比值动态参数D-f_G；

D-f_G计算方法如下：

式中，∑W_post为冬小麦开花期-成熟期期间不同时期累积的地上生物量(kg/hm²)；∑W_whole为采样时期对应的全部地上生物量(kg/hm²)；t为采样时间，W_t为t采样时间的干物质量的重量(kg/hm²)，W_a为开花期干物质量的重量(kg/hm²)，D-f_G,t表示t采样时间的比值参数；

A2、基于地面作物冠层高光谱数据构建的任意两个冠层高光谱窄波段光谱指数NDSI，建立NDSI与冬小麦D-f_G之间的线性模型；

A3、绘制并分析NDSI与冬小麦D-f_G间的拟合精度R²二维图；

A4、通过确定R²极大值区域和极大值区域重心，从而得到对冬小麦D-f_G敏感的波段中心；

A5、确定D-f_G估算最优波段组合；

A6、基于NDSI和D-f_G关系的D-f_G遥感估算模型；

A7、D-f_G的遥感估算；

A8、获得基于D-f_G和动态收获指数D-HI关系的动态收获指数估算模型；

A9、D-HI的遥感估算。

2.根据权利要求1所述的方法，其特征在于，所述步骤A4，根据R²极大值区域重心法获得敏感波段中心，从而确定D-f_G估算的敏感波段，即在冠层高光谱每个波段对应的窄波段光谱指数NDSI和D-f_G间相关性计算基础上，根据相关系数满足统计显著性要求的阈值确定极大值区域，在此基础上，计算相关系数极大值区域的重心，从而获得NDSI与D-f_G参数相关性较大的光谱波段中心和波段组合；具体过程如下：

首先，在绘制NDSI与冬小麦D-f_G间拟合R²二维图的基础上，确定NDSI与冬小麦D-f_G间相关性高的波段区域；其次，在该区域内寻找R²极大值点，并遍历该点8邻域内满足显著性条件的所有点，并将这些点的集合标记为R²极大值区域Ω；最后，计算R²极大值区域的重心，将其作为每个R²极大值区域的敏感波段中心；重心的计算公式(7)如下：

式中，f(u,v)为波段坐标为(u,v)的R²值，Ω为极大值区域，敏感波段中心坐标。

3.根据权利要求1所述的方法，其特征在于，所述步骤A4，对冬小麦D-f_G敏感的6个敏感波段中心包括λ(443nm，506nm)、λ(442nm，635nm)、λ(732nm，834nm)、λ(787nm，804nm)、λ(810nm，877nm)和λ(861nm，985nm)。

4.根据权利要求1所述的方法，其特征在于，所述步骤A6，基于NDSI和D-f_G关系的D-f_G遥感估算模型为：

D-f_G,t＝m×NDSI_i，j，t+n (5)

其中，i、j为分别为350nm～1000nm间高光谱波段，t为不同取样时间，m和n为拟合后得到线性方程中的拟合参数；根据该公式计算求出作物开花期至t时期累积地上生物量与t生育期地上生物量间比值参数D-f_G,t。