[发明专利]基于自动差分聚类算法的遥感图像分割方法

权利要求书

1.一种基于自动差分聚类算法的遥感图像分割方法，包括以下步骤：

(1)输入待分割图像I，分别提取待分割图像I的小波特征向量和纹理特征向量，并用小波特征向量和纹理特征向量表示待分割图像I的每一个像素点v；

(2)产生聚类数据：

用待分割图像I的膨胀变换减去其腐蚀变换，得到形态梯度图像I₁，计算形态梯度图像I₁的灰度值矩阵的平方积，得到反映图像边缘的浮点活动图像I₂，用分水岭方法初分割浮点活动图像I₂，得到不同的图像块，对每一个图像块的所有像素点的特征取平均值，获得代表初始聚类数据的每一块的特征值；

(3)随机抽取聚类数据初始化种群：

种群中每个个体编码长度为L+L×d，每个个体分为两个部分，前L位为标签位，后L×d位为L个待激活的聚类中心，d为聚类数据的维数，随机初始化每个个体的标签位，并随机抽取L个聚类数据作为待激活的聚类中心，令当前迭代次数t=1；

(4)根据每个个体的标签位激活相应个体的聚类中心；

判断个体标签T每一位值的大小是否大于0.5，如果大于0.5，则激活相应的聚类数据点，否则不予激活；所有被激活的聚类数据点组合成为该个体的聚类中心；

(5)根据激活的聚类中心，采用PBM有效性指标公式计算每个个体适应度值；

(6)用改进的差分进化方法对种群进行变异和交叉，并利用个体适应度值进行种群更新；

(7)对更新后的种群进行类别数振荡操作：

7a)对种群中的每个个体依照振荡规则进行类别数振荡操作，每个个体获得新的类别数；

7b)计算每个个体中L个数据点的密度，并依照密度大小进行排序；

7c)将新的类别数k_new与当前染色体类别数k_old进行比较，如果k_new＜k_old，则依照密度排名从该个体中选择密度较大的聚类数据点加入聚类中心，否则，从已有聚类中心中淘汰密度排名较小的聚类数据点。

(8)利用FCM方法对振荡操作后的聚类中心进行更新；

(9)判断种群迭代的次数是否达到使用者设定的最大进化代数T=100，或者相邻两代最优个体的适应度差值是否小于规定数值E=10^-4，如果是，则保存最优个体并执行步骤(10)，否则，返回步骤(4)；

(10)对获得的最优个体进行解码，求得每个聚类数据的类别标号，并用这些类别标号作为像素的灰度值，输出图像的分割结果。

2.根据权利要求1所述的基于自动差分聚类算法的遥感图像分割方法，其特征在于，所述步骤1中提取待分割图像I的小波特征向量，是采用小波分解方法获得，即对图像进行窗口大小为16×16的三层小波变换，得到由子带系数所构成的10维小波特征向量。

3.根据权利要求1所述的基于自动差分聚类算法的遥感图像分割方法，其特征在于，所述步骤1中提取待分割图像I的纹理特征向量，采用灰度共生矩阵方法获得，具体步骤如下：

第1步，将待处理图像量化为16个灰度级；

第2步，依次令两个像素点连线与横轴的方向夹角为0°、45°、90°和135°，按照下式分别计算四个方向的灰度共生矩阵：

P(i,j)＝#{(x₁,y₁),(x₂,y₂)∈M×N|f(x₁,y₁)＝r,f(x₂,y₂)＝s}

其中，P(i，j)为灰度共生矩阵在坐标(i，j)位置上的元素，#为集合{}的元素个数，(x₁，y₁)和(x₂，y₂)为距离等于1的两个像素点坐标，∈为集合中的属于符号，M×N为待处理图像的大小，|为概率论中的条件符号，r为(x₁，y₁)处像素点矢量化后的灰度值，s为(x₂，y₂)处像素点矢量化后的灰度值；

第3步，根据灰度共生矩阵分别选取该矩阵四个方向上的对比度、同质性和能量值，获得像素的12维纹理特征向量。