[发明专利]一种洪灾变化检测方法、装置和存储介质

权利要求书

1.一种洪灾变化检测方法，其特征在于，所述方法包括：

获取待检测地区的多源时序遥感图像序列，所述多源时序遥感图像序列包括多张光学遥感图像与多张雷达遥感图像；

对所述多源时序遥感图像序列中的每个遥感图像进行预处理；

根据预处理后的多源时序遥感图像序列，生成所述多源时序遥感图像序列中的每个遥感图像各自对应的洪灾检测结果；

根据所述每个遥感图像各自对应的洪灾检测结果，利用所述多源时序遥感图像序列中遥感图像间的时序互补特性，通过时序增强算法生成所述待检测地区的多源时序遥感图像序列对应的洪灾变化检测结果，所述检测结果包括洪水区域与非洪水区域。

2.根据权利要求1所述的方法，其特征在于，根据预处理后的多源时序遥感图像序列，生成所述多源时序遥感图像序列中的每个遥感图像各自对应的洪灾检测结果，包括：

基于第一算法，根据预处理后的多源时序遥感图像序列，确定所述每个遥感图像对应的初始检测结果；

基于第二算法，对所述每个遥感图像对应的初始检测结果进行迭代，直到达到标定迭代状态或迭代次数等于标定次数时，输出每个遥感图像各自对应的洪灾检测结果；

其中，所述第一算法用于基于所述每个遥感图像整体确定初始检测结果；所述第二算法用于提取所述每个遥感图像中的局部洪水特征。

3.根据权利要求2所述的方法，其特征在于，根据预处理后的多源时序遥感图像序列，确定所述每个遥感图像对应的初始检测结果，包括：

通过分级模糊C均值聚类算法，确定所述预处理后的多源时序遥感图像序列每个遥感图像对应的初始化设置的动态阈值矩阵；

根据所述每个遥感图像对应的初始化设置的动态阈值矩阵，确定每个遥感图像对应的初始空间聚类结果，将所述初始空间聚类结果作为所述每个遥感图像对应的初始检测结果。

4.根据权利要求3所述的方法，其特征在于，对所述每个遥感图像对应的初始检测结果进行迭代，直到达到标定迭代状态或迭代次数等于标定次数时，输出每个遥感图像各自对应的洪灾检测结果，包括：

根据所述每个遥感图像的初始检测结果中洪水区域的轮廓边界，更新所述每个遥感图像对应的动态阈值矩阵；

根据更新后的所述每个遥感图像对应的动态阈值矩阵，更新所述每个遥感图像对应的空间聚类结果；

若所述每个遥感图像的更新前后的空间聚类结果对应的稀疏马尔可夫随机场的能量函数的差值小于0或迭代次数等于标定次数时，停止迭代，并将本次迭代过程更新的空间聚类结果作为所述每个遥感图像各自对应的洪灾检测结果。

5.根据权利要求4所述的方法，其特征在于，根据更新后的所述每个遥感图像对应的动态阈值矩阵，更新所述每个遥感图像对应的空间聚类结果之后，所述方法还包括：

基于稀疏扩张邻域空间，计算所述每个遥感图像的更新前后的空间聚类结果对应的稀疏马尔可夫随机场的能量函数；

根据每个遥感图像的更新前后的空间聚类结果对应的稀疏马尔可夫随机场的能量函数，计算能量函数衰减速率；

根据所述能量函数衰减速率，对当前稀疏扩张邻域空间的窗口大小进行调节。

6.根据权利要求1所述的方法，其特征在于，根据所述每个遥感图像各自对应的洪灾检测结果，利用所述多源时序遥感图像序列中遥感图像间的时序互补特性，通过时序增强算法生成所述待检测地区的多源时序遥感图像序列对应的洪灾变化检测结果，包括：

对于所述多源时序遥感图像序列中任一遥感图像，根据排列在该遥感图像之前的多张遥感图像，确定历史洪灾信息；

计算所述历史洪灾信息与该遥感图像中每个连通区域的重叠比；

根据所述重叠比与该遥感图像中洪水区域的比例，确定该遥感图像的可靠的洪水区域，并确定该遥感图像对应的时序增强结果；

根据所述多源时序遥感图像序列中所有图像各自对应的时序增强结果，确定所述多源时序遥感图像序列对应的洪灾变化检测结果。