CN109697466A - 一种自适应区间型空间模糊c均值的地物分类方法 - Google Patents

Abstract

发明公开了一种自适应区间型空间模糊C均值的地物分类方法，属于数据挖掘与遥感图像处理交叉领域。本发明设计了一种新的自适应的区间模型，利用像素点的最大隶属度控制区间宽度变化，避免了初始参数的依赖性，提高了数据的可分性；定义了新的空间相关性度量方法，综合考量局部空间内像素点的隶属度与空间距离，在控制“噪声”的同时，保留了图像的细节信息；基于空间邻域信息模型，设计了新的区间型模糊C均值算法的类别中心与隶属度的计算方法。自适应区间型空间模糊C均值的地物分类方法具有更好的分类结果，对多光谱遥感影像数据具有良好的准确性与鲁棒性，相比于现有的模糊C均值的方法精度有了提升。

Description

一种自适应区间型空间模糊C均值的地物分类方法

技术领域

本发明属于数据挖掘、遥感图像处理交叉领域，是一种基于自适应区间建模与空间邻域信息的模糊C均值的遥感影像地物分类方法。

背景技术

遥感影像分类是遥感地物信息提取的关键手段，其最终目的是对遥感图像数据中的每一个像素点赋予唯一的类别标志。目前，遥感影像分类的主要问题是：①多光谱遥感影像中的像元信息往往是多种地表成分的复合，存在一定的模糊性，因此对于影像中细节的划分非常困难。②复杂的自然环境以及遥感光谱相互作用，使像素光谱测度的空间相关性越来越复杂，产生许多类属不确定性的“噪声”。多数分类方法对于“噪声”的敏感性强、鲁棒性差，而处理“噪声”较好的分类方法，往往会产生细节丢失的问题。由于以上的问题，导致了已有分类方法不适用于具有严重“噪声”与细节丰富的遥感影像分类。

发明内容

为解决现有技术中的上述问题，本发明提供了一种自适应区间型空间模糊C均值的地物分类方法进行遥感地物分类。本发明的分类方法结合区间建模方法与空间邻域信息，设计了一种新的自适应的区间模型，并定义了新的空间相关性度量方法。

发明公开了一种自适应区间型空间模糊C均值的地物分类方法，属于数据挖掘与遥感图像处理交叉领域。本发明设计了一种新的自适应的区间模型，利用像素点的最大隶属度控制区间宽度变化，避免了初始参数的依赖性，提高了数据的可分性；定义了新的空间相关性度量方法，综合考量局部空间内像素点的隶属度与空间距离，在控制“噪声”的同时，保留了图像的细节信息；基于空间邻域信息模型，设计了新的区间型模糊C均值算法的类别中心与隶属度的计算方法。技术思路如下：

1)初始化各分类参数，把多光谱遥感数据的每个像素点作为分类特征向量；

2)由各像素点的局部空间偏差确定各像素点的最大区间宽度，利用各像素点的最大隶属度作为调节因子，对待测数据进行动态区间计算；

3)结合隶属度与空间距离，计算各局部空间内的空间相关性；基于空间相关性模型，分别计算类别中心区间的上、下边界；

4)利用最大相异性度量方法，计算待测数据与类别中心的距离；结合空间邻域信息，计算隶属度；

5)重复2)-4)过程，直到满足终止条件，分类完成；

6)保持待测数据的隶属度不变，用最大隶属度规则对待测样本点进行类别归一化。

具体的技术方案如下所述：

一种自适应区间型空间模糊C均值的地物分类方法(本发明中简称AIVSFCM方法)，包括如下步骤：

1)按照需求确定类别数c，模糊加权指数m，迭代次数最大值T，阈值ε，初始化迭代次数t＝1；设定待测多光谱遥感影像数据为n维矢量空间数据集，即待测数据集X＝{x₁，x₂，…，x_i，…，x_N}，x_i∈Rⁿ，N是矢量特征的个数，R是实数集，随机初始化待测数据集X对类别质心V＝{v₁，v₂，…，v_c}，(v_i＝{v_i1，v_i2，…，v_in}，i＝1，2，…，c)的隶属度矩阵U＝{u₁，u₂，…，u_c}(u_i＝{u_i1，u_i2，…，u_iN}，i＝1，2，…，c)；

2)计算待测数据集的自适应区间值，计算公式如式(1)所示，

其中，u_i表示像素点x_i隶属于各个类别的隶属度矢量，max(u_i)表示隶属度的最大值，ω_i表示x_i所处局部空间(如3×3邻域)的偏差；

3)更新中心的上、下边界，计算公式分别如式(2)、(3)所示，

其中，N_i表示以x_i为中心的局部空间，||N_i||表示局部空间内像素点个数，u_ki表示像素点x_i隶属于第k个簇的隶属度，SC_ij(k)表示像素点x_i和x_j属于第k类别的空间相关性，计算公式如(4)所示，

其中，D_ij表示两像素点的空间距离；

4)更新隶属度u_ki，如式(5)所示，

其中，表示矢量区间数据与矢量区间中心的距离，计算公式如(6)所示，

5)如果或t＜T，终止迭代，否则t＝t+1，返回2)；

6)根据最大隶属度原则，对每个像素点P_i进行类别划分，如式(7)所示，

P_i＝arg_k∈[1，c]{max(u_ki)} (7)

优选，步骤3)的D_ij表示两像素点的空间距离为欧氏距离。

上述一种自适应区间型空间模糊C均值的地物分类方法，其中，多光谱遥感影像为同一传感器、多个波段的遥感影像。

本发明的创新在于设计了一种新的自适应的区间模型，利用像素点的最大隶属度控制区间宽度变化，避免了初始参数的依赖性，提高了数据的可分性；定义了新的空间相关性度量方法，综合考量局部空间内像素点的隶属度与空间距离，在控制“噪声”的同时，保留了图像的细节信息；基于空间邻域信息模型，设计了新的区间型模糊C均值算法的类别中心与隶属度的计算方法。分析如下：步骤2)新的区间建模方法将像素点最大隶属度作为控制区间宽度变化的参数，通过对像素点隶属度大小的调整控制像素点的区间宽度，避免了对初始参数的依赖性，即步骤2)的建模不需要对不同遥感图像数据进行特定的参数初始设定；步骤3)结合隶属度与空间距离设计了新的空间相关性度量模型，在控制敏感“噪声”的同时，保留了图像的细节信息；步骤4)结合空间邻域信息和计算隶属度设计了最大相异性度量方法，考虑了数据各维度间的各向异性，取不同维度之间的最大值，进一步提高了数据的可分离性，因此本专利的方法适用于多光谱的区间遥感图像数据，相比于多光谱遥感数据中单一的度量方式难以刻画光谱近似或混叠像素的最大分离性具有显著的进步。

本发明结合自适应区间建模与空间领域信息，通过新的自适应区间建模方法，解决遥感影像数据分类中数据的难分或不可分问题；引入空间邻域信息，不仅可以减少“噪声”的影响，而且较好的平衡了图像的细节保留问题，最终改善了分类结果。

本发明的遥感影像分类方法具有更好的分类结果，对多光谱遥感影像数据具有良好的准确性与鲁棒性，相比于现有的模糊C均值的方法精度有了提升。

附图说明

图1(a)是广东省珠海市横琴岛区域原图像SPOT5多光谱遥感影像假彩色合成图；

图1(b)是广东省珠海市横琴岛区域SPOT5多光谱遥感影像FCM方法的分类结果图；

图1(c)是广东省珠海市横琴岛区域SPOT5多光谱遥感影像ADFLICM方法的分类结果图；

图1(d)是广东省珠海市横琴岛区域SPOT5多光谱遥感影像IVFCM方法的分类结果图；

图1(e)是广东省珠海市横琴岛区域SPOT5多光谱遥感影像AIVFCM方法的分类结果图；

图1(f)是广东省珠海市横琴岛区域SPOT5多光谱遥感影像AIVSFCM方法的分类结果图。

其中，图1(a)-图1(f)中每幅图下方的三幅小图选取的是存在明显分类差异的区域放大图(分别标为区域1、区域2、区域3)。

图2(a)是北京市昌平区域原图像SPOT5多光谱遥感影像假彩色合成图；

图2(b)是北京市昌平区域SPOT5多光谱遥感影像FCM方法的分类结果图；

图2(c)是北京市昌平区域SPOT5多光谱遥感影像ADFLICM方法的分类结果图；

图2(d)是北京市昌平区域SPOT5多光谱遥感影像IVFCM方法的分类结果图；

图2(e)是北京市昌平区域SPOT5多光谱遥感影像AIVFCM方法的分类结果图；

图2(f)是北京市昌平区域SPOT5多光谱遥感影像AIVSFCM方法的分类结果图。

其中，图1(a)-图1(f)中每幅图下方的两幅小图选取的是存在明显分类差异的区域放大图(分别标为区域1、区域2)，图2(a)-图2(f)中每幅图下方的三幅小图选取的是存在明显分类差异的区域放大图(分别标为区域1、区域2、区域3)。

图3为图1(a)-图1(f)和图2(a)-图2(f)中地形标识说明。

具体实施方式

实施例

针对选自广东省珠海市横琴岛区域(北到磨刀门，南到二井国，西到滩尾角，东到大、小横琴岛中间区域)和北京市昌平区域(北到秦上路，南到温榆河，西到翠明湖，东到吕各庄)的影像数据，对地物类型复杂、光谱混叠的多光谱遥感影像分类，具体实施方式如下所述：

1)按照需求确定类别数c，模糊加权指数m，迭代次数最大值T，阈值ε，初始化迭代次数t＝1；设定待测多光谱遥感影像数据是一个n维矢量空间数据集，即待测数据集X＝{x₁，x₂，…，x_i，…，x_N}，x_i∈Rⁿ，N是矢量特征的个数。随机初始化待测数据集X对类别质心V＝{v₁，v₂，…，v_c}，(v_i＝{v_i1，v_i2，…，v_in}，i＝1，2，…，c)的隶属度矩阵U＝{u₁，u₂，…，u_c}(u_i＝{u_i1，u_i2，…，u_iN}，i＝1，2，…，c)；

2)计算待测数据集的自适应区间值，计算公式如式(1)所示，

其中，u_i表示像素点x_i隶属于各个类别的隶属度矢量，max(u_i)表示隶属度的最大值，ω_i表示x_i所处局部空间的偏差；

3)更新中心的上、下边界，计算公式分别如式(2)、(3)所示，

其中，N_i表示以x_i为中心的局部空间，||N_i||表示局部空间内像素点个数，u_ki表示像素点x_i隶属于第k个簇的隶属度，SC_ij(k)表示像素点x_i和x_j属于第k类别的空间相关性，计算公式如(4)所示，

其中，D_ij表示两像素点空间坐标的欧氏距离；

4)更新隶属度u_ki，如式(5)所示，

其中，表示矢量区间数据与矢量区间中心的距离，计算公式如(6)所示，

5)如果或t＜T，终止迭代，否则t＝t+1，返回2)；

6)根据最大隶属度原则，对每个像素点进行类别划分，如式(7)所示，

P_i＝arg_k∈[1，c]{max(u_ki)} (7)

本例中设置参数为：c＝5，m＝2，T＝300，ε＝10^(-5)。

仿真试验如下：

1)遥感影像的地物分类

实验数据选取了具有高阶模糊性、干扰性、地物复杂的广东省珠海市横琴岛区域与北京市昌平区域的多光谱遥感影像数据。

分类结果的评价计算方法和度量标准如下：

①总体分类精度(OA，Overall Accuracy)：指被正确分类的类别像元数与总的类别个数的比值，计算方法如式(8)所示。

其中，N表示像素总个数，x_ii表示混淆矩阵对角线元素，即各类别结果的正确个数。若P_e越大，则表明分类精度越高，反之亦然。

②Kappa系数(Cohen's kappa)。Kappa系数是一种比例，代表着分类与完全随机的分类产生错误减少的比例，用于衡量分类精度，计算方法如式(9)所示。

其中，N表示像素总个数，x_ii表示混淆矩阵对角线元素，即各类别结果的正确个数，x_i+、x_+i分别表示混淆矩阵各行、各列之和。若K越大，则表示分类精度越高，反之亦然。

参照的分类方法分别为：一型模糊C均值方法(FCM)、自适应模糊局部信息C均值方法(ADFLICM)、区间型模糊C均值方法(IVFCM)、自适应区间型模糊C均值方法(AIVFCM)、自适应区间型空间模糊C均值方法(AIVSFCM)。

仿真1：广东省珠海市横琴岛区域遥感影像分类

本组实验数据的高分辨率遥感影像选自广东省珠海市横琴岛的部分区域(北到磨刀门，南到二井国，西到滩尾角，东到大、小横琴岛中间区域)，地物类型如下表1所示，土地覆盖了林地、水域、草地、建筑用地、农业用地5个类别。多光谱遥感影像选自SPOT5号卫星拍摄的相应地区的分辨率为10m的多光谱图像，波谱范围为0.43～0.89μm，由源影像数据1、2、3三个波段合成的假彩色图像(400×400pixels)，如图1(a)所示。表2给出了不同方法分类结果的统计参数对比，参考图像为待分类的多光谱遥感影像。

FCM(图1(b))的分类结果中，难以区分光谱近似的地物类别，如区域1中的鱼塘错分为水域与裸地，区域2中的水库面积划分明显变小，水库边缘淹没在裸地中。ADFLICM(图1(c))出现了与FCM(图1(b))相同的错分问题，但是由于ADFLICM结合了局部空间信息，分类结果明显更加平滑紧凑，“噪声”更少。IVFCM(图1(d))使用区间建模的方法提高了数据的可分性，在区域1与区域2的水域划分中得到了较为理想的结果，边界较为清晰，但是静态的区间宽度并不能提高所有像元的可分性，因此在区域1中还有少许水域错分为裸地。AIVFCM(图1(e))则得到了区域1与区域2中更加清晰完整的划分结果，这是因为AIVFCM采用了动态的区间建模方法。AIVSFCM(图1(f))相较于AIVFCM(图1(e))，综合了区间建模与空间信息的优势，分类结果中类别间边界清晰、类内更加紧凑、“噪声”更少。

表1广东省珠海市横琴岛区域SPOT5遥感影像数据地物类别描述

表2实验方法聚类结果统计

仿真2：北京市昌平区域遥感影像分类

本组实验采用的数据是北京市昌平区域影像数据(北到秦上路，南到温榆河，西到翠明湖，东到吕各庄)，地物类型如下图表3所示，土地覆盖了林地、水域、草地、建筑用地、农业用地5个类别。多光谱遥感影像选自SPOT5号卫星拍摄的相应地区的分辨率为10m的多光谱图像，波谱范围为0.43～0.89μm，由源影像数据1、2、3三个波段合成的假彩色图像(800×600pixels)，如图2(a)所示。表4给出了不同方法分类结果的统计参数对比，参考图像为待分类的多光谱遥感影像。

FCM(图2(b))对于分布复杂、光谱混叠地物的分类性能较差，如区域1中草地被错分为林地，区域2中的受林地阴影影响将林地错分为水域，在区域3中出现了与区域1、区域2相同的错分问题，上方的实线箭头处受林地阴影影响严重。ADFLICM(图2(c))在区域1中草地的错分变少，受区域2、区域3中林地阴影影响较小，但是出现了由于过分考虑空间信息导致的细节消失问题，如下方实线箭头处。IVFCM(图2(d))静态的区间建模并没有得到更好的分类结果，如草地与林地错分现象严重。AIVFCM(图2(e))得到了更为准确的划分结果，但是结果中“噪声”较多。AIVSFCM(图2(f))得到了最佳的分类结果，更好的平衡了图像细节与类内异质性，如区域1中裸地清晰可见，错分为林地的像素点更少，区域2中受林地阴影影响更小，区域3中准确抑制了林地阴影的影响，下方的实线箭头处划分更加平滑且保留了细节信息。

表3北京市昌平区域SPOT5遥感影像数据地物类别描述

表4实验方法聚类结果统计

由仿真1和仿真2的实验结果，从目视效果(图1(a)-图1(f)和图2(a)-图2(f))和客观指标统计(表2和表4)两个方面分别对原始遥感影像和分类后图像进行对比和分析，可以得出如下结论：

(1)FCM对分布复杂、光谱混叠地物的分类性能较差，ADFLICM由于过分考虑空间邻域信息易出现细节消失问题与错分问题，IVFCM对于不同的遥感影像数据适应性不强，没有整体考虑对像素点的区间建模，AIVFCM易受“噪声”影响，分类效果不理想。

(2)目视来看，相比FCM、ADFLICM、IVFCM、AIVFCM、AIVSFCM的遥感影像分类方法，AIVSFCM得到分类后图像的类别与土地利用图最相似，这说明本发明的方法(AIVSFCM)的分类结果更精确，更好的平衡了图像细节与类内异质性的矛盾，类别紧凑，边界明显，具有一定的抗噪能力。

(3)由表2和表4可见，各方法的分类总体精度与KAPPA系数相近，AIVSFCM指标数值最高。

Claims (2)

1.一种自适应区间型空间模糊C均值的地物分类方法，其特征在于包括如下步骤：

1)按照需求确定类别数c，模糊加权指数m，迭代次数最大值T，阈值ε，初始化迭代次数t＝1；设定待测多光谱遥感影像数据为n维矢量空间数据集，即待测数据集X＝{x₁，x₂，…，x_i，…，x_N}，x_i∈Rⁿ，N是矢量特征的个数，R是实数集，随机初始化待测数据集X对类别质心V＝{v₁，v₂，…，v_c}，(v_i＝{v_i1，v_i2，...，v_in}，i＝1，2，…，c)的隶属度矩阵U＝{u₁，u₂，...，u_c}(u_i＝{u_i1，u_i2，...，u_iN}，i＝1，2，…，c)；

2)计算待测数据集的自适应区间值，计算公式如式(1)所示，

其中，u_i表示像素点x_i隶属于各个类别的隶属度矢量，max(u_i)表示隶属度的最大值，ω_i表示x_i所处局部空间的偏差；

3)更新中心的上、下边界，计算公式分别如式(2)、(3)所示，

其中，N_i表示以x_i为中心的局部空间，||N_i||表示局部空间内像素点个数，u_ki表示像素点x_i隶属于第k个簇的隶属度，SC_ij(k)表示像素点x_i和x_j属于第k类别的空间相关性，计算公式如(4)所示，

其中，D_ij表示两像素点的空间距离；

4)更新隶属度u_ki，如式(5)所示，

其中，表示矢量区间数据与矢量区间中心的距离，计算公式如(6)所示，

5)如果或t＜T，终止迭代，否则t＝t+1，返回2)；，

6)对每个像素点进行类别划分，如式(7)所示，

P_i＝arg_k∈[1，c]{max(u_ki)} (7)

2.根据权利要求1所述的方法，其特征在于D_ij表示两像素点的空间距离为欧氏距离。