[发明专利]一种基于直方图智能阈值和SRM辅助的遥感图像水体提取方法

权利要求书

1.一种基于直方图智能阈值和SRM辅助的遥感图像水体提取方法，其特征在于，包括以下步骤：

a.选择遥感图像中对水体光谱具有诊断性的波段图像A，并对图像A统计直方图；

b.用直方图中水体直方图波峰和其右侧其他地物直方图波峰间的波谷最低处作为区分和提取水体的阈值n；

c.采用阈值n，对图像A进行阈值处理，获得水体信息数据A1；

d.根据遥感图像A成像时的太阳高度角（α）、方位角（a）和DEM数据B，由公式（1）计算地形光照系数，得到地形光照阴影数据B1；

B1=cosθpcos(π/2-α)+sinθpsin(θπ/2-α)sin(a-θ0) （式1）

其中：

θ_p 为坡度角，θ₀坡向，均由DEM数据计算，a:太阳方位角；α：太阳高度角

e.对地形光照阴影数据B1进行直方图统计处理，采用b方法，以山体阴影区直方图波峰和其右侧山体迎光区直方图波峰间的波谷最低处作为区分和提取山体阴影的阈值m；

f.以阈值m为参数，对阴影数据B1进行阈值处理，得到其相应的地形阴影图斑数据B2；

g.将水体信息数据A1同地形阴影图斑数据B2做逻辑与运算，得到剔除山体阴影信息的水体信息数据A2。

2.根据权利要求1所述的一种基于直方图智能阈值和SRM辅助的遥感图像水体提取方法，其特征在于，选择遥感图像中对水体光谱具有诊断性的波段图像A，图像A为水体反射率最低的波段，通常为近红外波段，利用水体在近红外波段上具有的反射率较低，易与其他地物区分的特点，选取单一的红外波段影像显示，进行图像A的水体提取。

3.根据权利要求1所述的一种基于直方图智能阈值和SRM辅助的遥感图像水体提取方法，其特征在于，智能阈值处理是指自动选取直方图两波峰之间的波谷最低处的灰度值作为阈值。

4.根据权利要求1所述的一种基于直方图智能阈值和SRM辅助的遥感图像水体提取方法，其特征在于，水体信息数据A1的数据格式为：水体像元值为1，非水体像元素值为0，将图像中阈值部分标记为水体，其余部分则标记为背景，由此获得水体范围的二值图。

5.根据权利要求1所述的一种基于直方图智能阈值和SRM辅助的遥感图像水体提取方法，其特征在于，SRM（地形阴影）计算是指在图像A成像时的太阳高度角（α）和方位角（a）等信息因素参与下的地形阴影数据的辅助计算，用于地形阴影数据B2的计算。

6.根据权利要求1所述的一种基于直方图智能阈值和SRM辅助的遥感图像水体提取方法，其特征在于，地形阴影数据B2格式为：地形阴影区像元值为0，非阴影区像元值为1，水体信息数据A2格式同数据A1。

7.根据权利要求1所述的一种基于直方图智能阈值和SRM辅助的遥感图像水体提取方法，其特征在于，A1 and -B2指：水体信息数A1同地形阴影数据B2做逻辑与，即and运算，用于去除水体提取中山体阴影的影响。