基于非局部的三马尔可夫随机场SAR图像分割方法.pdf

(19)中华人民共和国国家知识产权局*CN102903102A*(12)发明专利申请(10)申请公布号CN102903102A(43)申请公布日2013.01.30(21)申请号201210334910.4(22)申请日2012.09.11(71)申请人西安电子科技大学地址710071陕西省西安市太白南路2号(72)发明人侯彪焦李成牛佳颖马文萍张向荣王爽(74)专利代理机构陕西电子工业专利中心61205代理人王品华朱红星(51)Int.Cl.G06T7/00(2006.01)权利要求书权利要求书3页3页说明书说明书66页页附图附图22页(54)发明名称基于非局部的三马尔可夫随机场SAR图像分割方法(57)摘要本发明公开了一种基于非局部的三马尔可夫随机场SAR图像分割方法，它属于图像处理技术领域，主要解决传统TMF方法应用于SAR图像分割易产生区域一致性差，边界凌乱的问题。其实现步骤为：(1)输入待分割图像；(2)利用FCM聚类初始化每个像素类标;(3)利用k-means初始化每个像素的场景类别，并利用非局部冗余信息对场景类别进行一次更新迭代；(4)计算图像的势能；(5)构建三马尔可夫场联合分布，并用吉布斯采样器对该分布函数采样，得到后验概率；(6)计算后验边缘概率，并逐点更新各像素点类标；(7)判断各像素类标变化率是否大于阈值，若变化率大于阈值，重复(4)-(6),否则输出分割结果。本发明收敛速度快，分割结果区域一致性好,保留信息完整，可用于SAR图像目标识别。CN10293ACN102903102A权利要求书1/3页1.一种基于非局部的三马尔可夫随机场SAR图像分割方法，包括如下步骤：1)输入待分割的SAR图像；2)利用FCM聚类的方法得到待分割图像每个像素点的初始类标；3)提取待分割图像的灰度共生矩阵，采用k-means聚类方法得到待分割图像每个像素点的场景类别，并利用图像非局部冗余信息对场景类别进行一次迭代更新，得到每个像素点新的场景类别；4)根据步骤2)得到的每个像素点的初始类标和步骤3)得到的每个像素点新的场景类别，利用下式计算势能W(x,u)：式中，x为像素点的类标，u为像素点的场景类别，s、t为一对相邻像素点，xs、xt为一对相邻像素点的类标值，us、ut为一对相邻像素点的场景类别，CH为图像中水平方向上的相邻像素点对的集合，CV为图像中垂直方向上的相邻像素点对的集合，a,b为场景类别的两种取值，和分别表示不同的权重参数，为类标的惩罚函数，为场景类别取a时的惩罚函数，为场景类别取b时的惩罚函数；5)根据步骤2)中得到的每个像素点的类标，计算各像素点的条件概率：0其中，p(ys|xs)服从统计模型G分布，ys为像素点s的灰度值，Γ为Gamma函数，n为SAR图像的等效视数，n通过图像的先验知识获得，α为形状参数，γ为尺度参数；6)根据上述得到的势能W(x,u)和各像素点的条件概率p(ys|xs)，更新各像素点的类标；6a)根据势能W(x,u)和各像素点的条件概率p(ys|xs)，使用统计概率公式计算三马尔可夫场联合概率分布；6b)利用三马尔可夫场联合概率分布，获得各像素点的后验边缘概率；6c)利用贝叶斯最大后验边缘概率准则逐点更新图像中各像素点的类标；7)将更新前后类标变化的像素点个数和图像中像素点总数的比率作为终止条件，如果比率大于输入的阈值10-9，或达到最大迭代次数100，返回步骤4)，否则将各像素点的最终类标作为最终分割结果输出。2.根据权利要求1所述的基于非局部的三马尔可夫随机场SAR图像分割方法，其中步骤3)所述的利用图像非局部冗余信息对场景类别进行一次迭代更新，按如下步骤进行：3a)提取待分割图像的灰度共生矩阵，采用k-means聚类估计每个像素点的初始场景类别；2CN102903102A权利要求书2/3页3b)根据步骤3a)得到的每个像素点s的初始场景类别，取邻域窗口m*m，m=11，分别计算窗口内每一类像素点对s的影响因子，取影响因子最大的一类作为s点更新后的场景类别，第k类像素点对s的影响因子的计算公式为：式中，t为邻域大窗口中属于第k类的像素点，v(Ns)代表以s为中心的邻域矢量，v(Nt)代表以t为中心的邻域矢量，为高斯核函数，Z(s)为归一化因子，其表达式为：其中，h是控制指数函数衰减速度的衰减因子，α为高斯核函数的标准差，α>0。3.根据权利要求1所述的基于非局部的三马尔可夫随机场SAR图像分割方法，其中步骤6a)所述的使用统计概率公式计算三马尔可夫场联合概率分布，按如下公式进行：式中，p(x,u,y)为类标x、场景类别u和灰度值y的三马尔可夫场联合概率分布，γ′为归一化常数，W(x,u)为势能，s为像素点，p(ys|xs)为条件概率，ys为像素点s的灰度值，xs为像素点s的类标。4.根