[发明专利]一种基于水体指数的亚像元级水体提取方法

权利要求书

1.一种基于水体指数的亚像元级水体提取方法，其特征在于，包括以下步骤：

1)水陆混合像元提取；

2)水陆端元光谱选取；

3)水体丰度估计；

所述的2)具体为：

利用遥感影像上邻近地物相似性原理，动态选取邻近纯净水体和邻近最优纯净陆地光谱作为水陆混合像元中水体端元和陆地端元的光谱；

所述的3)具体为：利用全约束线性混合像元分解的方法对水陆混合像元中水体百分比含量进行估计；

以亚像元级细小水体为提取目标，以Landsat OLI多光谱遥感影像为实验数据，对亚像元级水体提取方法展开探索，实现对亚像元级细小水体的准确提取；

所述的水陆混合像元提取具体为：

在对水陆混合像元进行水体百分数含量进行估计之前，首先要提取出含有水体的混合像元，在像元级水体提取方法中，利用水体指数对遥感影像进行波段计算时，水体和陆地在水体指数值图像直方图统计中会形成两个区分度明显的波峰，两个波峰之间有一段过渡带，在对像元级水体提取时，选取过渡带中心作为水体阈值进行水陆分割，水体指数的原理是根据水体和陆地的光谱特征对水体信息进行增强，同时对陆地信息进行抑制，以此在影像上形成区分度明显的两类像元；对于水体指数值图像直方图统计两个峰值周围的像元，其像元光谱特征明显，为纯净的水体或纯净陆地像元，而在两峰之间过渡带上的像元，由于其像元光谱特征混淆，水体指数难以实现对其进行增强或抑制的效果，通过实验发现，该过渡带像元为水陆混合像元，选择水陆区分性能最好的水体指数MNDWI来进行水陆混合像元的提取；

水体丰度估计具体为：

在线性混合像元分解模型中，混合像元光谱被认为是由各端元光谱根据端元在混合像元中所占比重线性加权组合而成，在本次亚像元级水体提取研究中，同样将混合像元光谱看作是由水体端元与陆地端元线性组合而成，同时需满足两个约束条件：(1)端元丰度之和为1，(2)端元丰度取值区间为0～1，混合像元分解模型定义如下：

ρ_m＝ρ_l×f_l+ρ_w×f_w+σ (2)

f_l+f_w＝1 (3)

0≤f_w,f_l≤1 (4)

其中ρ_l、ρ_w和ρ_m分别表示陆地端元、水体端元以及混合像元的反射率，f₁和f_w表示陆地端元和水体端元的丰度，σ表示残差，根据公式(2)(3)，水体丰度即为：

|σ|＝min＜0.1 (6)

在公式(5)中，ρ_m为已知量，ρ_w和ρ_l根据ρ_w和ρ_l根据水陆混合像元提取方法获得，当满足残差σ的一维范数最小且最小值小于0.1时，所求f_w为混合像元水体丰度值，若残差σ的一维范数最小但不小于0.1，则令混合像元水体丰度值为0。

2.根据权利要求1所述的一种基于水体指数的亚像元级水体提取方法，其特征在于，所述的水体指数MNDWI具体为：

其中ρ(Green)表示地物在绿光波段上的反射率，ρ(MIR)表示地物在中红外波段上的反射率。