[发明专利]基于模糊化高斯隶属函数的高分辨率遥感影像分割方法

说明书

技术领域

本发明属于图像处理领域，具体是一种基于模糊化高斯隶属函数的高分辨率遥感影像分割方法。

背景技术

影像分割是遥感影像处理中最为基本和关键的任务，一直以来是图像处理的热点和难点问题。高分辨率遥感影数据由于能够更加清晰细致的呈现地表覆盖信息，因此在地物目标的精准分割中具有巨大的潜力与优势；同时，高分辨率的特征也带来了新的分割问题：(1)不同类别地物分布曲线重叠范围大，使像素类属的不确定性增加；(2)同一类别区域局部特征数据差异性显著，从而引起建模的不确定性。

高分辨率遥感影像的上述两种不确定性特征可用“模糊化隶属函数”表征，该函数以待分割影像的原始隶属函数为基础，通过模糊化原始隶属函数参数的方式实现模糊化隶属函数的建模。相对于只能处理像素类属的不确定性的传统不确定性建模方法，如FCM方法，模糊化隶属函数可以同时刻画像素类属的不确定性及建模的不确定性，因此，其处理不确定性的能力更强。

合理应用模糊化隶属所提供的信息构建目标函数是实现高分辨率遥感影像正确划分的核心问题。应用模糊化隶属函数提供的信息建模神经网络模型，并将该模型作为目标函数是最常用的方法，该建模思想是将所有类别的模糊化隶属函数信息以加权和的形式联系起来，用权重符号判断信息在目标函数中是否发挥作用(权重为正表示该隶属函数信息被激活，权重为负和零则表示被抑制)，通过权重大小判断其发挥作用的大小(权重越大发挥的作用越大)，相对于传统建模方法，该建模方法使目标函数的质量得到很大提高。但是，上述建模目标函数的过程中没有考虑到像素间的空间相关性对分割结果的影响，故对噪声敏感。

发明内容

针对现有技术的不足，本发明提出一种基于模糊化高斯隶属函数的高分辨率遥感影像分割方法。

本发明的技术方案是：

一种基于模糊化高斯隶属函数的高分辨率遥感影像分割方法，包括以下步骤：

步骤1：读取待分割的高分辨率遥感影像；

步骤2：对待分割的高分辨率遥感影像中的每个地物类别进行监督采样提取训练样本，计算训练样本中每个像素的灰度值在对应的地物类别中出现的频率值；

步骤3：利用监督采样提取的训练样本及其在对应的地物类别中出现的频率值，在高分辨率遥感影像的灰度测度空间对不同地物类别建立高斯隶属函数模型；

步骤4：在灰度测度范围内，利用高斯隶属函数模型计算各个地物类别所有像素灰度值的高斯隶属函数值，并对各地物类别中的所有像素的灰度值的高斯隶属函数值进行归一化，使每个地物类别内的各高斯隶属函数值之和为1；

步骤5：高斯隶属函数模型参数模糊化，得到模糊化隶属函数：保持高斯隶属函数模型系数不变，初始化高斯隶属函数均值的调节因子或初始化高斯隶属函数标准差的调节因子，将高斯隶属函数均值或高斯隶属函数标准差进行区间取值得到模糊化高斯隶属函数，进而确定模糊化高斯隶属函数的上边界、下边界；

步骤6：以所有训练样本的高斯隶属函数值及模糊化高斯隶属函数值的上边界、下边界值作为输入，建立用于描述高分辨率遥感影像不同地物类别像素灰度值分布特征的线性神经网络模型，作为高分辨率遥感影像的目标函数，将该模型融入空间关系得到高分辨率遥感影像的目标函数矩阵；

步骤7：按最大隶属度原则划分高分辨率遥感影像的目标函数矩阵，实现高分辨率遥感影像分割；

步骤8：按照设定步长改变步骤5中的调节因子，即改变模糊化隶属函数均值或标准差的上边界、下边界，并重复步骤5至步骤7，对所有分割结果进行比较，取最优分割作为最终结果。

所述步骤3具体步骤如下：

步骤3.1：设定高斯隶属函数模型参数的初始值及上限、下限；高斯隶属函数模型参数，包括高斯隶属函数系数、高斯隶属函数均值、高斯隶属函数标准差；

步骤3.2：以已知训练样本在对应的地物类别中出现的频率值与对应的高斯隶属函数值的最小误差平方和为标准，设定迭代阈值；

步骤3.3：利用最小二乘法进行地物类别直方图拟合，建立每个地物类别的高斯隶属函数模型，循环迭代求解高斯隶属函数模型参数；

步骤3.4：当训练样本在对应的地物类别中出现的频率值与对应的高斯隶属函数值的误差平方和大于迭代阈值时，将当前的高斯隶属函数模型参数作为初始值重复步骤3.2～步骤3.4，循环迭代至所述误差平方和小于迭代阈值时迭代结束，此时的高斯隶属函数模型参数为最佳高斯隶属函数模型参数，此时的高斯隶属函数模型为最佳高斯隶属函数模型。

所述步骤6具体步骤如下：