[发明专利]一种基于环境要素的湖库稳定性分区方法

权利要求书

1.一种基于环境要素的湖库稳定性分区方法，其特征在于，包括：

获取多个年份包含完整目标湖库的区域遥感影像，计算每个年份区域遥感影像的地表反射率，得到每个年份包含区域地表反射率的遥感影像；

根据所述每个年份包含区域地表反射率的遥感影像，提取每个年份区域内目标湖库水面，生成每个年份区域内目标湖库的水面范围矢量数据；

根据所述每个年份区域内目标湖库的水面范围矢量数据，通过自然间断法，得到目标湖库多年岸线复杂度指数的分类栅格影像；

根据所述每个年份包含区域地表反射率的遥感影像，通过自然间断法，得到目标湖库多年水面平均温差的分类栅格影像；

对所述目标湖库多年岸线复杂度指数的分类栅格影像和所述目标湖库多年水面平均温差的分类栅格影像进行叠加，通过自然间断法得到所述目标湖库的稳定性分区结果。

2.根据权利要求1所述的湖库稳定性分区方法，其特征在于，所述多个年份之间间隔相等。

3.根据权利要求1所述的湖库稳定性分区方法，其特征在于，所述区域遥感影像在每个年份中的相同月份进行采集；若任一年份未采集到与其他年份相同月份的遥感影像，则该年份采用与所述月份相邻月份的遥感影像。

4.根据权利要求1所述的湖库稳定性分区方法，其特征在于，所述根据所述每个年份包含区域地表反射率的遥感影像，提取每个年份区域内目标湖库水面，生成每个年份区域内目标湖库的水面范围矢量数据，包括：

1)对任一年份包含区域地表反射率的遥感影像，分别计算每个栅格的改进的归一化差异水体提取MNDWI指数和自动水面提取AWEIsh指数，得到该年份分别对应MNDWI指数计算结果的栅格图像和对应AWEIsh指数计算结果的栅格图像；

2)根据步骤1)得到的栅格图像，分别统计MNDWI指数值和AWEIsh值在对应栅格图像中的分布情况，得到每个指数分别对应的双峰分布图；将每张双峰分布图中双峰间的中值作为该年份对应指数的阈值；

3)根据该年份对应AWEIsh指数计算结果的栅格图像，将该图像中每个栅格对应的AWEIsh指数值与该年份AWEIsh指数阈值进行比较：

若任一栅格对应的AWEIsh指数值小于该年份AWEIsh指数阈值，则判定该栅格为初始水面范围；否则，判定该栅格为非水面范围；

对该年份对应AWEIsh指数计算结果的栅格图像的所有栅格判定完毕后，保留所有判定为初始水面范围的栅格；

4)根据该年份对应MNDWI指数计算结果的栅格图像，获取步骤3)保留的所有栅格对应的MNDWI指数值并判定：

若步骤3)保留的任一栅格对应的MNDWI指数值小于该年份MNDWI指数阈值，则判定该栅格为水面范围；否则，判定该栅格为非水面范围；

判定完毕后，该年份所有判定为水面范围的栅格组成该年份区域内目标湖库的水面范围矢量数据；

5)重复步骤1)至4)，得到每个年份区域内目标湖库的水面范围矢量数据。