一种基于超像素与模糊C均值聚类SAR图像分割方法.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN110349160A(43)申请公布日2019.10.18(21)申请号201910555710.3(22)申请日2019.06.25(71)申请人电子科技大学地址611731四川省成都市高新区（西区）西源大道2006号(72)发明人陈彦陈云坪冉崇敬(74)专利代理机构成都行之专利代理事务所(普通合伙)51220代理人温利平(51)Int.Cl.G06T7/10(2017.01)G06K9/62(2006.01)权利要求书4页说明书9页附图4页(54)发明名称一种基于超像素与模糊C均值聚类SAR图像分割方法(57)摘要本发明公开了一种基于超像素与稀疏表示的模糊C均值聚类SAR图像分割方法，首先生成超像素，之后对图像中每一子区域提取灰度和纹理特征作为基础特征，并在稀疏表示理论基础上根据SAR图像不同地物类别散射特性差异提出稀疏自表示矩阵校正处理办法，以得到准确的判别特征，最终实现对SAR图像中相干斑噪声具有较强鲁棒性且运行效率高的图像分割处理；由于在超像素层面进行相应图像处理可在保留图像内部信息和边界信息的基础上，通过像素集合的整体性，减弱相干斑噪声的影响，同时通过整合相邻像素信息使提取的特征更为稳定。CN110349160ACN110349160A权利要求书1/4页1.一种基于超像素与稀疏表示的模糊C均值聚类SAR图像分割方法，其特征在于，包括以下步骤：(1)、获取SAR图像获取待观测区域内包含合适且丰富地物类型的单极化合成孔径雷达SAR图像；(2)、使用PILS算法生成超像素(2.1)、初始化聚类中心设SAR图像共计N个像素点，初始化k个聚类中心，每个聚类中心为一个超像素，相邻的聚类中心的距离为将每个聚类中心用三维特征表示为：[i,x,y]，其中，i表示聚类中心的像素强值，x与y表示聚类中心的横纵坐标值；(2.2)、聚类中心重定位以某一聚类中心为中心，找到n×n邻域内梯度最小位置处，将该聚类中心重定位至该位置处；(2.3)、图像聚类以某一像素点为中心，建立搜索域S×S；遍历整幅SAR图像，确定该像素点在搜索域S×S内距离最近的聚类中心，并将该像素点划分为该聚类中心的同一类；(2.4)、迭代更新分配完所有像素点的类别属性后，计算各类别的聚类中心的三维特征[i,x,y]的均值，并以该均值更新该类别的聚类中心，并迭代至收敛；(3)、超像素特征提取在搜索域S×S内计算灰度共生矩阵GLCM提取出超像素纹理特征，包括：角二阶矩、对比度、相关性、逆差矩、方差、和均值、和熵、和方差、熵、差分熵、差分方差、相关性信息测度和最大相关系数；提取超像素内的灰度直方图作为超像素灰度特征；将超像素纹理特征和灰度特征构成图像基本特征集X＝(x1,x2,…,xn)，xn表示第n个特征；(4)、创建稀疏自表示模型(4.1)、初始化Z＝G＝Ψ＝0，Z为稀疏自表示矩阵，G为对角矩阵，Ψ为拉格朗日乘子矩阵；初始化当前迭代次数t＝1；(t)(t)(4.2)、判断第t次迭代后是否满足||Z-G||∞＜η，如果满足，则进入步骤(4.3)，否则，输出第t次迭代后得到稀疏自表示矩阵Z；(4.3)、逐列更新G(t)：其中，j＝1,2,…,J，gj＝(g1j,g2j,…,g(j-1)j,0,…,gJj)，μ为惩罚系数，zj为稀疏自表示(t)(t)(t)矩阵Z的第j列，lj为矩阵Ψ的第j列；将gjj设置为0，以满足约束条件diag(G)＝0，因而在逐列更新时，当遇到G(t)的对角线元素时，不予更新；2CN110349160A权利要求书2/4页更新Z(t)：其中，X为图像基本特征集，I为单位矩阵，γ为设置参数；对更新后的Z(t)进行校正：(t)对Z中的元素zij按如下公式进行进行处理，得到元素zij的估计值将所有元素的估计值构成矩阵遍历计算每一超像素集合的后向散射均值，令(t)超像素i和j的后向散射值的均值差为Δσij，得到校正后Z的元素为：更新Ψ(t)：Ψ(t)＝Ψ(t-1)+μ(Z(t)-G(t))更新惩罚系数μ：μ＝min(εμ,1010)其中，ε为常数；(4.4)、待步骤(4.3)完成后，输出稀疏自表示矩阵Z(t)，并将迭代次数t自加1，返回步骤(4.2)；(5)、图像分割(5.1)、定义SSR_FCM算法的目标函数Jm：其中，m表示模糊指数，vk、分别表示矩阵X和Z的第k类聚类中心，β为用于控制稀疏判别特征在目标函数Jm中所占权重值，uki为引入的模糊隶属度参数，表示超像素i对于第k个类别的隶属度值；(5.2)、对目标函数Jm进行最小化：其中，ξ表示拉格朗日乘子的微分；(5.3)、求解最小化后的目标函数Jm，得到：3CN110349160A权利要求书3/4页(5.4)、代入约束可求解出：进而得到