[发明专利]一种基于时序遥感影像的水体变化自动识别方法

权利要求书

1.一种基于时序遥感影像的水体变化自动识别方法,其特征在于，包括以下步骤：

步骤S01：采集研究区历年的植被指数时序数据，并建立植被指数时序数据集；

步骤S02：根据得到的植被数据时序数据集，逐年提取植被丰度指标；

步骤S03：根据研究区最大化合成MODIS遥感影像波段反射率数据，通过缨帽变换建立亮度、绿度和湿度指标时序数据集；

步骤S04：根据得到的亮度、绿度和湿度指标时序数据集，建立干扰指数时序数据集；

步骤S05：计算干扰指数、植被丰度指标和湿度指标的变化趋势；

步骤S06：根据得到的干扰指数、植被丰度和湿度指标的变化趋势，获得研究区水体变化分布图；

所述步骤S03具体为：基于2001-2017年8天最大化合成MODIS遥感影像波段1-7反射率数据，通过缨帽变换计算亮度、绿度和湿度三个指标，建立2001-2017年亮度、绿度和湿度指标时序数据集，其中亮度B、绿度G和湿度W指标的计算公式分别为：

B=0.3956*Band1 + 0.4718*Band2 + 0.3354*Band3 + 0.3834*Band4 + 0.3946*Band5 + 0.3434* Band6+0.2964* Band7

G= -0.3399* Band1 + 0.5952*Band2 - 0.2129*Band3 - 0.2222*Band4 + 0.4617*Band5- 0.1037* Band6 - 0.4600* Band7

W=0.10839 *Band1+ 0.0912*Band2 + 0.5065*Band3 + 0.4040*Band4 - 0.2410*Band5- 0.4658* Band6 - 0.5306* Band7

其中波段Band1，Band2，Band3，Band4，Band5，Band6和Band7分别为MODIS影像波段1-7对应的波段反射率；

所述步骤S04具体为：

步骤S401:对亮度、绿度和湿度指标进行标准化处理;

亮度、绿度和湿度指标的标准化计算方法依次为：

B_r= (B-B_u)/ B_ǒ

G_r= (G-G_u)/ G_ǒ

W_r= (W-W_u)/ W_ǒ

其中，B_r、G_r和W_r分别是经过标准化后的亮度、绿度和湿度指标；B_u, G_u, W_u分别为研究区内亮度、绿度、湿度指标的平均值；B_ǒ, G_ǒ, W_ǒ分别为研究区内亮度、绿度、湿度指标的标准差；

步骤S402:利用标准化后的亮度、绿度和湿度指标，建立干扰指数DI；

DI=B_r-(G_r+W_r)

其中，DI为干扰指数，Br、Gr和Wr分别是经过标准化后的亮度、绿度和湿度指标；逐像元逐年逐期计算干扰指数，生成2001-2017年研究区干扰指数时序数据集。