[发明专利]一种基于时空加权的生产高时空分辨率NDVI的方法

权利要求书

1.一种基于时空加权的生产高时空分辨率NDVI的方法，用于将高空间分辨率和低时间分辨率的Landsat NDVI数据与低空间分辨率和高时间分辨率的MODIS NDVI数据融合，从基准时刻t₁获得预测时刻t_p的高时空分辨率的NDVI数据，其特征在于，所述方法包括如下步骤：

步骤A：以同一基准时刻t₁为起点，利用MODIS NDVI数据和Landsat NDVI数据，分别获得预测时刻t_p的每个像元在NDVI上的第一增量△NDVI_d和第二增量△NDVI_s；其中，

所述第一增量△NDVI_d指的是，从基准时刻t₁到预测时刻t_p，假定MODIS NDVI数据的每个像元的NDVI增量，等于同样面积下的Landsat NDVI数据的不同地物类别的像元的NDVI增量的平均值，则利用t₁时刻的MODIS NDVI数据和Landsat NDVI数据，以及t_p时刻的MODISNDVI数据，计算获得t₁时刻到t_p时刻的Landsat NDVI数据的每个像元的NDVI增量，该增量即为所述第一增量△NDVI_d；

所述第二增量△NDVI_s指的是，采用现有ENVI/IDL软件中集成的薄板样条方法工具，直接将t₁时刻与t_p时刻的低空间分辨率MODIS像元分别插值到高空间分辨率的Landsat像元上去，然后使用t_p时刻的插值结果与t₁时刻的插值结果相减，其结果即为每个像元的所述第二增量△NDVI_s；

步骤B，使用贝叶斯模型平均方法，将步骤A获得的每个高空间分辨率的像元在NDVI上的所述第一增量与所述第二增量进行权重组合计算，获得每个像元的NDVI综合增量；

步骤C，通过步骤B计算得到的t₁时刻到t_p时刻，每个像元的所述NDVI综合增量，加到t₁时刻对应的高空间分辨率NDVI数值上，即可得到t_p时刻的高空间分辨率NDVI数值。

2.如权利要求1所述的方法，其特征在于，所述步骤B中，每个高空间分辨率像元(x_i，y_j)的所述NDVI综合增量用公式表示为，

ΔNDVI_com(x_i,y_j)＝w_s*ΔNDVI_s(x_i,y_j)+w_d*ΔNDVI_d(x_i,y_j) (1)

其中，公式(1)中，w_d和w_s分别表示所述第一增量和所述第二增量的权重，△NDVI_com(x_i，y_j)表示所述第一增量与所述第二增量结合得到的所述NDVI综合增量。

3.如权利要求2所述的方法，其特征在于，所述公式(1)中的每个像元的所述权重w_d和w_s通过贝叶斯模型平均方法求解获得。

4.如权利要求2或3所述的方法，其特征在于，所述步骤B进一步包括对所述NDVI综合增量进行误差修正的步骤。