基于语义条件随机场模型的SAR图像分割方法技术

本发明专利技术公开了一种基于语义条件随机场模型的SAR图像分割方法。主要解决现有技术中不能保持图像细节信息的问题。其实现步骤为：1.根据SAR图像的区域图，将SAR图像划分为混合聚集结构地物子空间、结构区域子空间和匀质区域子空间；2.对混合聚集结构地物子空间采用词袋模型提取特征，用AP聚类的方法进行分割；3.构建语义条件随机场模型对结构区域子空间和匀质区域子空间进行分割；4.将混合聚集结构地物子空间、结构区域子空间和匀质区域子空间的分割结果合并，得到SAR图像的分割结果。本发明专利技术获得了SAR图像的良好分割效果，可用于SAR图像的语义分割。

SAR image segmentation method based on semantic conditional random field model

The invention discloses a SAR image segmentation method based on a semantic conditional random field model. The invention mainly solves the problem that the details of the image can not be maintained in the prior art. The method comprises the following steps: 1. according to the SAR image, SAR image is divided into mixed aggregation structure of object space, domain subspace and homogeneous region subspace; 2. pairs of hybrid aggregation structure of object subspace using the bag of words model feature extraction, using AP clustering method to construct semantic segmentation; 3. conditions random field model for regional structure subspace and subspace and homogeneous regions segmentation; 4. the results of mixed aggregation structure of object space, domain subspace and subspace and homogeneous regions segmentation with SAR image segmentation results. The invention obtains the good segmentation effect of the SAR image and can be used for the semantic segmentation of the SAR image.

【技术实现步骤摘要】
基于语义条件随机场模型的SAR图像分割方法
本专利技术属于图像处理
，特别涉及SAR图像分割方法，可用于图像分类、识别和检测。
技术介绍
随机场方法是SAR图像分割中一类比较受欢迎的方法。典型的随机场方法是马尔可夫随机场MRF模型，它是一种概率生成模型。在MRF模型中，后验概率等于似然概率和先验概率的乘积。似然概率描述了SAR图像的特征，通常用SAR图像的统计分布来表示，分布的选择主要根据SAR图像的特性。先验概率描述了图像的空间上下文信息，通常用Gibbs分布来表示。然而，MRF模型的假设中需要强烈的依赖关系，并且MRF模型的先验模型中没有考虑观测数据间的相互关系。针对MRF模型的上述不足，条件随机场CRF模型应用而生，其主要由一元势函数和二元势函数组成。它是一种判别性的模型，并且把后验概率直接定义为Gibbs分布。该CRF模型不仅捕获了单个像素点的信息和邻域像素的信息，而且还捕获了图像类标之间的相互作用和图像观测之间的相互作用。由于CRF模型具有比其他分割方法的优势，被广泛经用于SAR图像分割。但是由于原始的CRF模型没有考虑SAR图像本身的特性，针对此问题，张鹏等人提出了对原始CRF模型的一元势函数进行改进的一般CRF模型，其整合了SAR图像的纹理特征和SAR图像的统计特性。该一般CRF模型的二元势函数采用传统的多层逻辑斯蒂函数捕获图像上下文的信息。然而，该二元势函数只捕获了图像空间上下文中各向同性的关系，忽视了SAR图像本身的各向异性的关系，造成分割结果中细节信息的丢失，且分割结果不具有语义一致性，影响后续对SAR图像的分类、识别和检测。专利
技术实现思路
本专利技术的目的在于针对上述已有方法的不足，提出一种基于语义条件随机场模型的SAR图像分割方法，以提升SAR图像分割的效果。本专利技术的技术思路是：通过对一般的CRF模型进行改进，改善SAR图像分割的效果，即根据SAR图像的素描图和区域图组成的语义空间捕获SAR图像各向异性的关系，将语义空间的信息嵌入到一般的CRF模型中，构建语义条件随机场模型，其实现步骤如下：(1)根据SAR图像的区域图，将SAR图像划分为混合聚集结构地物子空间、结构区域子空间和匀质区域子空间；(2)对混合聚集结构地物子空间，采用词袋模型提取区域的特征，并用仿射传播AP聚类的方法对该地物子空间进行分割；(3)构建语义条件随机场模型：(3a)定义一元势函数为：其中，us(xs,ys)是一元势函数，Z＝{1,2,...,N}是整个SAR图像像素的集合，N是SAR图像像素的总数目；ys是SAR图像中第s个像素，s∈Z；xs是SAR图像中第s个像素的类标；fs(ys)是像素的类别特征函数，p(ys)是像素的类别特征概率，p(ys|xs)是似然概率；K是图像类别的总数目，l∈{1,2,...,K}，δ(xs，l)为第一指示函数，如果xs＝l，则δ(xs，l)＝1，如果xs≠l，则δ(xs，l)＝0；(3b)定义二元势函数为：其中，yt是SAR图像中第t个像素，t∈Z，且s≠t；xt是SAR图像中第t个像素的类标；l∈{1,2,...,K}，m∈{1,2,...,K}，δ(xt，m)为第二指示函数，如果xt＝m，则δ(xt，m)＝1，如果xt≠m，则δ(xt，m)＝0；ψst(ys,yt,φ(ys,yt))是混合核函数，φ(ys,yt)是像素ys和像素yt之间的语义函数，p(xs,xt|ψst(ys,yt,φ(ys,yt)))是先验概率；(3c)根据(3a)和(3b)，定义语义条件随机场模型的目标函数如下：p(xs|ys)＝p(ys)×p(ys|xs)×p(xs,xt|ψst(ys,yt,φ(ys,yt)))<3>其中，p(ys)是像素的类别特征概率，p(ys|xs)是似然概率，p(xs,xt|ψst(ys,yt,φ(ys,yt)))是先验概率；(4)采用(3c)得到的语义条件随机场模型的目标函数对结构区域子空间和匀质区域子空间进行分割，即对结构区域子空间和匀质区域子空间的每个像素，取式<3>的最大值，得到每个像素的类标为：(5)将混合聚集结构地物子空间、结构区域子空间和匀质区域子空间的分割结果进行合并，得到SAR图像的分割结果。本专利技术与现有技术相比具有如下优点：第一、本专利技术将SAR图像像素空间的信息和语义空间的信息相结合，进行SAR图像分割，能有效的完成了SAR图像分割的任务。第二、本专利技术通过建立语义条件随机场模型，捕获了SAR图像中各项异性的关系，不仅提高了分割结果的区域一致性，而且有效的保留了图像的细节信息。附图说明图1是本专利技术对SAR图像分割的实现流程图；图2是本专利技术中对SAR图像子空间的划分结果图；图3是用本专利技术与现有方法对Ku波段分辨率为1米的SAR图像分割结果图；图4是用本专利技术与现有方法对C波段分辨率为3米的SAR图像分割结果图。具体实施方式参照图1，本专利技术是根据SAR图像的区域图，将SAR图像划分为混合聚集结构地物子空间、结构区域子空间和匀质区域子空间；对混合聚集结构地物子空间采用词袋模型提取区域的特征，然后用AP聚类的方法得到该地物子空间的分割结果；对结构区域子空间和匀质区域子空间，构建语义条件随机场进行分割；语义条件随机场的一元势函数采用多项式逻辑斯蒂回归函数和SAR图像的统计特性表示；语义条件随机场的二元势函数采用基于混合核函数的多项式逻辑斯蒂回归函数表示；将混合聚集结构地物子空间、结构区域子空间和匀质区域子空间的分割结果合并，得到SAR图像的分割结果，具体实施步骤如下:步骤1，根据SAR图像的区域图，将SAR图像划分为混合聚集结构地物子空间、结构区域子空间和匀质区域子空间。(1.1)根据Jie-Wu和Fang-Liu等人于2014年发表在IEEETransactionsonGeoscienceandRemoteSensing杂志上的文章《LocalmaximalhomogenousregionsearchforSARspecklereductionwithsketch-basedgeometricalkernelfunction》中所提出的模型得到SAR图像素描图；(1.2)根据SAR图像素描图中素描线段的聚集度，将素描线划分为聚集的素描线和非聚集的素描线；(1.3)构造一个半径为最优聚集度区间上界的圆形基元，用该圆形基元对聚集的素描线进行膨胀，对膨胀后的线段集合由外向内进行腐蚀，在素描图上得到聚集区域；(1.4)对非聚集的素描线，构造大小为5×5的几何结构窗，得到结构区域；(1.5)将素描图中除去聚集区域和结构区域以外的部分作为不可素描区域；(1.6)将素描图中的聚集区域、结构区域和不可素描区域，映射到SAR图像上，将SAR图像划分为混合聚集结构地物子空间、结构区域子空间和匀质区域子空间，如图2所示，图2(a)为原始的SAR图像，图2(b)为SAR图像素描图，图2(c)中白色区域为混合聚集结构地物子空间，灰色区域为结构区域子空间，黑色区域为匀质区域子空间。步骤2，对混合聚集结构地物子空间采用词袋模型提取区域的特征，然后用仿射传播AP聚类的方法对该地物子空间进行分割。对混合聚集结构地物子空间中的每个像素ya，是混合聚集结构地物子空间像素本文档来自技高网...

【技术保护点】
基于语义条件随机场模型的SAR图像分割方法，其特征在于：包括如下步骤：(1)根据SAR图像的区域图，将SAR图像划分为混合聚集结构地物子空间、结构区域子空间和匀质区域子空间；(2)对混合聚集结构地物子空间，采用词袋模型提取区域的特征，并用仿射传播AP聚类的方法对该地物子空间进行分割；(3)构建语义条件随机场模型：(3a)定义一元势函数为：

【技术特征摘要】
1.基于语义条件随机场模型的SAR图像分割方法，其特征在于：包括如下步骤：(1)根据SAR图像的区域图，将SAR图像划分为混合聚集结构地物子空间、结构区域子空间和匀质区域子空间；(2)对混合聚集结构地物子空间，采用词袋模型提取区域的特征，并用仿射传播AP聚类的方法对该地物子空间进行分割；(3)构建语义条件随机场模型：(3a)定义一元势函数为：其中，us(xs,ys)是一元势函数，Z＝{1,2,...,N}是整个SAR图像像素的集合，N是SAR图像像素的总数目；ys是SAR图像中第s个像素，s∈Z；xs是SAR图像中第s个像素的类标；fs(ys)是像素的类别特征函数，p(ys)是像素的类别特征概率，p(ys|xs)是似然概率；K是图像类别的总数目，l∈{1,2,...,K}，δ(xs，l)为第一指示函数，如果xs＝l，则δ(xs，l)＝1，如果xs≠l，则δ(xs，l)＝0；(3b)定义二元势函数为：其中，yt是SAR图像中第t个像素，t∈Z，且s≠t；xt是SAR图像中第t个像素的类标；l∈{1,2,...,K}，m∈{1,2,...,K}，δ(xt，m)为第二指示函数，如果xt＝m，则δ(xt，m)＝1，如果xt≠m，则δ(xt，m)＝0；ψst(ys,yt,φ(ys,yt))是混合核函数，φ(ys,yt)是像素ys和像素yt之间的语义函数，p(xs,xt|ψst(ys,yt,φ(ys,yt)))是先验概率；(3c)根据(3a)和(3b)，定义语义条件随机场模型的目标函数如下：p(xs|ys)＝p(ys)×p(ys|xs)×p(xs,xt|ψst(ys,yt,φ(ys,yt)))<3>其中，p(ys)是像素的类别特征概率，p(ys|xs)是似然概率，p(xs,xt|ψst(ys,yt,φ(ys,yt)))是先验概率；(4)采用(3c)得到的语义条件随机场模型的目标函数对结构区域子空间和匀质区域子空间进行分割，即对结构区域子空间和匀质区域子空间的每个像素，取...

【专利技术属性】
技术研发人员：刘芳，段一平，李婷婷，焦李成，郝红侠，陈璞华，马晶晶，尚荣华，
申请(专利权)人：西安电子科技大学，
类型：发明
国别省市：陕西,61