一种改进的时序植被指数数据合成方法

摘要

本发明公开了一种改进的时序植被指数数据合成方法，包括：A、对地表反射率数据质量进行分级；B、计算各个像元的植被指数，包括NDVI和EVI；C、依据所述地表反射率数据质量评估结果，逐像元地统计清洁像元的个数N；D、依据所述清洁像元个数N的不同，选择不同的操作得到最终的合成值。采用本发明的方法，能够对现有CV‑MVC合成方法中存在的不足进行改进，提高其精细度和限制性。

主权项

一种改进的时序植被指数数据合成方法，其特征在于，包括如下步骤：A、对地表反射率数据质量进行分级；具体为：根据是否有云、大气气溶胶厚度、观测角度以及太阳高度，得到各个像元的质量分级Rank，依据分级的不同，分为Rank＝0～Rank＝7；B、计算各个像元的植被指数，包括NDVI和EVI；所述计算各个像元的植被指数，具体为：$NDVI=\frac{{\rho }_{nir}-{\rho }_{red}}{{\rho }_{nir}+{\rho }_{red}},$正常地区：或云、雪覆盖的地区：其中：G＝2.5，L＝1.0，C₁＝6，C₂＝‑7.5；所述ρ_nir为近红外波段反射率，ρ_red为红光波段反射率，ρ_blue为蓝光波段反射率；C、依据所述地表反射率数据质量评估结果，逐像元地统计清洁像元的个数N；所述的清洁像元为Rank＝0和Rank＝1的像元；D、依据所述清洁像元个数N的不同，选择不同的操作得到最终的合成值；步骤D所述依据清洁像元个数N的不同，选择不同的操作得到最终的合成值，具体为：当N≥2时，从当前点多天植被指数中选取清洁像元，得到一个数组VIARR；计算所述VIARR的最大值MaxVi，对VIARR按观测天顶角VZA升序排序，并依次将VIARR中的元素与MaxVi对比，当VIARR中第i个元素VIARR_i与MaxVi满足如下关系时，$\frac{\left(MaxVi-{\mathrm{VIARR}}_i\right)}{MaxVi}<{\Delta }_{\theta }$则选取VIARR_i作为最终的合成值；其中，Δ_θ为最大值容忍度，取值范围在0.1～0.3之间。